Structural Analysis and Conformational Dynamics of STN1 Gene Mutations Involved in Coat Plus Syndrome

The human CST complex (CTC1–STN1–TEN1) is associated with telomere functions including genome stability. We have systemically analyzed the sequence of STN and performed structure analysis to establish its association with the Coat Plus (CP) syndrome. Many deleterious non-synonymous SNPs have been identified and subjected for structure analysis to find their pathogenic association and aggregation propensity. A 100-ns all-atom molecular dynamics simulation of WT, R135T, and D157Y structures revealed significant conformational changes in the case of mutants. Changes in hydrogen bonds, secondary structure, and principal component analysis further support the structural basis of STN1 dysfunction in such mutations. Free energy landscape analysis revealed the presence of multiple energy minima, suggesting that R135T and D157Y mutations destabilize and alter the conformational dynamics of STN1 and thus may be associated with the CP syndrome. Our study provides a valuable direction to understand the molecular basis of CP syndrome and offer a newer therapeutics approach to address CP syndrome.

Introduction

Telomere is a complex of protein and nucleic acid, crucial for protecting chromosome ends from degradation, end-to-end chromosome fusion, and activation of DNA damage response (DDR) (de Lange, 2009). Abnormal length of a chromosome or deprotected chromosome ends compromise replication potential and thus genome stability (Bodnar et al., 1998). Telomere shortening or degradation are coupled with many human diseases, termed as telomeropathies (Armanios and Blackburn, 2012; Holohan et al., 2014). In addition, the CST complex (CTC1–STN1–TEN1) plays a vital role in the synthesis of C-strand of telomere (Chen and Lingner, 2013) and helps in genome-wide replication and recovery from replication fork stalling during replication stress (Stewart et al., 2012, 2018). The conserved CST complex interacts with the G-strand of the telomere, promotes C-strand synthesis, and suppresses telomerase-mediated elongation of telomere (Chen et al., 2012; Stewart et al., 2012; Feng et al., 2017, 2018). Interruption of CST complex directly affecting the synthesis of C-strand and telomere length leads to the formation of elongated 3′ overhang (Gu et al., 2012).

Regulation of C-strand synthesis and telomere length is primarily regulated by CTC1–STN1 and STN1–TEN1 complex formation (Feng et al., 2018). STN1 consists of an N-terminal OB-fold comprising of a β-barrel, consisting of two- to three-stranded β-sheet sandwiched by three α-helices (Bryan et al., 2013). The N-terminal OB-fold of STN1 forms a strong complex with the OB-fold of TEN1 protein (STN1–TEN1) and thus enables the CST complex for C-strand synthesis. The C-terminal domain of STN1 is composed of β-strands and 11 α-helices organized into two different winged helix-turn-helix (wHTH) motifs (Ganduri and Lue, 2017). Mutations of amino acid residues important for STN1–TEN1 complex formation results in elongated telomeres and several other chromosomal abnormalities including Coat Plus (CP) syndrome (Bryan et al., 2013; Simon et al., 2016).

Owing to the significance of the STN1–TEN1 complex in maintaining telomere homeostasis, mutations in the STN1 gene are reported to cause CP syndrome. CP is identified by calcification in the intracranial region, hematological abnormalities, and neurologic and retinal defects (Simon et al., 2016). Patients with CP often present shortened telomeres, indicating that telomerase malfunctioning is associated with the pathogenesis. To date, only two STN1 mutations (R135T and D157Y) have been reported that causes CP syndrome. However, the molecular basis of such pathogenesis remained elusive (Simon et al., 2016). In addition, in vitro mutational analysis in the STN1 and TEN1 gene depicts that mutation in the TEN1 gene, in particular R27Q, Y115A, and R119Q, shows a marked change in their dissociation constant. However, STN1 double mutants (D78A/I164A and D78A/M167A) show a complete loss of binding with TEN1 (Simon et al., 2016).

Herein, we have analyzed the complete mutational spectrum in the STN1 gene to identify the disease-causing mutations and subsequent pathogenic characterization based on their impact on structure and functions. To understand the molecular basis of CP syndrome, the structural and conformational changes in R135T and D157Y mutants were extensively studied at an atomic level using 100 ns molecular dynamics (MD) simulation. The results possibly unveil an understanding of R135T and D157Y mutations and their association with the CP syndrome.

Materials and Methods

Collection of Dataset

FASTA sequence of STN1 was retrieved from the UniProt database (UniProt ID: Q9H668). Distribution of SNPs was collected from Ensembl (Hubbard et al., 2002), dbSNP (Sherry et al., 2001), and OMIM (Amberger et al., 2008) databases. Functional annotation of each SNP was extracted from the dbSNP database. Structures of STN1 were obtained from the Protein Data Bank (PDB code: 4JOI and 4JQF) (Berman et al., 2000).

Prediction of Deleterious nsSNPs

