Receptor aggregation by intermembrane interactions: a Monte Carlo study

The lateral organization of receptors on cell surfaces is critically important to their function; many receptors transmit transmembrane signals when redistributed into clusters, while the response of others is potentiated by their aggregation. Cell-cell contact can play a crucial role in receptor aggregation, even when the bonds between receptors on one cell and ligands on the other are monovalent. Monte Carlo simulations on a two-membrane model were carried out to determine whether weak enthalpic interactions among receptors in one membrane, and among ligands in another, can work synergistically to give large-scale clustering when the two membranes are brought into contact. The simulations give support to such a clustering mechanism. In addition, because clustering is a cooperative process akin to a phase separation, individual receptors and ligands may undergo repeated binding and unbinding while in a clustered "phase," and a single ligand could interact with multiple different receptor partners. The results suggest a resolution of the dichotomy between serial triggering and aggregation models of T cell activation.     

The relationship of DNA base composition and individual protein composition in micro-organisms

Summary To investigate the dependence of protein composition on DNA base composition, a set of data on individual proteins with known amino acid compositions from a spectrum of bacterial species has been compiled. It is found that similar relationships of amino acid frequency to G + C content exist for these proteins as for the bulk proteins studied by Sueoka (1961). The data are analysed by linear and cubic regression, and a measure of the proportions of A + T-rich and G + C-rich codons in the underlying messenger RNAs is put forward. The theoretical limits on the G + C content of coding DNA are discussed, and inference are made about the various selective forces acting on DNAs of different G + C contents.     

Helix-forming tendencies of nonpolar amino acids predicted by Monte Carlo simulated annealing

Monte Carlo simulated annealing is applied to the study of the α-helix-forming tendencies of seven nonpolar amino acids, Ala, Leu, Met, Phe, Ile, Val, and Gly. Homooligomers of 10 amino acids are used and the helix tendency is calculated by folding α-helicies from completely random initial conformations. The results of the simulation imply that Met, Ala, and Leu are helix formers and that Val, Ile, and Gly are helix breakers, while Phe comes in between the two groups. The differences between helix formers and breakers turned out to be large in agreement with the recent experiments with short peptides. It is argued from the energy distributions of the obtained conformations that the helix tendency is small for the helix breakers because of steric hindrance of side chains. Homoglycine is shown to favor a random coil conformation. The β-strand tendencies of the same homooligomers are also considered, and they are shown to agree with the frequencies of amino acids in β-sheet from the protein data base. © 1994 Wiley-Liss, Inc.     

Monte Carlo simulation studies on the prediction of protein folding types from amino acid composition

A number of methods to predicting the folding type of a protein based on its amino acid composition have been developed during the past few years. In order to perform an objective and fair comparison of different prediction methods, a Monte Carlo simulation method was proposed to calculate the asymptotic limit of the prediction accuracy [Zhang and Chou (1992),Biophys. J.63, 1523–1529, referred to as simulation method I]. However, simulation method I was based on an oversimplified assumption, i.e., there are no correlations between the compositions of different amino acids. By taking into account such correlations, a new method, referred to as simulation method II, has been proposed to recalculate the objective accuracy of prediction for the least Euclidean distance method [Nakashimaet al. (1986),J. Biochem.99, 152–162] and the least Minkowski distance method [Chou (1989),Prediction in Protein Structure and the Principles of Protein Conformation, Plenum Press, New York, pp. 549–586], respectively. The results show that the prediction accuracy of the former is still better than that of the latter, as found by simulation method I; however, after incorporating the correlative effect, the objective prediction accuracies become lower for both methods. The reason for this phenomenon is discussed in detail. The simulation method and the idea developed in this paper can be applied to examine any other statistical prediction method, including the computersimulated neural network method.     

Activation of Bax by joint action of tBid and mitochondrial outer membrane: Monte Carlo simulations

