Ed the best score (Table 1). The site score takes into account parameters such as volume, density, solvent exposure, hydrophilic and hydrophobic nature of residues and donor to acceptor ratio and hence is a comprehensive representation of the possibility of it being a binding site.MDS Analyses of HtrA2 and HtrA2?Peptide ComplexesThe peptides GSAWFSF was chosen for MDS studies as it gave the best XP and E-model scores (Table 2). GQYYFV has been reported to be a well known activator of HtrA2 [19] and hence used as another representative peptide for simulation studies. Moreover, the two peptides were chosen such that one is a designed peptide (GQYYFV) while the other is a part of a well-known HtrA2 binding protein Pea-15 (GSAWFSF). In addition to this, GQYYFV with docking score lesser than GSAWFSF was chosen for MDS analysis to understand whether different affinity for the substrate results in similar movements in the protease. MDS analyses of HtrA2-GQYYFV and HtrA2GSAWFSF complexes demonstrated significant difference in conformation as well as dynamics when compared with unbound HtrA2. Visual inspection of the domain wise movements in peptide bound HtrA2 indicated large fluctuations in hinge/linker region (211?26) as shown in Figures 3a and b. Although theseAllosteric Regulation of HtrAFigure 1. Ribbon model of HtrA2 structures (PDB ID: 1LCY). a. Domain organization of HtrA2 298690-60-5 site protease which comprises N-terminal region (blue), protease domain denoted as PD (yellow) and PDZ domain (red) at C-terminal end. b. Structural alignment of loop refined (light magenta) and unrefined (light green) structures of HtrA2 protein with modelled N-terminal AVPS, loop L3 (residues 142?62) and hinge region (residues 211?25) ?built with Prime (Schrodinger 2011). On refinement, loop L3 and hinge region are reorganized so as to define new regions at the protease and PDZ domain interface. c. Selective binding pocket (SBP) on HtrA2. The energy minimised structure of HtrA2 after modelling flexible regions in the protein is represented as a ribbon model. The binding site designated as SBP, selected on the basis of the MedChemExpress SC1 SiteMap score and residue analyses, is located at the interface of PDZ and protease domain and shown as a multi-coloured mesh. doi:10.1371/journal.pone.0055416.gmovements were larger for GSAWFSF than GQYYFV bound complex, the movement pattern remained similar in these two peptides. Enhanced dynamic movement in the 18325633 former complex could be attributed to the peptide length (heptameric as comparedto hexameric in the latter). Domain wise RMSD analysis of these trajectories provided quantitative output of deviations with respect to time. The trajectory graphs (Figures 3c ) show that along the ?entire sequence, hinge region (211 2226) has RMSD of 2.5 A forTable 1. Putative binding sites in HtrA2 identified by SiteMap tool.Site Number from SiteMap Site 2 Site 1 Site 3 Site 4 SiteResidues present in the site K214, K215, N216,S217,S219, R226, R227, Y228, I229, G230,V231,M232,M233, L234, T235, L236, S237, S239, I240, E243, H256, K262, I264,Q289, N290, A291,E292, Y295,E 296, R299, S302 H65, D69, R71, A89, V90, P92, D95,T324 N48, H65, D169, S173,K191, M232, H261,L265 V192, F251 I33,L34,D35,R36,V73,RSite score 1.092716 0.957142 0.936056 0.807891 0.doi:10.1371/journal.pone.0055416.tAllosteric Regulation of HtrAFigure 2. Representative surface structures of peptide activator docked HtrA2. a. Peptide GSAWFSF -HtrA2 complex and b. Peptide GQYYFV-HtrA2 complex. The former.Ed the best score (Table 1). The site score takes into account parameters such as volume, density, solvent exposure, hydrophilic and hydrophobic nature of residues and donor to acceptor ratio and hence is a comprehensive representation of the possibility of it being a binding site.MDS Analyses of HtrA2 and HtrA2?Peptide ComplexesThe peptides GSAWFSF was chosen for MDS studies as it gave the best XP and E-model scores (Table 2). GQYYFV has been reported to be a well known activator of HtrA2 [19] and hence used as another representative peptide for simulation studies. Moreover, the two peptides were chosen such that one is a designed peptide (GQYYFV) while the other is a part of a well-known HtrA2 binding protein Pea-15 (GSAWFSF). In addition to this, GQYYFV with docking score lesser than GSAWFSF was chosen for MDS analysis to understand whether different affinity for the substrate results in similar movements in the protease. MDS analyses of HtrA2-GQYYFV and HtrA2GSAWFSF complexes demonstrated significant difference in conformation as well as dynamics when compared with unbound HtrA2. Visual inspection of the domain wise movements in peptide bound HtrA2 indicated large fluctuations in hinge/linker region (211?26) as shown in Figures 3a and b. Although theseAllosteric Regulation of HtrAFigure 1. Ribbon model of HtrA2 structures (PDB ID: 1LCY). a. Domain organization of HtrA2 protease which comprises N-terminal region (blue), protease domain denoted as PD (yellow) and PDZ domain (red) at C-terminal end. b. Structural alignment of loop refined (light magenta) and unrefined (light green) structures of HtrA2 protein with modelled N-terminal AVPS, loop L3 (residues 142?62) and hinge region (residues 211?25) ?built with Prime (Schrodinger 2011). On refinement, loop L3 and hinge region are reorganized so as to define new regions at the protease and PDZ domain interface. c. Selective binding pocket (SBP) on HtrA2. The energy minimised structure of HtrA2 after modelling flexible regions in the protein is represented as a ribbon model. The binding site designated as SBP, selected on the basis of the Sitemap score and residue analyses, is located at the interface of PDZ and protease domain and shown as a multi-coloured mesh. doi:10.1371/journal.pone.0055416.gmovements were larger for GSAWFSF than GQYYFV bound complex, the movement pattern remained similar in these two peptides. Enhanced dynamic movement in the 18325633 former complex could be attributed to the peptide length (heptameric as comparedto hexameric in the latter). Domain wise RMSD analysis of these trajectories provided quantitative output of deviations with respect to time. The trajectory graphs (Figures 3c ) show that along the ?entire sequence, hinge region (211 2226) has RMSD of 2.5 A forTable 1. Putative binding sites in HtrA2 identified by SiteMap tool.Site Number from SiteMap Site 2 Site 1 Site 3 Site 4 SiteResidues present in the site K214, K215, N216,S217,S219, R226, R227, Y228, I229, G230,V231,M232,M233, L234, T235, L236, S237, S239, I240, E243, H256, K262, I264,Q289, N290, A291,E292, Y295,E 296, R299, S302 H65, D69, R71, A89, V90, P92, D95,T324 N48, H65, D169, S173,K191, M232, H261,L265 V192, F251 I33,L34,D35,R36,V73,RSite score 1.092716 0.957142 0.936056 0.807891 0.doi:10.1371/journal.pone.0055416.tAllosteric Regulation of HtrAFigure 2. Representative surface structures of peptide activator docked HtrA2. a. Peptide GSAWFSF -HtrA2 complex and b. Peptide GQYYFV-HtrA2 complex. The former.