Computational and Experimental Approaches to Reveal the Effects of Single Nucleotide Polymorphisms with Respect to Disease Diagnostics

DNA mutations are the cause of many human diseases and they are the reason for natural differences among individuals by affecting the structure, function, interactions, and other properties of DNA and expressed proteins. The ability to predict whether a given mutation is disease-causing or harmless is of great importance for the early detection of patients with a high risk of developing a particular disease and would pave the way for personalized medicine and diagnostics. Here we review existing methods and techniques to study and predict the effects of DNA mutations from three different perspectives: in silico, in vitro and in vivo. It is emphasized that the problem is complicated and successful detection of a pathogenic mutation frequently requires a combination of several methods and a knowledge of the biological phenomena associated with the corresponding macromolecules

Predicting residue contacts using pragmatic correlated mutations method: reducing the false positives

BACKGROUND: Predicting residues' contacts using primary amino acid sequence alone is an important task that can guide 3D structure modeling and can verify the quality of the predicted 3D structures. The correlated mutations (CM) method serves as the most promising approach and it has been used to predict amino acids pairs that are distant in the primary sequence but form contacts in the native 3D structure of homologous proteins. RESULTS: Here we report a new implementation of the CM method with an added set of selection rules (filters). The parameters of the algorithm were optimized against fifteen high resolution crystal structures with optimization criterion that maximized the confidentiality of the predictions. The optimization resulted in a true positive ratio (TPR) of 0.08 for the CM without filters and a TPR of 0.14 for the CM with filters. The protocol was further benchmarked against 65 high resolution structures that were not included in the optimization test. The benchmarking resulted in a TPR of 0.07 for the CM without filters and to a TPR of 0.09 for the CM with filters. CONCLUSION: Thus, the inclusion of selection rules resulted to an overall improvement of 30%. In addition, the pair-wise comparison of TPR for each protein without and with filters resulted in an average improvement of 1.7. The methodology was implemented into a web server that is freely available to the public. The purpose of this implementation is to provide the 3D structure predictors with a tool that can help with ranking alternative models by satisfying the largest number of predicted contacts, as well as it can provide a confidence score for contacts in cases where structure is known

In Silico and In Vitro Investigations of the Mutability of Disease-Causing Missense Mutation Sites in Spermine Synthase

Spermine synthase (SMS) is a key enzyme controlling the concentration of spermidine and spermine in the cell. The importance of SMS is manifested by the fact that single missense mutations were found to cause Snyder-Robinson Syndrome (SRS). At the same time, currently there are no non-synonymous single nucleoside polymorphisms, nsSNPs (harmless mutations), found in SMS, which may imply that the SMS does not tolerate amino acid substitutions, i.e. is not mutable.To investigate the mutability of the SMS, we carried out in silico analysis and in vitro experiments of the effects of amino acid substitutions at the missense mutation sites (G56, V132 and I150) that have been shown to cause SRS. Our investigation showed that the mutation sites have different degree of mutability depending on their structural micro-environment and involvement in the function and structural integrity of the SMS. It was found that the I150 site does not tolerate any mutation, while V132, despite its key position at the interface of SMS dimer, is quite mutable. The G56 site is in the middle of the spectra, but still quite sensitive to charge residue replacement.The performed analysis showed that mutability depends on the detail of the structural and functional factors and cannot be predicted based on conservation of wild type properties alone. Also, harmless nsSNPs can be expected to occur even at sites at which missense mutations were found to cause diseases

ZC4H2, an XLID gene, is required for the generation of a specific subset of CNS interneurons

Miles-Carpenter syndrome (MCS) was described in 1991 as an XLID syndrome with fingertip arches and contractures and mapped to proximal Xq. Patients had microcephaly, short stature, mild spasticity, thoracic scoliosis, hyperextendable MCP joints, rocker-bottom feet, hyperextended elbows and knees. A mutation, p.L66H, in ZC4H2, was identified in a XLID re-sequencing project. Additional screening of linked families and next generation sequencing of XLID families identified three ZC4H2 mutations: p.R18K, p.R213W and p.V75in15aa. The families shared some relevant clinical features. In silico modeling of the mutant proteins indicated all alterations would destabilize the protein. Knockout mutations in zc4h2 were created in zebrafish and homozygous mutant larvae exhibited abnormal swimming, increased twitching, defective eye movement and pectoral fin contractures. Because several of the behavioral defects were consistent with hyperactivity, we examined the underlying neuronal defects and found that sensory neurons and motoneurons appeared normal. However, we observed a striking reduction in GABAergic interneurons. Analysis of cell-type-specific markers showed a specific loss of V2 interneurons in the brain and spinal cord, likely arising from mis-specification of neural progenitors. Injected human wt ZC4H2 rescued the mutant phenotype. Mutant zebrafish injected with human p.L66H or p.R213W mRNA failed to be rescued, while the p.R18K mRNA was able to rescue the interneuron defect. Our findings clearly support ZC4H2 as a novel XLID gene with a required function in interneuron development. Loss of function of ZC4H2 thus likely results in altered connectivity of many brain and spinal circuit

ZC4H2, an XLID gene, is required for the generation of a specific subset of CNS interneurons