Deleterious or damaging nsSNPs in the STN1 gene were predicted by utilizing Sorting Intolerant from Tolerant (SIFT) (Kumar et al., 2009), PolyPhen 2.0 (Adzhubei et al., 2013), and PROVEAN (Choi and Chan, 2015) web servers. SIFT predictions are based on the sequence homology and it differentiates nsSNPs as tolerant (neutral) or intolerant (disease) on the basis of a predicted score (deleterious if a score is <0.05 and neutral if a score >0.05). PolyPhen 2.0 calculates the impact of point mutations on the structure of protein as well as its effects on phenotype. A detailed description of methods for deleterious nsSNP prediction is provided in our earlier communication (Amir et al., 2019).

Prediction of Destabilizing nsSNPs

Protein stability is represented by the change in the Gibbs free energy (ΔG) upon folding. In this study, we have utilized five stability prediction tools, DUET (Pires et al., 2014), SDM2 (Pandurangan et al., 2017), mCSM (Pires et al., 2013), CUPSAT (Parthiban et al., 2006), and STRUM (Quan et al., 2016). DUET is an integrated computational method used for calculating the impact of mutations on the stability of the protein. SDM2 differentiates the competence of amino acid residues between the WT and mutant proteins, whereas the mCSM algorithm establishes effects of mutations depending on the structural signature. CUPSAT uses protein structure environment-specific potential and torsion angle potential to find the differences in stability of WT and mutant proteins. It requires a PDB structure and the position of the substituted residue. The output entails information about its structural features, sites of mutation, and inclusive information about the change in the stability of protein. STRUM calculates the impact of mutations using conservation score obtained from alignment of the multiple-threading template.

Prediction of Pathogenic nsSNPs

MutPred2 is a computational server used to predict the molecular basis of disease-associated substitutions. It employs many attributes such as protein structure, its function, and evolution. MutPred2 utilizes three different servers, PSI-BLAST (Altschul et al., 1997), SIFT (Ng and Henikoff, 2003), and Pfam (Finn et al., 2015) together with some structural disorder prediction methods such as MARCOIL (Delorenzi and Speed, 2002), TMHMM (Krogh et al., 2001), and DisProt (Sickmeier et al., 2006). Thus, MutPred combines the score of these computational servers and then provides prediction results. PhD-SNP is a single-sequence SVM-based tool that distinguishes disease-associated mutations depending on the local environment of substitution (Capriotti et al., 2006).

Aggregation Propensity Analysis

Solubility is an important feature required for the function of a protein. Predicting solubility of a particular protein is critical for protein engineering and normal functioning. SODA is a newly developed aggregation prediction method (Paladin et al., 2017). It predicts the change in protein solubility based on a number of physicochemical features. SODA utilizes the propensity of amino acid sequence to aggregate, intrinsic disorder, secondary structure, and hydrophobicity to calculate the differences in protein solubility (Paladin et al., 2017).

Structure Refinement

The atomic coordinates of the STN1 structure were retrieved from the STN1–TEN1 complex available in the Protein Data Bank (PDB ID: 4JOI). The STN1–TEN1 complex comprises of the N-terminal domain of STN1 (residues 24–184) and full-length TEN1 (residue 2–124). The structure of the N-terminal domain of STN1 has a missing loop (residues 92–110) that was modeled using MODELER 9.20 imbedded in PyMOL plugin PyMod 2.0. Further, the WT N-terminal domain of STN1 was used to create R135T and D157Y mutations using mutagenesis plugin embedded in PyMOL (www.pymol.org). Mutant structures and the WT structure of STN1 were energy-minimized using SPDBV to remove high energy configurations by changing their coordinate geometries in such a way as to release internal constraints and reduce the total potential energy.

Molecular dynamics simulations

To in-depth understand the effects of mutations on the STN1 structure, all-atom MD simulations were carried out for 100 ns under explicit water solvent conditions for WT and its mutants (R135T and D157Y). MD simulations were performed at 300 K at the molecular mechanics level implemented in the GROMACS 5.1.2 software package using the GROMOS96 43a1 force field. All systems were soaked in a cubic box of water molecules with a dimension of 10 Å using the gmx editconf module for setting boundary conditions and gmx solvate module for solvation. Further, the systems were subsequently immersed in a box having a simple point charge (SPC16) water model. Na⁺ and Cl⁻ ions were aided further in the systems for neutralizing and preserving a physiological concentration (0.15 M) using the gmx genion module. All the systems were minimized using 1,500 steps of steepest descent. All systems were equilibrated at a constant temperature, 300 K, by utilizing the two-step ensemble process (NVT and NPT) for 100 ps. Initially, the Berendsen thermostat with no pressure coupling was employed for the NVT (i.e., constant number of particles, volume, and temperature) canonical ensemble, and then we use the Parrinello–Rahman method pressure of 1 bar (P) for the NPT ensemble (i.e., constant particle number, pressure, and temperature). The final simulations were performed for each system for 100 ns where leap-frog integrator was applied for the time evolution of trajectories. The details of MD simulations have been described elsewhere (Gulzar et al., 2018; Naqvi et al., 2018).

Analysis of MD Trajectories