The mitochondrial pathway of apoptosis proceeds when molecules, such as cytochrome c, sequestered between the outer and inner mitochondrial membranes are released to the cytosol by mitochondrial outer membrane (MOM) permeabilization. Bax, a member of the Bcl-2 protein family, plays a pivotal role in mitochondrion-mediated apoptosis. In response to apoptotic stimuli, Bax integrates into the MOM, where it mediates the release of cytochrome c from the intermembrane space into the cytosol, leading to caspase activation and cell death. The pro-death action of Bax is regulated by interactions with both other prosurvival proteins, such as tBid, and the MOM, but the exact mechanisms remain largely unclear. Here, the mechanisms of integration of Bax into a model membrane mimicking the MOM were studied by Monte Carlo simulations preceded by a computer prediction of the docking of tBid with Bax. A novel model of Bax activation by tBid was predicted by the simulations. In this model, tBid binds to Bax at an interaction site formed by Bax helices α1, α2, α3 and α5 leading, due to interaction of the positively charged N-terminal fragment of tBid with anionic lipid headgroups, to Bax reorientation such that a hydrogen-bonded pair of residues, Asp98 and Ser184, is brought into close proximity with negatively charged lipid headgroups. The interaction with these headgroups destabilizes the hydrogen bond which results in the release of helix α9 from the Bax-binding groove, its insertion into the membrane, followed by insertion into the membrane of the α5–α6 helical hairpin. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Predicting membrane protein type by functional domain composition and pseudo-amino acid composition

Given the sequence of a protein, how can we predict whether it is a membrane protein or non-membrane protein? If it is, what membrane protein type it belongs to? Since these questions are closely relevant to the function of an uncharacterized protein, their importance is self-evident. Particularly, with the explosion of protein sequences entering into databanks and the relatively much slower progress in using biochemical experiments to determine their functions, it is highly desired to develop an automated method that can be used to give a fast answers to these questions. By hybridizing the functional domain (FunD) and pseudo-amino acid composition (PseAA), a new strategy called FunD-PseAA predictor was introduced. To test the power of the predictor, a highly non-homologous data set was constructed where none of proteins has 25% sequence identity to any other. The overall success rates obtained with the FunD-PseAA predictor on such a data set by the jackknife cross-validation test was 85% for the case in identifying membrane protein and non-membrane protein, and 91% in identifying the membrane protein type among the following 5 categories: (1) type-1 membrane protein, (2) type-2 membrane protein, (3) multipass transmembrane protein, (4) lipid chain-anchored membrane protein, and (5) GPI-anchored membrane protein. These rates are much higher than those obtained by the other methods on the same stringent data set, indicating that the FunD-PseAA predictor may become a useful high throughput tool in bioinformatics and proteomics.     

Melting of individual lipid components in binary lipid mixtures studied by FTIR spectroscopy, DSC and Monte Carlo simulations

Monte Carlo (MC) simulations, Differential Scanning Calorimetry (DSC) and Fourier Transform InfraRed (FTIR) spectroscopy were used to study the melting behavior of individual lipid components in two-component membranes made of DMPC and DSPC. We employed Monte Carlo simulations based on parameters obtained from DSC profiles to simulate the melting of the different lipids as a function of temperature. The simulations show good agreement with the FTIR data recorded for deuterated and non-deuterated lipids, which demonstrates that the information on the differential melting of the individual components is already contained in the calorimetric profiles. In mixtures, both lipids melt over a wide temperature range. As expected, the lipid melting events of the lipid with the lower melting temperature occur on average at lower temperatures. The simulations also yield information on the lateral distribution of the lipids that is neither directly contained in the DSC nor in the FTIR data. In the phase coexistence region, liquid disordered domains are typically richer in the lower-melting-temperature lipid species.     

Monte Carlo simulations of a protein molecule with and without hydration energy calculated by the hydration-shell model

