Investigation on the mechanism of crystallization of soluble protein in the presence of nonionic surfactant

The mechanism of crystallization of soluble, globular protein (lysozyme) in the presence of nonionic surfactant C8E4 (tetraoxyethylene glycol monooctyl ether) was examined using both static and dynamic light scattering. The interprotein interaction was found to be attractive in solution conditions that yielded crystals and repulsive in the noncrystallizing solution conditions. The validity of the second virial coefficient as a criterion for predicting protein crystallization could be established even in the presence of nonionic surfactants. Our experiments indicate that the origin of the change in interactions can be attributed to the adsorption of nonionic surfactant monomers on soluble proteins, which is generally assumed to be the case with only membrane proteins. This adsorption screens the hydrophobic attractive force and enhances the hydration and electrostatic repulsive forces between protein molecules. Thus at low surfactant concentration, the effective protein-protein interaction remains repulsive. Large surfactant concentrations promote protein crystallization, possibly due to the attractive depletion force caused by the intervening free surfactant micelles.     

Protein specificity of charged sequences in polyanions and heparins

Long-range electrostatic interactions are generally assigned a subordinate role in the high-affinity binding of proteins by glycosaminoglycans, the most highly charged biopolyelectrolytes. The discovery of high and low sulfation domains in heparan sulfates, however, suggests selectivity via complementarity of their linear sulfation patterns with protein charge patterns. We examined how charge sequences in anionic/nonionic copolymers affect their binding to a protein with prominent charge anisotropy. Experiments and united-atom Monte Carlo simulations, together with Delphi electrostatic modeling for the protein, confirm strongest binding when polyanion sequences allow for optimization of repulsive and attractive electrostatics. Simulations also importantly identified retention of considerable polyion conformational freedom, even for strong binding. The selective affinity for heparins of high and low charge density found for this protein is consistent with nonspecific binding to distinctly different protein charge domains. These findings suggest a more nuanced view of specificity than previously proposed for heparinoid-binding proteins.     

Quantitative characterization of the lateral distribution of membrane proteins within the lipid bilayer.

The dependence of the lateral distribution of membrane proteins on the size, protein/lipoid molar ratio, and the magnitude of the interaction potentials has been investigated by computer modeling protein-lipid distributions with Monte Carlo calculations. These results have allowed us to develop a quantitative characterization of the distribution of membrane proteins and to correlate these distributions with experimental observables. The topological arrangement of protein domains, protein plus annular lipid domains, and free lipid domains is described in terms of radial distribution, pair connectedness, and cluster distribution functions. The radial distribution functions are used to measure the distribution of intermolecular distances between protein molecules, whereas the pair connectedness functions are used to estimate the physical extension of compositional domains. It is shown that, at characteristic protein/lipid molar ratios, previously isolated domains become connected, forming domain networks that extend over the entire membrane surface. These changes in the lateral connectivity of compositional domains are paralleled by changes in the calculated lateral diffusion coefficients and might have important implications for the regulation of diffusion controlled processes within the membrane.     

Domain formation in a fluid mixed lipid bilayer modulated through binding of the C2 protein motif

The role and mechanism of formation of lipid domains in a functional membrane have generally received limited attention. Our approach, based on the hypothesis that thermodynamic coupling between lipid-lipid and protein-lipid interactions can lead to domain formation, uses a combination of an experimental lipid bilayer model system and Monte Carlo computer simulations of a simple model of that system. The experimental system is a fluid bilayer composed of a binary mixture of phosphatidylcholine (PC) and phosphatidylserine (PS), containing 4% of a pyrene-labeled anionic phospholipid. Addition of the C2 protein motif (a structural domain found in proteins implicated in eukaryotic signal transduction and cellular trafficking processes) to the bilayer first increases and then decreases the excimer/monomer ratio of the pyrene fluorescence. We interpret this to mean that protein binding induces anionic lipid domain formation until the anionic lipid becomes saturated with protein. Monte Carlo simulations were performed on a lattice representing the lipid bilayer to which proteins were added. The important parameters are an unlike lipid-lipid interaction term and an experimentally derived preferential protein-lipid interaction term. The simulations support the experimental conclusion and indicate the existence of a maximum in PS domain size as a function of protein concentration. Thus, lipid-protein coupling is a possible mechanism for both lipid and protein clustering on a fluid bilayer. Such domains could be precursors of larger lipid-protein clusters ('rafts'), which could be important in various biological processes such as signal transduction at the level of the cell membrane.     