All the trajectory files were analyzed using trajectory analysis module embedded in the GROMACS simulation package and Visual Molecular Dynamics (VMD) software. The trajectory files were analyzed by using gmx confirms, gmx rmsd, gmx rmsf , gmx gyrate, gmx sasa, gmx hbond, gmx covar, gmx anaeig, gmx energy, gmx do_dssp, and gmx sham GROMACS utilities to extract the graph of root-mean-square deviation (RMSD), root-mean-square fluctuations (RMSFs), radius of gyration (R_g), solvent accessible surface area (SASA), hydrogen bond, principal component, secondary structure, etc. (Syed et al., 2018). All the graphs were plotted using the Qt Grace Visualization tool.

Principal Component and Free Energy Landscape Analysis

Principal component analysis (PCA) was done to study the collective motions of WT and mutant STN1. This method employs the calculation and diagonalization of the covariance matrix. The covariance matrix is calculated as:

$\begin{array}{l}{C}_{ij}=<(x_i-<x_i>)(x_j-<x_j>)>\end{array}$

where x_i/x_j is the coordinate of the ith/jth atom of the systems, and < —> represents an ensemble average.

Free energy landscape (FEL) can be used to understand the stability, folding and function of the protein. The FEL can be constructed as:

$\begin{array}{l}\Delta G(X)=-K_{\text{B}}TlnP(X)\end{array}$

where K_B and T are the Boltzmann constant and absolute temperature, respectively, and P(X) is the probability distribution of the molecular system along the PCs.

Results and Discussion

This study provides the structural and functional consequences of nsSNPs in the STN1 gene. In addition, disease-causing or pathogenic spectrum, aggregation behavior, and conservation score were screened using advanced computational methods. Finally, the atomistic levels of two pathogenic mutations (R135T and D157Y) causing CP syndrome have been analyzed in detail using all-atom MD simulation approach.

Prediction of Deleterious and Destabilizing nsSNPs in STN1 Gene

For prediction of deleterious and destabilizing nsSNPs, we have cross-checked information present in the dbSNP and UniProt databases, removed invalid mutations based on wrong amino acid position and alignment, and merged or removed data with other nsSNPs in dbSNP. As a result, a sum of 243 nsSNPs in the STN1 gene was considered for analysis (Tables S1, S2). Distribution of deleterious nsSNPs in STN1 gene using SIFT, PolyPhen 2.0, and PROVEAN predicated that 97 (~40%), 98 (40%), and 67 (27%) nsSNPs were found deleterious, respectively (Figure 1 and Table S1). Similarly, predictions of destabilizing nsSNPs employing CUPSAT 163 (67%), SDM 170 (70%), DUET 183 (75%), mCSM 219 (90%), and STRUM 232 (95%) nsSNPs depict a destabilizing effect (Figure 1 and Table S2). The nsSNPs that depict a deleterious effect by at least two different tools and the destabilizing behavior by at least three different methods were collected for phenotypic spectrum analysis. A total of 76 nsSNPs were taken for this analysis.

FIGURE 1

Figure 1. Prediction of deleterious and destabilizing mutations using different computational methods. Numbers of neutral and deleterious mutations are indicated in green and red, respectively.

Analysis of Molecular Phenotype

Molecular phenotype analysis complements the characterization of disease-associated nsSNPs in the human genome. All 76 deleterious and destabilizing nsSNPs were subjected to MutPred 2.0 and PhD-SNP computational servers to find the relation of these deleterious and destabilizing nsSNPs with the disease phenotype. Prediction of disease phenotype using MutPred 2.0 and PhD-SNP shows that 39 (51%) and 41 (53%) nsSNPs are linked with disease phenotype, respectively (Table S3). Further, about 30 (39%) mutations identified as disease-causing from both prediction methods were further used for aggregation propensity scanning.

Aggregation Propensity Analysis

The solubility of protein is one of the important characteristics that primarily belong to the concentration, conformation, and location of the protein. It plays a critical role in protein turnover in the cells (Ross and Poirier, 2004; Sami et al., 2017). Protein aggregation often tends to the progression of a range of pathologies such as Alzheimer's (Thal et al., 2015) and Parkinson's (Tan et al., 2009). Analysis of aggregation behavior using the SODA web server aids in molecular identification of a disease or pathogenic nsSNPs to the conformational level. Of 30 pathogenic mutations identified using MutPred and PhD-SNP prediction servers, 10 mutations (G45D, G51V, C88Y, R135T, D157Y, P158S, R166G, Y174C, G278R, and C312F) show a decrease in the solubility score (Table 1). Results depict that majority of aggregate-forming mutations lie in the N-terminal domain of STN1 in comparison to the C-terminal domain.

TABLE 1

Table 1. Predictions of aggregation propensities of pathogenic mutations in the STN1 gene.

Sequence Conservation Analysis

Analysis of conserved amino acid residues in protein provides a better understanding of the significance of a particular amino acid residue and probably its localized evolution. Sequence conservation analyses in the STN1 protein have been performed using the ConSurf tool, which depicts that most of the amino acid residues in the STN1 protein are highly conserved (Figure 2). In particular, the stretches of highly conserved amino acid residues such as 23–30, 42–45, 56–58, 60–65, 75–86, 88–90, 153–155, 173–175, 178–180, 228–230, and 321–323 were found. From this conservation analysis, it is very clear that the N-terminal domain of STN1 is considerably more conserved than the C-terminal domain. The stretches of variable residues in the STN1 protein had the following ranges: 2–10, 90–120, 180–207, and 230–260. Variable stretches are usually long and mainly concentrated in the middle of protein length compared to conserved stretches, which are usually short and scattered throughout the length of the STN1 protein.

