Distance matrix-based approach to protein structure prediction

Much structural information is encoded in the internal distances; a distance matrix-based approach can be used to predict protein structure and dynamics, and for structural refinement. Our approach is based on the square distance matrix D = [r _ij ² ] containing all square distances between residues in proteins. This distance matrix contains more information than the contact matrix C, that has elements of either 0 or 1 depending on whether the distance r _ij is greater or less than a cutoff value r _cutoff. We have performed spectral decomposition of the distance matrices $ {\mathbf{D}} = \sum {\lambda_{k} {\mathbf{v}}_{k} {\mathbf{v}}_{k}^{T} } Much structural information is encoded in the internal distances; a distance matrix-based approach can be used to predict protein structure and dynamics, and for structural refinement. Our approach is based on the square distance matrix D = [r _ij²] containing all square distances between residues in proteins. This distance matrix contains more information than the contact matrix C, that has elements of either 0 or 1 depending on whether the distance r _ij is greater or less than a cutoff value r _cutoff. We have performed spectral decomposition of the distance matrices , in terms of eigenvalues and the corresponding eigenvectors and found that it contains at most five nonzero terms. A dominant eigenvector is proportional to r ²—the square distance of points from the center of mass, with the next three being the principal components of the system of points. By predicting r ² from the sequence we can approximate a distance matrix of a protein with an expected RMSD value of about 7.3 ?, and by combining it with the prediction of the first principal component we can improve this approximation to 4.0 ?. We can also explain the role of hydrophobic interactions for the protein structure, because r is highly correlated with the hydrophobic profile of the sequence. Moreover, r is highly correlated with several sequence profiles which are useful in protein structure prediction, such as contact number, the residue-wise contact order (RWCO) or mean square fluctuations (i.e. crystallographic temperature factors). We have also shown that the next three components are related to spatial directionality of the secondary structure elements, and they may be also predicted from the sequence, improving overall structure prediction. We have also shown that the large number of available HIV-1 protease structures provides a remarkable sampling of conformations, which can be viewed as direct structural information about the dynamics. After structure matching, we apply principal component analysis (PCA) to obtain the important apparent motions for both bound and unbound structures. There are significant similarities between the first few key motions and the first few low-frequency normal modes calculated from a static representative structure with an elastic network model (ENM) that is based on the contact matrix C (related to D), strongly suggesting that the variations among the observed structures and the corresponding conformational changes are facilitated by the low-frequency, global motions intrinsic to the structure. Similarities are also found when the approach is applied to an NMR ensemble, as well as to atomic molecular dynamics (MD) trajectories. Thus, a sufficiently large number of experimental structures can directly provide important information about protein dynamics, but ENM can also provide a similar sampling of conformations. Finally, we use distance constraints from databases of known protein structures for structure refinement. We use the distributions of distances of various types in known protein structures to obtain the most probable ranges or the mean-force potentials for the distances. We then impose these constraints on structures to be refined or include the mean-force potentials directly in the energy minimization so that more plausible structural models can be built. This approach has been successfully used by us in 2006 in the CASPR structure refinement ().     

优化的绿色荧光蛋白在超嗜热嗜酸古菌冰岛硫化叶菌中的表达与应用

【目的】开发可用于在极端嗜热嗜酸模式泉古菌冰岛硫化叶菌(Sulfolobus islandicus)中进行高效表达的eCGP123(enhanced consensus green protein variant 123)荧光蛋白,并用作S.islandicus的细胞内蛋白定位工具。【方法】绿色荧光蛋白突变体eCGP123具有极高的热稳定性、耐酸性和可逆的荧光特性等。本研究主要对eCGP123的基因根据S.islandicus密码子偏好性进行优化与合成,在大肠杆菌(Escherichia coli)中表达并研究其蛋白性质;通过在eCGP123的C末端分别融合具有不同细胞内定位的蛋白(包括E.coli来源的Fts Z和S.islandicus来源的Ups E、PCNA1和SiRe_1200等),构建eCGP123及其融合蛋白的表达菌株,用激光共聚焦显微镜分析eCGP123及其融合蛋白在E.coli和S.islandicus活细胞中的亚细胞定位。【结果】我们确认了在E.coli中表达并纯化密码子优化后的e CGP123具有与野生型绿色荧光蛋白相同的吸光值和较高的热稳定性。细胞学分析显示细胞分裂相关蛋白FtsZ和Si Re_1200分别主要定位于E.coli和S.islandicus分裂细胞的中间;鞭毛组分蛋白Ups E呈点状均匀分布,可能定位于细胞膜上;DNA复制滑动夹亚基PCNA1呈区域性点状分布,暗示了DNA复制区域的位置。蛋白的亚细胞定位与预期结果基本吻合。【结论】绿色荧光蛋白e CGP123可以作为报告蛋白,应用于S.islandicus细胞的蛋白定位分析中,可作为该模式菌株中功能基因研究的重要工具,但需要进一步优化条件。     