Monte Carlo simulations of a small protein, carmbin, were carried out with and without hydration energy. The methodology presented here is characterized, as compared with the other similar simulations of proteins in solution, by two points: (1) protein conformations are treated in fixed geometry so that dihedral angles are independent variables rather than cartesian coordinates of atoms; and (2) instead of treating water molecules explicitly in the calculation, hydration energy is incorporated in the conformational energy function in the form of g _i A _i, whereA _i is the accessible surface area of an atomic groupi in a given conformation, andg _i is the free energy of hydration per unit surface area of the atomic group (i.e., hydration-shell model). Reality of this model was tested by carrying out Monte Carlo simulations for the two kinds of starting conformations, native and unfolded ones, and in the two kinds of systems,in vacuo and solution. In the simulations starting from the native conformation, the differences between the mean propertiesin vacuo and solution simulations are not very large, but their fluctuations around the mean conformation during the simulation are relatively smaller in solution thanin vacuo. On the other hand, in the simulations starting from the unfolded conformation, the molecule fluctuates much more largely in solution thanin vacuo, and the effects of taking into account the hydration energy are pronounced very much. The results suggest that the method presented in this paper is useful for the simulations of proteins in solution.     

Inversion of terrestrial ecosystem model parameter values against eddy covariance measurements by Monte Carlo sampling

Effective measures to counter the rising levels of carbon dioxide in the Earth's atmosphere require that we better understand the functioning of the global carbon cycle. Uncertainties about, in particular, the terrestrial carbon cycle's response to climate change remain high. We use a well‐known stochastic inversion technique originally developed in nuclear physics, the Metropolis algorithm, to determine the full probability density functions (PDFs) of parameters of a terrestrial ecosystem model. By thus assimilating half‐hourly eddy covariance measurements of CO₂ and water fluxes, we can substantially reduce the uncertainty of approximately five model parameters, depending on prior uncertainties. Further analysis of the posterior PDF shows that almost all parameters are nearly Gaussian distributed, and reveals some distinct groups of parameters that are constrained together. We show that after assimilating only 7 days of measurements, uncertainties for net carbon uptake over 2 years for the forest site can be substantially reduced, with the median estimate in excellent agreement with measurements.     

Numerical and Monte Carlo simulations of phenolic polymerizations catalyzed by peroxidase

Numerical and Monte Carlo simulations of horseradish peroxidase-catalyzed phenolic polymerizations have been performed. Kinetic constants for the simulations were fit to data from the oxidation and polymerization of bisphenol A. Simulations of peroxidase-catalyzed phenolic polymerization were run as a function of enzyme concentration and radical transfer and radical coupling rate constants. Predictions were performed with respect to conversion vs. time and number average molecular weight and polydispersity vs. conversion. It is shown that the enzymatic polymerization of phenols can be optimized with respect to high molecular weights by employing low enzyme concentrations and phenols with low radical coupling rate constants coupled with relatively high radical transfer rate constants. Such phenols may be identified by using linear free energy relationships that relate radical reactivity to electron donating/withdrawing potential of the phenolic substituent. (c) 1993 John Wiley & Sons, Inc.     

A New Monte Carlo Method for Direct Calculation of the Critical Size and the Formation Work of a Microdrop

New Monte Carlo procedures in open ensembles are proposed. Non-stationary Markov chain procedure in the μl;pT - ensemble provides a direct estimation for the critical size of a condensation nucleus at given p and T. A stationary procedure in the μlpT ensemble with two allowed particle numbers n and n + 1 provides the direct way to calculate the chemical potential and Gibbs free energy of a cluster; in the grand canonical (μlVT) ensemble the same approach gives μl and the Helmholtz free energy. The same procedures are readily applicable to periodic systems representing bulk phases.     

Conformational analysis and dynamics of mannobiosides and mannotriosides using Monte Carlo/stochastic dynamics simulations

The conformation and dynamics of alpha-(1-->2)-mannobioside, alpha-(1-->6)-mannobioside, and of the trisaccharide alpha-Man-(1-->2)-alpha-Man-(1-->6)-alpha- Man-OMe were studied using Monte Carlo/stochastic dynamics (MC/SD) simulations, the AMBER* force field, and the GB/SA implicit water solvation model. The results are in agreement with available experimental data.     

Size and shape of bovine interphotoreceptor retinoid-binding protein by electron microscopy and hydrodynamic analysis