FIGURE 2

Figure 2. Sequence conservation analysis in STN1 protein using ConSurf. The conservation score from highly variable to conserve is represented in a scale of 1–9.

Association of STN1 Domains With Disease-Causing Mutations

Analysis of pathogenic/disease-causing mutations in different domains of the STN1 gene entails relative information about a particular domain to be disease-causing or neutral. Distributions of pathogenic mutations were analyzed in all three domains (OB1, wHTH1, and wHTH2) of STN1. We have calculated the percentage of pathogenic mutations in each domain by dividing total mutation occurring in a particular domain to the pathogenic mutations found in the same domain. The nsSNPs in OB1, wHTH1, and wHTH2 domains have a 16%, 5%, and 14% chance of being pathogenic, respectively (Figure 3). These results depict that the N- and C-terminal domain of STN1 (OB1 and wHTH2) are possibly more prone to pathogenesis. These interpretations are further complemented by a high conservation score of the N- and C-terminal domain of STN1 in contrast to the wHTH1 domain. How does a particular mutation disrupt the structure of STN1 and its function? To answer this, we have performed 100-ns MD simulation of R135T and D157Y, to understand the molecular mechanism of the disease.

FIGURE 3

Figure 3. Association of pathogenic mutations in different domains of STN1 protein.

Mutations Induced Conformational Stability, Flexibility, and Dynamic Changes

MD simulation was employed to investigate the disruptive effects of R135T and D157Y (run in duplicate) mutations on the conformational stability of STN1 through 100-ns MD simulations. The initial and final structures of WT, R135T, and D157Y during simulation runs are shown in Figure S1. We did not find any significant structural changes in WT except for the changes in the loop connecting strand β3 and β4, which become more flexible after 100 ns of simulation. Interestingly, loss of the C-terminal helix structure was observed in the case of D157Y. The final structures obtained at 100-ns simulation were superimposed and are shown in Figure S2. The figure clearly shows that there is significant change upon D157Y mutation as the RMSD between WT and D157Y structures is 3.51 as compared to R135T with an RMSD of 1.96. Overall changes in the STN1 stability upon mutations were investigated by RMSD calculation. It depicts the differences in backbone flexibility of C^α atoms of WT and mutants. The average RMSD values for WT, R135T, and D157Y mutations were calculated as 0.64, 0.62, and 0.60 nm, respectively (Table 2). These results indicate no significant differences in the RMSD values of mutants (Figure 4A). Moreover, we did not find any significant changes in any of the structure parameters of two simulation runs in D157Y (Figure S3). For consistency, the probability distribution function (PDF) of WT conformers was confined within a range of 0.5–0.8 nm, whereas mutations R135T and D157Y conformers lie within a range of 0.4–0.8 and 0.3–0.8 nm, respectively (Figure 4B). These results also indicated no remarkable differences in the values of PDF of the WT and mutants.

TABLE 2

Table 2. The calculated parameters for all the systems obtained after 100-ns MD simulations.

FIGURE 4

Figure 4. Conformational changes in STN1 protein and its mutant. (A) Graphical representation of backbone conformation. (B) Probability distribution functions of RMSD. (C) Radius of gyration. (D) Probability distribution functions of the compactness of wild-type and STN1 mutants. Wild-type (black), R135T (green), and D157Y (red).

Compactness of WT and STN1 mutants is denoted by R_g, which represents the compactness of the protein structure and is associated with stability (Lobanov et al., 2008). Differences in R_g values between WT and mutants are shown in Figure 4C. Average R_g values of WT, R135T, and D157Y mutations were calculated as 1.59, 1.61, and 1.83 nm, respectively (Table 2). The mutation D157Y shows a significant increase in R_g, suggesting a loss in compactness. PDF analysis further suggested a higher probability of R_g value of D157Y compared to R135T and WT (Figure 4D).

We have also investigated the hydrophobic core region of WT and both mutants by calculating the change in solvent accessibility surface area (SASA). A significant increment in average SASA has been observed in the D157Y mutant (Figure 5A and Table 2). However, R135T mutation does not show any significant change in the SASA value compared to WT as Arg135 is located on the surface of the protein. D157Y mutation brings larger changes in the SASA value because Asp157 is partially buried and replacement of Asp157 by Tyr renders this residue more accessible to solvent. The increase in average SASA value and PDF suggested that D157Y has a large surface exposed to solvent, and this might cause the exposure of hydrophobic residues and subsequently unfolding of the protein (Figure 5B).

FIGURE 5

Figure 5. (A) Graphical representation of SASA. (B) Probability distribution of SASA. (C) Residual fluctuations. (D) Intra-protein hydrogen bond analyses of wild-type (black), R135T (green), and D157Y (red) mutants of STN1.