Calculation of resonance energy transfer in crowded biological membranes.

Analytical and numerical models were developed to describe fluorescence resonance energy transfer (RET) in crowded biological membranes. It was assumed that fluorescent donors were linked to membrane proteins and that acceptors were linked to membrane lipids. No restrictions were placed on the location of the donor within the protein or the partitioning of acceptors between the two leaflets of the bilayer; however, acceptors were excluded from the area occupied by proteins. Analytical equations were derived that give the average quantum yield of a donor at low protein concentrations. Monte Carlo simulations were used to generate protein and lipid distributions that were linked numerically with RET equations to determine the average quantum yield and the distribution of donor fluorescence lifetimes at high protein concentrations, up to 50% area fraction. The Monte Carlo results show such crowding always reduces the quantum yield, probably because crowding increases acceptor concentrations near donor-bearing proteins; the magnitude of the reduction increases monotonically with protein concentration. The Monte Carlo results also show that the distribution of fluorescence lifetimes can differ markedly, even for systems possessing the same average lifetime. The dependence of energy transfer on acceptor concentration, protein radius, donor position within the protein, and the fraction of acceptors in each leaflet was also examined. The model and results are directly applicable to the analysis of RET data obtained from biological membranes; their application should result in a more complete and accurate determination of the structures of membrane components.     

Effect of membrane microheterogeneity and domain size on fluorescence resonance energy transfer

Studies of multicomponent membranes suggest lateral inhomogeneity in the form of membrane domains, but the size of small (nanoscale) domains in situ cannot be determined with current techniques. In this article, we present a model that enables extraction of membrane domain size from time-resolved fluorescence resonance energy transfer (FRET) data. We expand upon a classic approach to the infinite phase separation limit and formulate a model that accounts for the presence of disklike domains of finite dimensions within a two-dimensional infinite planar bilayer. The model was tested against off-lattice Monte Carlo calculations of a model membrane in the liquid-disordered (l(d)) and liquid-ordered (l(o)) coexistence regime. Simulated domain size was varied from 5 to 50 nm, and two fluorophores, preferentially partitioning into opposite phases, were randomly mixed to obtain the simulated time-resolved FRET data. The Monte Carlo data show clear differences in the efficiency of energy transfer as a function of domain size. The model fit of the data yielded good agreement for the domain size, especially in cases where the domain diameter is <20 nm. Thus, data analysis using the proposed model enables measurement of nanoscale membrane domains using time-resolved FRET.     

The More the Tubular: Dynamic Bundling of Actin Filaments for Membrane Tube Formation

Tubular protrusions are a common feature of living cells, arising from polymerization of stiff protein filaments against a comparably soft membrane. Although this process involves many accessory proteins in cells, in vitro experiments indicate that similar tube-like structures can emerge without them, through spontaneous bundling of filaments mediated by the membrane. Using theory and simulation of physical models, we have elaborated how nonequilibrium fluctuations in growth kinetics and membrane shape can yield such protrusions. Enabled by a new grand canonical Monte Carlo method for membrane simulation, our work reveals a cascade of dynamical transitions from individually polymerizing filaments to highly cooperatively growing bundles as a dynamical bottleneck to tube formation. Filament network organization as well as adhesion points to the membrane, which bias filament bending and constrain membrane height fluctuations, screen the effective attractive interactions between filaments, significantly delaying bundling and tube formation.     

Protein shape and crowding drive domain formation and curvature in biological membranes