The new understanding of Arabidopsis thaliana proteins associated with salinity

Abstract

Salt stress is a major abiotic stress limiting the productivity and the geographical distribution of many plant species. Arabidopsis thaliana is an excellent model with rich genetic resources for modern plant biology research. To comprehensively and representatively understand salt-response mechanisms in A. thaliana, we applied the first attempt to use the most data (252 of 10,469 reviewed A. thaliana protein) from public protein database for displaying the enriched protein domains, Kyoto Encyclopedia of Genes and Genomes pathways, molecular functions, and cell localizations involved in salt-response. The data were analyzed by Database for Annotation Visualization and Integrated Discovery. Our results indicated salt-response proteins cross-talked not only with drought and temperature stress as previously reported but also with further stresses such as bacterium, light, metal ion, radiation, and wounding stress. Multiple cellular localizations under salt stress indicated proteins were versatile. In addition, 27 proteins have the characteristics with response to multiple stresses and localization in multiple places. We called it the ‘space-stress’ double cross-talk effects, which indicated that A. thaliana proteins dealt with salt stress and other stresses in a reciprocal economical way. An enriched bioinformatics analysis of the large data could provide clues and basis for the development of salt-response potential biomarkers for plant growth and crop productivity.     

Theoretical Model of the Three-dimensional Structure of a Disease Resistance Gene Homolog Encoding Resistance Protein in Vigna mungo

Abstract

Plant disease resistance (R) genes, the key players of innate immunity system in plants encode ‘R’ proteins. ‘R’ protein recognizes product of avirulance gene from the pathogen and activate downstream signaling responses leading to disease resistance. No three dimensional (3D) structural information of any ‘R’ proteins is available as yet. We have reported a ‘R’ gene homolog, the 'VMYR1′, encoding ‘R’ protein in Vigna mungo. Here, we describe the homology modeling of the 'VMYR1′ protein. The model was created by using the 3D structure of an ATP-binding cassette transporter protein from Vibrio cholerae as a template. The strategy for homology modeling was based on the high structural conservation in the superfamily of P-loop containing nucleoside triphosphate hydrolase in which target and template proteins belong. This is the first report of theoretical model structure of any ‘R’ proteins.     

Distance restraints from crosslinking mass spectrometry: Mining a molecular dynamics simulation database to evaluate lysine–lysine distances

Integrative structural biology attempts to model the structures of protein complexes that are challenging or intractable by classical structural methods (due to size, dynamics, or heterogeneity) by combining computational structural modeling with data from experimental methods. One such experimental method is chemical crosslinking mass spectrometry (XL‐MS), in which protein complexes are crosslinked and characterized using liquid chromatography‐mass spectrometry to pinpoint specific amino acid residues in close structural proximity. The commonly used lysine‐reactive N‐hydroxysuccinimide ester reagents disuccinimidylsuberate (DSS) and bis(sulfosuccinimidyl)suberate (BS³) have a linker arm that is 11.4 Å long when fully extended, allowing C_α (alpha carbon of protein backbone) atoms of crosslinked lysine residues to be up to ～24 Å apart. However, XL‐MS studies on proteins of known structure frequently report crosslinks that exceed this distance. Typically, a tolerance of ～3 Å is added to the theoretical maximum to account for this observation, with limited justification for the chosen value. We used the Dynameomics database, a repository of high‐quality molecular dynamics simulations of 807 proteins representative of diverse protein folds, to investigate the relationship between lysine–lysine distances in experimental starting structures and in simulation ensembles. We conclude that for DSS/BS³, a distance constraint of 26–30 Å between C_α atoms is appropriate. This analysis provides a theoretical basis for the widespread practice of adding a tolerance to the crosslinker length when comparing XL‐MS results to structures or in modeling. We also discuss the comparison of XL‐MS results to MD simulations and known structures as a means to test and validate experimental XL‐MS methods.     

Discrete analyses of protein dynamics

Abstract

