in the present world , the prime requisite is to maintain global sulphur balance in the environment including wadden sea and other associated ocean belts .
this therefore would make a better and safe environment to dwell in . in the foremost stage ,
formation of thiosulphate ( s2o3 ) occurs from sulphide ( s ) in aerobic situations . at this step ,
the end product is such that sulphur has  + 2  oxidation state but it needs to be transformed to its elemental form ( s ) . in this phase , microorganisms play a crucial role to completely oxidize thiosulphate ( the most abundant sulphur anion ) to maintain the proper sulphur cycle throughout the environment .
. moreover , deep ocean bed surviving microorganisms use sulphides and such varied sulphur compounds for their chemolithoautotrophy [ 2 , 3 ] .
one such microorganism is an epsilon - proteobacterium sulfurimonas denitrificans [ 2 , 3 ] .
sox multienzyme complex  within their genes coding for the necessary proteins or few of them follow the pathway for  reverse sulphate reduction  .
the  sox multienzyme complex organization  encodes a cluster of genes ( sox operon , sulphur oxidizing operon ) and is the most widely predominant one [ 47 ] .
nevertheless , it is worth mentioning that the epsilon - proteobacteria sulfurimonas denitrificans ( s. denitrificans ) retain a sox multienzyme complex .
intriguingly , the explicit participation of soxyz complex , soxxa complex , soxcd complex , and the soxb protein entirely accomplishes the thiosulphate oxidation mechanism [ 6 , 9 ] . amongst them , the molybdenum - cofactor preserving protein is soxc , while soxd protein is a cytochrome protein of type c. it comprises duo - heme clusters : one is well known as d1 and the other one is known as d2 . altogether , soxcd aids in the performance of sulfane dehydrogenase . first and foremost , it is essential to study the interactive complex of soxcd protein from s. denitrificans .
further after that , the interactive nature amongst soxy and soxz proteins is necessary to be explored .
after that , the soxyzcd complex formation and the residual participation for thiosulphate oxidation are highly paramount to be examined .
this is , moreover , based upon the wet - laboratory experiments and their documentation , which states soxyz to be crucial and to serve as an important interacting partner for soxcd to complete the thiosulphate oxidation process [ 8 , 9 ] . in the presence of soxyzcd , sox system is proficient to oxidize elemental sulphur , sulphide , and explicitly sulfane sulphur ( -s- ) in thiosulphate [ 4 , 6 , 7 , 9 ] . within the complex ,
further soxcd aids in recycling the thiosulphate - bound protein soxy from soxyz complex via releasing of the sulphate group from soxy protein , hydrolytically [ 4 , 6 , 7 , 9 ] .
formerly , the modeling of soxc and soxd protein from sulfurimonas denitrificans was performed via simulation .
but hitherto , there was no such molecular and residual disclosure for the analysis of the entire pattern of interaction in soxyzcd with thiosulphate for the oxidation of thiosulphate in wadden sea and deep ocean belts .
in addition to that , even previously the analysis of the conformational switches in soxyz that benefits the recycling process along with the thermodynamic stability in soxyz for the fulfilment of the entire phenomena in wadden sea and associated ocean belts was not dealt with .
so , the 3d protein structures of soxc and soxd proteins from s. denitrificans were developed through comparative modeling technique . additionally , soxy and soxz proteins from s. denitrificans were modeled efficiently using fold recognition technique .
the individual modeled proteins were docked amongst themselves to obtain soxyz and soxcd complex , first .
the models were individually optimized and simulated for a stable interactive investigation . furthermore , the simulated soxcd protein of s. denitrificans was picked to dock with the simulated soxyz protein from the same microorganism .
moreover , to the final soxyz complex , thiosulphate was docked via protein - ligand docking procedure .
the detailed interactive residues indulged in the interaction of the quadracomplex ; soxyzcd was further studied .
the stability and conformational alterations in soxyz protein upon soxcd interaction for its recycling were additionally examined .
therefore , this novel present scenario provides a way out to examine the fields of binding , stability , conformational switches , and the arrangement of interactions in soxyzcd protein from s. denitrificans ( epsilon - proteobacteria ) .
this might cover up a massive part for thiosulphate oxidation in dark ocean beds for the chemolithoautotrophy of microorganisms .
in our previous investigation , we analyzed the residual interaction pattern between two essential proteins , soxc and soxd , from sulfurimonas denitrificans for global sulphur oxidation .
for that purpose , we modeled the 3d structures of soxc and soxd protein and further disclosed the residue - level binding strength between two modeled proteins .
the comparative study was examined for the stability of soxc and soxd protein after the duo - protein was simulated . in this present study , first we have opted for an improved process for the preparation of the homology models of the two paramount proteins .
the individual amino acid sequences of soxc and soxd from s. denitrificans were obtained from uniprotkb ( accession numbers q30nu7 and q30nu8 for soxc and soxd , resp . ) .
the sequences were substantiated from ncbi protein sequence database too . for template selection , templates were obtained by swiss - model template searching procedure . with the usage of varied force fields , the explorations of templates for the individual proteins ( soxc and soxd ) were performed .
the same templates were also obtained by the utilization of psi - blast with reference to pdb .
psi - blast helps to bring the distantly related species into consideration and is also acknowledged to be more sensitive when compared to pblast .
the respective templates for building the respective comparative ( homology ) models were from the x - ray crystal structure of(i)paracoccus pantotrophus for soxc ( pdb code : 2xts_a with 87.00% query coverage sharing 43.78% sequence identity , e - value : 2e  114),(ii)paracoccus pantotrophus for soxd ( pdb code : 2xts_b with 51.00% query coverage and sharing 37.82% sequence identity , e - value : 8e  38 ) .
paracoccus pantotrophus for soxc ( pdb code : 2xts_a with 87.00% query coverage sharing 43.78% sequence identity , e - value : 2e  114 ) , paracoccus pantotrophus for soxd ( pdb code : 2xts_b with 51.00% query coverage and sharing 37.82% sequence identity , e - value : 8e  38 ) .
soxc and soxd protein models were constructed utilizing homology modeling with the aid of swiss model software .
the modeled proteins were selected from the gmqe ( global model quality estimation ) score .
furthermore , the next step for preferring the best modeled protein is to examine the sequence identity of the target and template sequence .
the gmqe values are 0.70 and 0.38 for soxc and soxd , respectively . additionally , the qmean score was also supportive with 3.56 and 3.45 for soxc and soxd , respectively .
the individual homology models were superimposed upon their relevant crystal templates to observe the root mean square deviations ( rmsd ) with the aid of pymol .
dalilite server ( analysis tool web services from the embl - ebi , 2013 ) from embl packages also helped to validate the results .
0.068   and 0.079   were the rmsd values with respect to c atoms for soxc and soxd , respectively .
this implies that the models were potentially akin to their crystal template structures . for the purpose of modeling the soxy and soxz proteins ,
fold recognition technique was utilized by operating phyre2 , rather than utilizing traditional homology modeling technique . to overcome the certain specific problems , fold recognition technique ( threading )
fold recognition is a technique for modeling a specific protein that demonstrates the comparable fold with respect to the proteins of well - known tertiary structures .
so , phyre2 was operated here for the execution of fold recognition . in this method , first , the query sequence is used to perform psi - blast to classify remotely related amino acid sequences .
the usage of hidden markov model ( hmm ) was to produce a familiar sequence profile .
template i d from fold library , chain , pdb i d , confidence score , and query coverage were c4uwqk , chain k , 4uwq_k , 100% , and 77% , respectively , for soxy protein .
template i d from fold library , chain , pdb i d , confidence score , and coverage were d1v8ha1 , chain a , 1v8h_a , 100% , and 99% , respectively , for soxz protein . during homology modeling , when the sequence identity between the protein to be modeled and its competent template remains within the range of 30% to 50% , then certain disparities in loop conformations and fallacies might occur in the side - chain packing [ 1517 ] .
so , optimization and remodeling of the loop regions are essential to accomplish the proper orientation of - angles .
modloop via modeller was operated for the establishment of the loop optimization . for the refinement and energy optimization of the single loop optimized protein structures
(  soxy , soxz , soxc , and soxz  ) , an efficient atomic - level algorithm , modrefiner , was followed .
the algorithm executes via a combination of two different force fields : physics - based and knowledge - based . in terms of the proper h - bonding , topology of backbone , and positioning of side chains , the proteins were well optimized by drawing them closer to their native states , where the highest state of interaction exists .
the stereochemical properties of the four homology modeled structures were substantiated by saves server . to analyze the residue - level profile of individual amino acid from the respective protein 3d models ,
examination of ramachandran plots from procheck   web server revealed that most of the residues occupied the most preferred regions with zilch amino acids in the disallowed regions [ 21 , 22 ] .
first , soxyz and soxcd complexes were formed which was followed by the docking of the duo - protein complexes amongst themselves to form the soxyzcd complex . for this purpose , at each case , the entirely automated program cluspro 2.0 was operated [ 23 , 24 ] . for the individual duo - complexes ( soxyz and soxcd ) , both soxy and soxc , being longer , served to be the receptors while the other two modeled proteins were uploaded as ligands .
but for docking the duo - complexes , soxcd served as receptor and soxyz as ligand .
the unstructured residues were confiscated by the execution of the advanced option of cluspro 2.0 . in each of the cases , amongst all the 10 docked structures , the finest cluster sized model was opted for for further investigations .
docking outcomes from zdock and gramm - x were shown to be inclusive of the results [ 25 , 26 ] .
the modeled and docked complexes were energy optimized using charmm force field to diminish the net overall energy .
it was also performed to alter the respective conformations of the dihedral angles to allow them to vary concurrently .
the steepest descent technique preceded by conjugate gradient mechanism was executed through operation of discovery studio software .
the phenomena continued till the modeled structures individually attained a final rms gradient of 0.0001 . to attain a stable conformation for the individual protein complexes with an overall optimization of the net energy , two - step molecular dynamics ( md ) simulation
initially , by template information of modeled structures and based on the principle of md simulation , fragment - guided molecular dynamics ( fg - md ) was operated .
it first disassembles ( fragments ) the provided docked protein structure and further assembles each fragment to accomplish the energy optimized protein modeled structure .
the physical based potentials , except for amber99 , explicitly carried the decoys away from the native protein state .
fg - md was extensively tested with experimentally validated x - ray crystal structures and demonstrated noteworthy potential in the protein structure refinements at the atomic level .
the complete simulation process comprises mainly three steps : ( i ) identification and analysis of the fragment structure , ( ii ) utilization of the simulated annealing technique for md simulation , and ( iii ) selection of the final model .
thus , every time the structure is brought closer to its native state with a higher standard of accuracy .
further , these reconstructed optimized complexes from its fragments were simulated via the operation of the chiron energy minimization tool .
charmm force field is pursued for the optimization of the whole docked complex at a time .
further , a high temperature exchange rate for the complexes is utilized with subjection to simulations by molecular dynamics ( md ) .
thus , finally , not only were the steric clashes purged rapidly but also minimal disconcertion of the backbone of the protein was achieved . with the aid of protein interaction calculator ( pic ) ,
several interactions like disulphide bonds , main - chain side - chain interactions , and hydrophobic , ionic , aromatic - aromatic , aromatic - sulphur , and cation- interactions were examined .
analysis of the interactive docked complexes was also carried out using the analyze binding site tool from discovery studio ( ds ) package .
this is because the ionic interaction predominantly plays a pivotal role in strengthening the interaction . for the examination of the stability of soxyz complex after participation of soxcd complex ( for the recycling of soxyz ) ,
free energy of folding for soxyz complex was assessed and compared ( before and after the interaction of soxcd complex ) .
a descent in the free energy of folding implies better folds and a strong interactive protein complex .
increment in the percentage of residues adopting turn - like and helical conformations led to stronger interaction with a better conformational stability in the protein [ 33 , 34 ] . furthermore , ascent in the residues involved in turn - like regions leads to firmer backbone flexibility which additionally aids the protein to interact with its partner proteins .
moreover , it has been previously investigated with evidence that proteins exhibiting few 310 helices accompanied by the increased percentage for pure -helices are conformationally steadier than the proteins adopting only pure -helices .
so , the secondary structure distribution in the soxyz protein complex was evaluated and compared utilizing dssp method , before and after its interaction with soxcd protein complex .
the stability and strength of interaction of the soxyz complex were thus again perceived through the conformational switches .
calculation , analysis , and comparison of the solvent accessible surface area for the essential interacting amino acids residues from soxyz protein complex were performed .
the specific amino acid residues that take part in the interaction pattern were taken into consideration from the preinteractive complex of soxyz with soxcd and after the tetracomplex formation .
it was observed to rationalize the aforementioned energy estimation outcomes . as the interaction grows stronger , the area available for solvent will keep descending .
the literature suggests that more negative electrostatic potential value apprehends a better interaction of the protein complex of interest and thereby leads to an overall complex stability .
the pictorial representation of the electrostatic mapping upon the protein surface was illustrated by the assistance of pymol ( molecular viewer using python language ) .
so , the optimized complex structure and the preoptimized tetracomplex structure were further utilized for the electrostatic potential mapping upon the surface of the protein .
electrostatic surface potential was evaluated via vacuum electrostatics in kt / e units , with the support of pymol . to examine and rationalize all the estimated outcomes through statistical significance
, p value was evaluated through the paired t - test in each of the cases . in this t - test
, the only assumption usually accepted is that the standard deviations ( sds ) are not approximately akin to both measurements for the individual cases ( before and after soxcd interaction ) .
analysis of the interaction amongst the essential proteins upon the involvement of the sulphur anion , thiosulphate , was also essential .
thiosulphate was obtained from pubchem and was docked as a ligand with the receptor ( final stable soxyzcd protein complex ) .
autodock vina was operated for the purpose . not only the condition and declaration of the search space but also the specification of the binding region was provided .
further again , the molecular surface of the protein - ligand complex was computed with the utilization of segmentation algorithm using patchdock .
further , the best suited and finest model was opted for , depending upon the relevant data for interaction energy , that is , gibbs free energy ( g ) .
the soxyzcd - thiosulphate model having the maximum and best suited (  ) g value is chosen for having the highest strong binding affinity .
the obtained final complex of soxyzcd - thiosulphate was analyzed with the aid of pymol plugins and the interacting residues and binding patterns were thereby examined .
thus , the complete model having four interactive proteins and a ligand ( thiosulphate ) was investigated with the establishment of the interactive residues for the purpose of recycling the thiosulphate - bound soxy protein from the soxyz complex with the aid of soxcd .
the tertiary structure for the homology modeled soxc protein from sulfurimonas denitrificans was analogous to its template from paracoccus pantotrophus ( pdb code : 2xts ; a chain ) .
the 391-residue - long protein comprises mainly 18 sets of -sheets and 10 sets of -helices interspersed with coil regions .
14.3% of the amino acid residues contributed to the respective 10 -helices , whereas 29.2% of the residues were responsible for the adoption of -sheet conformation in the protein .
the modeled soxc structure is well illustrated in figure 1 with red , cyan , and magenta shades representing -sheets , -helices , and interspersing coils , respectively .
the tertiary structure for the homology modeled soxd protein from sulfurimonas denitrificans was also analogous to its template from paracoccus pantotrophus ( pdb code : 2xts ; b chain ) .
only two sets of antiparallel -sheets were observed in the 3d structure with the involvement of 2.8% of the amino acid residues .
soxd structure is illustrated in figure 1 with -helices and coils in yellow and pink shades , respectively .
the functional tertiary structure for the homology modeled soxy protein from sulfurimonas denitrificans was 122 amino acids long .
soxy protein comprises mainly 8 sets of -sheets and 2 sets of -helix conformations at the n - terminal end ( residues : 220 and 5961 ) .
the model is well presented in the interactive picture of figure 1 with -helix , -sheets , and turn - like regions in marine - blue shade .
the functional tertiary structure for the homology modeled soxz protein from sulfurimonas denitrificans was 97 amino acid residues long .
the -sheets occupy 47.82% of the entire structure . the model is well presented in the interactive picture of figure 1 with -sheets and -turns in lime - green shade . the exploration of the interactive residues with the binding patterns of the soxyzcd tetracomplex using pic web server and discovery studio platform from accelrys software displays the four proteins to participate strongly via 9 sets of ionic interactions .
several h - bonding interactions were perceived indulging the main chains and the side chains of the four individual proteins .  
table 1 tabulates the ionic interactions in the optimized and simulated tetracomplex ( soxyzcd ) with a pictorial representation through figure 2 , showing few paramount residues for the firmer interaction .
the polar negatively charged asp2 and glu31 from soxc protein participated strongly in the interaction pattern .
again the polar positively charged lys34 and arg93 were observed to actively participate in the formation of two ionic interactions with glu8 and asp57 from soxz .
lys4 , a polar positively charged residue from soxz , further formed two ionic interactions with glu44 and asp35 from the soxd protein .
this was additionally supported by the interaction of polar negatively charged asp5 ( adjacently positioned to lys4 that formed two ionic interactions ) from soxz with lys50 from soxd for a better stability to the tetracomplex .
furthermore , lys61 , arg67 , and asp37 from soxd participated strongly in building up a strong electrostatic interaction in the soxyzcd tetracomplex .
altogether , these interactions helped to accommodate the partner proteins to interact firmly . to analyze and compare the stability of soxyz complex after interaction with soxcd complex , free energy of folding
vadar 2.0   aided to disclose a decline in the value from 144.71  kcal / mol ( in soxyz complex ) to 577.31  kcal / mol for the soxyzcd complex ( table 2 ) .
therefore , the reduction in the free energy of folding indicates spontaneous interaction with better folds in the soxyz protein complex after interaction with soxcd complex . on comparison and analysis of the conformational switches in the soxyz complex before and after interaction with soxcd complex , it was disclosed that both soxy and soxz underwent abrupt conformational alterations upon the participation of soxcd in recycling soxy protein . for soxy
captivatingly , the presence of few 310 helices in soxy protein ( after soxcd interaction ) accompanied by the increased percentage of residues adopting pure -helices additionally led to a firmer and steadier conformation in the protein ( figure 3 ) .
the increase in percentage of pure -helices with the presence of 310 helices and -sheet conformations in the soxy protein ( after soxcd interaction ) makes soxy exhibit a steadier and more interactive conformation [ 34 , 35 ] .
again , the increment in the -turn conformations in soxz from 21.73% to 24.00% after soxcd interaction ( figure 4 ) affirms even soxz to exhibit a stronger conformation and better stability .
thereby , the entire conformational shift in soxyz protein complex benefits the recycling process of soxy protein by soxcd . the pictorial representation with tabulation of the conformational switches is illustrated in figures 4 and 5 , respectively , for soxy and soxz proteins .
herein , from table 2 , the net area available for solvent accessibility from the interacting residues in soxyz duo - protein complex in soxyzcd complex was perceived to get decreased from 486.92   to 153.92   , after interaction .
there was even a decline in the overall value for net area of solvent accessibility ( considering the summation of the accessible area values from all the residues from soxy and soxz ) from soxyz protein . in kt / e units , electrostatic surface potential displayed the postoptimized and simulated soxyzcd protein complex to exhibit better and firmer interaction ( figure 5 ) .
the blue region symbolizes the electropositive regions , while the red shaded ones are for the electronegative areas . a vivid illustration from figure 5
provides the notion that , from the tetraprotein complex , soxyz complex ( after optimization and simulation ) exhibited an enhanced stability .
the alteration in the electrostatic potential on the soxyz surface was from 48.710 to 56.295 .
the respective outcomes for the alterations in soxyz complex upon soxcd interaction were observed to be statistically significant with an evaluated p value of less than 5% in each of the calculations .
the p values for the alterations in the free energy of folding and net area for solvent accessibility values were 0.021342 and 0.01218 , respectively .
furthermore , again the conformational switches in soxy were observed to be statistically significant with a p value of 0.011327 and 0.012138 for the ascent in the helical conformations and for the decrease in the -sheet conformations ( to allow few residues to adopt 310 helices ) .
p value of 0.021231 for the presence in the residues undergoing 310 helices , which accompanied the pure -helical conformations , also is statistically significant . additionally ,
even the conformational shifts in soxz ( having no helices ) were also perceived to be statistically significant with a p value of 0.013219 and 0.005136 for decrement in -sheets ( to allow few residues to exhibit an increase in the turn - like regions ) and increment in turn - like regions , respectively . along with the strong and stable ionic interactions within the tetraprotein complex ( soxyzcd ) , thiosulphate binding interactions were therefore analyzed . for the interaction between soxy and thiosulphate ( to get recycled ) , mainly his65 , his80 , and asp78 from soxy formed h - bonds as depicted in figure 6 .
strong hydrogen bonds were accomplished by these aforementioned residues from soxy protein with the thiosulphate molecule .
formation of a cavity in soxy by the positively charged environment of soxz helps to accommodate the thiosulphate ligand within the protein strongly .
previously , other molecular level studies involving manifold molecular and computational probes were examined with several species of proteobacteria like allochromatium vinosum , p. pantotrophus , pseudaminobacter salicylatoxidans ( kct001 ) , and many more [ 57 , 41 , 42 ] .
in addition to that , molecular level studies were also explored for the eradication of environmental pollutants like phenol .
literature study documents the species ( mainly proteobacteria ) to perform thiosulphate oxidation , but a computational investigation into the oxidation of sulphur after the formation of the tetracomplex protein ( soxyzcd ) from s. denitrificans ( which is responsible for performing thiosulphate oxidation on a major part on this earth including deep ocean belts ) was nowhere dealt with .
so this novel molecular analysis with the epsilon - proteobacteria sulfurimonas denitrificans was essential to be explored .
figure 1 in supplementary material available online at http://dx.doi.org/10.1155/2016/8683713 . to summarize , the four relevant energy optimized sox proteins , 
soxy , soxz , soxc , and soxd ,  were modeled discretely after the satisfaction of their necessary stereochemical properties .
firstly , the necessary duo - complexes ( soxyz and soxcd ) were formed which was further followed by the generation of the tetracomplex soxyzcd .
the finally docked best modeled tetracomplex of soxyzcd was energy optimized and underwent two - step molecular dynamics simulation further . to examine and explore the residual participation and the binding pattern for the steady simulated soxyzcd complex , the predominant strengthening ionic interactions were disclosed .
nine sets of ionic interactions in the entire soxyzcd complex made the complex a firmly interactive one .
the polar negatively charged asp2 and glu31 from soxc protein participated strongly in the interaction pattern .
several other paramount residues helped in the proper accommodation of soxcd protein complex through pocket formation in soxyz complex .
furthermore , the thiosulphate molecule was observed to interact with soxy protein efficiently with strong h - bonds by the participation of asp78 , his65 , and his80 amino acids .
overall interactions between the four proteins and the thiosulphate molecule ( with soxy ) aided in dragging thiosulphate - bound soxy protein for its recycling from its respective duo - complex ( soxyz ) .
thus , these residual contribution and binding patterns help in the performance of sulphur oxidation mechanism potentially . from the biological point of view
, it was observed that soxz protein interacted strongly with more ionic interactions to participate with soxcd complex .
this joint effort therefore aided in recycling the soxy protein , only to which the thiosulphate molecule binds , for the recycling of the soxy - thiosulphate complex .
evaluations from g values in free energy of folding , net area available for solvent molecules , and electrostatic potential upon soxyz protein surface ( an abrupt increase in the negative value ) displayed that the sole soxyz protein turned out to be more stable and firmer after interaction with soxcd complex .
the conformational transitions in soxyz complex also help to analyze the improved conformational stability and steady interaction of soxyz complex in the presence of soxcd . in soxy ,
the statistically significant increase in the pure -helices accompanied by 3-ten helices and -sheets presents the structure to become a more stable and interactive one .
likewise , soxz ( having no helical regions ) also experienced conformational alterations with increased -turn conformations .
this led to the better backbone flexibility in the proteins permitting them to easily interact with their partner proteins .
this further affirms the plausible recycling of soxy protein from the soxyz complex , by soxcd complex .
so , herein , the deep seas residing microorganism ( s. denitrificans ) was explored for the first time .
the formation of soxyz and soxcd and further soxyzcd complex formation as well as their residual disclosure with importance to their respective positions were disclosed . additionally , the conformational shifts and thermodynamic stability in soxyz protein complex for recycling soxy - thiosulphate for the progression of sulphur oxidation mechanism served to analyze the core basis for rendering a toxic - free environment .
finally , the outcomes of this present study focused mainly on investigating the conformational transitions and stability of soxyz complex upon the interaction of soxcd protein complex from the epsilon - proteobacteria s. denitrificans .
in addition to that , it finally also dealt with the interaction and cooperative participation of thiosulphate for recycling soxy - thiosulphate complex from the soxyz complex , with the aid of soxcd protein complex .
the 3d homology modeled functional protein structures of soxy , soxz , soxc , and soxd were analyzed and demonstrated , followed by optimization .
the relevant necessary interactions in the optimized and simulated soxyzcd protein complex were explored from the respective positions of the pertinent proteins .
predominantly , 9 sets of interactions for this recycling phenomenon for thiosulphate oxidation were served by the ionic interactions from mainly soxy , soxc , soxd , and soxz .
the stability and stronger conformation for soxyz protein upon interaction with soxcd protein complex were explored further with corroboration of statistical significance .
additionally , in the presence of soxcd , soxyz protein complex also interacts firmly and steadily .
further , on the interaction with thiosulphate , four strong h - bonding interactions were observed to strengthen the protein - ligand interaction . from this stable simulated tetracomplex formation
, soxy - thiosulphate gets susceptible to being dragged out for getting recycled by soxcd protein complex .
though there exist former investigations on the cooperation of soxyz and soxcd proteins from other proteobacteria species for thiosulphate oxidation , a comparative in silico examination into the interactive pattern , binding residues , and conformational switches for the recycling of soxy - thiosulphate to oxidize thiosulphate ( from the wadden sea and associated ocean belts ) by soxcd complex of epsilon - proteobacteria ( s. denitrificans ) was not delved into hitherto .
thus , this present novel research presents a rational outlook to explore and analyze the soxcd complex participation with soxyz protein complex from sulfurimonas denitrificans , an epsilon - proteobacterium .
this aids in maintaining a systemic sulphur cycle in the biota , even giving importance to the ocean belts .
the study further instigates the computational analysis of the participation of other proteins from the same microorganism for the help in thiosulphate oxidation .
it is also exciting to investigate any mutational alterations that might affect the interactive pattern among the proteins and thereby might alter the sequential cellular mechanism .
the computational investigation for the sustainability of the global sulphur balance through essential microorganisms can also be provoked in future .
 thiosulphate oxidation ( an essential mechanism ) serves to maintain the global sulphur cycle . 
 earlier experimental and computational studies dealt with environmental thiosulphate oxidation but none dealt with thiosulphate oxidation from deep ocean belts . 
 wet - laboratory experimental research shows that epsilon - proteobacteria sulfurimonas denitrificans possess sox ( sulphur - oxidizing ) operon and perform thiosulphate oxidation efficiently underneath the oceans . from this specific sox operon , 
 soxcd complex recycles the thiosulphate - bound soxy from soxyz complex to balance the environmental sulphur cycle . 
 so , four chief proteins were variedly modeled and relevant simulated interactive structures were obtained . 
 the final simulated tetraprotein complex ( soxyzcd ) from docked soxyz and soxcd complexes was disclosed to be a highly interactive one with predominant ionic residues . 
 free energy of folding , solvent accessibility , and conformational shifts ( coil - like conformation to helices and sheets ) were observed in soxyz complex after interacting with soxcd . 
 the stability of the complex ( soxyzcd ) after simulation was also observed through the electrostatic surface potential values . 
 these evaluations were rationalized via biostatistics . 
 this aids soxcd for recycling soxy along with thiosulphate , which remains interconnected by four h - bonds with soxy . 
 therefore , this novel exploration is endowed with the detailed molecular viewpoint for maintaining the sulphur cycle ( globally ) including the ocean belts . 
1. Introduction
2. Methodology
3. Results and Discussion
4. Conclusion and Future Scope