Residual flexibility of WT and mutants was estimated and is shown in Figure 5C. It is evident that the D157Y mutation reduces the fluctuations of Asp78, which is implicated in salt bridge formation with the residue of TEN1 (Bryan et al., 2013). Disruption in salt bridge formation between Asp78 of STN1 and Arg27 of TEN1 results in loss of about 5.5-fold affinity between both proteins. A significant increase in fluctuations of residues from 94 to 124 was observed. These residues, which are highly non-conserved and belong to the α2 helix, lie between strands β5 and β6 of STN1 (Figure 2). The exact function of the α2 helix is not known possibly due to partial visibility of the electron density map (Bryan et al., 2013). Moreover, we observed the increase in the fluctuations of residues 156–174 belonging to the α3 helix of STN1, which is implicated in hydrophobic interaction with residues of the α2 helix of the TEN1 protein. An interruption in the formation of a hydrophobic interaction of these residues reduces the interaction of STN1 and TEN1 with binding affinity ranges from 2.5- to 5.5-fold (Bryan et al., 2013). Thus, Asp78 forms a salt bridge with Arg27 of TEN1. This salt bridge plays a critical role for proper functioning of CST complex. We have analyzed the change in deviation, compactness, and SASA of Asp78 during the course of simulation (Figure S4).

Hydrogen bond (HB) analysis in WT and STN1 mutants is essential to understand the stability and flexibility of the protein. Figure 5D illustrates the number of intramolecular HB network in the WT and mutant STN1 protein during the course of simulation. The number of HB decreases in R135T and D157Y mutants as compared to the WT (Table 2). WT STN1 has an average of 114 HB, whereas R135T and D157Y mutants have an average of 105 and 106 HBs, respectively. This reduction in HB reflects a loss of structure and compactness of the STN1 protein. Overall, D157Y mutation disrupts the structure, stability, and dynamics of STN1, which leads to interruption of STN1 and TEN1 interaction and thus results in telomere signal free ends and elongated and fragile telomeres, acquiring phenotypes related to telomere length deregulation and dysfunctional CST complex.

Secondary Structure Analysis

The propensity of secondary structural content is a crucial component to study the structural behavior of the protein. We have investigated the changes in secondary structure in WT and R135T and D157Y mutants as shown in Figure 6. The mutation R135T does not induce any significant change in the secondary structure content (Figure 6B), while a large decrease in overall structure was noticed in the case of the D157Y mutant (Figure 6C). The large changes in secondary structure content in the D157Y mutant are consistent with the results obtained from other conformational analysis. Our MD analysis has clearly indicated that the D157Y mutation disrupts the conformational stability, dynamics, and flexibility of STN1, which may be a possible cause of CP syndrome.

FIGURE 6

Figure 6. Secondary structure changes during the course of 100-ns MD simulation in (A) wild-type, (B) R135T, and (C) D157Y mutants.

Principal Component and Free Energy Landscape Analysis

Generally, proteins carry out their specific functions through collective atomic motions. Hence, a collective atomic motion of a particular protein is used as a parameter to understand the stability of protein. PCA is used to investigate the global motions of protein into a few principal motions, characterized by eigenvectors and eigenvalues. Figure 7 shows the conformational sampling of WT and mutant SNT1 in the essential subspace by projecting the C^α atom, showing the tertiary conformations along eigenvectors 1 and 2. From the projection of PC1 and PC2 of WT and mutant, we could envisage a cluster of stable states in WT and mutants. Results clearly depict that the mutants cover a wide range of phase spaces as compared to WT, especially D157Y. The large increase in overall motion in D157Y might be considered as the prime cause in the impairment of protein function.

FIGURE 7

Figure 7. Projection of C^α atoms in essential subspace along the first two eigenvectors of wild-type (black), R135T (green), and D157Y (red) mutant of STN1.

To distinguish the conformational states of WT and mutant STN1, the Gibbs free energy landscape (FEL) was calculated using the first two principal components as reaction coordinates. Using PCA, Helmoltz free energy change is calculated and the FELs thus obtained from the simulations are plotted, as shown in Figure 8. The FEL can provide remarkable information about the different conformational states accessible to the protein in the simulation. Figure 8 shows the FEL of (A) WT and its mutants (B) R135T and (C) D157Y along the two principal components. As can be seen from the figure, WT STN1 displayed a well-defined single large global energy minima basin associated with its conformational state. R135T, in addition to native energy minima, explores wide global energy minima through a transition state. This energy minimum corresponds to a structure very much similar to WT STN1 with some loss of irregular secondary structure like bend and turn.

FIGURE 8

Figure 8. Free energy landscape analysis of (A) wild-type, (B) R135T, and (C) D157Y mutants of STN1 along with the structural snapshots present in the stable minima of each FEL.

Interestingly, in D157Y, native basin splits into two relative energy minima, suggesting the destabilization of native state. The presence of multiple energy minima in the conformational space achieved in D157Y indicates the significant destabilization of the protein. Moreover, the structural snapshots revealed the presence of a structure with significant loss of C-terminal helices α3 along with minor loss of β2 and β3 residues with the overall increase in coil structure.