Protein structures are highly dynamic macromolecules. This dynamics is often analysed through experimental and/or computational methods only for an isolated or a limited number of proteins. Here, we explore large-scale protein dynamics simulation to observe dynamics of local protein conformations using different perspectives. We analysed molecular dynamics to investigate protein flexibility locally, using classical approaches such as RMSf, solvent accessibility, but also innovative approaches such as local entropy. First, we focussed on classical secondary structures and analysed specifically how β-strand, β–turns, and bends evolve during molecular simulations. We underlined interesting specific bias between β–turns and bends, which are considered as the same category, while their dynamics show differences. Second, we used a structural alphabet that is able to approximate every part of the protein structures conformations, namely protein blocks (PBs) to analyse (i) how each initial local protein conformations evolve during dynamics and (ii) if some exchange can exist among these PBs. Interestingly, the results are largely complex than simple regular/rigid and coil/flexible exchange. Abbreviations N_eq number of equivalent

PB Protein Blocks

PDB Protein DataBank

RMSf root mean square fluctuations

Communicated by Ramaswamy H. Sarma     

Predicting the Structure of the Flavodoxin from Eschericia coli by Homology Modeling,Distance Geometry and Molecular Dynamics

Abstract

As part of our on-going development of a method, based upon distance geometry calculations, for predicting the structures of proteins from the known structures of their homologues, we have predicted the structure of the 176 residue Flavodoxin from Escherichia coli. This prediction was based upon the crystal structures of the homologous Flavodoxins from Anacystis nidulans, Chondrus crispus, Desulfovibrio vulgaris and Clostridium beijerinckii, whose sequence identities with Escherichia coli were 44%, 33%, 23% and 16%, respectively. A total of 13,043 distance constraints among the alpha-carbons of the Escherichia coli structure were derived from the sequence alignments with the known structures, together with 8,893 distance constraints among backbone and sidechain atoms of adjacent residues, 978 between the alpha-carbons and selected atoms of the flavin mononucleotide cofactor, 116 constraints to enforce conserved hydrogen bonds, and 452 constraints on the torsion angles in conserved residues. An ensemble of ten random Escherichia coli structures was computed from these constraints, with an average root mean square coordinate deviation (RMSD) among the alpha carbons of 0.85 Ångstroms (excluding the first 1 and last 6 residues, which have no corresponding residues in any of the homologues and hence were unconstrained); the corresponding average heavy-atom RMSD was 1.60 Å.

Since the distance geometry calculations were performed without hydrogen atoms, protons were added to the resulting structures and these structures embedded in a 50 × 50 × 40 Å solvent box with periodic boundary conditions. They were then subjected to a 20 picosecond dynamical simulated annealing procedure, starting at 300 K and gradually reduced to 10K, in which all the distance and torsion angle constraints were maintained by means of harmonic restraint functions. This was followed up by 1000 iterations of unrestrained conjugate gradients minimization. The goal of this energy refinement procedure was not to drastically modify the structures in an attempt at a priori prediction, but merely to improve upon the predictions obtained from the geometric constraints, particularly with regard to their local backbone and sidechain conformations and their hydrogen bonds. The resulting structures differed from the respective starting structures by an average of 1.52 Å in their heavy atom RMSD's, while the average RMSD among the heavy atoms of residues 2-170 increased slightly to 1.66 Å. We hope these structures will be good enough to enable the phase problem to be solved for the crystallographic data that is now being collected on this protein.     

The uniqueness of protein sequences.o-Uniqueness and infrequent peptides

Concepts of the uniqueness of the amino acid sequences of proteins were defined in a prior report (Saroff, H. A. and F. A. Kutyna. 1981. “The Uniqueness of Protein Sequences: A Monte Carlo Analysis.”Bull. math. Biol. 43, 619–639), which presented a detailed discussion ofi-uniqueness, i.e. the tendency of small peptides to be repeated within an amino acid sequence of a protein. We now report on the quantitative analysis ofo-uniqueness, which evaluates the tendency of small peptides to be repeated amongst different proteins, usually of a single species. A detailed analysis of theo-uniqueness of several proteins is presented to illustrate the method and the range of values encountered. Uniqueness data on sequences of human proteins in a data bank of sequences containing about 32,500 amino acids are made available in the form of a microfiche. Analysis of biologically active subsequences such as the angiotensins and the enkephalins suggest a tendency of the subsequences contributing to the property ofo-uniqueness to cluster in portions of the parent protein sequence which are biologically active. This property may provide a general method for predicting biologically active areas of proteins. Current data may already be adequate to permit useful predictions, and the rapidly accumulating and interrelated new data on nucleic acid and protein sequences will further enhance the power ofo-uniqueness analysis.     

