Identification of a Novel Class of Farnesylation Targets by Structure-Based Modeling of Binding Specificity
[['Protein prenylation is a post-translational modification in which a prenyl group ( farnesyl or geranylgeranyl ) is attached to the protein via a thioether bond to a cysteine at or near the carboxy terminus of the protein ( reviewed in [1] , [2] ) .', 'Protein farnesyltransferase ( FTase ) and geranylgeranyltransferase type I ( GGTase-I ) are also called CaaX prenyltransferases , due to their ability to catalyze modification of peptides and substrate proteins bearing the carboxy terminal ( C’ ) Cys-aliphatic-aliphatic-variable amino acid ( Ca1a2X ) motif [3] .', 'Upon binding of the substrate and the C-terminal Ca1a2X motif , the catalytic zinc ion of FTase coordinates the thiol side chain of the cysteine and catalyzes the covalent attachment of the lipid anchor to this residue .', 'A detailed view of this mechanism has been obtained by a series of structures solved at different stages of the reaction [4] .', 'After the covalent attachment of the isoprenoid in the cytoplasm , substrate proteins can undergo further processing , resulting in a C’ structure that is able to serve as a specific recognition motif in certain protein-protein interactions [5] and to direct the modified protein towards incorporation into cellular membranes [6] .', 'A wide range of proteins involved in diverse cellular functions require this post-translational modification for their action [2] .', 'While numerous proteins have been experimentally shown to undergo farnesylation in vivo [7] , [8] , [9] , it is likely that many FTase substrates remain to be discovered .', 'There is a wide interest in the mapping of FTase targets in the genome , in part due to the therapeutic potential of FTase inhibitors against cancer [10] , [11] , [12] , as well as parasitic infection [13] , [14] .', 'Identification of new targets might lead to novel therapeutic approaches [15] .', 'Moreover , the elucidation of cellular FTase targets might shed light on the function of various proteins , as well as on the cellular network of interactions .', 'Computational approaches have predicted FTase targets based on sequence features of known targets [7] , [8] .', 'These methods show good performance in terms of sensitivity , i . e . known targets are correctly identified .', 'Thus , prenylation is mainly defined by the last four residues of the protein , although additional weaker sequence constraints have also been identified upstream in the sequence [16] .', 'Other approaches were based on manual inspection and derived from structural features [9] .', 'Substrate specificity has also been examined using peptide libraries .', 'A comprehensive study by Hougland et al . on the farnesylation of a large synthetic peptide library has allowed a detailed characterization of FTase specificity [17] .', 'In addition to compiling a large and clean dataset of peptides that contains both efficient substrates and non-substrates for FTase , this study discovered a third group of sequences that are farnesylated only under single-turnover ( STO ) conditions ( [E]>[S] ) .', 'Analysis of peptide substrates has also demonstrated that reactivity depends on synergy between the side chains at the a2 and X positions [18] .', 'These findings indicate that FTase substrate recognition is more complex than the simple Ca1a2X motif model , and that non-canonical sequences can serve as substrates .', 'A large number of structures have been determined for FTase and FTase-substrate peptide complexes [19] .', 'The peptide binding pocket is well-characterized , although a structure of the ternary FTase•farnesyl diphosphate ( FPP ) •peptide in an active conformation has not been determined [9] .', 'The Ca1a2X cysteine sulfur atom ( prior to the product formation ) coordinates the catalytic Zn2+ ion together with side chains ( D297 , C299 and H362 ) of the FTase β-subunit .', 'The a1 side chain points out of the binding pocket and faces the solvent , while the a2 side chain is buried within the binding pocket and interacts both with the farnesyl chain of FPP and the residues lining the pocket .', 'The C’ X position interacts with residues mostly from the FTase β-subunit and is considered the main determinant for the specificity between FTase and GGTase-I 9 .', 'Finally , two highly conserved hydrogen bonds are formed:', '1 ) between the C-terminal carboxylate group and the side chain of FTase Q167α and', '2 ) between the a2 backbone carbonyl oxygen and the side chain of FTase R202β ( Figure 1 ) .', 'Despite this detailed structural information , only a handful of different peptide sequences have been solved in complex with FTase .', 'We previously developed a scheme for modeling the structures of peptide-protein complexes ( Rosetta FlexPepDock [20] ) , which is incorporated within the Rosetta modeling suite framework [21] .', 'This protocol is the starting point for the development of a structure-based scheme for the prediction of peptide binding specificity ( FlexPepBind ) .', 'Specifically , to refine FlexPepBind for the prediction of FTase binding peptides , we have incorporated constraints derived from the conserved features in solved FTase structures and adapted the energy function to distinguish between reacting and non-reacting tetrapeptides ( based on an underlying assumption that tetrapeptides that bind will react , while those that do not bind will not react ) .', 'We trained and tested this protocol against the recent dataset published by Hougland et al . [17] .', 'Validation of the protocol against several independent sets showed accurate prediction of peptides that could be farnesylated , both under multiple turnover ( MTO ) and single turnover ( STO ) conditions .', 'Evaluation of all possible Cxxx peptides identified a previously uncharacterized class of farnesylation targets that contain an acidic C-terminal residue .', 'The 13 peptides predicted to bind with best affinity were experimentally shown to indeed undergo farnesylation in vitro .', 'Finally , a genomic scan for novel FTase targets revealed 77 novel putative FTase targets previously undetected by sequence-based approaches .', 'Among these , 13 out of 16 selected novel putative farnesylation targets were indeed farnesylated by FTase in an in vitro experimental validation .', 'FTase-peptide binding is a model system for our approach to peptide-protein binding specificity prediction and design .', 'Our protocol can easily be adapted to additional peptide-protein interactions where both experimental structure and affinity data are available , thereby providing a mechanism to identify targets not detectable by sequence conservation only .'], ['Recently Hougland et al . performed a large-scale study , in which they characterized a TKCxxx peptide library for reactivity with rat protein farnesyltransferase ( rat FTase ) [17] .', 'Out of an initial library of 213 sequences , 77 peptides are farnesylated under multiple turnover ( MTO ) conditions , and 51 sequences are not farnesylated under any conditions .', 'Interestingly , the remaining 85 sequences are farnesylated under single turnover ( STO ) conditions but not under MTO conditions .', 'We set out to use FlexPepBind and the structural data available for FTase to discriminate MTO sequences from non-reactive ( NON ) peptide sequences , using the 77 MTO and 51 NON peptide sequences as our training set ( 128 peptides in total; Dataset S1A ) .', 'Towards this end , we used the high resolution structure of human FTase in complex with a peptide derived from the carboxy terminus of Rap2a and a farnesyl diphosphate ( FPP ) analog ( PDB: 1tn6 [9] ) to create a starting model .', 'The bound peptide was truncated to include only the terminal Ca1a2X motif .', 'Different peptide sequences were then threaded onto the peptide backbone and used as starting structures .', 'Initially , we modeled peptide-FTase complex structures for different peptide sequences by applying the Rosetta FlexPepDock protocol to the threaded starting models .', 'This protocol was developed previously in our lab for the modeling and refinement of peptide-protein complex structures to high resolution [20] .', 'Our simulations included three constraints , namely the conservation of the 2 structurally conserved hydrogen bonds ( C’ carboxylate - FTase Q167α; a2 backbone carbonyl oxygen - FTase R202β ) and the location of the cysteine sulfur atom coordinating the Zn2+ ion ( Figure 1 , see Methods for more details ) .', 'For each simulation , the energy of the best scoring Cxxx peptide was extracted ( see Methods for further details ) .', 'Figure 2A shows the Receiver Operating Characteristic ( ROC ) plot for the ability of the peptide energy to discriminate between MTO sequences and non-substrate sequences .', 'The plot shows very good discrimination with an Area Under the ROC Curve ( AUC ) value of 0 . 915 on our training set .', 'These results demonstrate that a structure-based evaluation of the peptide energy can distinguish very well between farnesylated and non-farnesylated peptide sequences .', 'Since the known constraints restrict the simulation to a closely defined region in the binding site , we reasoned that a simpler and faster protocol might be able to model the peptides with similar accuracy .', 'Our simplified protocol therefore includes only a minimization using the Rosetta energy function [21] , [22] under constraints to retain the 2 structurally conserved hydrogen bonds and the cysteine sulfur atom location coordinating the Zn2+ ion ( see above and Methods for more details ) .', 'This protocol yielded similar results with an AUC value of 0 . 875 on the training set .', 'A peptide energy threshold of -0 . 4 ( i . e . sequences with energy below/above -0 . 4 are predicted to be binders/non-binders and therefore farnesylated/non-farnesylated , respectively ) corresponds to a 69% True Positive Rate ( TPR ) and 8% False Positive Rate ( FPR ) .', 'A more stringent threshold of -1 . 1 energy units corresponds to a 44% TPR and 2% FPR ( Figure 2A ) .', 'With the two protocols exhibiting similar performance , we decided to proceed further using the fast minimization protocol .', '( Performance on the training set using additional sampling and scoring schemes is summarized in Table S1 . )', 'To assess FlexPepBind using the selected thresholds , we evaluated performance on three independent test sets ( Dataset S1B-D online ) .', 'Using FlexPepBind , we modeled all of the 8000 possible Cxxx sequences and scored them according to our protocol .', 'The thresholds for the discrimination of MTO/NON predict that 1349 ( 17%; stringent threshold\u200a=\u200a−1 . 1 ) and 2309 ( 29%; threshold\u200a=\u200a−0 . 4 ) of all tetramer peptide sequences could be possible substrates ( see Figure 3 ) .', 'This set of putative farnesylation targets suggest a much more versatile binding motif than previously accepted ( see Figure 4 ) : while position a2 of the Ca1a2X motif is still prominently aliphatic ( ILE/VAL/LEU/PHE ) , positions a1 and X are less restricted than previously reported ( compare Figure 4C to Figures 4A&B ) .', 'In particular , we identify within this set a novel class of farnesylation targets that contain an acidic residue at the C-terminus ( 238/1349 putative targets; ∼20%; see Figure 4D ) .', 'Figure 4C indicates that the minimization-based protocol tends to miss larger residues at the C-terminal X position .', 'Indeed , assessment of the prediction accuracy for this position on the training set shows that only 1/8 CxxF and 0/3 CxxW sequences are correctly predicted with the chosen protocol ( CxxM peptides are predicted with higher accuracy: 10/14 ) .', 'Using the FlexPepDock based protocol , performance increases to: 6/8 CxxF; 2/3 CxxW and 11/14 CxxM , demonstrating that CxxF peptides are indeed rescued by the additional backbone flexibility .', 'Therefore , it might be advisable to use the FlexPepDock based protocol for peptides that contain a bulky C-terminal side chain .', 'We compared our predictions to the PrePS [7] prediction of prenylation targets on the initial training set of peptides .', 'Regarding the discrimination of MTO substrates from non-active peptides , PrePS results are comparable to FlexPepBind ( AUC of 0 . 92 , with a threshold corresponding to 60% TPR for 2% FPR ) .', 'However , the performance for STO peptides is significantly better for our structure-based approach: while FlexPepBind recovers 47% and 32% of the STOs with the loose and stringent thresholds concordantly , PrePS predicts only 14% of these sequences as substrates .', 'Since our retrospective studies indicated that our approach can very accurately retrieve actual farnesylation targets , we were interested in testing it prospectively – could novel targets be indeed identified ?', 'We selected the 13 best scoring peptides ( i . e . predicted tightest binders ) , yet previously uncharacterized for experimental validation .', 'These are mostly ‘non-canonical’ peptides , including 5 peptides with an acidic C-terminal residue .', 'Indeed , PrePS [7] predicts only 2 out of the top-scorers to be FTase substrates .', 'In vitro farnesylation assays indicate that all of these peptides indeed undergo farnesylation catalyzed by FTase: 10 under MTO conditions and 3 under STO conditions ( Table 1A ) .', 'These results demonstrate the robustness of our protocol and its exceptional accuracy .', 'Importantly , they confirm the novel class of farnesylation substrates that contain a negatively charged C-terminal residue ( Figure 4D ) .', "Structural investigation of this novel class of substrates suggests that the negatively charged C' side-chain is stabilized by FTase residue His 149βwhile accepting a hydrogen bond from Trp102β ( GLU ) and creating an additional hydrogen bond with the side-chain of Ser99β ( GLU & ASP ) ( see Figure S1 ) .", 'Additional polar interactions with water molecules are possible but were not explicitly modeled .', 'Equipped with a score that can predict both known and novel FTase targets , we set out to scan the human genome for proteins that may undergo farnesylation .', 'Our protocol was developed based on experimental assays on rat FTase ( and the structure of human FTase [9] ) .', 'Since rat and human FTases show very high sequence identity ( 92% and 96% for subunits α and β respectively ) , and none of the sequence differences are located at or near the peptide binding site , we are confident that our prediction scheme can be applied to human farnesylation as well .', 'We identified 756 unique proteins in human SwissProt [24] that contain the Cxxx motif at their carboxy terminus .', '167 and 309 of these protein sequences obtained scores lower than the −1 . 1 and −0 . 4 threshold values , respectively , indicating that these proteins might be farnesylated by FTase .', 'We focused on the group of 167 proteins detected with the more stringent threshold .', 'Could these proteins indeed be FTase substrates ?', 'Several indications support our predictions: First , amongst the 167 candidates , 42 contain a Cxxx motif of a known FTase substrate .', 'Secondly , the Gene Ontology ( GO ) [25] cellular compartment annotation for most of these 167 proteins is Membrane related ( see Figure S2; see Methods for more details ) .', 'This supports their association with membranes , possibly by farnesylation ( albeit this localization annotation might have been inferred from sequence similarity ) .', 'Furthermore , peptide library studies have demonstrated FTase-catalyzed farnesylation ( under STO or MTO conditions ) of 50 of these Cxxx motifs ( representing 66 human proteins ) [17] .', 'Finally , analysis of the putative target proteins with the PrePS server predicts that most of them ( 90/167 ) are indeed FTase targets , while the other 77 are not predicted to be farnesylated ( see Figure S3 ) .', 'To further characterize the latter , we proceeded with in vitro experimental validation of selected sequences .', 'Among these 77 proteins ( containing 72 unique Cxxx motifs ) , 39 motifs had not yet been tested for in vitro farnesylation .', 'The second set chosen for experimental validation consisted of 16 top-scoring peptides selected from these 39 motifs .', 'Of the 16 tested peptides , 9 and 4 peptides are farnesylated in vitro under MTO and STO conditions , respectively , while only 3 were not farnesylated by FTase ( Table 1B ) .', 'None of the 16 sequences in this second set are predicted to serve as farnesylation targets by PrePS .', 'Interestingly , for 9 of these 16 sequences , PrePS predicts that the upstream context of the motif is suitable for farnesylation .', 'In these cases , the PrePS negative prediction is based on the sequence of the Cxxx motif .', 'This suggests that improved characterization of the contribution of the 4 C-terminal residues to farnesylation can identify more farnesylation targets .', 'Finally , for 8 of these 16 sequences , PrePS would predict farnesylation of the Cxxx motifs in the background of the favorable H-Ras upstream sequence .', 'The balance between the upstream signal and the C-terminal Cxxx motif is therefore an interesting subject for future research .', 'Most of the proteins identified by this study as novel FTase substrates have not been well characterized to date .', 'Consequently , in vivo experiments that evaluate the cellular localization and prenylation status of these proteins , in conjunction with the in vitro farnesylation demonstrated in this study , will advance their functional characterization .'], ['The protocol that we developed essentially estimates the binding affinity of FTase for Cxxx peptides , using a training set of reactive peptides , rather than predicting the farnesylation activity of these sequences .', 'This has several implications and limitations .', 'Remarkably , the ability to discriminate peptides that undergo MTO reaction from non-active peptides according to binding energy suggests that the non-active peptides may bind weakly or not at all to FTase ( see Figure 3 ) .', 'This finding is supported by results from an in vitro inhibition experiment in which none of the tested non-active peptides inhibited FTase-catalyzed farnesylation of a known substrate [17] .', 'In turn , the members of the small class of FlexPepBind false positive peptides may bind to FTase with high affinity but still not be farnesylated .', 'These false positive peptides could therefore serve as FTase inhibitors and represent an interesting set to characterize in future work .', 'Previous studies have shown that the sequence immediately upstream of the conserved cysteine residue may also play a role in substrate selectivity [16] .', 'These sequences modulate peptide affinity and reactivity with FTase , i . e . a high-affinity terminal tetramer sequence does not necessarily ensure farnesylation of the protein .', 'For half of the proteins tested in the study , the PrePS [7] program predicts favorable upstream sequences .', 'This result coupled with the high-affinity -Cxxx motif predicted by FlexPepBind ( see Results and Table 1B ) increases the confidence that the human proteins containing the said Cxxx motif could be farnesylated in vivo .', 'In turn , a favorable upstream sequence might compensate for a weak C-terminal signal .', 'Our future work will therefore further characterize the balance between these two signals in determining farnesylation .', 'One puzzling aspect of FTase substrate recognition is the large number of peptides that exhibit single turnover activity .', 'The single turnover rate constant , kfarn , reflects all of the rate constants up to but not including the release of the farnesylated product [4] , [26] , [27] , [28] .', 'Therefore , the STO peptides bind to FTase and are readily farnesylated , but the product dissociates very slowly so multiple turnover activity is very slow .', 'Consistent with this , FlexPepBind achieves an AUC value of 0 . 776 in the discrimination between STO and non-active peptides on the training set , indicating that STO peptides have higher affinity for FTase than the non-active peptides ( see Figure 3 ) .', 'Our protocol thus identifies STO peptides much better than sequence-based methods ( see Results and Hougland et al . [17] ) .', 'What then discriminates between MTO and STO peptides ?', 'Hougland et al . postulated that the farnesylated STO peptides might bind more tightly to FTase than farnesylated MTO peptides , and as a consequence FPP-catalyzed product dissociation is slow [17] .', 'However , binding energy , as approximated by our approach , seems to be a poor discriminator between MTO and STO peptides ( AUC value of 0 . 625 on the training set – Dataset S1B ) .', 'That is , estimation of the binding affinity of peptides in the context of static conformations of the protein cannot explain the difference in reactivity .', 'Furthermore , application of this approach to models of MTO and STO peptides at different stages of the reaction sequence ( pre-farnesylation , post-farnesylation with the farnesyl group in the exit groove ) was not able to account for this difference as well .', 'Hence , rather than binding affinity , a parameter related to the dynamics of product dissociation might dictate turnover .', 'We therefore conclude that a dynamical approach , such as molecular dynamics , will be required to explain the mechanism that distinguishes STO from MTO peptides .', 'Past in vitro peptide farnesylation experiments with FTase have measured kcat/KMpeptide under MTO conditions and kfarn rate constants under STO conditions [17] .', 'The estimated reactivity of MTO and STO peptides ( see Methods ) measured in this work falls within the range of previously measured activity [17] .', 'Therefore , these peptides have comparable reactivity to other substrates , including peptides that correspond to proteins that are farnesylated in vivo .', 'Measured under MTO conditions , the kinetic parameter kcat/KMpeptide is termed the specificity constant and best reflects the reactivity of an enzyme in the presence of multiple substrates , as observed in vivo [29] .', 'In a cell , the reactivity of a protein substrate with FTase depends on the value of kcat/KMpeptide as well as on the concentration of the substrate within the cytosol .', 'Although a protein substrate with a higher value of kcat/KMpeptide is more likely to be farnesylated in vivo , it is unclear what level of in vitro activity corresponds to a true FTase substrate in vivo .', 'Furthermore , in vivo the optimal levels of farnesylation of a given substrate may vary and a low fraction of modification may still be biologically relevant .', 'Additionally , a substrate must be localized to the proper cellular locale in order for modification to occur and the C-terminus of the protein must be structurally available .', 'Peptide library studies and this work have aided in determining potential FTase substrates and have also identified already known substrates , but more work is needed to characterize the reactivity of these substrates in vivo .', 'As for the STO-only peptides , these substrates are readily farnesylated but the product does not dissociate rapidly .', 'One possibility is that these proteins function as FTase inhibitors and consequently play a regulatory role within the cell [17] .', 'However , both FPP and peptides have been implicated in catalyzing product dissociation of farnesylated STO peptides [17] , [30] , [31] and therefore it is possible that other cellular components could activate product dissociation allowing rapid farnesylation of these proteins in vivo .', 'Therefore , competition or synergy among different FTase substrates could play an important functional role for modification and localization of proteins .', 'Improved identification of STO peptides using the structure-based FlexPepBind approach presented here will expand our understanding of regulatory aspects of this reaction .', 'In addition , the overlap in substrate preference of FTase and GGTase-I [3] indicates that modulation of the type of prenyl modification ( e . g . changes in relative enzyme availability or magnesium concentration ) might be functionally important as well .', 'Our future focus on structure-based characterization of GGTase-I specificity will allow an improved investigation of this regulatory feature , complementary to sequence-based studies [7] , [8] .', 'Scanning the human genome for putative FTase targets using our structure-based approach revealed many putative , not yet detected , farnesylated proteins .', 'These new farnesylation substrates may provide novel disease targets for farnesyltransferase inhibitors .', 'Moreover , the prediction that these proteins are farnesylated might shed light on their function .', "As an example , the putative proteins Q8NA34 , A6NHS1 , and P0C7P2 ( UniProt identifiers [24] ) all contain C' sequences strongly predicted to serve as farnesylation targets suggesting that the proteins are membrane localized .", 'Additionally , our method also predicts FTase substrates that have recently been identified from in vivo experiments .', 'For example , Kho et al . used a tagging-via-substrate proteomic approach to discover novel farnesylation targets [32] .', 'They found a total of 18 farnesylated proteins: 13 are well known , and of the remaining 5 our approach predicts 4 to be farnesylated , including one hypothetical protein .', 'Furthermore , it was recently found that pathogens can hijack the host farnesylation machinery to their own advantage , for example , anchoring effector proteins to the membrane of Legionella-containing vacuoles [33] , [34] , [35] .', 'Thus , in addition to the identification of putative new farnesylation targets in the human genome , FlexPepBind can be used to scan pathogen genomes for farnesylation as well .', '13/16 motifs derived from human proteins tested for in vitro farnesylation indeed undergo the reaction .', 'Will this also happen in vivo ?', 'In the following we compile additional available details on these targets that might help answer this question .', 'One way to assess the in vivo relevance of the observed in vitro ability to undergo farnesylation of the C-terminus of a protein is to look for homologous proteins that also undergo farnesylation .', 'Such information can easily be retrieved from PRENbase [8] .', 'A search in this database revealed that Kinesin-like protein KIF21B variant ( Q2UVF0; CFLT ) maps to a cluster of 9 highly similar eukaryotic sequences ( E-val<e-20 ) that are all predicted to undergo farnesylation by PrePS .', 'Similarly , Ankyrin repeat and BTB/POZ domain-containing protein BTBD11 ( A6QL63-3; CWLS ) maps to a cluster of 25 sequences of related proteins in PRENbase .', 'Zinc finger protein 64 homolog ( Q9NTW7-3; CYVA ) also contains a number of conserved homologs in PRENbase , however in this specific isoform the target cysteine is part of the Zinc-finger structural motif , and therefore it might not readily be farnesylated .', 'Another interesting putative farnesylation target that we have identified is the short isoform of Intersectin-2 protein ( Q9NZM3-3; CCLS ) .', 'This protein is involved in clathrin-mediated endocytosis [36] , [37] , and farnesylation could be a mechanism for regulation and localization to the membrane , similar to the prenylation of Rho GTPases for endocytosis [38] .', 'In particular , the long isoform of intersectin-2 contains additional domains [39] , including a PH domain known to bind phosphoinositides [40] , and a C2 domain known to be involved in Ca-dependent and independent binding of phospholipids [41] .', 'Consequently , in the short isoform that lacks these domains , farnesylation might indeed be used as an alternative way to achieve membrane proximity and attachment .', "While the localization of some Rho GAP proteins ( e . g . p190 [42] ) is regulated by phosphorylation , the short isoform of Rho GTPase-activating protein ( GAP ) 19 ( Q14CB8-5; CSLI ) exposes a new C' motif that may target it to the membrane ( while keeping the Rho GAP domain intact ) .", 'The same goes for MAPKAP1 isoform 6 ( Q9BPZ7-6; CKLA ) , a subunit of mTORC2 .', "While the full length protein was shown to contain a functional PH and Ras binding domains [43] , the truncated isoform reveals a C' putative farnesylation motif instead .", 'Thus , for all but three MTO sequences we could gather additional information that supports actual in vivo farnesylation .', 'We further discuss alternative splicing as a regulatory mechanism below .', 'Four motifs were found to undergo in vitro farnesylation under STO conditions .', "The Homeobox protein ESX1 ( Q8N693; CPFF ) is cleaved into an N' and C' domain; while the N' enters the nucleus , the C' domain is localized to the cytoplasm where it inhibits cyclin degradation[44] .", 'A search for homologues in PRENbase produced a cluster with 2 sequences predicted to undergo farnesylation by PrePS .', 'While the latter could support actual farnesylation of this protein , in this case this modification would serve for purposes other than membrane association , such as the interaction with new partners [5] .', "Isoform 2 of the integral membrane protein solute carrier family 7 member 13 ( Q8TCU3-2; CHFH ) is missing an intracellular domain , and therefore places its C' in proximity to the membrane .", 'Here farnesylation could play a role in targeting this transmembrane protein to a specific membrane compartment [45] , resulting in different membrane distributions for alternative spliced isoforms .', 'Decaprenyl-diphosphate synthase subunit 1 isoform ( Q5T2R2-2; CTTE ) is a nuclear encoded mitochondrial protein .', 'If indeed farnesylated , this would be a first example where an isoform of a mitochondrial protein is farnesylated in the cytosol .', 'Finally , the proton-coupled amino acid transporter 1 ( Q7Z2H8; CAFI ) is likely not a farnesylation target , since mutation of the target cysteine to alanine did not affect its function [46] .', 'As discussed above , the biological role of farnesylation under STO conditions is not yet clear; furthermore , if these proteins are farnesylated in vivo , the function is likely more complex than localization to the membrane .', 'For the three motifs that were not farnesylated under in vitro conditions , additional information about the cognate proteins indeed suggests that the C-terminal cysteines are likely not farnesylated in vivo .', 'The target cysteines of Growth/differentiation factor 15 ( Q99988; CHCI ) and the extracellular C-type lectin domain family 2 member D isoform ( Q9UHP7-3; CLFE ) are part of a conserved disulfide bridge and therefore most likely not farnesylated in vivo .', 'In this study , we chose peptide motifs for in vitro experimental characterization based on their predicted ability to bind FTase and their novelty ( i . e . not predicted by PrePS , and not yet experimentally tested ) .', 'While our post-hoc literature analysis reinforces some of the predictions , other targets will apparently undergo farnesylation only in vitro .', 'The latter represent an interesting set of proteins that allow the investigation of additional factors that regulate the actual farnesylation in vivo , and that therefore distinguish between the ability of a protein to undergo farnesylation in vitro and in vivo .', 'In any case , future in vivo validation is required for all putative targets to unequivocally define their functional importance in the cell .', 'Approximately half of the proteins strongly predicted by FlexPepBind to undergo farnesylation ( 86/167 ) appear in alternative splicing isoforms ( according to Swissprot [24]; the actual number of isoforms is expected to be higher , as more experimental data accumulate from large scale sequencing efforts ) .', 'Among these 86 proteins , most ( 61 ) contain the Cxxx motif only in some of the isoforms .', 'This may present a second layer of regulation for the localization of such proteins , in which a protein can reside in different cellular compartments as a function of the isoform expressed at a given time or a given tissue and therefore perform different functions .', 'This form of regulation may be a consequence of the irreversible nature of farnesylation .', 'On the other hand , farnesylation can be maintained despite alternative splicing .', 'For example , in Rab28 the two reported isoforms ( hRab28S , hRab28L ) differ only by a 95-bp insertion within the coding region [47] .', 'This insertion generates two alternative sequences in the 30 C-terminal amino acids , which strikingly both contain a high-affinity farnesylation motif ( CSVQ – L isoform , CAVQ – S isoform ) 
['Introduction', 'Results', 'Discussion', 'Methods']
['Farnesylation is an important post-translational modification catalyzed by farnesyltransferase ( FTase ) .', 'Until recently it was believed that a C-terminal CaaX motif is required for farnesylation , but recent experiments have revealed larger substrate diversity .', 'In this study , we propose a general structural modeling scheme to account for peptide binding specificity and recapitulate the experimentally derived selectivity profile of FTase in vitro .', 'In addition to highly accurate recovery of known FTase targets , we also identify a range of novel potential targets in the human genome , including a new substrate class with an acidic C-terminal residue ( CxxD/E ) .', 'In vitro experiments verified farnesylation of 26/29 tested peptides , including both novel human targets , as well as peptides predicted to tightly bind FTase .', 'This study extends the putative range of biological farnesylation substrates .', 'Moreover , it suggests that the ability of a peptide to bind FTase is a main determinant for the farnesylation reaction .', 'Finally , simple adaptation of our approach can contribute to more accurate and complete elucidation of peptide-mediated interactions and modifications in the cell .']
['Linear sequence motifs serve as recognition sites for protein-protein interactions as well as for post-translational modifications .', 'One such motif is the CaaX box located at protein C-termini that serves as prenylation site .', 'This prenylation is critical for many signal transduction related proteins and it is thus an important goal to uncover the range of prenylated proteins .', 'Due to poor generalization ability , sequence based computational methods can only go so far in predicting novel targets .', 'In this study , we introduce a novel structure based modeling approach that allows both recovery of known farnesylation substrates , as well as detection of a new class of farnesylation targets .', 'We demonstrate high accuracy in retrospective discrimination between substrates and non-substrates of farnesyltransferase ( FTase ) .', 'More importantly , in a prospective study , in vitro experiments validate that 26/29 predicted peptides indeed undergo farnesylation .', 'These novel peptides were derived either from actual human proteins , or predicted to bind particularly well to FTase .', 'Other than the discovery of putative novel farnesylation targets in the human genome , as well as possible inhibitors , we provide insights into the main determinants of farnesylation .', 'Our approach could be easily extended to additional peptide-protein interactions and help the elucidation of the cellular peptide-protein interaction network .']
['biochemistry', 'enzyme structure', 'biochemical simulations', 'enzymes', 'protein structure', 'biology', 'computational biology', 'macromolecular structure analysis']