Individual molecules of interphotoreceptor retinoid-binding protein (IRBP), a protein likely to be important in the visual cycle, were visualized by means of electron microscopy. IRBP was coated with a very thin layer of tungsten and photographed by dark-field imaging. IRBP is seen to be a flexible, elongated molecule about 24 nm in length by 3-4 nm in width (statistical modes). These dimensions agree very well with those calculated from the frictional ratio obtained from sedimentation data. Approximately half of these rod-shaped IRBP molecules are straight, and half are bent in the middle, usually with an angle of 60-90 degrees between the two arms. A representation of IRBP as a bendable string of beads yields calculations of dimensions and of hydrodynamic parameters consistent with the electron microscopic and sedimentation data; the sedimentation coefficients derived from this representation are nearly insensitive to molecular bending. When IRBP is bound to saturating amounts of its endogenous ligands, all-trans- or 11-cis-retinol, its sedimentation behavior is unchanged, and the same types of particles are visualized by electron microscopy as with the free protein; however, a greater proportion of the molecules are bent. Deglycosylation of IRBP (with peptide:N-glycosidase F) results in a somewhat smaller molecule that retains its rod-like shape, as shown by gel filtration and sedimentation data. The results indicate that IRBP is an elongated molecule and suggest that a structural change may occur upon ligand binding.     

Mean residence time of molecules diffusing in a cell bounded by a semi-permeable membrane: Monte Carlo simulations and an expression relating membrane transition probability to permeability

 The rapid exchange of water across erythrocyte membranes is readily measured using an NMR method that entails doping a suspension of cells with a moderately high concentration of Mn²⁺ and measuring the rate of transverse relaxation of the nuclear magnetisation. Analysis of the data yields an estimate of the rate constant for membrane transport, from which the membrane permeability can be determined. It is assumed in the analysis that the efflux rate of the water is solely a function of the rate of membrane permeation and that the time it takes for intracellular water molecules to diffuse to the membrane is relatively insignificant. The limits of this assumption were explored by using random-walk simulations of diffusion in cells modelled as parallel planes, spheres, and biconcave discs. The rate of membrane transport was specified in terms of a transition probability but it was not initially clear what the relationship should be between this parameter and the diffusional membrane permeability P _d. This relationship was derived and used to show that the mean residence time for a water molecule is determined by P _d when the diffusion coefficient is above a certain threshold value; it is determined by the distance to the membrane below that value. Received: 7 January 2000 / Revised version: 4 April 2000 / Accepted: 4 April 2000     

Role of membrane sialic acid and glycophorin protein in thorium induced aggregation and hemolysis of human erythrocytes

Thorium-232 (²³²Th), a natural radionuclide from the actinide family, is abundantly present in monazite and other ores. It is used as one of the prime fuel materials in nuclear industry and may pose an exposure risk to nuclear workers and members of the public. Human erythrocytes, as a classical cellular membrane model, were coincubated with ²³²Th in order to elucidate whether this naturally occurring important radionuclide produced perturbations to cell membrane. Present study revealed that erythrocytes underwent aggregation or lysis depending on the ratio of ²³²Th to cell. Scanning electron micrographs showed that erythrocytes transformed into equinocytes and/or spherocytes after ²³²Th treatment. Further examination of erythrocyte by atomic force microscopy suggested significant increase in surface roughness after ²³²Th treatment. Experiments on neuraminidase treated and/or anti-GpA antibody blocked erythrocytes suggested significant role of membrane sialic acid and glycophorin A (GpA) protein in aggregation or hemolytic effects of ²³²Th. Further results showed that ²³²Th caused hemolysis by colloid osmotic mechanism, as evidenced by potassium efflux, osmotic protection and osmotic fragility studies. Osmoprotection experiments indicated that hemolysis get elicited through the formation of membrane pores of ∼2.0 nm in size. Hemolysis studies in presence of inhibitors (TEA, bumetanide, DIDS and amiloride) revealed the role of K⁺ channel, Na⁺/K⁺/2Cl⁻ channel, Cl⁻/HCO₃⁻ anion exchanger and Na⁺/H⁺ antiporter in ²³²Th induced erythrolysis. Presence of non-diffusible cation (N-methyl d-glucasamine) or anion (gluconate) in erythrocyte suspending medium further confirm the role of Na⁺ and Cl⁻ influx in hemolytic effect of ²³²Th. These findings provide significant insight in structural, biochemical and osmotic toxic effects of ²³²Th on human erythrocytes.     