Structural homology of spinach acyl carrier protein and Escherichia coli acyl carrier protein based on NMR data

Acyl carrier proteins (ACPs) from spinach and from Escherichia coli have been used to demonstrate the utility of proton NMR for comparison of homologous structures. The structure of E. coli ACP had been previously determined and modeled as a rapid equilibrium among multiple conformational forms (Kim and Prestegard, Biochemistry 28:8792–8797, 1989). Spinach ACP showed two slowly exchanging forms and could be manipulated into one form for structural study. Here we compare this single form to postulated multiple forms of E. coli ACP using the limited amount of NOE data available for the spinach protein. A number of long-range NOE contacts were present between homologous residues in both spinach and E. coli ACP, suggesting tertiary structural homology. To allow a more definitive structural comparison, a method was developed to use spinach ACP NOE constraints to search for regions of structural divergence from two postulated forms of E. coli ACP. The homologous regions of the two protein sequences were aligned, additional distance constraints were extracted from the E. coli structure, and these were mapped onto the spinach sequence. These distance constraints were combined with experimental NOE constraints and a distance geometry simulated annealing protocol was used to test for compatibility of the constraints. All of the experimental spinach NOE constraints could be successfully combined with the E. coli data, confirming the general hypothesis of structural homology. A better fit was obtained with one form, suggesting a preferential stabilization of that form in the spinach case. Proteins 27:131–143 © 1997 Wiley-Liss, Inc.     

Intrinsic nature of the three-dimensional structure of proteins as determined by distance geometry with good sampling properties

Summary A protocol for distance geometry calculation is shown to have excellent sampling properties in the determination of three-dimensional structures of proteins from nuclear magnetic resonance (NMR) data. This protocol uses a simulated annealing optimization employing mass-weighted molecular dynamics in four-dimensional space (Havel, T.F. (1991) Prog. Biophys. Mol. Biol., 56, 43–78). It attains an extremely large radius of convergence, allowing a random coil conformation to be used as the initial estimate for the succeeding optimization process. Computations are performed with four systems of simulated distance data as tests of the protocol, using an unconstrained l-alanine 30mer and three different types of proteins, bovine pancreatic trypsin inhibitor, the -amylase inhibitor Tendamistat, and the N-terminal domain of the 434-repressor. The test of the unconstrained polypeptide confirms that the sampled conformational space is that of the statistical random coil. In the larger and more complicated systems of the three proteins, the protocol gives complete convergence of the optimization without any trace of initial structure dependence. As a result of an exhaustive conformational sampling by the protocol, the intrinsic nature of the structures generated with distance restraints derived from NMR data has been revealed. When the sampled structures are compared with the corresponding X-ray structures, we find that the averages of the sampled structures always show a certain pattern of discrepancy from the X-ray structure. This discrepancy is due to the short distance nature of the distance restraints, and correlates with the characteristic shape of the protein molecule.Abbreviations r.m.s.d. root-mean-square deviation - MD molecular dynamics - NMR nuclear magnetic resonance - NOE nuclear Overhauser enhancement - BPTI bovine pancreatic trypsin inhibitor     

蛋白质结构确定领域中的几种矩阵填充算法的对比评估

目的 目前,如何从核磁共振(nuclear magnetic resonance,NMR)光谱实验中准确地确定蛋白质的三维结构是生物物理学中的一个热门课题,因为蛋白质是生物体的重要组成成分,了解蛋白质的空间结构对研究其功能至关重要,然而由于实验数据的严重缺乏使其成为一个很大的挑战。方法 在本文中,通过恢复距离矩阵的矩阵填充(matrix completion,MC)算法来解决蛋白质结构确定问题。首先,初始距离矩阵模型被建立,由于实验数据的缺乏,此时的初始距离矩阵为不完整矩阵,随后通过MC算法恢复初始距离矩阵的缺失数据,从而获得整个蛋白质三维结构。为了进一步测试算法的性能,本文选取了4种不同拓扑结构的蛋白质和6种现有的MC算法进行了测试,探究了算法在不同的采样率以及不同程度噪声的情况下算法的恢复效果。结果 通过分析均方根偏差(root-mean-square deviation,RMSD)和计算时间这两个重要指标的平均值及标准差评估了算法的性能,结果显示当采样率和噪声因子控制在一定范围内时,RMSD值和标准差都能达到很小的值。另外本文更加具体地比较了不同算法的特点和优势,在精确采样情况下...     

