The NSP6 protein of SARS-CoV-2 is composed of 290 amino acid sequences; accession number YP_009742613 was obtained from the UniProtKB and NCBI databases.
Prediction of secondary structure and solvent accessibility
The prediction of the secondary structure of the SARS-CoV-2 NSP6 protein was investigated using a consensus result that was generated from the eleven servers Hhpred, pep2d, Jpred4, Porter, PredictProtein, Psipred, SOPMA, Raptorx, SSPRO, Yaspin and APSSP2 (Table 1). The predicted secondary structure of the SARS-CoV-2 NSP6 protein from the Psipred server may be composed of 77.93% α-helices, 2.76% β-strands and 19.31% coils. In addition, the predict protein server predicted that the secondary structure of the NSP6 protein of SARS-CoV-2 consists of 68% alpha helices, 19% coils and 13% beta strands. The Helices were analyzed to be 8-35, 41-85, 88-100, 107-130, 136-149, 158-172, 176-196, 216-236, 248-256 and 265-272; Coils, 1-7, 36-40, 86, 87, 101-103, 131-135, 155-157, 173-175, 197, 198, 204-206, 237-241, 247, 257-264, 273, 274 and 278-281; and strands, 104-106, 150-154, 199-203, 207-215, 242-246, 275-277 and 282-289. Additionally, the PredictProtein server could predict solvent accessibility; 76.2% and 23.79% of the proteins were mostly buried intermediately and exposed, respectively. The NSP6 protein of SARS-CoV-2 is found to be a transmembrane protein with a transmembrane helix of 56.5%, an extracellular helix of 16.5% and a cytoplasmic helix of 27%.
Prediction and refinement of tertiary structure
To construct a tertiary structure model, we built many models, refined them and evaluated the full models to select the model with the highest quality.
Construction of the initial model using target-template alignment
Several servers, such as the GalaxyWEB, SWISS-MODEL and LOMETS, are based on the structuring of aligned regions of proteins, while the I-TASSER, Robetta and Phyre2 servers are designed to predict the structures of proteins that have low similarity with known structured proteins. We used AlphaFold which is a server that uses deep learning techniques and neural networks to predict the distance and torsion of proteins. Its data depends on training on known schemes of PDB structures with known amino acid sequences and MSA features. Additionally, the AlphaFold model is based on a template-free modeling algorithm that can predict new unstructured proteins. In addition, the Robetta, C-Quark, LOMETS and I-TASSER servers can predict the five full models for the query according to the confidence score (C-score), which was -1.90 for the best model and a TM-score of 0.45±0.15. Additionally, GalaxyWEB predicts five full models for the query arranged according to the HHsearch results rescoring system. In addition, the Phyre2 server was able to construct eight models of small regions of the query, as was the SWISS-MODEL server.
Reduced-level structure assembly and refinement simulations
To improve the quality of the tertiary structure prediction process, we needed to refine the predicted models of the NSP6 protein of SARS-CoV-2. As a result, we used the 3D-REFINE, ModRefiner, Refold3, DeepRefiner, GalaxyWEB Refine and GalaxyWEB Refine2 servers to refine the models by refining the backbone structure, side chain positioning and hydrogen bonds. The model refinement results revealed that the RMSD decreased to 0.25 and 0.3 A and the TM-score increased to 0.9952 for the AlphaFold model of the NSP6 protein of SARS-CoV-2, which indicates an increase in the quality of the predicted structure compared to that of the initially predicted structure.
Model evaluation and selection
Fig 1 (a, b) shows the best predicted three-dimensional structure of the NSP6 protein of SARS-CoV-2 using the AlphaFold server, which was generated by Autodock Vina software different tools of view ribbon view, solvent-accessible surface and atomic representation of the protein. Fig 1 (c) shows the binding sites of SARS-CoV-2 NSP6 using Autodock Vina software. Finding the closest superimposed model is an important process in which we can construct a structural framework for obtaining additional information in protein design, interpreting functional data, determining the molecular function of the protein and designing and targeting drugs. I-TASSER server was used to predict the closest superimposed experimentally structured models to the NSP6 protein of SARS-CoV-2. I-TASSER server revealed that the experimental structure 5aexJ template was the closest predicted structural model of the SARS-CoV-2 NSP6 protein.
Model evaluation and selection of the best model
To determine the best model of the NSP6 protein with good topology, all the predicted and refined models were evaluated, Table 2. The model predicted by AlphaFold and then refined by DeepRefiner presented the best results, with an overall quality of 99.64%, according to the ERRAT server compared to all the other prediction models. The Ramachandran plot generated with the PROCHECK server revealed 94.4% of the residues in the core regions in addition to 5.6% in the allowed regions and no residues in the disallowed regions. The SWISSMODEL structure Assessment server supported the results of the Ramachandran plot generated by the PROCHECK server, which included 98.35% of the residues in the favored regions and 0.41% of the residues in the outgroup regions. Additionally, according to the PROSA server, the model had a molprobity value of 2.61 and a clash score of 139.01, the mean Z score was -1.33 and the Z score was -4.12. These Z scores clarified that the modeled protein was within the range of X-ray-solved protein structures, which indicated that it was a good model. The PROQ server showed that the LG score was 8.807, which indicates an extremely good model according to the evaluation rules and the maximum substitution value was -0.544.
The absence of an experimental structure for SARS-CoV-2 NSP6 has been a significant impediment to targeted antiviral development. This study was designed to bridge this gap by implementing a rigorous, multi-step computational pipeline to generate a high-confidence tertiary structure model. Our approach was predicated on the rationale that a consensus from multiple, independent prediction and refinement methods, followed by stringent validation, would yield a more reliable model than any single method alone. The integration of multiple computational platforms has also been successfully adopted in comparative genomic investigations, where combining complementary bioinformatics approaches improves the robustness and biological interpretation of computational analyses
(Liu et al., 2025).
Our initial analysis confirmed that NSP6 is a highly helical transmembrane protein, a finding consistent with its known biological role in anchoring replication complexes to endoplasmic reticulum membranes
(Abdelkader et al., 2022; Lu et al., 2020). This secondary structure profile provided a crucial sanity check for the subsequent tertiary models, as a plausible model would need to reflect this fundamental characteristic.
The cornerstone of our strategy was the generation of a diverse model set using eight distinct prediction servers, encompassing both template-based and
ab initio approaches. The initial models displayed considerable variation, underscoring the challenge of predicting transmembrane protein structures and the necessity for a comparative approach. The subsequent refinement step proved critical, significantly enhancing model quality by resolving steric clashes and optimizing hydrogen-bonding networks. The marked improvement in metrics such as RMSD and TM-score post-refinement highlights that this is not a superfluous step but an essential one for achieving atomic-level realism
(Gordon et al., 2020; Zhang et al., 2020).
The comprehensive validation phase, utilizing a battery of assessment tools, provided clear and convergent evidence for the superiority of the model generated by AlphaFold and refined by DeepRefiner. This model consistently ranked highest across multiple independent quality metrics. Its exceptional ERRAT score (99.64%) signifies a high-quality non-bonded atomic interaction profile, while its Ramachandran plot (94.4% core, 0% disallowed) confirms the stereochemical excellence of its backbone dihedral angles
(Zhou et al., 2011). Furthermore, the PROSA Z-score of -1.33 definitively places this model within the range of scores typically observed for experimentally determined X-ray structures of similar size (
Wiederstein and Sippl, 2007). While other servers like Robetta also produced good models, the AlphaFold-DeepRefiner combination demonstrated superior and consistent performance across all validation parameters.
It is noteworthy that AlphaFold, a deep learning system, outperformed traditional homology modeling and threading-based servers. This aligns with the broader revolution in protein structure prediction, where neural networks trained on known structures can now accurately predict distances and angles, even for proteins with few homologous templates
(Guo et al., 2020; van der Hoeven et al., 2016). Our study provides a concrete example of this superiority for a challenging viral transmembrane protein.
A potential limitation of any computational study is the lack of experimental validation. However, our multi-pronged validation strategy, demonstrating convergence across multiple independent quality metrics, provides a high degree of confidence in the selected model. Collectively, these findings support the reliability and structural plausibility of the predicted NSP6 model. Importantly, the availability of a high-confidence NSP6 structure may facilitate future structure-based drug discovery approaches, including molecular docking, virtual screening and molecular dynamics simulations aimed at identifying potential inhibitors targeting NSP6-mediated host-virus interactions. Similar structure-guided in silico approaches have been successfully applied for identifying and evaluating potential therapeutic compounds against multiple biological targets, highlighting the value of reliable structural models in rational drug discovery (
Veni et al., 2022). Future molecular dynamics simulations in membrane-mimetic environments may further validate the structural stability of the predicted NSP6 structure.