Folding, curvature, and domain formation are characteristics of many biological membranes. Yet the mechanisms that drive both curvature and the formation of specialized domains enriched in particular protein complexes are unknown. For this reason, studies in membranes whose shape and organization are known under physiological conditions are of great value. We therefore conducted atomic force microscopy and polarized spectroscopy experiments on membranes of the photosynthetic bacterium Rhodobacter sphaeroides. These membranes are densely populated with peripheral light harvesting (LH2) complexes, physically and functionally connected to dimeric reaction center-light harvesting (RC-LH1-PufX) complexes. Here, we show that even when converting the dimeric RC-LH1-PufX complex into RC-LH1 monomers by deleting the gene encoding PufX, both the appearance of protein domains and the associated membrane curvature are retained. This suggests that a general mechanism may govern membrane organization and shape. Monte Carlo simulations of a membrane model accounting for crowding and protein geometry alone confirm that these features are sufficient to induce domain formation and membrane curvature. Our results suggest that coexisting ordered and fluid domains of like proteins can arise solely from asymmetries in protein size and shape, without the need to invoke specific interactions. Functionally, coexisting domains of different fluidity are of enormous importance to allow for diffusive processes to occur in crowded conditions.     

Membrane-Mediated Interaction between Strongly Anisotropic Protein Scaffolds

Specialized proteins serve as scaffolds sculpting strongly curved membranes of intracellular organelles. Effective membrane shaping requires segregation of these proteins into domains and is, therefore, critically dependent on the protein-protein interaction. Interactions mediated by membrane elastic deformations have been extensively analyzed within approximations of large inter-protein distances, small extents of the protein-mediated membrane bending and small deviations of the protein shapes from isotropic spherical segments. At the same time, important classes of the realistic membrane-shaping proteins have strongly elongated shapes with large and highly anisotropic curvature. Here we investigated, computationally, the membrane mediated interaction between proteins or protein oligomers representing membrane scaffolds with strongly anisotropic curvature, and addressed, quantitatively, a specific case of the scaffold geometrical parameters characterizing BAR domains, which are crucial for membrane shaping in endocytosis. In addition to the previously analyzed contributions to the interaction, we considered a repulsive force stemming from the entropy of the scaffold orientation. We computed this interaction to be of the same order of magnitude as the well-known attractive force related to the entropy of membrane undulations. We demonstrated the scaffold shape anisotropy to cause a mutual aligning of the scaffolds and to generate a strong attractive interaction bringing the scaffolds close to each other to equilibrium distances much smaller than the scaffold size. We computed the energy of interaction between scaffolds of a realistic geometry to constitute tens of k_BT, which guarantees a robust segregation of the scaffolds into domains.     

Three separable domains regulate GTP-dependent association of H-ras with the plasma membrane

The microlocalization of Ras proteins to different microdomains of the plasma membrane is critical for signaling specificity. Here we examine the complex membrane interactions of H-ras with a combination of FRAP on live cells to measure membrane affinity and electron microscopy of intact plasma membrane sheets to spatially map microdomains. We show that three separable forces operate on H-ras at the plasma membrane. The lipid anchor, comprising a processed CAAX motif and two palmitic acid residues, generates one attractive force that provides a high-affinity interaction with lipid rafts. The adjacent hypervariable linker domain provides a second attractive force but for nonraft plasma membrane microdomains. Operating against the attractive interaction of the lipid anchor for lipid rafts is a repulsive force generated by the N-terminal catalytic domain that increases when H-ras is GTP loaded. These observations lead directly to a novel mechanism that explains how H-ras lateral segregation is regulated by activation state: GTP loading decreases H-ras affinity for lipid rafts and allows the hypervariable linker domain to target to nonraft microdomains, the primary site of H-ras signaling.     

Lateral interactions among membrane proteins. Implications for the organization of gap junctions.

We have studied the relationship between interprotein forces and the lateral distribution of proteins in disordered mouse liver gap junctions. Data on protein positions are obtained from freeze-fracture electron micrographs. Short-ranged correlations in observed positions are characteristic of interacting particles in a fluid state. An analysis derived from statistical mechanics allows the determination of the magnitude and functional form of interprotein forces. We find that jap junction proteins are mutually repulsive, in a manner consistent with electrostatics and excluded volume. This dictates that long-ranged protein aggregation into jap junction plaques cannot arise solely from interparticle interactions. An alternative is the balance of lateral pressures between the junction and the surrounding glycocalyx. This idea is quantified into a model. Junctional pressure arises from protein-protein interactions and is computed from a pressure equation based on the force and a radial distribution function describing order. The pressure from the glycocalyx is assumed to arise from mixing, electrostatic, and elastic interactions of sugar residues, and is described with terms from Flory-Krigbaum and McMillan-Mayer theories. The results of this modeling are in reasonable agreement with available experimental data.     