Tandem photoaffinity labeling of a target protein using a linker with biotin,alkyne and benzophenone groups and a bioactive small molecule with an azide group

A novel linker containing biotin, alkyne and benzophenone groups (1) was synthesized to identify target proteins using a small molecule probe. This small molecule probe contains an azide group (azide probe) that reacts with an alkyne in 1 via an azide–alkyne Huisgen cycloaddition. Cross-linking of benzophenone to the target protein formed a covalently bound complex consisting of the azide probe and the target protein via 1. The biotin was utilized via biotin–avidin binding to identify the cross-linked complex. To evaluate the effectiveness of 1, it was applied in a model system using an allene oxide synthase (AOS) from the model moss Physcomitrella patens (PpAOS1) and an AOS inhibitor that contained azide group (3). The cross-linked complex consisting of PpAOS1, 1 and 3 was resolved via SDS–PAGE and visualized using a chemiluminescent system. The method that was developed in this study enables the effective identification of target proteins.     

A computer-based protocol for semiautomated assignments and 3D structure determination of proteins

Summary A strategy is presented for the semiautomated assignment and 3D structure determination of proteins from heteronuclear multidimensional Nuclear Magnetic Resonance (NMR) data. This approach involves the computer-based assignment of the NMR signals, identification of distance restraints from nuclear Overhauser effects, and generation of 3D structures by using the NMR-derived restraints. The protocol is described in detail and illustrated on a resonance assignment and structure determination of the FK506 binding protein (FKBP, 107 amino acids) complexed to the immunosuppressant, ascomycin. The 3D structures produced from this automated protocol attained backbone and heavy atom rmsd of 1.17 and 1.69 Å, respectively. Although more highly resolved structures of the complex have been obtained by standard interpretation of NMR data (Meadows et al. (1993) Biochemistry, 32, 754–765), the structures generated with this automated protocol required minimal manual intervention during the spectral assignment and 3D structure calculations stages. Thus, the protocol may yield an approximate order of magnitude reduction in the time required for the generation of 3D structures of proteins from NMR data.     

沈阳地区北虫草野生与市售菌株人工培育子实体的蛋白质组学比较

【背景】北虫草作为冬虫夏草的代用品,具有与冬虫夏草类似的药理活性,其富含的蛋白质和氨基酸通常作为衡量真菌营养价值的重要指标,从中分离纯化具有潜在临床应用价值的蛋白质或多肽,已成为一个研究热点。【目的】检测沈阳北虫草野生与市售菌株人工培育子实体的蛋白质组成,分析相同培育条件下获得的蛋白种类、数量及其功能的差异,为深入研究鉴定沈阳地区北虫草药用蛋白和针对性驯化提供了蛋白质组学数据基础。【方法】采集沈阳棋盘山野生北虫草菌株,与市售人工栽培北虫草菌株同期分别经组织分离、液体发酵后培育获得子实体,通过蛋白提取、胰酶酶解后,采用非标定量技术液相色谱-质谱联用方法,对野生和市售来源培育的子实体样本进行定量蛋白组的研究。【结果】共鉴定到9 233条特异性肽段和1 923个蛋白,其中含有1 163个可定量蛋白,野生来源培育子实体有214个蛋白表达发生上调,181个蛋白表达发生下调,对这些差异蛋白进行功能富集分析发现,其主要参与能量生产/转换、氨基酸转运/代谢、抗氧化功能。在相同的营养摄取条件下,野生来源培育菌种在各个能量代谢、氨基酸代谢功能中的相关蛋白表达量高于市售来源培育的菌种。野生来源培育菌种的一种抗氧化重要蛋白(Gene Name:ISF_02112)表达量远远高于(Fold Change9)市售来源培育菌种。同时与抗氧化和代谢功能相关的差异蛋白有22个。【结论】沈阳地区北虫草野生菌株经适当人工培育会保留部分优良的生物学特性,2种来源菌株培育的子实体具有丰富及优异抗氧化功能的蛋白,子实体蛋白的抗氧化能力与其整体代谢能力相关。本研究结果为深入研究鉴定北虫草药用蛋白和针对性驯化提供蛋白质组学数据基础。     

An integrated in silico approach to understand protein–protein interactions: human meprin-β with fetuin-A

Abstract

