In silico functional, structural and pathogenicity analysis of missense single nucleotide polymorphisms in human MCM6 gene

Single nucleotide polymorphisms (SNPs) are one of the most common determinants and potential biomarkers of human disease pathogenesis. SNPs could alter amino acid residues, leading to the loss of structural and functional integrity of the encoded protein. In humans, members of the minichromosome maintenance (MCM) family play a vital role in cell proliferation and have a significant impact on tumorigenesis. Among the MCM members, the molecular mechanism of how missense SNPs of minichromosome maintenance complex component 6 (MCM6) contribute to DNA replication and tumor pathogenesis is underexplored and needs to be elucidated. Hence, a series of sequence and structure-based computational tools were utilized to determine how mutations affect the corresponding MCM6 protein. From the dbSNP database, among 15,009 SNPs in the MCM6 gene, 642 missense SNPs (4.28%), 291 synonymous SNPs (1.94%), and 12,500 intron SNPs (83.28%) were observed. Out of the 642 missense SNPs, 33 were found to be deleterious during the SIFT analysis. Among these, 11 missense SNPs (I123S, R207C, R222C, L449F, V456M, D463G, H556Y, R602H, R633W, R658C, and P815T) were found as deleterious, probably damaging, affective and disease-associated. Then, I123S, R207C, R222C, V456M, D463G, R602H, R633W, and R658C missense SNPs were found to be highly harmful. Six missense SNPs (I123S, R207C, V456M, D463G, R602H, and R633W) had the potential to destabilize the corresponding protein as predicted by DynaMut2. Interestingly, five high-risk mutations (I123S, V456M, D463G, R602H, and R633W) were distributed in two domains (PF00493 and PF14551). During molecular dynamics simulations analysis, consistent fluctuation in RMSD and RMSF values, high Rg and hydrogen bonds in mutant proteins compared to wild-type revealed that these mutations might alter the protein structure and stability of the corresponding protein. Hence, the results from the analyses guide the exploration of the mechanism by which these missense SNPs of the MCM6 gene alter the structural integrity and functional properties of the protein, which could guide the identification of ways to minimize the harmful effects of these mutations in humans.


Materials and methods
The schematic diagram of the in-silico analyses conducted in the study is presented in Fig. 1.

