- Research article
- Open Access
Molecular modeling and molecular dynamics simulations based structural analysis of the SG2NA protein variants
BMC Research Notes volume 7, Article number: 446 (2014)
SG2NA is a member of the striatin sub-family of WD-40 repeat proteins. Striatin family members have been associated with diverse physiological functions. SG2NA has also been shown to have roles in cell cycle progression, signal transduction etc. They have been known to interact with a number of proteins including Caveolin and Calmodulin and also propagate the formation of a multimeric protein unit called striatin-interacting phosphatase and kinase. As a pre-requisite for such interaction ability, these proteins are known to be unstable and primarily disordered in their arrangement. Earlier we had identified that it has multiple isoforms (namely 35, 78, 87 kDa based on its molecular weight) which are generated by alternative splicing. However, detailed structural information of SG2NA is still eluding the researchers.
This study was aimed towards three-dimensional molecular modeling and characterization of SG2NA protein and its isoforms. One structure out of five was selected for each variant having the least value for C score. Out of these, m35 kDa with a C score value of −3.21was the most poorly determined structure in comparison to m78 kDa and m87 kDa variants with C scores of −1.16 and −1.97 respectively. Further evaluation resulted in about 61.6% residues of m35 kDa, 76.6% residues of m78 kDa and 72.1% residues of m87 kDa falling in the favorable regions of Ramchandran Plot. Molecular dynamics simulations were also carried out to obtain biologically relevant structural models and compared with previous atomic coordinates. N-terminal region of all variants was found to be highly disordered.
This study provides first-hand detailed information to understand the structural conformation of SG2NA protein variants (m35 kDa, m78 kDa and m87 kDa). The WD-40 repeat domain was found to constitute antiparallel strands of β-sheets arranged circularly. This study elucidates the crucial structural features of SG2NA proteins which are involved in various protein-protein interactions and also reveals the extent of disorder present in the SG2NA structure crucial for excessive interaction and multimeric protein complexes. The study also potentiates the role of computational approaches for preliminary examination of unknown proteins in the absence of experimental information.
SG2NA (Striatin-3) is a WD-40 repeat protein, originally identified as a nuclear auto-antigen with augmented expression during S to G2 phase of cell cycle . Subsequent studies established its structural and functional relatedness to striatin and zinedin, thereby constituting a three member sub-family of WD-40 superfamily [2–4]. In addition to WD-40 repeats at the carboxylic terminal, they also contain a caveolin binding motif, a coiled-coiled structure and a calmodulin binding domain located in the same order at the amino terminus . As expected from these domain organizations, Striatin family members have been shown to interact with a diverse group of proteins including the regulatory subunit (C subunit) of serine/threonine protein phosphatase PP2A , caveolin , tumor suppressor protein APC  and estrogen receptor . Recently, a large multiprotein complex named Striatin-interacting phosphatase and kinase (STRIPAK) containing Mob3, STRIP1, STRIP2, the cerebral cavernous malformation 3 (CCM3) protein, and Ste20 kinase has been described . In agreement with the existence of diverse interacting partners, Striatin family members have been associated with functions as diverse as non genomic estrogen signalling , cytokinesis , spinogenesis , Golgi polarization  and vesicular trafficking . Members of the striatin family are conserved in metazoan evolution but are absent in plants and prokaryotes. Also, SG2NA is the first among the three to evolve from prokaryotes . Noticeably, Drosophila has only one Striatin homologue i.e., CKA; which acts as a platform for organizing the components of JNK signaling and the transcription factor AP-1, suggesting similar but distinct functions of the three members in higher metazoans [3, 14].
Although SG2NA has been shown to have roles in cell cycle progression, non-genomic estrogen signaling and signal transduction; the structural basis of such diverse functions is not known. Earlier we had shown that it has multiple isoforms generated by alternative splicing . Variants of SG2NA localize in cytosol, mitochondria, plasma membrane, lysosomes and nucleus (unpublished results). It is thus likely that variants of SG2NA perform multiple functions in a context specific manner. Further, due to the presence of WD-40 repeats, it is anticipated that SG2NAs have scaffolding function, but the interactomes are not known.
Since most of the functional properties of the proteins intricately depend upon their structures, the requirement of the 3D structures is indispensable for understanding the biological functions. Three of the mouse SG2NA variants (35 kDa, 78 kDa, and 87 kDa) have been cloned and purified. They consist of 301, 712 and 796 amino acid residues respectively. Since among homologous proteins, structure is more conserved than the sequences, the three-dimensional (3D) structure prediction of SG2NA variants can help us in the predicting their functions. The crystallographic structure of SG2NA has not been reported yet and its active site residues have also not been identified. In the present study we carried out a molecular modeling study of SG2NA proteins using iterative threading assembly refinement (I-TASSER) to obtain their 3D structures. These structures were further energy minimized and corrected by molecular dynamics (MD) simulations to arrive at a better model. The simulations of SG2NA and its isoforms were performed for a long time duration of 10 ns and the stabilized structures so obtained were further analyzed.
Retrieval of sequence data for SG2NA protein variants
The amino acid sequences of SG2NA protein variants (35 kDa, 78 kDa, and 87 kDa) were retrieved from UNIPROT (IDs: B3STQ1 , B2RQS1  and Q9ERG2  respectively). These SG2NA sequences were used to screen similar sequences with known 3D structures against Protein Data Bank (PDB) using NCBI’s BLAST program .
Primary structural analysis of SG2NA protein variants (35 kDa, 78 kDa and 87 kDa)
Expasy’s Prot Param server  was used to study the physiochemical characters of SG2NA protein variants such as theoretical isoelectric point (pI), molecular weight, molecular formula, total number of positive and negative residues, instability index , extinction coefficient , aliphatic index  and grand average hydropathy (GRAVY) [21, 22]. CYS_REC tool has been applied for the anticipation of disulphide bond (S-S) pattern. The instability index provides the estimate of a protein’s stability in vitro and is based on the difference in occurrence of certain dipeptides in the unstable and stable proteins and calculated by assigning weights of unstability to 400 different dipeptides (DIWV) and using the equation [Eq 1]
where i = 1, L is the length of the sequence and DIWV(x[i] x[i + 1]) stands for the instability weight value for the dipeptide starting from positon i. The extinction coefficient predicts the amount of light a protein can absorb at a certain wavelength. It can be calculated for any protein given the molar extinction coefficients of Trp, Tyr and Cys are provided, using the following equation [Eq 2].
where (for proteins in water measured at 280 nm): Ext(Tyr) = 1490; Ext(Trp) = 5500; Ext(Cys) =125. The aliphatic index can be regarded as a positive factor in increased thermostability of globular proteins by measuring the relative volume occupied by aliphatic side chains (alanine, valine, leucine and isoleucine). It is calculated using the formula given in Eq 3.
where X(Ala), X(Val), X(Ile), and X(Leu) are mole percent (100 X mole fraction) of alanine, valine, isoleucine, and leucine. The coefficients a and b are the relative volume of valine side chain (a = 2.9) and of Leu/Ile side chains (b = 3.9) to the side chain of alanine.
The GRAVY score is the ratio of the sum of hydropathy values of all amino acids and the number of residues in the sequence.
Secondary structural analysis of SG2NA protein variants
Secondary structural analysis was done using PSIPRED view web server . It is a new and highly accurate secondary structure prediction method. It was employed to enumerate the secondary structural features of SG2NA proteins Two feed-forward neural networks integrated by PSIPRED, perform an analysis on output received from PSI-BLAST (position specific Iterated BLAST). A putative secondary structure is obtained with respect to each residue associated with a confidence value for the prediction.
Ab Initio modeling, Structure refinement and identification of functional site
Pertaining to lack of suitable structural template for SG2NA and its isoform, we resorted to ab initio modeling of SG2NA using I-TASSER [24, 25], a web based, full chain protein prediction server. I-TASSER follows a comprehensive protocol for protein modeling. It first scans a representative PDB structure library to search for possible folds by threading the target sequence using Profile-profile threading alignment (PPA). The threading aligned regions are excised and threading unaligned regions are modeled using ab initio modeling before reassembling the complete model together. The knowledge of a protein’s conformational space is crucial and to search that it uses replica exchange Monte-Carlo simulations. Models having conformational similarities are clustered by SPICKER. By averaging the coordinates of structures in a cluster, a cluster centroid is obtained. To rule out steric clashes and refine the models further, fragment assembly simulation is done starting from the cluster centroid obtained from the previous simulation. After retrieval of spatial constraints from these clusters, TM-align is used to search for PDB structures to guide second simulation. After clustering the structure decoys, the lowest energy structure in each cluster is chosen. Simulations are run for full chain as well as the domains. Finally, the complete model is obtained by docking the model of domains together by Metropolis Monte Carlo simulation. This docked complex is compared to the I-TASSER full chain model and arrives at the closest one with least steric clashes.
The C score is defined by
where M is the multiplicity of structures in the SPICKER cluster; Mtot is the total number of the I-TASSER structure decoys used in the clustering; RMSD is the average RMSD of the decoys to the cluster centroid; Z(i) is the highest Z score (the energy to mean in the unit of standard deviation) of the templates by the ith PPA threading program and Z0(i) is a program-specified Z-score cutoff for distinguishing between good and bad templates, i.e. Z0(1) = 7.0, Z0(2) = 8.5, Z0(3) = 8.0, Z0(4) = 10.5.
The TM score is defined to assess the topological similarity of two proteins
where d i is the distance of the ith pair of residues between two structures after an optimal superposition, and L is the protein length. TM score lies in between 0 to 1 with higher values indicating better models. Statistically, a TM-score ≤ 0.17 corresponds to a similarity between two randomly selected structures from the PDB library while a TM-score > 0.5 corresponds approximately to two structures of the similar topology. Structure based protein functional site of SG2NA and of three isoforms were predicted by Q-site Finder .
Structure Validation of the models obtained by I-TASSER
I-TASSER server provided five structures for each of the protein. Then PSVS web server was used to assess the qualities of the models. The models were evaluated by Verifiy3D, Prosall, PROCHECK  and Molprobity programs to check the correctness of the overall fold/structure, errors over localized regions and stereochemical parameters such as bond lengths and angles. Structural validation of target proteins model were done to determine stereochemical aspects along with main chain and side chain parameters with comprehensive analysis. This shows the various residues of SG2NA falling under allowed, favoured and disallowed regions of a Ramchandran plot obtained from PROCHECK and Richardson Lab’s MolProbity . The post-simulated structures were also assessed for their quality using the same server.
Molecular Dynamics simulations of SG2NA variants (35, 78 and 87)
The AMBER v.10 package  was used to prepare the proteins for the Molecular Dynamics (MD) simulations. The coordinates of protein atoms were optimized in the protein preparation wizard of Schrodinger’s Maestro suite where hydrogens were added; water molecules and UDP were removed, and the complex structure was minimized using the OPLS2001 force field. The 1-μm simulation used the CHARM27 forcefield , and the simple point charge model for water. The CHARM27 forcefield was applied to the system using the VIPARR utility. The default Desmond relaxation was performed before simulation, and molecular dynamics were run at constant temperature (300 K) and pressure (1 bar). The simulation was performed by using the program Desmond, compiled by SB Grid on an optimized 64-node Linux-based InfiniBand cluster.
Comparative split-domain modeling using QUARK and MUSTER algorithms
The SG2NA isoforms m78 and m87 were split into three domains namely, caveolin-binding and putative coiled-coil (res. 1–140), calmodulin binding domain (res. 141–300) and the WD 40 repeat domain (301–712 residues in m78 and 301–797 residues in case of m87 respectively) owing to long amino acid sequences. As the first two domains were small in size, they were modeled using ab initio server QUARK  and the last one was submitted to template based server MUSTER . QUARK is a template-free prediction tool which accepts sequences of less than 200 residues. It splits the sequences into smaller fragments up to of 1–20 residues and obtains a number of fragment structures from unrelated experiments for each position. Using replica-exchange Monte Carlo simulations the fragments are assembled based on knowledge based force fields to yield complete structural models. On the other hand, MUSTER is template-based modeling software. It stands for MUlti-Source ThreadER and is an extension to the simple threading algorithm named, PPA. The algorithm uses dynamic programming methods and includes single-body features like (1) sequence profiles (2) secondary structure prediction (3) depth-dependent structure profiles (4) solvent accessibility (5) backbone dihedral torsion angles and (6) hydrophobic scoring matrix. After selection of appropriate templates the models are built using the inbuilt Modeller algorithm.
Disorder prediction of SG2NA proteins
The disorder prediction for SG2NA proteins was done using a web-based server named IUPred  which is based on the assumption that relatively stable proteins can form inter-residue contacts that stabilize the entropy loss during folding while the intrinsically disordered proteins (IDPs) of which Striatin is also a part, have certain sequence regions that are unable to form sufficient inter-contacts. In this method, the pairwise interaction energies of a few selected globular proteins can be estimated using a quadratic expression set for each residue in terms of whether it contributes to the ordered or disordered regions. This contribution depends upon its own chemical composition as well as its sequential interaction partners as a whole. It resulted in a disorder prediction graph along with the disorder probability for each residue position.
Results & discussion
SG2NAs belong to the Striatin family constituted of two other members i.e., Striatin and Zinedin. Although the expression of mammalian Striatin is largely restricted to the central nervous system (CNS), that of SG2NA is ubiquitous [6, 13]. All three members share a number of protein-protein interaction domains like WD-repeat (except 35 kDa SG2NA), Ca2+-calmodulin binding and the putative coiled-coil structure, suggesting that they are involved in Ca2+ dependent signal transduction. Thus, the knowledge of their structure is imperative for a better understanding of their functions. The three variants of the SG2NA viz., 35 kDa, 78 kDa and 87 kDa differ in their size and sequences. In agreement, structures with different topology were observed.
Primary structural analysis of SG2NA protein
The primary structural features of SG2NAs (35 kDa, 78 kDa, and 87 kDa) are described and compared in Table 1. The calculated isoelectric points (pI) for all three isoforms are 5.90, 5.06 and 5.12 respectively, suggesting the presence of more negatively charged residues. Amidst the three isoforms under study, the 87 kDa possesses higher extinction coefficient, indicating the presence of Cys, Trp and Tyr in abundance. Higher number of these residues aid to the quantitative study of protein-protein and protein-ligand interactions in solution. The disulphide (S-S) bonding pattern depicts three cysteine residues in 78 kDa variant at positions 544, 546 and 670 and two cysteine residues in 87 kDa variant at positions 556 and 613. The 35 kDa variant appears to have no disulphide bond. As the presence of a disulphide bond increases the enthalpy of the folded state by stabilizing local interactions, these results suggest that SG2NA has high enthalpy at the folded state. Although, at a wide range of temperatures SG2NAs were found to be stable (due to the higher aliphatic index [AI]), an instability index of more than 40 suggests instability and thus predicts that the variants of SG2NA are thermally unstable . The value of aliphatic index for SG2NAs (35 kDa, 78 kDa, and 87 kDa) were 69.44, 78.40 and 77.60 respectively while that of instability index were 60.06, 46.0 and 45.15. GRAVY (Grand average hydropathy) values of SG2NAs were −0.869, −0.470 and −0.532, suggesting a hydrophilicity pattern with better interaction with water. Also, the GRAVY value of 35 kDa (−0.869) was lowest among all three variants, indicating a better solubility as compared to 78 and 87 kDas.
Secondary structure prediction of SG2NA proteins
Currently, most of the algorithms for the anticipation of secondary structure of a protein are based on machine learning techniques. Among the available methods to date, PSIPRED-view achieved highest level of accuracy, an average Q3 score of 76.5% . It was used for the prediction and analysis of secondary structure of SG2NA proteins and is shown in Table 2. Only α helices were found in the predicted structure of 35 kDa SG2NA, It reveals the unfavored structural property in non-polar solvent. The amount of alpha helix, beta sheets, and coils in 78 SG2NA are 13.5%, 23.7% and 62.78% respectively while that of in 87 kDa are 13.94%, 20.22%, and 65.82% respectively. These results can help in experimental verification of a predicted folding motif as it may be gained by measurements of protein secondary structural elements of which the motif is composed. Another insight into the secondary structure of these isoforms was provided by the PSVS validation suite and has been summarized in Table 2 along with PSI-PRED results to compare the accuracy of the servers and reach a consensus. According to PSVS, 35 kDa SG2NA is constituted of only α-helices which amount to 26.24% of the complete protein. In case of the 78 kDa isoform, amount of β-sheets exceeded α-helices standing at 25.28% and 10.81% respectively. For the 87 kDa isoform, α-helices take up 7.66% and β-sheets take up 28.39% of the complete structure.
Prediction of three dimensional (3D) structures of 35 kDa, 78 kDa and 87 kDa SG2NA
Three dimensional structure of SG2NA proteins were predicted and compared. The comparative protein structure analysis for SG2NA variants is still untouched and unavailable. The tertiary structure prediction was performed by I-TASSER server using the best aligned template obtained by searching against Protein Data Bank database. The search for optimal templates is inbuilt in the I-TASSER algorithm. Templates that were chosen are listed: 1wp1A, the OprM lipoprotein of Pseudomonas aeruginosa for 35 kDa and 2ymuA, WD-40 repeat protein of Nostoc punctiforme for 78 kDa and 87 kDa. The template was selected to model and analyze 3D structure because a high level of sequence identity should guarantee a more accurate alignment between the target sequence and template whose structure is known. The final template chosen is used to model the initial structure. This final template is chosen after aligning the query sequence with a number of other protein sequences whose structures are already known. Out of five generated models of the target sequence, the best ones have been chosen owing to the criteria of good alignment with the chosen template and are measured by C-Score, TM score and RMSD values as provided in Table 3. The developed 3D model of SG2NA proteins were deposited to the PMDB database. The diagram showing the post simulated topological structures of the three variants has been provided in Figure 1.
The C score refers to the confidence in the quality of the predicted structure and is based on threading template alignments and convergence parameters involved in simulation. It lies between −5 to 2, where a higher C score signifies absolute quality of the model and vice versa. The isoform 35 kDa with a C score of −3.21 is a poor model in comparison to the 78 kDa and 87 kDa models having a score of −1.16 and −1.97 respectively. To measure the structural similarity between two structures consequential in determining the quality of the predicted model, TM score and RMSD are calculated. But, owing to the error-prone state of the RMSD calculations, we rely on the TM score in which the smaller distances between structures are weighted stronger than the larger distances thus making it free of local errors. A TM value of more than 0.5 ensures a model of correct topology and a value of less than 0.17 refers to random similarity and hence should be discarded. The isoform 35 kDa with a TM score of 0.36 ± 0.12 is a poor model in comparison to the 78 kDa and 87 kDa models having a score of 0.57 ± 0.0.15 and 0.48 ± 0.15 respectively.
As discussed before, I-TASSER works by excising continuous fragments from the template alignments and reassembles those using Monte Carlo simulations. These low temperature replica structures also called decoys are clustered using SPICKER. The number of such decoys in a unit space of a SPICKER cluster is known as the cluster density, a high value representing how often a structure occurs in the simulation trajectory, signifies good quality of a model. The total number of structure decoys employed to model a structure has been provided in Table 3. The isoform 78 kDa with the highest value of cluster density at 0.0943 amongst its counterparts, 35 kDa and 87 kDa with cluster density of 0.0498 and 0.0404 respectively, signifies a model of better quality.
Structure validation of the modeled SG2NA isoforms using PSVS
The quality of the predicted structures of SG2NA were further assessed and confirmed by Verifiy3-D, Prosall, PROCHECK, and Molprobity which are a part of the PSVS server. The scores (from −1 to +1) were added and plotted for individual residues. The residues falling in the area where the orange line crosses 0.0 have low prediction accuracy and less stable conformation whereas, most of the residues fall above 0.15-0.4 so we can say that the model is of good quality. The stereochemical quality and accuracy of the predicted model of SG2NA were evaluated after the refinement process using Ramachandran Map calculation with the PROCHECK program. The Ramachandran plot has been shown a tight clustering of phi ~ −50 and psi ~ −50. This plot is a way to visualize dihedral angles φ against ψ of amino acid residues in protein structure. It shows the possible conformations of φ and ψ angles for a polypeptide. In Ramachandran plot, the white areas correspond to conformations where atoms in the polypeptide come closer than the sum of their Vander Waals radii. These regions are sterically disallowed for all amino acids except glycine, which is unique in that it lacks a side chain. The red regions correspond to the allowed regions namely the alpha-helical and beta-sheet conformations where there are no steric clashes. The yellow areas show the partially allowed regions of left handed helix wherein the atoms are allowed to come a little closer together. Glycine is restrained by triangles and other residues are represented by squares.
The Ramchandran plot analysis obtained by PROCHECK and Richardson Lab’s MolProbity have been summarized in Table 4 and the plots are provided in Figure 2A. The residues falling in most favored regions of the Ramchandran plot obtained by PROCHECK amounted to 61.6%, 76.6% and 72.1% for m35kDa, m78kDa and m87kDa respectively while for the plot obtained by MolProbity it amounted to 69.9%, 80.7% and 78.5% respectively. Similarly, residues falling in disallowed regions amounted to 2.5%, 2.7% and 2.4% in PROCHECK while it stood at 9.7%, 7.6% and 8.7% obtained by MolProbity for the three SG2NA variants respectively. This provides an insight into the correctness of the modeled structures, m78kDa being the best modeled structure both in terms of PROCHECK and MolProbity. A low number of residues falling in the favored regions signify errors in the modeled structure of m35kDa protein.
Molecular dynamics simulations of SG2NA protein structures
Molecular dynamics simulations were performed on Desmond utility provided in the Schrodinger’s Maestro suite to ascertain the absolute quality of the models obtained to refine them by subjecting to 10000 picosecond (ps) long simulations. The final structures were obtained by taking the average of structures pertaining to that part of the trajectory where it becomes stable. A comparative trajectory graph has been provided in Figure 3 for the three SG2NA variants: m35kda, m78kDa and m87kDa. These post-simulated structures were again subjected to validation by the PSVS suite and the results have been summarized in Table 4 for comparison with the pre-simulated structures. As can be gathered from the results, the number of residues falling under most favored regions and allowed regions in the Ramchandran plot (both from PROCHECK and MolProbity) have significantly increased for all the three variants especially for m35kDa. These results are obtained due to the removal of unfavorable contacts from the protein model during MD simulation. This may, however, affect the global quality of protein; a limitation that is faced due to unavailability of more accurate modeling techniques. Similarly, residues falling in disallowed regions of the Ramchandran plot decreased by a good margin signifying an overall improvement in the quality of these models. Post MD Ramchandran Plot analysis has been provided in Figure 2B. Figure 4 summarizes the extent of variation in the topology of these models occurring after molecular dynamics simulations were performed.
Comparison of structures obtained after split-domain modeling of m78 and m87 isoforms
As the SG2NA sequence contains a number of conserved domains also mentioned before, it was split into three domain regions for m78 and m87 isoforms and their structures were predicted separately using QUARK (for sequences less than 200 residues) and MUSTER as I-TASSER does not accept very short stretches of protein sequences. The three split domain regions comprised of (1) Caveolin + putative coiled-coil domain (1–140 residues) (2) Calmodulin binding domain (141–300 residues) (3) WD-40 domain. The amino acid sequences from isoforms m78 and m87 are same for caveolin and calmodulin binding domains and were submitted to QUARK. The WD-40 region is different with 301–712 residues in the m78 isoform and 301–797 residues in the m87 isoform and thus was submitted in MUSTER.
The separately modeled region containing Caveolin binding domain and putative coiled-coil domain (from m78 and m87 isoforms) is significantly different in the topology when compared to the complete models. As represented in Figure 5, it consists of a continuous α-helix hairpin containing caveolin binding domain (in blue) and coiled coil region (in orange). The region from 74 – 100 residues is more or less similar in all the three structures especially the caveolin binding region and thus overlaps significantly with the 3 isoforms. However, the coiled coil region predominately consists of α-helix and is completely different than all the three isoforms. The calmodulin binding region shows a completely different structure when modeled exclusively consisting majorly of helical regions whereas only loop regions are present in the complete m35 model and β-sheets are present in the same region of m78 and m87 isoforms. The possible reason for this observation can be the various stereochemical constraints under which a complete protein is modeled. In the absence of other protein domains, the algorithm just takes the single fragment into perspective and puts stereochemical constraints occurring within that part. This may not result in credible results as the true geometry is determined along with the complete protein. Finally, the WD-40 repeat domain, when compared with m78 and m87 isoforms, was found to be more or less similar exuding the same β-propeller topology except the residues from 301–500 which is topologically different even when consisting of β-sheets predominantly (Figure 6). The WD-40 repeat domain was modeled using the template based server MUSTER. The top template selected for modeling m78 isoform was the crystal structure of human TBL1XR1 WD-40 repeat protein [PDB ID: 4LG9] with a Z score of 19.63 and the top scoring template for modeling m87 isoform was highly repetitive propeller structure of WD-40 repeat protein from Nostoc punctiforme [PDB ID: 2YMU] with a Z score of 17.846. The Z score is given by
Where R’score is the normalized score obtained by raw alignment score (Rscore) divided by partial alignment length (Lpartial i.e. excluding query ending gap) or Rscore divided by full alignment length (Lfull i.e. including query and template ending gaps). If the template chosen using Lpartial is more similar to the query, it is chosen otherwise the template chosen by Lfull is chosen. If the Z-score is more than 7.5, the chosen template is said to be a good template. The SG2NA proteins are known by their WD-40 repeat propeller like domain which plays a key role in its functions.
Analysis of disorder prediction for SG2NA protein variants (35, 78 and 87 kDa)
As the SG2NA protein belongs to the three-member superfamily of WD-40 repeat containing proteins supposedly having scaffolding functions and is also known to be involved in a number of protein interactions and multimeric assembly formations. In order to exhibit such functionality with high speed of interaction but less binding strength, these proteins are ought to be disordered in their basic arrangement. Such proteins are known to have crucial functional roles in signaling pathways, cell-cycle regulation and transcription. Similarly, SG2NA is known to be involved in estrogen signaling, vesicular trafficking, spinogenesis etc. The disorder prediction results for the three protein variants have been graphically depicted in Figure 7. A disorder prediction probability value of more than 0.5 signifies disordered state of a particular amino acid residue. SG2NA variant 35 kDa, constituting just the N-terminal region was predicted to be highly disordered for most of the residue positions. This region contains the Caveolin and Calmodulin-binding domains and the putative coiled-coil structure which are expected to be highly unstable. Also, the 3D predicted model was also found to be very unstable during MD simulations further validating our disorder prediction findings. SG2NA variant 78 kDa and 87 kDa were also found to be highly disordered at the N-terminal region but became ordered for further residue positions up to the carboxylic end.
In this study, we have determined the structural specifications of the SG2NA protein that constitutes the WD-40 repeat (present at the carboxylic terminal) superfamily. The three SG2NA protein variants (35 kDa, 78kDA and 87 kDa) have been modeled using the I-TASSER server and validated using the PSVS server. The obtained models provide crucial insight into the topological specifications that are responsible for their functionality. Striatin-3 (35 kDa) is the N-terminal sequence region and thus only contains the caveolin and calmodulin-binding domains. It was found to be the most hydrophilic amongst all during initial analysis. The C-score that signifies overall quality of the model was found to be very low for this variant indicating poor model prediction. The other two variants Striatin-3 (78 and 87 kDa) have slight difference in the sequence size and thus in the overall structure. Along with a small helical caveolin-binding domain (residues 71–79) and a loop region signifying calmodulin-binding domain (residues 166–183), it contains β-sheet hairpins as the WD-40 repeat region (residues 353–753) arranged in a circular fashion in both the variants. The presence of WD-40 repeat domain in Striatin-3 is touted to be responsible for scaffolding function. A putative coiled-coil region (residues 86–132) is also present in all the three variants. The three protein models were subjected to MD simulations which resulted in relatively different folds and patterns. The SG2NA variant 35 kDa was the most variable structure amongst all during simulations. As crucial as the MD simulation can be for bringing accuracy in protein modeling process, it is accompanied with its own limitations. The present inability of computational methods for modeling native-like folds for larger proteins poses a limitation for this study as well. However, as we try to mimic biological protein folding by carrying out long MD simulations and then report the average structure of energetically stable frames, it is almost crucial to carry out MD simulation procedure. The final structures were also validated based on Ramchandran analysis with most of the residues falling in sterically allowed regions. The Ramchandran analysis scores improve significantly for MD simulated proteins as a result of removal of all unfavorable contacts in the structure. This, however, imposes another limitation to in silico analysis in terms of compromising the global quality of the protein model. Making as much use of available techniques, we attempted to identify the protein conformations closest to biologically relevant conformations using the available methods. Finally, owing to the multimeric organization and multi-protein interaction property of SG2NA proteins, the extent of presence of disordered regions was estimated. These results establish and propagate the structural knowledge along with certain physicochemical facts regarding SG2NA proteins that were unknown till now.
SG2NA proteins are known to play vital role as diverse as cell signaling, cytokinesis, cell cycle progression, spinogenesis etc. SG2NA, owing to its disordered nature, is difficult to isolate, thus posing strong hindrances in experimental structural analysis. In view of these, computational approaches are a viable option to obtain 3D structure and to understand further the functional relationships. The analysis of the sequences revealed that 78 and 87 kDa proteins are mixture of alpha, beta and coil but containing more than 50% of random coils which make them highly flexible and nothing can be concluded about the 35 kDa protein due to low confidence score of the structure. On splitting the two isoforms m78 and m87 into three domains, it could be deduced that the caveolin binding and calmodulin binding domains are highly disordered and their structures are difficult to predict using in silico methods as no similar structures were found. The characteristic WD-40 repeat domain has been modeled in other organisms and bears sequence similarity with the mouse SG2NA sequence, thus, a relatively reliable structure could be obtained for this domain. Structure prediction and functional analysis of SG2NA will give an insight to the location of these proteins along with site of utilization for protein-protein interaction. This study paves way for further attempts to elucidate the experimental structures and thus validate the accuracy of the predictions. The structural knowledge of these proteins is important for deducing its mechanism of action while interacting with substrates like caveolin and calmodulin. The structural information is a pre-requisite for understanding and assigning functional role to the protein.
Muro Y, Chan EK, Landberg G, Tan EM: A cell-cycle nuclear autoantigen containing WD-40 motifs expressed mainly in S and G2 phase cells. Biochem Biophys Res Commun. 1995, 207 (3): 1029-1037. 10.1006/bbrc.1995.1288.
Castets F, Rakitina T, Gaillard S, Moqrich A, Mattei MG, Monneron A: Zinedin, SG2NA, and striatin are calmodulin-binding, WD repeat proteins principally expressed in the brain. J Biol Chem. 2000, 275 (26): 19970-19977. 10.1074/jbc.M909782199.
Benoist M, Gaillard S, Castets F: The striatin family: a new signaling platform in dendritic spines. J Physiol Paris. 2006, 99 (2–3): 146-153.
Sanghamitra M, Talukder I, Singarapu N, Sindhu KV, Kateriya S, Goswami SK: WD-40 repeat protein SG2NA has multiple splice variants with tissue restricted and growth responsive properties. Gene. 2008, 420 (1): 48-56. 10.1016/j.gene.2008.04.016.
Gordon J, Hwang J, Carrier KJ, Jones CA, Kern QL, Moreno CS, Karas RH, Pallas DC: Protein phosphatase 2a (PP2A) binds within the oligomerization domain of striatin and regulates the phosphorylation and activation of the mammalian Ste20-Like kinase Mst3. BMC Biochem. 2011, 12: 54-10.1186/1471-2091-12-54.
Gaillard S, Bartoli M, Castets F, Monneron A: Striatin, a calmodulin-dependent scaffolding protein, directly binds caveolin-1. FEBS Lett. 2001, 508 (1): 49-52. 10.1016/S0014-5793(01)03020-4.
Breitman M, Zilberberg A, Caspi M, Rosin-Arbesfeld R: The armadillo repeat domain of the APC tumor suppressor protein interacts with Striatin family members. Biochim Biophys Acta. 2008, 1783 (10): 1792-1802. 10.1016/j.bbamcr.2008.04.017.
Bernelot Moens SJ, Schnitzler GR, Nickerson M, Guo H, Ueda K, Lu Q, Aronovitz MJ, Nickerson H, Baur WE, Hansen U, Iyer LK, Karas RH: Rapid estrogen receptor signaling is essential for the protective effects of estrogen against vascular injury. Circulation. 2012, 126 (16): 1993-2004. 10.1161/CIRCULATIONAHA.112.124529.
Hyodo T, Ito S, Hasegawa H, Asano E, Maeda M, Urano T, Takahashi M, Hamaguchi M, Senga T: Misshapen-like kinase 1 (MINK1) is a novel component of striatin-interacting phosphatase and kinase (STRIPAK) and is required for the completion of cytokinesis. J Biol Chem. 2012, 287 (30): 25019-25029. 10.1074/jbc.M112.372342.
Chen YK, Chen CY, Hu HT, Hsueh YP: CTTNBP2, but not CTTNBP2NL, regulates dendritic spinogenesis and synaptic distribution of the striatin-PP2A complex. Mol Biol Cell. 2012, 23 (22): 4383-4392. 10.1091/mbc.E12-05-0365.
Kean MJ, Ceccarelli DF, Goudreault M, Sanches M, Tate S, Larsen B, Gibson LC, Derry WB, Scott IC, Pelletier L, Baillie GS, Sicheri F, Gingras AC: Structure-function analysis of core STRIPAK Proteins: a signaling complex implicated in Golgi polarization. J Biol Chem. 2011, 286 (28): 25065-25075. 10.1074/jbc.M110.214486.
Bailly YJ, Castets F: Phocein: A potential actor in vesicular trafficking at Purkinje cell dendritic spines. Cerebellum. 2007, 6 (4): 1-9.
Castets F, Bartoli M, Barnier JV, Baillat G, Salin P, Moqrich A, Bourgeois JP, Denizot F, Rougon G, Calothy G, Monneron A: A novel calmodulin-binding protein, belonging to the WD-repeat family, is localized in dendrites of a subset of CNS neurons. J Cell Biol. 1996, 134 (4): 1051-1062. 10.1083/jcb.134.4.1051.
Chen HW, Marinissen MJ, Oh SW, Chen X, Melnick M, Perrimon N, Gutkind JS, Hou SX: CKA, a novel multidomain protein, regulates the JUN N-terminal kinase signal transduction pathway in Drosophila. Mol Cell Biol. 2002, 22 (6): 1792-1803. 10.1128/MCB.22.6.1792-1803.2002.
Church DM, Goodstadt L, Hillier LW, Zody MC, Goldstein S, She X, Bult CJ, Agarwala R, Cherry JL, DiCuccio M, Hlavina W, Kapustin Y, Meric P, Maglott D, Birtle Z, Marques AC, Graves T, Zhou S, Teague B, Potamousis K, Churas C, Place M, Herschleb J, Runnheim R, Forrest D, Amos-Landgraf J, Schwartz DC, Cheng Z, Lindblad-Toh K, Eichler EE, et al: Lineage-specific biology revealed by a finished genome assembly of the mouse. PLoS Biol. 2009, 7 (5): e1000112-10.1371/journal.pbio.1000112.
Gasteiger E, Hoogland C, Gattiker A, Duvaud S, Wilkins MR, Appel RD, Bairoch A: The Proteomics Protocols Handbook. Edited by: Walker JM. 2005, Totowa, NJ: The Proteomics Protocols Handbook, Humana Press, 571-607.
Guruprasad K, Reddy BV, Pandit MW: Correlation between stability of a protein and its dipeptide composition: a novel approach for predicting in vivo stability of a protein from its primary sequence. Protein Eng. 1990, 4 (2): 155-161. 10.1093/protein/4.2.155.
Ikai A: Thermostability and aliphatic index of globular proteins. J Biochem. 1980, 88 (6): 1895-1898.
Jeong JH, Kim SW, Yoon SM, Park JK, Lee JS: Characterization of the conserved region of the mxaF gene that encodes the large subunit of methanol dehydrogenase from a marine methylotrophic bacterium. Mol Cells. 2002, 13 (3): 369-376.
Pelton JT, McLean LR: Spectroscopic methods for analysis of protein secondary structure. Anal Biochem. 2000, 277 (2): 167-176. 10.1006/abio.1999.4320.
Kyte J, Doolittle RF: A simple method for displaying the hydropathic character of a protein. J Mol Biol. 1982, 157 (1): 105-132. 10.1016/0022-2836(82)90515-0.
McGuffin LJ, Bryson K, Jones DT: The PSIPRED protein structure prediction server. Bioinform Appl Note. 1999, 16 (4): 404-405.
Roy A, Kucukural A, Zhang Y: I-TASSER: a unified platform for automated protein structure and function prediction. Nat Protoc. 2010, 5 (4): 725-738. 10.1038/nprot.2010.5.
Zhang Y: I-TASSER server for protein 3D structure prediction. BMC Bioinformatics. 2008, 9: 40-10.1186/1471-2105-9-40.
Laurie ATR, Jackson RM: Q-SiteFinder: an energy-based method for the prediction of protein–ligand binding sites. Bioinformatics. 2005, 21 (9): 1908-1916. 10.1093/bioinformatics/bti315.
Vaguine AA, Richelle J, Wodak SJ: SFCHECK: a unified set of procedures for evaluating the quality of macromolecular structure-factor data and their agreement with the atomic model. Acta Crystallogr D Biol Crystallogr. 1999, 55 (Pt 1): 191-205.
Case DA, Darden T, Cheatham TE, Simmerling CL, Wang J, Duke RE, Luo R, Crowley MWR, Zhang W, Merz KM, Wang B, Hayik S, Roitberg A, Seabra G, Kolossvary I, Wong KF, Paesani F, Vanicek J, Wu X, Brozell SR, Steinbrecher T, Gohlke H, Yang L, Tan C, Mongan J, Hornak V, Cui G, Mathews DH, Seetin MG, et al: AMBER 10. 2008, San Francisco: University of California
MacKerell AD, Banavali N, Foloppe N: Development and current status of the CHARMM force field for nucleic acids. Biopolymers. 2000, 56 (4): 257-265. 10.1002/1097-0282(2000)56:4<257::AID-BIP10029>3.0.CO;2-W.
Xu D, Zhang Y: Ab initio protein structure assembly using continuous structure fragments and optimized knowledge-based force field. Proteins. 2012, 80 (7): 1715-1735.
Wu S, Zhang Y: MUSTER: Improving protein sequence profile-profile alignments by using multiple sources of structure information. Proteins. 2008, 72 (2): 547-556. 10.1002/prot.21945.
Dosztányi Z, Csizmok V, Tompa P, Simon I: IUPred: web server for the prediction of intrinsically unstructured regions of proteins based on estimated energy content. Bioinformatics. 2005, 21 (16): 3433-3434. 10.1093/bioinformatics/bti541.
Wedemeyer WJ, Welker E, Narayan M, Scheraga HA: Disulfide bonds and protein folding. Biochemistry. 2000, 39 (23): 7032-10.1021/bi005111p.
Jones DT: Protein secondary structure prediction based on position-specific scoring matrices. J Mol Biol. 1999, 292 (2): 195-202. 10.1006/jmbi.1999.3091.
Authors are thankful to Jawaharlal Nehru University for usage of all computational facilities and for providing funds for publication of this article. AG is grateful to University Grants Commission, India for the Faculty Recharge Position.
The authors declare that they have no competing interests.
SS, AG and SKG designed the methods and experimental setup. SS and CT carried out the implementation of the various methods. SS, CT, AG and SKG wrote the manuscript. All authors have read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Soni, S., Tyagi, C., Grover, A. et al. Molecular modeling and molecular dynamics simulations based structural analysis of the SG2NA protein variants. BMC Res Notes 7, 446 (2014). https://0-doi-org.brum.beds.ac.uk/10.1186/1756-0500-7-446
- Molecular modeling
- Molecular dynamics simulations
- Disorder prediction