Human meprin-β, a zinc metalloprotease belonging to the astacin family, have been found to be associated with many pathological conditions like inflammatory bowel disease, fibrosis and neurodegenerative disease. The inhibition of meprin-β by various inhibitors, both macromolecular and small molecules, is crucial in the control of several diseases. Human fetuin-A, a negative acute phase protein involved in inflammatory disease, has recently been identified as an endogenous inhibitor for meprin-β. In this computational study, an integrated in silico approach was performed using existing structural information of meprin-β coupled with ab initio modelling of human fetuin-A to predict a rational model of the complex through protein–protein docking. Further, the models were optimized and validated to generate an ensemble of conformations through extensive molecular dynamics simulation. Virtual alanine scanning mutagenesis was explored to identify hotspot residues on both proteins significant for protein–protein interaction (PPI). The results of the study provide structural insight into PPI between meprin-β and fetuin-A which can be useful in designing molecules to modulate meprin-β activity.

Communicated by Ramaswamy H. Sarma     

A Study of Conformational Equilibria in Chain Molecules. I. Liquid n-Butane

Abstract

We have performed molecular dynamics simulations for liquid n-butane in order to understand liquid structures in terms of both inter- and intra-molecular interactions. Each n-butane molecule consists of four sites interacting with LJ potential and only a dihedral angle is taken into account as the internal degree of freedom. The population of gauche conformations with respect to the ideal gas state is found to increase in the liquid state. To investigate how the intermolecular interaction affects the dihedral angle distribution, we also adopt the repulsive LJ potential (RLJ) model. It is found that the nearest neighbor packing of the methyl and/or methylene groups can be approximately represented by using only the repulsive interaction. From the dihedral angle distribution, however, the rate of the shift of RLJ model to gauche is larger than that of LJ model and the attractive force also plays a significant role in the conformational equilibrium.     

利用蛋白质组学研究新德里金属b-内酰胺酶1对鲍曼不动杆菌的影响

【目的】比较临床分离的亲缘关系近的多药耐药鲍曼不动杆菌MDR-ZJ06(blaNDM-1–)和ABC3229(blaNDM-1+)的差异蛋白质组,以期发现新德里金属?-内酰胺酶1(New Delhimetallo-β-lactamase-1,NDM-1)对鲍曼不动杆菌生长代谢的影响。【方法】利用2-DE联合MALDI-TOF MS/MS技术鉴定差异表达蛋白,并在GO注析的基础上,对差异蛋白进行通路分析、功能分类和富集分析,并作出蛋白与蛋白相互作用网络。【结果】发现ABC3299相对于MDR-ZJ06有51个差异表达蛋白,其中11个蛋白表达上调,40个蛋白表达下调,并且这些差异蛋白主要涉及降低碳代谢、氨基酸代谢、脂肪酸代谢和细胞壁合成,增加铁离子转运系统形成。【结论】这个结果揭示了NDM-1可能是通过减缓细菌自身的代谢,增加自身铁的摄取使细菌机体系统地抵抗抗生素从而达到耐药。     

Effect of integral proteins upon bilayer permeability to small ions

We have modelled the effect of integral proteins upon the permeability of a lipid bilayer membrane to small ions. Our intention is to predict the temperature-and protein concentration-dependence of the relative ionic diffusion coefficient if it is assumed that the ions diffuse entirely through the lipid regions of the membrane, and do not cross the bilayer either by moving through the protein or along its hydrophobic surface. We used models of lipid-protein bilayers that have been used successfully before, together with a modification of the (protein-free) permeability mechanism of Kanchisa, Tsong, Cruzeiro-Hansson and Mouritsen (KTCM). We used Monte Carlo techniques to simulate the thermodynamics and permeability of the membrane because such a calculation, which depends upon spatial and time correlations, cannot be done using a mean field approximation. We took the protein to be analogous to bacteriorhodopsin and calculated the relative diffusion coefficients, R, of small ions through a DMPC bilayer containing such proteins. We found that although there is only a small change in R when T > T _m , R can change by an order of magnitude, depending upon the protein concentration, for T < T _m . We showed that regions of largest lipid area fluctuations and, hence, compressibility occur (i) at the interfaces between regions (which may be very small and not macroscopic phases) of extended-chain and excited-chain lipids, thus showing the similarity between this model and the KTCM model, and (ii) in regions where the temperature is close to T _m but the lipids are prevented from becoming ordered because of the proximity of proteins which try to keep them in their excited-chain states. We have plotted lipid fluctuation and ionic permeability maps to show which regions of the membrane display the largest permeability, and we show how these regions change as a function of temperature. Offprint requests to D. A. Pink