Conclusions

Here, we have investigated the effect of deleterious or destabilizing mutations in the STN1 protein that are presumably associated with the CP syndrome. Mutational landscape analysis suggested that about one-third of nsSNPs in the STN1 gene are deleterious, and a high rate of pathogenic occupancy was found in the OB1 domain in contrast to the wHTH1 and wHTH2 domains. Backbone conformation, residual flexibility, compactness, solvent accessibility, hydrogen bond, secondary structure, protein motion, and FEL analyses clearly indicated the significant structural and functional loss in the STN1 protein due to D157Y mutation. However, mutation R135T behaves similar to the WT. In conclusion, our results provide an in-depth understanding of destabilization and loss of conformational dynamics of STN1 mediated by D157Y mutation. The results may thus be further exploited to understand the cause of CP syndrome and development of novel strategies for the therapeutic management of CP.

Author Contributions

MA and MH: conceptualization and project administration. MA, TM, VK, and RD: methodology. MA, MR, AH, MFA, PA, and RD: software. TM, AH, MFA, AI, and MH: validation. TM, VK, MR, and FA: formal analysis. MA, TM, VK, and MH: investigation and writing—original draft preparation. MFA, AH, and MH: resources. MA, MR, and AH: data curation. MFA, AI, and FA: writing—review and editing. TM: visualization. MH: supervision. MFA and MH: funding acquisition.

Funding

This work is supported by the Council of Scientific and Industrial Research, India for the award of Senior Research Fellowship to MA (09/466(0197) 2k18 EMR-7).

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

MA acknowledges the Council of Scientific and Industrial Research, India, for the award of Senior Research Fellowship. MH is thankful to the Indian Council of Medical Research, India, and Department of Science and Technology, Government of India, for financial support. MFA acknowledges the generous support from the Deanship of Scientific Research at King Saud University, Riyadh, Kingdom of Saudi Arabia (Grant No. RGP-150). FA is thankful to the Indian National Science Academy for the award of Senior Scientist Position.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2019.00041/full#supplementary-material

References

Adzhubei, I., Jordan, D. M., and Sunyaev, S. R. (2013). Predicting functional effect of human missense mutations using PolyPhen-2. Curr. Protoc. Hum. Genet. 76, 7.20. 21–27.20. 41. doi: 10.1002/0471142905.hg0720s76

CrossRef Full Text | Google Scholar

Altschul, S. F., Madden, T. L., Schäffer, A. A., Zhang, J., Zhang, Z., Miller, W., et al. (1997). Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25, 3389–3402. doi: 10.1093/nar/25.17.3389

PubMed Abstract | CrossRef Full Text | Google Scholar

Amberger, J., Bocchini, C. A., Scott, A. F., and Hamosh, A. (2008). McKusick's online Mendelian inheritance in man (OMIM^®). Nucleic Acids Res. 37 (suppl_1), D793–D796. doi: 10.1093/nar/gkn665

CrossRef Full Text | Google Scholar

Amir, M., Kumar, V., Mohammad, T., Dohare, R., Hussain, A., Rehman, M. T., et al. (2019). Investigation of deleterious effects of nsSNPs in the POT1 gene: a structural genomics-based approach to understand the mechanism of cancer development. J. Cell. Biochem. 120, 10281–10294. doi: 10.1002/jcb.28312

PubMed Abstract | CrossRef Full Text | Google Scholar

Armanios, M., and Blackburn, E. H. (2012). The telomere syndromes. Nat. Rev. Genet. 13:693. doi: 10.1038/nrg3246

PubMed Abstract | CrossRef Full Text | Google Scholar

Berman, H. M., Westbrook, J., Feng, Z., Gilliland, G., Bhat, T. N., Weissig, H., et al. (2000). The protein data bank. Nucleic Acids Res. 28, 235–242. doi: 10.1093/nar/28.1.235

CrossRef Full Text | Google Scholar

Bodnar, A. G., Ouellette, M., Frolkis, M., Holt, S. E., Chiu, C.-P., Morin, G. B., et al. (1998). Extension of life-span by introduction of telomerase into normal human cells. Science 279, 349–352. doi: 10.1126/science.279.5349.349

PubMed Abstract | CrossRef Full Text | Google Scholar

Bryan, C., Rice, C., Harkisheimer, M., Schultz, D. C., and Skordalakes, E. (2013). Structure of the human telomeric Stn1–Ten1 capping complex. PLoS ONE 8:e66756. doi: 10.1371/journal.pone.0066756

PubMed Abstract | CrossRef Full Text | Google Scholar

Capriotti, E., Calabrese, R., and Casadio, R. (2006). Predicting the insurgence of human genetic diseases associated to single point protein mutations with support vector machines and evolutionary information. Bioinformatics 22, 2729–2734. doi: 10.1093/bioinformatics/btl423

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, L.-Y., and Lingner, J. (2013). CST for the grand finale of telomere replication. Nucleus 4, 277–282. doi: 10.4161/nucl.25701

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, L.-Y., Redon, S., and Lingner, J. (2012). The human CST complex is a terminator of telomerase activity. Nature 488:540. doi: 10.1038/nature11269