Miles-Carpenter syndrome (MCS) was described in 1991 as an XLID syndrome with fingertip arches and contractures and mapped to proximal Xq. Patients had microcephaly, short stature, mild spasticity, thoracic scoliosis, hyperextendable MCP joints, rocker-bottom feet, hyperextended elbows and knees. A mutation, p.L66H, in ZC4H2, was identified in a XLID resequencing project. Additional screening of linked families and next generation sequencing of XLID families identified three ZC4H2 mutations: p.R18K, p.R213W and p.V75in15aa. The families shared some relevant clinical features. In silico modeling of the mutant proteins indicated all alterations would destabilize the protein. Knockout mutations in zc4h2 were created in zebrafish and homozygous mutant larvae exhibited abnormal swimming, increased twitching, defective eye movement and pectoral fin contractures. Because several of the behavioral defects were consistent with hyperactivity, we examined the underlying neuronal defects and found that sensory neurons and motoneurons appeared normal. However, we observed a striking reduction in GABAergic interneurons. Analysis of cell-type-specificmarkers showed a specific loss of V2 interneurons in the brain and spinal cord, likely arising from mis-specification of neural progenitors. Injected human wt ZC4H2 rescued the mutant phenotype. Mutant zebrafish injectedwith human p.L66H or p.R213W mRNA failed to be rescued, while the p.R18K mRNA was able to rescue the interneuron defect. Our findings clearly support ZC4H2 as a novel XLID gene with a required function in interneuron development. Loss of function of ZC4H2 thus likely results in altered connectivity ofmany brain and spinal circuits

Structural, Dynamical, and Energetical Consequences of Rett Syndrome Mutation R133C in MeCP2

Rett Syndrome (RTT) is a progressive neurodevelopmental disease affecting females. RTT is caused by mutations in the MECP2 gene and various amino acid substitutions have been identified clinically in different domains of the multifunctional MeCP2 protein encoded by this gene. The R133C variant in the methylated-CpG-binding domain (MBD) of MeCP2 is the second most common disease-causing mutation in the MBD. Comparative molecular dynamics simulations of R133C mutant and wild-type MBD have been performed to understand the impact of the mutation on structure, dynamics, and interactions of the protein and subsequently understand the disease mechanism. Two salt bridges within the protein and two critical hydrogen bonds between the protein and DNA are lost upon the R133C mutation. The mutation was found to weaken the interaction with DNA and also cause loss of helicity within the 141-144 region. The structural, dynamical, and energetical consequences of R133C mutation were investigated in detail at the atomic resolution. Several important implications of this have been shown regarding protein stability and hydration dynamics as well as its interaction with DNA. The results are in agreement with previous experimental studies and further provide atomic level understanding of the molecular origin of RTT associated with R133C variant

Combining conformational flexibility and continuum electrostatics for calculating pK(a)s in proteins.

Protein stability and function relies on residues being in their appropriate ionization states at physiological pH. In situ residue pK(a)s also provides a sensitive measure of the local protein environment. Multiconformation continuum electrostatics (MCCE) combines continuum electrostatics and molecular mechanics force fields in Monte Carlo sampling to simultaneously calculate side chain ionization and conformation. The response of protein to charges is incorporated both in the protein dielectric constant (epsilon(prot)) of four and by explicit conformational changes. The pK(a) of 166 residues in 12 proteins was determined. The root mean square error is 0.83 pH units, and >90% have errors of <1 pH units whereas only 3% have errors >2 pH units. Similar results are found with crystal and solution structures, showing that the method's explicit conformational sampling reduces sensitivity to the initial structure. The outcome also changes little with protein dielectric constant (epsilon(prot) 4-20). Multiconformation continuum electrostatics titrations show coupling of conformational flexibility and changes in ionization state. Examples are provided where ionizable side chain position (protein G), Asn orientation (lysozyme), His tautomer distribution (RNase A), and phosphate ion binding (RNase A and H) change with pH. Disallowing these motions changes the calculated pK(a)

Mutations in the KDM5C ARID Domain and Their Plausible Association with Syndromic Claes-Jensen-Type Disease

Mutations in KDM5C gene are linked to X-linked mental retardation, the syndromic Claes-Jensen-type disease. This study focuses on non-synonymous mutations in the KDM5C ARID domain and evaluates the effects of two disease-associated missense mutations (A77T and D87G) and three not-yet-classified missense mutations (R108W, N142S, and R179H). We predict the ARID domain’s folding and binding free energy changes due to mutations, and also study the effects of mutations on protein dynamics. Our computational results indicate that A77T and D87G mutants have minimal effect on the KDM5C ARID domain stability and DNA binding. In parallel, the change in the free energy unfolding caused by the mutants A77T and D87G were experimentally measured by urea-induced unfolding experiments and were shown to be similar to the in silico predictions. The evolutionary conservation analysis shows that the disease-associated mutations are located in a highly-conserved part of the ARID structure (N-terminal domain), indicating their importance for the KDM5C function. N-terminal residues’ high conservation suggests that either the ARID domain utilizes the N-terminal to interact with other KDM5C domains or the N-terminal is involved in some yet unknown function. The analysis indicates that, among the non-classified mutations, R108W is possibly a disease-associated mutation, while N142S and R179H are probably harmless