Why do proteins divide into domains? Insights from lattice model simulations

It is known that larger globular proteins are built from domains, relatively independent structural units. A domain size seems to be limited, and a single domain consists of from few tens to a couple of hundred amino acids. Based on Monte Carlo simulations of a reduced protein model restricted to the face centered simple cubic lattice, with a minimal set of short-range and long-range interactions, we have shown that some model sequences upon the folding transition spontaneously divide into separate domains. The observed domain sizes closely correspond to the sizes of real protein domains. Short chains with a proper sequence pattern of the hydrophobic and polar residues undergo a two-state folding transition to the structurally ordered globular state, while similar longer sequences follow a multistate transition. Homopolymeric (uniformly hydrophobic) chains and random heteropolymers undergo a continuous collapse transition into a single globule, and the globular state is much less ordered. Thus, the factors responsible for the multidomain structure of proteins are sufficiently long polypeptide chain and characteristic, protein-like, sequence patterns. These findings provide some hints for the analysis of real sequences aimed at prediction of the domain structure of large proteins.     

Statistical thermodynamics of membrane bending-mediated protein-protein attractions

Highly wedge-shaped integral membrane proteins, or membrane-adsorbed proteins can induce long-ranged deformations. The strain in the surrounding bilayer creates relatively long-ranged forces that contribute to interactions with nearby proteins. In contrast, to direct short-ranged interactions such as van der Waal's, hydrophobic, or electrostatic interactions, both local membrane Gaussian curvature and protein ellipticity can induce forces acting at distances of up to a few times their typical radii. These forces can be attractive or repulsive, depending on the proteins' shape, height, contact angle with the bilayer, and a pre-existing local membrane curvature. Although interaction energies are not pairwise additive, for sufficiently low protein density, thermodynamic properties depend only upon pair interactions. Here, we compute pair interaction potentials and entropic contributions to the two-dimensional osmotic pressure of a collection of noncircular proteins. For flat membranes, bending rigidities of approximately 100k(B)T, moderate ellipticities, and large contact angle proteins, we find thermally averaged attractive interactions of order k(B)T. These interactions may play an important role in the intermediate stages of protein aggregation. Numerous biological processes where membrane bending-mediated interactions may be relevant are cited, and possible experiments are discussed.     

High-resolution imaging of antibodies by tapping-mode atomic force microscopy: attractive and repulsive tip-sample interaction regimes

A force microscope operated with an amplitude modulation feedback (usually known as tapping-mode atomic force microscope) has two tip-sample interaction regimes, attractive and repulsive. We have studied the performance of those regimes to imaging single antibody molecules. The attractive interaction regime allows determination of the basic morphologies of the antibodies on the support. More importantly, this regime is able to resolve the characteristic Y-shaped domain structure of antibodies and the hinge region between domains. Imaging in the repulsive interaction regime is associated with the irreversible deformation of the molecules. This causes a significant loss in resolution and contrast. Two major physical differences distinguish the repulsive interaction regime from the attractive interaction regime: the existence of tip-sample contact and the strength of the forces involved.     

Molecular simulation of the effect of cholesterol on lipid-mediated protein-protein interactions

Experiments and molecular simulations have shown that the hydrophobic mismatch between proteins and membranes contributes significantly to lipid-mediated protein-protein interactions. In this article, we discuss the effect of cholesterol on lipid-mediated protein-protein interactions as function of hydrophobic mismatch, protein diameter and protein cluster size, lipid tail length, and temperature. To do so, we study a mesoscopic model of a hydrated bilayer containing lipids and cholesterol in which proteins are embedded, with a hybrid dissipative particle dynamics-Monte Carlo method. We propose a mechanism by which cholesterol affects protein interactions: protein-induced, cholesterol-enriched, or cholesterol-depleted lipid shells surrounding the proteins affect the lipid-mediated protein-protein interactions. Our calculations of the potential of mean force between proteins and protein clusters show that the addition of cholesterol dramatically reduces repulsive lipid-mediated interactions between proteins (protein clusters) with positive mismatch, but does not affect attractive interactions between proteins with negative mismatch. Cholesterol has only a modest effect on the repulsive interactions between proteins with different mismatch.     