PubMed Abstract | CrossRef Full Text | Google Scholar

Choi, Y., and Chan, A. P. (2015). PROVEAN web server: a tool to predict the functional effect of amino acid substitutions and indels. Bioinformatics 31, 2745–2747. doi: 10.1093/bioinformatics/btv195

PubMed Abstract | CrossRef Full Text | Google Scholar

de Lange, T. (2009). How telomeres solve the end-protection problem. Science 326, 948–952. doi: 10.1126/science.1170633

PubMed Abstract | CrossRef Full Text | Google Scholar

Delorenzi, M., and Speed, T. (2002). An HMM model for coiled-coil domains and a comparison with PSSM-based predictions. Bioinformatics 18, 617–625. doi: 10.1093/bioinformatics/18.4.617

PubMed Abstract | CrossRef Full Text | Google Scholar

Feng, X., Hsu, S.-J., Bhattacharjee, A., Wang, Y., Diao, J., and Price, C. M. (2018). CTC1–STN1 terminates telomerase while STN1–TEN1 enables C-strand synthesis during telomere replication in colon cancer cells. Nat. Commun. 9, 2827. doi: 10.1038/s41467-018-05154-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Feng, X., Hsu, S.-J., Kasbek, C., Chaiken, M., and Price, C. M. (2017). CTC1-mediated C-strand fill-in is an essential step in telomere length maintenance. Nucleic Acids Res. 45, 4281–4293. doi: 10.1093/nar/gkx125

PubMed Abstract | CrossRef Full Text | Google Scholar

Finn, R. D., Coggill, P., Eberhardt, R. Y., Eddy, S. R., Mistry, J., Mitchell, A. L., et al. (2015). The Pfam protein families database: towards a more sustainable future. Nucleic Acids Res. 44 , D279–D285. doi: 10.1093/nar/gkv1344

PubMed Abstract | CrossRef Full Text | Google Scholar

Ganduri, S., and Lue, N. F. (2017). STN1 POLA2 interaction provides a basis for primase-pol α stimulation by human STN1. Nucleic Acids Res. 45, 9455–9466. doi: 10.1093/nar/gkx621

PubMed Abstract | CrossRef Full Text | Google Scholar

Gu, P., Min, J. N., Wang, Y., Huang, C., Peng, T., Chai, W., et al. (2012). CTC1 deletion results in defective telomere replication, leading to catastrophic telomere loss and stem cell exhaustion. EMBO J. 31, 2309–2321. doi: 10.1038/emboj.2012.96

PubMed Abstract | CrossRef Full Text | Google Scholar