Protein sequence and missense SNPs retrieval
National Center for Biotechnology Information (NCBI) (https:// www.ncbi.nlm.nih.gov/) and NCBI dbSNP (https:// www.ncbi.nlm.nih.gov/ snp/) databases were used to collect the protein sequence (FASTA format) and SNPs of MCM6 gene respectively.The missense SNPs were further analyzed using different software, as this type of mutation generates protein variants and induces crucial structural alterations that could decrease binding affinity and impair the protein function 31 .

Deleterious missense SNPs prediction using SIFT
Sorting Intolerant from Tolerant (SIFT) 32 (https:// sift.bii.a-star.edu.sg/) is a bioinformatic web server used to detect deleterious missense SNPs from tolerated SNPs.This is a homology-based sequence analysis that calculates the normalized probabilities for all possible substitutions from the alignment.In SIFT prediction, missense SNPs having score of > 0.05 regarded as 'tolerated' and less than or equal of that regarded as deleterious.

Functional effect prediction of missense SNPs using SNAP
Screening for non-acceptable polymorphisms (SNAP) 34 (http:// www.rostl ab.org/ servi ces/ SNAP) is a bioinformatics web server used to evaluate the functional effects of a single amino acid substitution in proteins using the neural network method.It predicts the changes that occur due to the missense SNPs on the secondary structure and compares the solvent accessibility of the native and mutated proteins to distinguish them into effect or neutral.The FASTA sequence of the native MCM6 protein was used as the input.

Disease association prediction of missense SNPs
PhD-SNP The predictor of human deleterious single nucleotide polymorphism (PhD-SNP) 35 (http:// snps.biofo ld.org/ phdsnp/ phd-snp.html) is a web server.It is based on Support Vector Machine (SVM) that is optimized to predict disease-related or neutral variants.FASTA sequences of the corresponding proteins and residue changes were submitted as inputs in the PhD-SNP server.

PANTHER
The protein analysis through evolutionary relationship (PANTHER) 36 (https:// panth erdb.org/) based web server was performed to evaluate the effect of the specific amino acid substitution in the biological function of the corresponding protein in the organism.Based on Hidden Markov Models (HMM), this server estimates the probability of how SNPs variants affect the structure of proteins based on their evolutionary origin.

SNPs&GO
The SNPs&GO 37 server (http:// snps-and-go.bioco mp.unibo.it/ snps-andgo/) also utilizes an SVM-based method that precisely predicts if the variants are disease-associated or not.This method calculates the score and evaluates the association of each mutated variant with human diseases.If the score of missense SNPs was ≥ 0.5, it was considered to be involved in the disease, while a score of < 0.5 was considered to have a neutral effect.

Conservation analysis
The ConSurf 38 (http:// consu rf.tau.ac.il/) web server detected a highly conserved functional network of the query protein.This tool creates a phylogenetic tree between homologous sequences to calculate the evolutionary conservation of the amino acids in a protein molecule.The I-Mutant2.0 39 web server (https:// foldi ng.biofo ld.org/ cgi-bin/i-mutan t2.0.cgi) was used to estimate the potential effects of missense SNPs on the structural reliability of the protein and free energy change DDG (Delta Delta G).This is an SVM-based prediction of changes in protein stability upon mutations in the corresponding protein.

Gene of
Here, stability increases when DDG is > 0 kcal/mol and decreases when DDG is < 0 kcal/mol.

MUpro
The MUpro 40 server (http:// mupro.prote omics.ics.uci.edu/) was used to predict the energy change and how mutations affect protein stability using both SVM and Neural Networks methods.A decrease in protein stability was predicted if the confidence score was < 0 while an increase in protein stability was predicted for a score of > 0.
MEDUSA MEDUSA 41 (https:// www.dsimb.inserm.fr/ MEDUSA/) web server was used to predict the flexibility of the corresponding protein.This provides a clear visualization of the prediction results.It predicts two, three and five classes of flexibility by using amino acid sequences.Following the evolutionary origin and physicochemical properties, the server categorized the flexibility class of each amino acid in the spatial arrangement of the protein.
The amino acid sequence was put onto the server in FASTA format to obtain the results.

Protein three-dimensional modeling
The Protein Homology/analogy Recognition Engine V 2.0 (Phyre2) 42 web server (http:// www.sbg.bio.ic.ac.uk/ phyre2) was used to generate the three-dimensional (3D) structure of representative MCM6 and other mutant proteins.The FASTA sequences of the wild-type (WT) and other MCM6 mutant proteins were used to generate 3D structures 43 .The PyMOL 44 software was used to visualize the homology models.

Prediction of harmful mutations using MutPred2
MutPred2 45 (http:// mutpr ed2.mutdb.org/) is a web server that explains the reasons for diseases at the molecular level based on amino acid submissions.It predicts the molecular cause of a disease using a general probability score based on the gain/loss of 14 different structural and functional properties.This score represents the probability that an amino acid substitution is associated with a disease, and the top 5 property scores are provided, where p represents the p-value that certain structural and functional properties are impacted.

Prediction of structural effects of MCM6 mutants using Project Hope server
Project Hope 46 server (http:// www.cmbi.ru.nl/ hope/) was used to calculate the structural and functional effects of point mutations.This investigation provides 3D structural visualization of mutated proteins and provides the results using the UniProt and DAS prediction servers.Here, the protein sequence, wild-type, and new amino acids were used as inputs and the output resulted in text, graphics, and animation format.

Structure-based analysis of mutations using DynaMut2
DynaMut2 47 (https:// biosig.lab.uq.edu.au/ dynam ut2/) was used to evaluate the mutations in protein stability and dynamics using the normal mode analysis (NMA) method.The predicted Gibbs free energy (ΔΔG) values of mutants less than zero (0) were classified as destabilizing, whereas those greater than 0 were classified as stabilizing.

Visualization of selected mutations using mutation3D server
The mutation3D 48 server (http:// mutat ion3d.org/) is a functional prediction and visualization tool for studying the spatial arrangement of amino acid substitutions (AAS) in protein models and structures.This server was used to identify the clusters of amino acid substitutions using the 3D clustering method.It is also useful for clustering other kinds of mutational data, or simply as a tool to quickly assess the relative locations of amino acids in proteins.Additionally, it can be employed to cluster other types of mutational data or as a tool to quickly assess the relative locations of amino acids in proteins.

Molecular dynamics simulations analysis
To evaluate the structural stability of the mutant protein, a 50 ns molecular dynamics (MD) simulation was performed using the "Desmond v6.3 Program" in Schrodinger 2020-3 under the Linux framework 49 .The simulation was performed following the three-site transferrable intermolecular potential (TIP3P) water model 50 .An orthorhombic box shape with a 10 Å distance from the center was used to maintain a specific volume, and Na + and Cl -were added to neutralize the whole system with a salt concentration of 0.15 M.An OPLS3e force field was applied 51 .The protein structure system was further minimized using a natural time and pressure (NPT) ensemble at a constant pressure of 1,01,325 Pascal's and a temperature of 300 K. To evaluate the stability and dynamic characteristics of the protein, RMSD (root means square deviation), RMSF (root means square fluctuation), Rg (radius of gyration), and hydrogen bonds were analyzed.

Gene-gene and protein-protein interaction networks
GeneMANIA Gene-gene interaction network was used to understand the disease phenomenon.The GeneMANIA www.nature.com/scientificreports/new genes in a pathway or complex.The human MCM6 protein sequence was used as input in GeneMANIA.
The analyzed results were based on genetic interactions, pathways, co-expression, co-localization, and shared protein domain similarity.

STRING
Search Tool for the Retrieval of Interacting Genes (STRING) 53 tool (https:// string-db.org/) was used to identify the protein-protein interaction (PPI) of the MCM6 protein with other proteins in the human genome.The PPI network showed correlations between proteins.The PPI network and functional analysis indicated that protein sets were enriched in the target network of the MCM6 protein.

Protein sequence and missense SNPs retrieval
The nsSNPs and sequence of the human MCM6 gene were retrieved from the NCBI database.A total of 15,009 SNPs were identified for the MCM6 gene.The automated computation resulted in 642 missense SNPs (4.28%), 291 synonymous SNPs (1.94%), and 12,500 intron SNPs (83.28%).Then, missense SNPs were further analyzed to identify the most deleterious variants.

Deleterious missense SNPs prediction using SIFT
Among 642 missense SNPs, 33 SNPs were predicted to be deleterious with a tolerance index of ≤ 0.05 (Table 1 and Supplementary Table S1).

Damaging missense SNPs prediction using PolyPhen-2
Based on PolyPhen-2 analysis, 27 and 23 missense SNPs were observed as probably damaging with high confidence in HumDiv and HumVar analyses respectively.Subsequently, 23 were overlapped in both HumDiv and HumVar analyses and were considered for downstream experiments (Table 1).

Functional effect prediction of missense SNPs using SNAP
Analysis of the 23 missense SNPs using SNAP program revealed that all inputted missense SNPs showed a significant effect.However, none of these SNPs was found to be neutral in this analysis (Table 1).www.nature.com/scientificreports/

Conservation analysis
Specific positions of amino acids are crucial for the correct function of a protein.The ConSurf tool was used to determine the conservation score of the MCM6 protein.This program identified highly conserved structural and functional amino acid regions essential for biological functions.The analysis revealed that residues R222C, L449F, D463G, H556Y, R602H and P815T with a conservation score of 9, R207C, V456M, and R633W with a conservation score of 8, and I123S and R658C with a conservation score of 5 had potential biological functions.Interestingly, among those missense SNPs, R222C, R222C, D463G, R602H, and R633W were functional and the rest were buried (Fig. 2).

Stability and flexibility prediction of missense SNPs on MCM6 protein
I-Mutant2.0 and MUpro servers were used to estimate the stability of 11 mutant proteins based on the free energy change value (Delta Delta G, DDG), and the confidence score respectively.Among these, 8 SNPs (I123S, R207C, R222C, V456M, D463G, R602H, R633W, and R658C) guided proteins were found to be the most unstable considering DDG and confidence score (Table 2).
The MEDUSA web server predicted and visualized the flexibility of corresponding proteins with dynamic properties.Based on the three-class flexibility prediction of MCM6 protein by MEDUSA (0 = rigid, 2 = flexible), the positions of all 8 SNPs were rigid, except D463G and R658C which were flexible.The amino acid sequence positions R207C and V456M had a confidence score > 0.5, while the position I123S had the confidence score of 0.5-0.6 and positions R222C, D463G, R602H, R633W and R658C had a confidence score < 0.5 (Fig. 3A).

Protein three-dimensional modeling
The Phyre2 homology-based modeling tool provided the 3D structure of the MCM6 WT and 8 mutant MCM6 proteins (Fig. 3B) 54 .Both WT and mutant proteins showed 100% confidence.The WT and mutant proteins showed 83% coverage of the corresponding proteins.The percentage of the alpha-helix in the MCM6 WT was 29%, whereas it varied from 29 to 31% in mutant proteins.Similarly, the beta strand in the MCM6 WT was 18%, whereas in the mutant protein it varied from 17 to 18%.Similarly, the disorder percentage varied from 22 to 24% (Fig. 3B and Supplementary Table S3).

Prediction of harmful mutations by MutPred2
The MutPred2 tool provides the structural and functional effects of a specific protein based on its physiochemical properties.The eight potentially destabilized missense SNPs caused significant variations in the structural and functional properties of the corresponding MCM6 mutated proteins.Interestingly, these potential SNPs significantly promoted disease pathogenesis by altering various aspects of the protein, including helices and interfaces, DNA binding sites, allosteric sites, catalytic sites, pyrrolidone carboxylic acid, methylation, and transmembrane proteins.Finally, all 8 missense SNPs (I123S, R207C, R222C, V456M, D463G, R602H, R633W, and R658C) resulted in highly harmful mutated proteins, as indicated by the MutPred2 general score (Table 3).

Prediction of structural effects of MCM6 mutants using Project Hope server
The Project HOPE server showed how mutation affects the structural variation of a protein in terms of size, charge, hydrophobicity, and spatial structure compared to the wild type.Among the 8 predicted mutants, the I123S, R207C, R222C, D463G, R602H, and R658C mutant residues were smaller than the wild residues, which caused an empty space in the core of the wild protein and a loss of hydrophobic interactions.However, the V456M and R633W mutants were found to be larger in size than the wild residues, resulting in their localization on the surface of the wild protein.Therefore, mutations of residues can interrupt the inter/intra-molecular interactions of the protein (Table 4).

Structure-based analysis of mutations using DynaMut2
The DynaMut2 server was used for mutation verification and its related effects on the spike protein structure and dynamics.Out of eight mutations, six (I123S, R207C, V456M, D463G, R602H, and R633W) mutations were found to be responsible for destabilizing the protein, whereas the remaining mutants were found to stabilize the protein structure (Table 5).

Radius of gyration (R g )
The Rg quantifies the distribution of atoms around a protein axis, and serves as a crucial metric for forecasting macromolecular structural behavior and evaluating alterations in protein compactness.Here, the complex stability of the wild type and mutated proteins was assessed by analyzing their Rg values throughout a 50 ns simulation period.

Hydrogen bonds
Hydrogen bonds are pivotal for ensuring the binding stability of the corresponding protein.The number of hydrogen bonds can define the protein characteristics, structural stability and ability to bind with other molecules.www.nature.com/scientificreports/Therefore, the number of hydrogen bonds in the wild protein and five mutated proteins (I123S, V456M, D463G, R602H, and R633W) (Fig. 5D).All the proteins formed multiple hydrogen bonds ranging from 670 to 790 in 50 ns simulation time.Higher specificity and less flexibility of hydrogen bonds in wild protein compared to the mutated protein might be due to the higher structural instability of the mutated protein compared to WT (Fig. 5D).

Gene-gene and protein-protein interaction networks
GeneMANIA GeneMANIA constructed a composite gene-gene functional interaction network for the MCM6 gene (Fig. 5).
The MCM6 gene was found to be associated with 20 other genes which play vital roles in various functions.Among these 20 genes, the most important were the MCM2, MCM4, CDC45, MCM7, and CDT1 genes (Fig. 6 and Supplementary Table S4).

STRING
The STRING database was used to partially describe the functional relationships and interaction networks of MCM6 gene.The analysis revealed that MCM6 gene was associated with 10 other genes.A significant correlation    49 .The simulation was performed following the three-site transferrable intermolecular potential (TIP3P) water model 50 .An orthorhombic box shape with a 10 Å distance from the center was used to maintain a specific volume, and Na + and Cl -were added to neutralize the whole system with a salt concentration of 0.15 M.An OPLS3e force field was applied 51 .
rigidity of the mutant protein (Fig. 3A).Mutation-mediated overall flexibility decreases and rigidity increases might affect the binding affinity of the mutant proteins 59 .Amino acid substitution significantly alters the 3D structure and function of the corresponding protein 60 .Alpha helix and beta strands are the structural elements of a protein where the former one represents intramolecular hydrogen (H) bonding and the latter consists of beta sheets.Alpha helix and beta strand significantly differ in mutation tolerance 61 .Alpha helix has potential to accumulate more mutations due to the higher numbers of inter-residue contacts and mutations to residues in β-strands reduce the volume of the amino acid 62 .The 29% alpha-helix in the WT and 29-31% in the mutant, along with 22-24 disordered percentage reveals that the mutant proteins lacked a fixed ordered 3D structure (Fig. 3B and Supplementary Table S3).Mutations affect the helix-helix interactions and disordered regions of corresponding protein by changing their properties 63 .Hence, mutations at these positions might alter the structural and functional properties of the mutant protein.
MutPred2 mediated a significantly higher general score of > 0.75 in all 8 missense SNPs (I123S, R207C, R222C, V456M, D463G, R602H, R633W, and R658C) reveal that these mutations significantly alter the structural and functional properties of the corresponding proteins (Table 3).Since V456M and R633W mutants were larger and the rest were smaller than the WT, these mutations might significantly alter the functional properties of the proteins (Table 4).Mutation of the residue interrupts inter/intra-molecular interactions of the protein which influence the function, characteristics or reactivity of the mutant protein 64 .Therefore, mutations in the residues at the above-mentioned eight positions might alter the function, characteristics or reactivity of the WT protein.Significant destabilization potential in six mutations (I123S, R207C, V456M, D463G, R602H, and R633W) as analyzed by DynaMut2 might alter the structure of the mutant protein leading to loss of function (Table 5).The distribution of five (I123S, V456M, D463G, R602H, and R633W) mutations in MCM (PF00493) and MCM_N (PF14551) domains might be due to the significant contribution of these mutations to the functional alteration of the protein and progression of disease pathogenesis (Fig. 4).Domains are functionally active sites in a protein structure and mutations at these sites may have a tremendous effect on their activity 65 .Hence, the corresponding proteins of these five mutations might have a harmful effect at supra-optimal level in the human genome.The molecular dynamics simulation evaluates how missense SNPs affect the stability, residual fluctuation, and compactness of the protein at different levels.This guide identifies novel mutations in the corresponding protein in an effective way.The protein structure is stabilized when the RMSD and RMSF values of a protein are within 1-3 Å 66 .The consistent fluctuation of the I123S, V456M, D463G, R602H, and R633W mutant in RMSD suggests that these five mutations cause unstable structures of the corresponding protein (Fig. 5A).The RMSF evaluates the mean fluctuation of WT and mutant structures to determine the compactness of the protein.The higher fluctuation in mutants compared to WT was due to the structural instability of the corresponding protein (Fig. 5B).Compared to the WT, a higher Rg value represents the disassociation of the respective protein.However, the increased number of hydrogen bonds causes the structural unsteadiness 67,68 .Hence, consistent fluctuations in RMSD and RMSF values, high Rg and hydrogen bonds in mutant proteins compared to WT might guide exploring the mechanism by which these missense SNPs alter the structure and function of the native MCM6 protein.In the GeneMANIA mediated gene-gene functional interaction network, MCM6 significantly interacted with 20 others genes.Among these, MCM2, MCM4, CDC45 (Cell division cycle 45), MCM7, and CDT1 (Chromatin licensing and DNA replication factor 1) genes interacted more significantly than the other genes (Fig. 6 and Supplementary Table S4).MCM2 regulates cell cycle and DNA replication-related pathways 8 , and while MCM4 acts as the replicative helicase and is required for DNA replication and genome stability 69 .CDC45 is essential for the establishment of an initiation complex at DNA origins 70 , MCM7 is responsible for markedly increased DNA synthesis, cell proliferation and an increased cell invasion in prostate cancer 71 , and CDT1 provides instructions for making a protein that is important in the copying of a cell's DNA before the cell divides 72 .Hence, missense SNPs in MCM6 might alter DNA replication and cell proliferation by interacting with other 20 genes [73][74][75][76][77][78][79][80][81][82][83][84] and resulting in serious health hazards even cancer.
In northern European populations, the MCM6 rs4988235 SNP (commonly referred to as LCT-13910 C/T) is highly correlated with lactase persistence 88,89 .The rs3754686 SNP in MCM6 gene occurs more frequently globally 90 .Numerous studies have been conducted on the lactase persistence-associated genetic variants of the MCM6 gene, including rs145946881, rs869051967, rs41380347, rs4988235 and rs41525747 91 .In children, rs1057031 contributes the most to the development of metabolically unhealthy obesity 92 .The situation is different in the Arabian Peninsula and East Africa, where four different mutations have been found to be associated with lactose persistence.These mutations include rs41525747, rs41380347, rs820486563, and rs145946881, all of which cluster in the MCM6 gene 25,93,94 .
Although different SNP in the MCM6 gene that significantly progress to different disease pathogenies have already been reported, the molecular mechanism of how our identified most significant missense SNPs (I123S, V456M, D463G, R602H, and R633W) cause disease pathogenesis has not yet been discovered.The distribution of these five missense SNPs in the PF00493 and PF14551 domains, which are involved in DNA replication, cell division and cell proliferation reveals that these missense SNPs in the MCM6 might alter their function.Mutated proteins might significantly contribute to the pathogenesis as MCM6 consistently interacted with different genes in the pathways involved in DNA replication and cell division as predicted by GeneMANIA and STRING.The results of the MD simulations also support the findings.Integrating the results of all our analyses on how missense SNPs of the MCM6 gene alter its structural integrity and functional properties, programming based synthetic genetic circuit enabled personalized drugs could be innovated for individuals with missense SNPs of the MCM6 gene in the mentioned positions.This requires a deep analysis of the mentioned missense SNPs along with the integration and application of synthetic biology, machine learning and artificial intelligence under in silico, in vitro, and in vivo conditions.

Conclusion
Here, we identified 642 (4.28%) missense SNPs from 15,009 SNPs for the MCM6 gene.Then, a series of precise bioinformatics analyses were performed to identify the deleterious, probably damaging, effective and diseaseassociated, highly harmful and destabilizing nsSNPs that can alter the structure and function of the MCM6 protein.After a series of analyses, 11 missense common SNPs (I123S, R207C, R222C, L449F, V456M, D463G, H556Y, R602H, R633W, R658C and P815T) were found to be deleterious, probably damaging, affective, and associated with diseases.Subsequently, 8 missense SNPs were found to be highly harmful and significantly contribute to disease pathogenesis.Finally, five mutations (I123S, V456M, D463G, R602H, and R633W) were found to be more harmful since those are located in two domains.Consistent fluctuations in RMSD and RMSF value and high Rg and hydrogen bond in mutant proteins compared to WT during MD simulations reveal that these mutations might alter the protein structure and stability of the WT protein and may have a significant contribution to disease pathogenesis.Considering the impact of these missense SNPs along with their interacting pathways, personalized medicine could be developed to mitigate the harmful effects of these missense SNPs on the diseased individual/population.

Figure 1 .
Figure 1.Schematic representation of the pipeline for in silico analysis of missense SNPs in the MCM6 gene.

Figure 2 .
Figure 2. Analysis of evolutionary conserved amino acids residues of MCM6 by ConSurf.Color coding bar showing conservation score.

Figure 3 .
Figure 3. (A) Flexibility prediction of missense SNPs on MCM6 protein, (B) Three-dimensional modeling of wild type MCM6 protein and other 8 mutants.The protein chains are coloured as blue at the N-terminus and red at the C-terminus.The secondary structures in the model are highlighted in green (for α-helices) and yellow (for β-strands).

Table 4 .
Prediction of how mutation affects MCM6 protein structure using Project Hope server.Grey color indicates the chain of protein, Green colored atoms are the wild amino acid residues, and red are the mutated amino acid residues.

Figure 5 .
Figure 5. Analysis of RMSD, RMSF, Rg, and H bonds of native and five mutants of MCM6 over 50 ns simulation.(A) RMSD values of Cα atoms of native and mutant structures.(B) RMSF values of Cα atoms of native mutant structures.(C) Calculation of Rg, which is represented as a time-dependent change during the simulation.(D) Calculation of H-bonds represented as a time-depend.The color scheme is as follows: native (blue), I123S mutant (red), V456M mutant (green), D463G mutant (violet), R602H mutant (yellow), and R633W mutant (orange).A 50 ns molecular dynamics (MD) simulation was performed using the "Desmond v6.3 Program" in Schrodinger 2020-3 under the Linux framework49 .The simulation was performed following the three-site transferrable intermolecular potential (TIP3P) water model50 .An orthorhombic box shape with a 10 Å distance from the center was used to maintain a specific volume, and Na + and Cl -were added to neutralize the whole system with a salt concentration of 0.15 M.An OPLS3e force field was applied51 .

Figure 7 .
Figure 7. Interaction of MCM6 protein with other proteins analyzed by STRING tool.

Table 2 .
Characterization of the effect of missense SNPs on protein stability.

Table 3 .
Prediction of pathogenicity of missense SNPs in MCM6 protein as predicted by MutPred2.

Table 5 .
Stability prediction of six most significant mutant proteins using DynaMut2.