Interaction Between Equally Charged Membrane Surfaces Mediated by Positively and Negatively Charged Macro-Ions

In biological systems, charged membrane surfaces are surrounded by charged molecules such as electrolyte ions and proteins. Our recent experiments in the systems of giant phospholipid vesicles indicated that some of the blood plasma proteins (macro-ions) may promote adhesion between equally charged membrane surfaces. In this work, theory was put forward to describe an IgG antibody-mediated attractive interaction between negatively charged membrane surfaces which was observed in experiments on giant phospholipid vesicles with cardiolipin-containing membranes. The attractive interactions between negatively charged membrane surfaces in the presence of negatively and positively charged spherical macro-ions are explained using functional density theory and Monte Carlo simulations. Both, the rigorous solution of the variational problem within the functional density theory and the Monte Carlo simulations show that spatial and orientational ordering of macro-ions may give rise to an attractive interaction between negatively charged membrane surfaces. It is also shown that the distinctive spatial distribution of the charge within the macro-ions (proteins) is essential in this process.     

Suppression of insulin aggregation by heparin

The aggregation of insulin near its isoelectric point and at low ionic strength was suppressed in the presence of heparin. To understand this effect, we used turbidimetry and stopped-flow to study the pH- and ionic strength ( I)-dependence of the aggregation of heparin-free insulin. The results supported the role of interprotein electrostatic interactions, contrary to the commonly held view that such forces are minimized at pH = pI. Electrostatic modeling of insulin (DelPhi) revealed that attractive interactions arise from the marked charge anisotropy of insulin near pI. We show how screening of the interprotein attractions by added salt lead to maximum aggregation near I = 0.01 M, corresponding to a Debye length nearly equal to the diameter of the insulin dimer, consistent with a dipole-like protein charge distribution. This analysis is also consistent with suppression of aggregation by heparin, a strong polyanion that by binding to the positive domain of one protein, inhibits its interaction with the negative domain of another.     

革兰氏阳性菌蛋白结构域特征分析

结构域是蛋白质序列中具有独特功能的区域,这些区域影响着蛋白质的功能,因此研究结构域的特征对于了解蛋白质功能很有帮助。构建革兰氏阳性菌蛋白质4个亚细胞位置数据集,对该数据集中的蛋白质进行结构域的搜索和功能分析,找到了革兰氏阳性菌的细胞壁、细胞质、细胞膜和细胞外四个蛋白质区域的结构域。分析这四个位置结构域的功能并在PDBsum数据库中找到了这些结构域的二级结构和三级结构图,利用这些特征信息可以更深入的了解革兰氏阳性菌蛋白质的结构和功能。     

Adaptive coarse-grained Monte Carlo simulation of reaction and diffusion dynamics in heterogeneous plasma membranes

Background  

An adaptive coarse-grained (kinetic) Monte Carlo (ACGMC) simulation framework is applied to reaction and diffusion dynamics in inhomogeneous domains. The presented model is relevant to the diffusion and dimerization dynamics of epidermal growth factor receptor (EGFR) in the presence of plasma membrane heterogeneity and specifically receptor clustering. We perform simulations representing EGFR cluster dissipation in heterogeneous plasma membranes consisting of higher density clusters of receptors surrounded by low population areas using the ACGMC method. We further investigate the effect of key parameters on the cluster lifetime.     

Simulation methods for protein structure fluctuations

Three numerical techniques for generating thermally accessible configurations of globular proteins are considered; these techniques are the molecular dynamics method, the Metropolis Monte Carlo method, and a modified Monte Carlo method which takes account of the forces acting on the protein atoms. The molecular dynamics method is shown to be more efficient than either of the Monte Carlo methods. Because it may be necessary to use Monte Carlo methods in certain important types of sampling problems, the behavior of these methods is examined in some detail. It is found that an acceptance ratio close to 1/6 yields optimum efficiency for the Metropolis method, in contrast to what is often assumed. This result, together with the overall inefficiency of the Monte Carlo methods, appears to arise from the anisotropic forces acting on the protein atoms due to their covalent bonding. Possible ways of improving the Monte Carlo methods are suggested.