Gulzar, M., Ali, S., Khan, F. I., Khan, P., Taneja, P., and Hassan, M. I. (2018). Binding mechanism of caffeic acid and simvastatin to the integrin linked kinase for therapeutic implications: a comparative docking and MD simulation studies. J. Biomol. Struct. Dyn. doi: 10.1080/07391102.2018.1546621. [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

Holohan, B., Wright, W. E., and Shay, J. W. (2014). Telomeropathies: an emerging spectrum disorder. J. Cell Biol. 205, 289–299. doi: 10.1083/jcb.201401012

PubMed Abstract | CrossRef Full Text | Google Scholar

Hubbard, T., Barker, D., Birney, E., Cameron, G., Chen, Y., Clark, L., et al. (2002). The Ensembl genome database project. Nucleic Acids Res. 30, 38–41. doi: 10.1093/nar/30.1.38

PubMed Abstract | CrossRef Full Text | Google Scholar

Krogh, A., Larsson, B., Von Heijne, G., and Sonnhammer, E. L. (2001). Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. J. Mol. Biol. 305, 567–580. doi: 10.1006/jmbi.2000.4315

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, P., Henikoff, S., and Ng, P. C. (2009). Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm. Nat. Protoc. 4:1073. doi: 10.1038/nprot.2009.86

PubMed Abstract | CrossRef Full Text | Google Scholar

Lobanov, M. Y., Bogatyreva, N., and Galzitskaya, O. (2008). Radius of gyration as an indicator of protein structure compactness. Mol. Biol. 42, 623–628. doi: 10.1134/S0026893308040195

PubMed Abstract | CrossRef Full Text | Google Scholar

Naqvi, A. A., Mohammad, T., Hasan, G. M., and Hassan, M. (2018). Advancements in Docking and Molecular Dynamics Simulations Towards Ligand-receptor Interactions and Structure-function Relationships. Curr. Top. Med. Chem. 20, 1755–1768. doi: 10.2174/1568026618666181025114157

CrossRef Full Text | Google Scholar

Ng, P. C., and Henikoff, S. (2003). SIFT: predicting amino acid changes that affect protein function. Nucleic Acids Res. 31, 3812–3814. doi: 10.1093/nar/gkg509

PubMed Abstract | CrossRef Full Text | Google Scholar

Paladin, L., Piovesan, D., and Tosatto, S. C. (2017). SODA: Prediction of protein solubility from disorder and aggregation propensity. Nucleic Acids Res. 45 , W236–W240. doi: 10.1093/nar/gkx412

PubMed Abstract | CrossRef Full Text | Google Scholar

Pandurangan, A. P., Ochoa-Montaño, B., Ascher, D. B., and Blundell, T. L. (2017). SDM: a server for predicting effects of mutations on protein stability. Nucleic Acids Res. 45, W229–W235. doi: 10.1093/nar/gkx439

PubMed Abstract | CrossRef Full Text | Google Scholar

Parthiban, V., Gromiha, M. M., and Schomburg, D. (2006). CUPSAT: prediction of protein stability upon point mutations. Nucleic Acids Res. 34 (suppl_2), W239–W242. doi: 10.1093/nar/gkl190

PubMed Abstract | CrossRef Full Text | Google Scholar

Pires, D. E., Ascher, D. B., and Blundell, T. L. (2013). mCSM: predicting the effects of mutations in proteins using graph-based signatures. Bioinformatics 30, 335–342. doi: 10.1093/bioinformatics/btt691

PubMed Abstract | CrossRef Full Text | Google Scholar

Pires, D. E., Ascher, D. B., and Blundell, T. L. (2014). DUET: A server for predicting effects of mutations on protein stability using an integrated computational approach. Nucleic Acids Res. 42 , W314–W319. doi: 10.1093/nar/gku411

PubMed Abstract | CrossRef Full Text | Google Scholar

Quan, L., Lv, Q., and Zhang, Y. (2016). STRUM: Structure-based prediction of protein stability changes upon single-point mutation. Bioinformatics 32, 2936–2946. doi: 10.1093/bioinformatics/btw361

PubMed Abstract | CrossRef Full Text | Google Scholar

Ross, C. A., and Poirier, M. A. (2004). Protein aggregation and neurodegenerative disease. Nat. Med. 10, S10–S17. doi: 10.1038/nm1066

PubMed Abstract | CrossRef Full Text | Google Scholar

Sami, N., Rahman, S., Kumar, V., Zaidi, S., Islam, A., Ali, S., et al. (2017). Protein aggregation, misfolding and consequential human neurodegenerative diseases. Int. J. Neurosci. 127, 1047–1057. doi: 10.1080/00207454.2017.1286339

PubMed Abstract | CrossRef Full Text | Google Scholar

Sherry, S. T., Ward, M.-H., Kholodov, M., Baker, J., Phan, L., Smigielski, E. M., et al. (2001). dbSNP: the NCBI database of genetic variation. Nucleic Acids Res. 29, 308–311. doi: 10.1093/nar/29.1.308

PubMed Abstract | CrossRef Full Text | Google Scholar

Sickmeier, M., Hamilton, J. A., LeGall, T., Vacic, V., Cortese, M. S., Tantos, A., et al. (2006). DisProt: The database of disordered proteins. Nucleic Acids Res. 35 (suppl_1), D786–D793. doi: 10.1093/nar/gkl893

PubMed Abstract | CrossRef Full Text | Google Scholar

Simon, A. J., Lev, A., Zhang, Y., Weiss, B., Rylova, A., Eyal, E., et al. (2016). Mutations in STN1 cause Coats plus syndrome and are associated with genomic and telomere defects. J. Exp. Med. 213, 1429–1440. doi: 10.1084/jem.20151618

PubMed Abstract | CrossRef Full Text | Google Scholar

Stewart, J. A., Wang, F., Chaiken, M. F., Kasbek, C., Chastain, P. D., Wright, W. E., et al. (2012). Human CST promotes telomere duplex replication and general replication restart after fork stalling. EMBO J. 31, 3537–3549. doi: 10.1038/emboj.2012.215

PubMed Abstract | CrossRef Full Text | Google Scholar

Stewart, J. A., Wang, Y., Ackerson, S. M., and Schuck, P. L. (2018). Emerging roles of CST in maintaining genome stability and human disease. Front. Biosci. 23, 1564–1586. doi: 10.2741/4661

PubMed Abstract | CrossRef Full Text | Google Scholar

Syed, S. B., Khan, F. I., Khan, S. H., Srivastava, S., Hasan, G. M., Lobb, K. A., et al. (2018). Mechanistic insights into the urea-induced denaturation of kinase domain of human integrin linked kinase. Int. J. Biol. Macromol. 111, 208–218. doi: 10.1016/j.ijbiomac.2017.12.164

PubMed Abstract | CrossRef Full Text | Google Scholar

Tan, J. M., Wong, E. S., and Lim, K.-L. (2009). Protein misfolding and aggregation in Parkinson's disease. Antioxid. Redox Signal. 11, 2119–2134. doi: 10.1089/ars.2009.2490

PubMed Abstract | CrossRef Full Text | Google Scholar

Thal, D. R., Walter, J., Saido, T. C., and Fändrich, M. (2015). Neuropathology and biochemistry of Aβ and its aggregates in Alzheimer's disease. Acta Neuropathol. 129, 167–182. doi: 10.1007/s00401-014-1375-y

PubMed Abstract | CrossRef Full Text | Google Scholar