Measurement of the glucose permeation rate across phospholipid bilayers using small unilamellar vesicles Effect of membrane composition and temperature

Small unilamellar vesicles were used to measure the permeability of saturated phosphatidylcholine bilayers to glucose. The presented method circumvents most of the common restrictions of classical permeability experiments. Increasing the fatty acid chain length of the lipids reduced the permeation rate significantly. Raising the temperature above that of the lipid phase transition drastically increased membrane permeability. Arrhenius plots demonstrated the activation energy to be independent of membrane composition and the phase-state of the lipids. The permeation process is discussed in terms of a constant energy to disrupt all hydrogen bonds between permeant and aqueous solvent prior to penetrating the membrane. The magnitude of the permeability coefficient is partly determined by a unfavourable change in entropy of activation on crossing the water/lipid interface. All results indicate that the penetration of the dehydrated permeant into the hydrophobic barrier is the rate-limiting step in the permeation of glucose.     

FRET analysis of domain formation and properties in complex membrane systems

The application of Förster Resonance Energy Transfer (FRET) to the detection and characterization of phase separation in lipid bilayers (both in model systems and in cell membranes) is reviewed. Models describing the rate and efficiency of FRET for both uniform probe distribution and phase separation, and recently reported methods for detection of membrane heterogeneity and determination of phase boundaries, probe partition coefficients and domain size, are presented and critically discussed. Selected recent applications of FRET to one-phase lipid systems, gel/fluid phase separation, liquid ordered/liquid disordered phase separation (lipid rafts), complex systems containing ceramide and cell membranes are presented to illustrate the wealth of information that can be inferred from carefully designed FRET studies of membrane domains.     

Micro-phase separation and configuration of ABC triblock copolymer in ultra-thin film by Monte Carlo simulation

A Monte Carlo simulation using the bond fluctuation and cavity diffusion algorithms was adopted to investigate the micro-phase separation of ABC triblock copolymer in ultra-thin film on simple cubic lattice. Simulations reveal that the morphologies of ABC copolymer films are dependent on not only the volume fraction of the middle block B (f _B) but also on the ratio of interaction between different kinds of blocks (?_(AC)/?_(AB)). As for the molecular orientation, the copolymers stretch parallel to the flat surface at lower f _B, but tend to align perpendicularly along z direction at higher f _B. Furthermore, the chain configuration was discussed in detail. Smaller ?_(AC)/?_(AB) is beneficial to the formation of a “loop” configuration, whereas, larger ?_(AC)/?_(AB) would result in a “bridge” configuration of ABC triblock copolymer chains. The formation of micro-phase structures was illustrated intuitively by the molecular orientation and the chain configuration.     

Structural analysis of membrane protein complexes by single particle electron microscopy

Single particle electron microscopy (EM) is an increasingly important tool for the structural analysis of macromolecular complexes. The main advantage of the technique over other methods is that it is not necessary to precede the analysis with the growth of crystals of the sample. This advantage is particularly important for membrane proteins and large protein complexes where generating crystals is often the main barrier to structure determination. Therefore, single particle EM can be employed with great utility in the study of large membrane protein complexes. Although the construction of atomic resolution models by single particle EM is possible in theory, currently the highest resolution maps are still limited to approximately 7-10A resolution and 15-30 A resolution is more typical. However, by combining single particle EM maps with high-resolution models of subunits or subcomplexes from X-ray crystallography and NMR spectroscopy it is possible to build up an atomic model of a macromolecular assembly. Image analysis procedures are almost identical for micrographs of soluble protein complexes and detergent solubilized membrane protein complexes. However, electron microscopists attempting to prepare specimens of a membrane protein complex for imaging may find that these complexes require different handling than soluble protein complexes. This paper seeks to explain how high-quality specimen grids of membrane protein complexes may be prepared to allow for the determination of their structure by EM and image analysis.