membrane proteins ( mps ) constitute about 30% of all the proteins encoded in the currently known genomes , and play critical roles in cell signaling , ion transport , and cell - cell communications , as well as assist the folding of other mps ( 1 ) . because of these biological significance
, mps represent the most important class of drug targets  about 50% of current molecular targets are membrane - bound ( 2 ) .
however , only about 2% ( 518 of 25,176 ) of the 3d structures deposited in the protein data bank ( pdb ; ref .
3 ) are for mps . and the number of high - resolution structures ( from x - ray diffraction and more recently from nmr ) remains even smaller , largely because of the difficulties in crystallizing mps .
recently , some new ideas and experimental approaches have been introduced in the area of mp crystallization ( 4 ) , all of which exploit the spontaneous self - assembling properties of lipids and detergent as vesicles ( vesicle - fusion method ) , discoidal micelles ( bicelle method ) , and liquid crystals or mesophases ( in meso or cubic - phase method ) . despite these promising new methods ,
the current gap between need and supply of mp 3d structures makes prediction algorithms important and essential .
mps come in a variety of sizes and shapes , though the available 3d structure principles are far less diverse than those of the globular proteins . from a structural point of view
one is the -helix bundle protein , in which one or several -helices span the membrane ; and the other is -barrel protein , in which eight or more antiparallel tm -strands form a closed barrel 5 . , 6 .. two recent examples 7 .
since jhnig and edholm in 1992 presented one of the first methods using secondary structure prediction to build suitable model structures as initial conformations for molecular dynamic studies ( 9 ) , several groups have tried computational approaches to elucidate mp structures . in 1993 , milik and skolnick presented a method based on the combination of a hydropathy scale for the prediction of trans - bilayer fragments with dynamic monte carlo simulation techniques ( 10 ) . in 1994 , taylor et al .
adapted some programs originally developed for the prediction of globular protein structures to derive a method for the prediction of integral mp structures ( 11 ) .
each step in the method is fully automated , from the initial sequence data bank searches to the final construction of 3d models .
consequently , estimates for prediction accuracy are perhaps overly optimistic . here , we summarize recent attempts within the field of computational biology and bioinformatics to predict an mp s structure .
most current methods of theoretical mp structure prediction do not actually deal with predicting the 3d structure , but rather try to predict the most likely topology of the protein , that is to say , the in / out location of the n and c termini relative to the membrane , and the number and position of transmembrane ( tm ) segments . a high - quality model of secondary structure and topology
is a prerequisite for experimental structure - function studies , and can be a starting point for attempts to model the 3d structure before molecular dynamics or simulated annealing simulations . in recent years , various accurate methods
because the number of high - resolution structures of -barrel proteins is less than that of the -helix proteins , the neural network has been more frequently adopted in the -strand topology prediction .
the details of some methods based on hidden markov models ( hmms ) are listed in table 2 .
many secondary structure prediction methods are based on statistical methods , physicochemical methods , sequence pattern maching , and evolutionary conservation ( 12 ) .
the main methods for identifing tm helices are on the basis of their hydrophobicity and known minimum length ( at least 15 residues ; ref .
membrane propensities were defined by a statistical analysis carried out on a set of 640 tm helices , belonging to 133 mps extracted from swiss - prot ( 14 ) that have experimentally defined topologies .
the five widely used prediction methods for predicting the topology of -helix bundle mps are tmhmm ( 15 ) , hmmtop ( 16 ) , memsat ( 17 ) , phdhtm ( 18 ) , and toppred ( 19 ) .
tmhmm , hmmtop , and memsat are all based on hmms with 5~7 types of structural states .
toppred was the first topology prediction method that combined hydrophobicity analysis and the positive - inside rule .
generally , these sequence - based methods for predicting the number and approximate location of tm helices within mps have about 85% accuracy . in 2003 , karin meln et al . tried to construct useful reliability scores for these methods ( 20 ) .
21 ) algorithm can provide a solution to the problem that non - transmembrane query sequences may give false positive hits ( 20%-30% ) in the prediction process .
the upgraded and modified version of the das - prediction method , das - tmfilter algorithm , has been distributed ( 22 ) .
the new algorithm is designed to make distinction between protein sequences with and without tm helices at a reasonably low rate of false positive prediction ( ~1 among 100 unrelated queries ) while the high efficiency of the original algorithm locating tm segments in queries is preserved ( sensitivity of ~95% among documented proteins with helical tm regions ) . in 2003 ,
xia and colleagues presented a new approach , conpred_elite ( 23 ) , that can predict the whole topology with accuracies of 98% for prokaryotic and 95% for eukaryotic proteins as they reported .
besides the tm helix , another tm segments type is -barrel , which consists of several tm strands . unlike -helical mps , there are no simple low - resolution experiments that yield large amounts of data for -barrel mps .
in fact , the overall hydrophobicity for -barrel mps is similar to that of soluble proteins ( 13 ) .
gromiha and colleagues combined amino acid preferences for -strands with the surrounding hydrophobicity of the respective residues to predict -strands ( 24 ) .
diederichs and colleagues proposed to use a neural network to predict the topology of the bacterial outer membrane -strand proteins and to locate residues along the axes of the pores ( 25 ) .
jacoboni and colleagues applied a method combining neural networks and dynamic programming to predict the location of membrane strands ( 26 ) .
the authors estimated that their system correctly predicted about 93% of all known membrane strands .
more recently , martelli et al . developed a sequence - profile - based hmm model that can predict the topology of -barrel mps cyclicing with 6 types of states ( 27 ) .
they reported that the accuracy of per residue of the model was about 83% . lately the following protocol starting from secondary structure prediction and tm segments topology prediction are often used .
secondary structure prediction followed by tm segments identification along with prediction of loops connecting the segments , and molecular dynamics or simulated annealing simulations , may be finally used to refine these primal models . during the last refinement step
, the protein is often inserted into a water / lipid bilayer / water or a water / n - octane / water environment to take into account the presence of the cell membrane .
charmm , gromos , amber , and cvff - insight are some widely used force fields in molecular dynamics calculation .
the slow dynamics of lipid molecules in the bilayer might bring the difficulties in equilibrating the system ( 28 ) .
for globular proteins , the major successful methods for structure prediction include homology modeling , threading , and ab initio folding . along with lucubrating the mechanism of mp folding and increasing the number of high - resolution mp structures , these methods will been applied to the direct prediction of whole mp 3d structures .
the question of how the controlled integration of an mp into the lipid bilayer takes place is still not fully worked out , and there are certainly aspects of mp structures that will probably not be fully appreciated until this step has been accomplished .
some pursuers educed the viewpoint that the prediction of mp structures from amino acid sequences was , in large measure , a problem of physicochemistry ( 29 ) . physical influences that shape mp structures include interactions of the polypeptide chains with water , bilayer hydrocarbon core , bilayer interfaces , and cofactors .
studies on the mechanism of insertion and folding of mps into membranes are relatively rare and have been mostly performed with two model proteins : bacteriorhodopsin ( br ; ref .
while br is a representative -helical bundle protein , ompa belongs to the class of -barrel protein .
homology modeling constructs structures ( targets ) that are homologous to other protein(s ) whose 3d structure is known ( templates ) .
because few high - resolution mp 3d structures are available to be used as templates , and the modeling can be unreliable when the sequence identity between the template and target proteins falls below 20%-30% , the applicability of homology modeling is limited .
an ab initio method was presented ( 32 ) , whose knowledge - based technique added a membrane potential to the energy terms ( pairwise , solvation , steric , and hydrogen bonding ) .
the method is based on the assembly of supersecondary structural fragments taken from a library of highly resolved protein structures using a standard simulated annealing algorithm .
results obtained by applying the method to small mps of known 3d structures showed that the method is able to predict , at a reasonable accurate level , both the helix topology and the conformations of these proteins .
the structure prediction of membrane proteins still remains an interesting scientific problem . because of the physical difference between mps and gps ( globular proteins )
current segment accuracy of reported algorithms are pretty high ( above 90% ) , while the overall accuracy are still around 50%-60% , which gives birth to hand - raising methods to combine the reports from several other algorithms .
the lack of both high - resolution and low - resolution experimental data of mp structures makes the algorithm development and their evaluation difficult , but the fact that most mp sequences are used as space blocks to get through the membrane bilayer that has predefined thickness makes the structure prediction of mps simple on functional aspects .
new algorithms will emerge and reported algorithms will be refined to give a better answer to this problem .
 there is a large gap between the number of membrane protein ( mp ) sequences and that of their decoded 3d structures , especially high - resolution structures , due to difficulties in crystal preparation of mps . 
 however , detailed knowledge of the 3d structure is required for the fundamental understanding of the function of an mp and the interactions between the protein and its inhibitors or activators . in this paper , 
 some computational approaches that have been used to predict mp structures are discussed and compared . 