the protein complex models of rep1:ggtase - ii ( protein data bank [ pdb ] i d 1ltx ) ( fig . 1 ) and that of rab7:rep1 ( pdb i d 1vg9 ) ( fig .
the software yet another scientific artificial reality application ( yasara ) dynamics and structure ( yasara biosciences gmbh , vienna , austria ) ( v 10.12.1 ) was used for generating the molecular model of ternary complex of rab7::rep1::ggtase - ii and also for its molecular dynamic simulation study .
the online q - sitefinder programme ( university of leeds , leeds , uk ) was used for identification of putative dbss over the ternary complex and pymol ( schrdinger , usa ) ( an open source software ) was used for the molecular visualisation . for phylogenetic analysis
, 60 human rab proteins were obtained from uniprot ( http://www.uniprot.org ) and the softwares used for the study were clustalw and phylip ( university of washington , seattle , usa ) .
the structure files of protein complex models of rep1::ggtase - ii ( pdb i d 1ltx ) ( fig . 1 ) and that of rab7::rep1 ( pdb i d 1vg9 ) ( fig .
2 ) were downloaded from pdb and the two complexes were joined using yasara by selecting the strands of rep1 molecule in both the complexes , a single complex was generated .
since the rep1 molecule was repeated , one of the rep1 molecules in the combining complexes was deleted by individually picking up the peptide chain .
the ternary complex of rab7::rep1::ggtase - ii thus obtained was subjected to solvation option of yasara to add water molecules . to study the solvation effects on side chains of the protein complex
a null model pka giving the lowest root - mean - square deviation ( rmsd ) without any shift in value was taken , where pka values were set as 3.22 , 4.09 , 6.20 , 10.8 and 10.76 for asp , glu , his , tyr and lys residues , respectively .
structure of protein complexes of rep1 and ggtase - ii ( pdb i d 1ltx ) .
the constructed and solvated system of ternary complex of rab7::ggtase - ii::rep was subjected to energy minimization without any constraints using steepest descent method followed by simulated annealing method .
particle mesh ewald method was employed to calculate the electrostatic interactions with a cut - off of 10  .
simulation snapshots were taken every 10,000 simulation steps , i.e. , with a timestep of 2.5 fs , i.e. , 10,000 * 12.5=25 ps .
the constructed ternary complex of rab7::ggtase - ii::rep was initially subjected to energy minimisation by amber force field , using dynamics or structure module of yasara .
then , after clicking onto simulation  force field , yamber99 force field was selected in yasara dynamics .
with the selection being made , simulation was run with command - options  macro and movie  play macro and double - clicking the standard macro  em_run.mcr. the molecular dynamics simulation was then carried out in yasara dynamics or structure , where a project directory was created and the structure of ternary complex of rab7::ggtase - ii::rep was stored after energy minimisation . the molecular dynamics simulation was carried out after selecting the saved target molecule and following the command - options  macro and movie  play macro and double - click the standard macro 
md_run.mcr. the simulation performance was increased ( when needed ) by increasing the simulation steps per screen update using 
the data obtained after simulations were analysed for trajectory projection . this was done by using 
this macro created a self - explanatory table and was saved in project directory with  tab - file extension .
the table was then imported in the data visualisation program ( ms - excel ) that included energies and rmsds from the starting structure during simulation .
the standard analysis macro also calculated the time average structure , whose atoms had the b - factors ( calculated from the root mean square fluctuations ) during simulation and the minimum energy structure was saved as 
energymin.sce in the project directory . to find the putative dbss of rab::rep::ggtase - ii complex ,
q - sitefinder programme was used in the present study . the q - sitefinder web site ( www.modelling.leeds.ac.uk/qsitefinder )
was visited and the pdb file of the ternary complex was uploaded and job was submitted .
the q - sitefinder programme uses several separate procedures to perform ligand binding site prediction such as ligandseek , liggrid and probe - clustering programme .
the probe clusters were ranked according to their total interaction energies , with the most favourable being identified as the first predicted binding site .
the clustering programme also calculated site volume , and identified the protein atoms that were within a defined range of cluster sites .
the programme generated the 10 best predicted sites ( interconnected dots ) and any ligands ( coloured by atoms ) were displayed as result .
predicted binding site were listed in the display sites ( atoms forming the binding sites on the protein ) window according to the likelihood ( energy ) of the cluster of probes ( site 1 , site 2 ) .
all together 60 different human rab proteins , relevant to the present study , were analysed for evolutionary relationship by multiple sequence alignment ( msa ) of their amino acid sequences .
the amino - acid sequences of these proteins were obtained from uniprot ( www.uniprot.org ) .
these sequences were then arranged in fasta format for msa by clustalw . the output file ( .phy )
was imported to the phylip for constructing the phylogenetic tree using neighbour - joining method .
the structure files of protein complex models of rep1::ggtase - ii ( pdb i d 1ltx ) ( fig . 1 ) and that of rab7::rep1 ( pdb i d 1vg9 ) ( fig .
2 ) were downloaded from pdb and the two complexes were joined using yasara by selecting the strands of rep1 molecule in both the complexes , a single complex was generated .
since the rep1 molecule was repeated , one of the rep1 molecules in the combining complexes was deleted by individually picking up the peptide chain .
the ternary complex of rab7::rep1::ggtase - ii thus obtained was subjected to solvation option of yasara to add water molecules . to study the solvation effects on side chains of the protein complex , electrostatic interactions were screened .
a null model pka giving the lowest root - mean - square deviation ( rmsd ) without any shift in value was taken , where pka values were set as 3.22 , 4.09 , 6.20 , 10.8 and 10.76 for asp , glu , his , tyr and lys residues , respectively .
structure of protein complexes of rep1 and ggtase - ii ( pdb i d 1ltx ) .
the constructed and solvated system of ternary complex of rab7::ggtase - ii::rep was subjected to energy minimization without any constraints using steepest descent method followed by simulated annealing method .
particle mesh ewald method was employed to calculate the electrostatic interactions with a cut - off of 10  .
simulation snapshots were taken every 10,000 simulation steps , i.e. , with a timestep of 2.5 fs , i.e. , 10,000 * 12.5=25 ps .
the constructed ternary complex of rab7::ggtase - ii::rep was initially subjected to energy minimisation by amber force field , using dynamics or structure module of yasara .
then , after clicking onto simulation  force field , yamber99 force field was selected in yasara dynamics . with the selection being made , simulation was run with command - options  macro and movie  play macro and double - clicking the standard macro  em_run.mcr. the molecular dynamics simulation was then carried out in yasara dynamics or structure , where a project directory was created and the structure of ternary complex of rab7::ggtase - ii::rep was stored after energy minimisation . the molecular dynamics simulation was carried out after selecting the saved target molecule and following the command - options  macro and movie  play macro and double - click the standard macro 
md_run.mcr. the simulation performance was increased ( when needed ) by increasing the simulation steps per screen update using 
the data obtained after simulations were analysed for trajectory projection . this was done by using  macro and movie option and playing the standard analysis macro .
this macro created a self - explanatory table and was saved in project directory with  tab - file extension .
the table was then imported in the data visualisation program ( ms - excel ) that included energies and rmsds from the starting structure during simulation .
the standard analysis macro also calculated the time average structure , whose atoms had the b - factors ( calculated from the root mean square fluctuations ) during simulation and the minimum energy structure was saved as 
to find the putative dbss of rab::rep::ggtase - ii complex , q - sitefinder programme was used in the present study . the q - sitefinder web site ( www.modelling.leeds.ac.uk/qsitefinder )
was visited and the pdb file of the ternary complex was uploaded and job was submitted .
the q - sitefinder programme uses several separate procedures to perform ligand binding site prediction such as ligandseek , liggrid and probe - clustering programme .
the probe clusters were ranked according to their total interaction energies , with the most favourable being identified as the first predicted binding site .
the clustering programme also calculated site volume , and identified the protein atoms that were within a defined range of cluster sites .
the programme generated the 10 best predicted sites ( interconnected dots ) and any ligands ( coloured by atoms ) were displayed as result .
predicted binding site were listed in the display sites ( atoms forming the binding sites on the protein ) window according to the likelihood ( energy ) of the cluster of probes ( site 1 , site 2 ) .
all together 60 different human rab proteins , relevant to the present study , were analysed for evolutionary relationship by multiple sequence alignment ( msa ) of their amino acid sequences .
the amino - acid sequences of these proteins were obtained from uniprot ( www.uniprot.org ) .
these sequences were then arranged in fasta format for msa by clustalw . the output file ( .phy )
was imported to the phylip for constructing the phylogenetic tree using neighbour - joining method .
the molecular model of the ternary complex of rab7::rep1::ggtase - ii ( fig . 3)was obtained using yasara edit option and was subjected to solvation . while studying the solvation effects on the side chains of protein complex ,
it was observed that the predicted pka range of lys , asp , glu , tyr and his were 10.46 - 11.84 , 2.82 - 5.09 , 3.74 - 5.08 , 10.8 and 6.34 - 7.42 , respectively .
the results thus obtained were in good accordance of assumed null model and was thus subjected to further energy minimisation programme and force field calculations .
structure of rab7::rep1::ggtase - ii ternary complex . after running the energy minimisation program , the ternary complex was subjected to simulation up to 500 ps ( 0.5 ns ) .
the different snapshots of simulating steps were saved with respect to set time - intervals ( fig .
thus , generated from these snapshots , reveals a very clear image of rab7 , ggtase - ii and rep1 interaction .
the trajectory was obtained for overall energy simulation of the ternary complex of rab7:ggtase - ii : rep for 500 ps and it revealed that overall energy stabilised after a peak of -6363401.4 kj / mol at 25 ps and tended to remain in plateau phase further for rest of the period ( fig .
this reflected that simulation was achieved with stable energy for rest of the period ( 50 - 700 ps ) for the said ternary complex .
almost , the similar trajectory was obtained for the plots of different energy contributors against simulation run time .
the contribution due to steric parameters like bond , angle , dihedral angle and planarity was found maximum at 25 ps with the values 828405 , 341330.665 , 357593.243 and 2679.375 kj / mol , respectively , which stabilised further to a stationary phase for rest of the period ( 50 - 700 ps ) , similar to the trajectory pattern of overall energy curve ( fig .
the contribution due to columbian charge and van der waal interactions was also in accordance to the above results and was maximum at 25 ps with the values of - 9042365.5 and 1193953.512 kj / mol ( fig . 5f and g ) .
the plots tended to stabilise further and attain the plateau phase for rest of the period ( 50 - 700 ps ) .
these trajectory patterns supports and validates the simulation profile of the ternary complex of rab7::rep1::ggtase - ii .
( a ) total energy versus time , ( b ) bond energy versus time , ( c ) angular energy versus time , ( d ) dihedral angular energy versus time , ( e ) energy of planarity versus time , ( f ) energy due to columbian charge versus time and ( g ) van der waal interaction energy versus time .
the simulation model investigation of the ternary complex of rab7::rep1::ggtase - ii suggests that the c - terminus of rab7 has to be in completely extended conformation during prenylation to reach the active site of rabggtase - ii .
the results obtained in the present study are in good accordance with previous investigations and has provided a good picture of molecular dynamics of ternary complex of rab7::rep1::ggtase - ii and the molecular events during interaction of these three proteins .
all together of about 60 - 70 putative dbss were generated by q - sitefinder programme .
to predict the dbs precisely , these sites were comparatively analysed using an inhouse perl algorithm .
the comparative analysis resulted in construction of comprehensive probe based on clustering method of energy criteria .
the probe was found to conjugate fit with the rab7::rep1 : : ggtase - ii ternary complex ( fig .
6b ) and thus was a good primary indication of the correctness of the generated comprehensive probe .
( a ) the probe generated from q - sitefinder s results analysis and ( b ) ternary complex of rab7::rep1::ggtase - ii with the consensus probe .
the comprehensive probe was used to integrate the putative dbss being generated by different q - sitefinder programme , as mentioned earlier . after an iterative process of selection and omission of different predicted dbss based on statistical parameters , a total of 10 major pockets were identified as putative dbss on rab7::rep1::ggtase - ii complex ( fig . 7 ) .
these pockets reflected the energy of interactions of different probes being generated by the different used algorithms , viz . ;
carbohydrate binding propensity , drug - like compound binding propensity , interacting pocket size threshold of radius 3  in three - dimension .
these 10 major putative dbss thus identified were then mapped onto the rab7::rep1::ggtase - ii ternary complex ( fig .
7 ) which revealed that a total of two major pockets of dbs ( one in the -barrel and one in -helical portion ) were located on the rep protein of rab7::rep1::ggtase - ii ternary complex . of these two dbs ,
one was found to be located within the major cavity produced at the interaction point of rab7-rep complex with ggtase - ii .
the other eight major putative dbss identified were found to be located on ggtase portion of the ternary complex .
the -strand of ggtase was found to have two major pockets of putative dbss while the -strand possessed the other six pockets of putative dbss .
all these putative dbs mapped completely over the rab7::rep1::ggtase - ii ternary complex , which was an indirect evidence of completeness of the computational calculation .
these results have paved the path for further investigation into assigning rigidity to these putative dbss and develop pharmacophoric patterns for the same .
these pharmacophore can be further used for virtual screening of ligand databases ( viz . ,
pubchem , zinc ) and are being presently investigated in our lab , apart from generating combinatorial libraries .
the drug binding sites mapped on the ternary complex of rab7::rep1::ggtase - ii . after performing msa of the selected 60 different human rab proteins and constructing the phylogenetic tree by neighbour - joining algorithm of phylip of the same , it was observed that human rab7a , rab7b , rab9a and rab9b share the common ancestor and are the nearest neighbours during the course of evolution .
so , the molecular dynamics studies of rab7::ggtase - ii::rep ternary complex will also provide a good platform for further investigation of related human rab proteins such as rab7a , rab7b , rab9a and rab9b .
rab7 and rab9 are related to each other being localised to late endosomes and thus the phylogenetic analysis in the present study is coherent with the earlier results.[2931 ] to conclude , we can say that the present study has produced a good insight into the molecular events occurring during the ternary complex formation of rab7::rep1::ggtase - ii .
the study also draws its significance because the protein studied is a membrane protein whose structural analysis is difficult through experimental way .
the present in silico study was a primary and initial attempt in this direction and has paved pathway for further in silico studies like identifying the potential drug targets by searching the putative dbss , generating pharmacophoric pattern , searching or constructing suitable ligand , docking studies .
 the structure - function correlation of membrane proteins have been a difficult task , particularly in context to transient protein complexes . 
 the molecular simulation of ternary complex of rab7::rep1::ggtase - ii was carried out to understand the basic structural events occurring during the prenylation event of rab proteins , using the software yasara . 
 the study suggested that the c - terminus of rab7 has to be in completely extended conformation during prenylation to reach the active site of rabggtase - ii . 
 also , attempt was made to find putative drug binding sites on the ternary complex of rab7::rep1::ggtase - ii using q - sitefinder programme . 
 the comprehensive consensus probe generated by the program revealed a total of 10 major pockets as putative drug binding sites on rab7::rep : : ggtase - ii ternary complex . 
 these pockets were found on rep protein and ggtase protein subunits . 
 the rab7 was found to be devoid of any putative drug binding sites in the ternary complex . 
 the phylogenetic analysis of 60 rab proteins of human was carried out using phylip and study indicated the close phylogenetic relationship between rab7 and rab9 proteins of human and hence with further in silico study , the present observations can be extrapolated to rab9 proteins . 
 the study paves a good platform for further experimental verifications of the findings and other in silico studies like identifying the potential drug targets by searching the putative drug binding sites , generating pharmacophoric pattern , searching or constructing suitable ligand and docking studies . 
MATERIALS AND METHODS
Construction of model of ternary complex of Rab7::REP1::GGTase-II:
Simulation of ternary complex of Rab7::REP1::GGTase-II:
Analysis of simulation trajectory of ternary complex of Rab7::GGTase-II::REP:
Identification of putative DBS of Rab7::REP1::GGTase-II:
Phylogenetic analysis of different Rab proteins found in humans:
RESULTS AND DISCUSSION