Impact of Rett Syndrome Mutations on MeCP2 MBD Stability

Nuclease-free precise genome editing corrects MECP2 mutations associated with Rett syndrome

Rett syndrome is an acquired progressive neurodevelopmental disorder caused by de novo mutations in the X-linked MECP2 gene which encodes a pleiotropic protein that functions as a global transcriptional regulator and a chromatin modifier. Rett syndrome predominantly affects heterozygous females while affected male hemizygotes rarely survive. Gene therapy of Rett syndrome has proven challenging due to a requirement for stringent regulation of expression with either over- or under-expression being toxic. Ectopic expression of MECP2 in conjunction with regulatory miRNA target sequences has achieved some success, but the durability of this approach remains unknown. Here we evaluated a nuclease-free homologous recombination (HR)-based genome editing strategy to correct mutations in the MECP2 gene. The stem cell-derived AAVHSCs have previously been shown to mediate seamless and precise HR-based genome editing. We tested the ability of HR-based genome editing to correct pathogenic mutations in Exons 3 and 4 of the MECP2 gene and restore the wild type sequence while preserving all native genomic regulatory elements associated with MECP2 expression, thus potentially addressing a significant issue in gene therapy for Rett syndrome. Moreover, since the mutations are edited directly at the level of the genome, the corrections are expected to be durable with progeny cells inheriting the edited gene. The AAVHSC MECP2 editing vector was designed to be fully homologous to the target MECP2 region and to insert a promoterless Venus reporter at the end of Exon 4. Evaluation of AAVHSC editing in a panel of Rett cell lines bearing mutations in Exons 3 and 4 demonstrated successful correction and rescue of expression of the edited MECP2 gene. Sequence analysis of edited Rett cells revealed successful and accurate correction of mutations in both Exons 3 and 4 and permitted mapping of HR crossover events. Successful correction was observed only when the mutations were flanked at both the 5' and 3' ends by crossover events, but not when both crossovers occurred either exclusively upstream or downstream of the mutation. Importantly, we concluded that pathogenic mutations were successfully corrected in every Rett line analyzed, demonstrating the therapeutic potential of HR-based genome editing.

MeCP2 ubiquitination and sumoylation, in search of a function†

MeCP2 (Methyl CpG binding protein 2) is an intrinsically disordered protein that binds to methylated genome regions. The protein is a critical transcriptional regulator of the brain, and its mutations account for 95% of Rett syndrome (RTT) cases. Early studies of this neurodevelopmental disorder revealed a close connection with dysregulations of the ubiquitin system (UbS), notably as related to UBE3A, a ubiquitin ligase involved in the proteasome-mediated degradation of proteins. MeCP2 undergoes numerous post-translational modifications (PTMs), including ubiquitination and sumoylation, which, in addition to the potential functional outcomes of their monomeric forms in gene regulation and synaptic plasticity, in their polymeric organization, these modifications play a critical role in proteasomal degradation. UbS-mediated proteasomal degradation is crucial in maintaining MeCP2 homeostasis for proper function and is involved in decreasing MeCP2 in some RTT-causing mutations. However, regardless of all these connections to UbS, the molecular details involved in the signaling of MeCP2 for its targeting by the ubiquitin-proteasome system (UPS) and the functional roles of monomeric MeCP2 ubiquitination and sumoylation remain largely unexplored and are the focus of this review.

A Severity Comparison between Italian and Israeli Rett Syndrome Cohorts

Rett syndrome (RTT) is a neurodevelopmental disorder marked by profound cognitive, communication, and motor impairments. Despite identified genotype/phenotype connections, the extent of clinical severity varies even among individuals sharing the same genetic mutation. Diverse sociocultural environments, such as the level of inclusivity of the scholar system, the time spent with family, and the intensity of the rehabilitative intervention provided, might influence their development diversely. This study examines the severity of RTT in people in Italy and Israel, countries with distinct contradictory approaches to caring for those with intricate disabilities, across two age groups. Data from 136 Italian and 59 Israeli girls and women with RTT were assessed and divided into two age categories: above and below 12 years. The RARS, a standardized RTT-specific clinical severity tool, was administered. Despite no differences in age and genetic characteristics, the Italian group showed better scores in the RARS motor and disease-related characteristics areas in both age groups. Moreover, the young Italian participants gathered better total RARS scores and emotional and behavioral characteristics area scores. Furthermore, the young group showed significantly less scoliosis, foot problems, and epilepsy than the older group. These findings endorse the inclusion of girls with RTT in the regular schooling system for a limited daily period, investing in high activity levels within the home and community environments, and suggest continuously providing the person with daily occasions of physical activity and socialization.

Quantitative investigation of the effects of DNA modifications and protein mutations on MeCP2-MBD-DNA interactions

DNA methylation as an important epigenetic marker, has gained attention for the significance of three oxidative modifications (hydroxymethyl-C (hmC), formyl-C (fC), and carboxyl-C (caC)). Mutations occurring in the methyl-CpG-binding domain (MBD) of MeCP2 result in Rett. However, uncertainties persist regarding DNA modification and MBD mutation-induced interaction changes. Here, molecular dynamics simulations were used to investigate the underlying mechanisms behind changes due to different modifications of DNA and MBD mutations. Alanine scanning combined with the interaction entropy method was employed to accurately evaluate the binding free energy. The results show that MBD has the strongest binding ability for mCDNA, followed by caC, hmC, and fCDNA, with the weakest binding ability observed for CDNA. Further analysis revealed that mC modification induces DNA bending, causing residues R91 and R162 closer to the DNA. This proximity enhances van der Waals and electrostatic interactions. Conversely, the caC/hmC and fC modifications lead to two loop regions (near K112 and K130) closer to DNA, respectively. Furthermore, DNA modifications promote the formation of stable hydrogen bond networks, however mutations in the MBD significantly reduce the binding free energy. This study provides detailed insight into the effects of DNA modifications and MBD mutations on binding ability. It emphasizes the necessity for research and development of targeted Rett compounds that induce conformational compatibility between MBD and DNA, enhancing the stability and strength of their interactions.

MeCP2 heterochromatin organization is modulated by arginine methylation and serine phosphorylation

Rett syndrome is a human intellectual disability disorder that is associated with mutations in the X-linked MECP2 gene. The epigenetic reader MeCP2 binds to methylated cytosines on the DNA and regulates chromatin organization. We have shown previously that MECP2 Rett syndrome missense mutations are impaired in chromatin binding and heterochromatin reorganization. Here, we performed a proteomics analysis of post-translational modifications of MeCP2 isolated from adult mouse brain. We show that MeCP2 carries various post-translational modifications, among them phosphorylation on S80 and S421, which lead to minor changes in either heterochromatin binding kinetics or clustering. We found that MeCP2 is (di)methylated on several arginines and that this modification alters heterochromatin organization. Interestingly, we identified the Rett syndrome mutation site R106 as a dimethylation site. In addition, co-expression of protein arginine methyltransferases (PRMT)1 and PRMT6 lead to a decrease of heterochromatin clustering. Altogether, we identified and validated novel modifications of MeCP2 in the brain and show that these can modulate its ability to bind as well as reorganize heterochromatin, which may play a role in the pathology of Rett syndrome.

Electrostatics in Computational Biophysics and Its Implications for Disease Effects

This review outlines the role of electrostatics in computational molecular biophysics and its implication in altering wild-type characteristics of biological macromolecules, and thus the contribution of electrostatics to disease mechanisms. The work is not intended to review existing computational approaches or to propose further developments. Instead, it summarizes the outcomes of relevant studies and provides a generalized classification of major mechanisms that involve electrostatic effects in both wild-type and mutant biological macromolecules. It emphasizes the complex role of electrostatics in molecular biophysics, such that the long range of electrostatic interactions causes them to dominate all other forces at distances larger than several Angstroms, while at the same time, the alteration of short-range wild-type electrostatic pairwise interactions can have pronounced effects as well. Because of this dual nature of electrostatic interactions, being dominant at long-range and being very specific at short-range, their implications for wild-type structure and function are quite pronounced. Therefore, any disruption of the complex electrostatic network of interactions may abolish wild-type functionality and could be the dominant factor contributing to pathogenicity. However, we also outline that due to the plasticity of biological macromolecules, the effect of amino acid mutation may be reduced, and thus a charge deletion or insertion may not necessarily be deleterious.

Genetics behind Cerebral Disease with Ocular Comorbidity: Finding Parallels between the Brain and Eye Molecular Pathology

Cerebral visual impairments (CVIs) is an umbrella term that categorizes miscellaneous visual defects with parallel genetic brain disorders. While the manifestations of CVIs are diverse and ambiguous, molecular diagnostics stand out as a powerful approach for understanding pathomechanisms in CVIs. Nevertheless, the characterization of CVI disease cohorts has been fragmented and lacks integration. By revisiting the genome-wide and phenome-wide association studies (GWAS and PheWAS), we clustered a handful of renowned CVIs into five ontology groups, namely ciliopathies (Joubert syndrome, Bardet–Biedl syndrome, Alstrom syndrome), demyelination diseases (multiple sclerosis, Alexander disease, Pelizaeus–Merzbacher disease), transcriptional deregulation diseases (Mowat–Wilson disease, Pitt–Hopkins disease, Rett syndrome, Cockayne syndrome, X-linked alpha-thalassaemia mental retardation), compromised peroxisome disorders (Zellweger spectrum disorder, Refsum disease), and channelopathies (neuromyelitis optica spectrum disorder), and reviewed several mutation hotspots currently found to be associated with the CVIs. Moreover, we discussed the common manifestations in the brain and the eye, and collated animal study findings to discuss plausible gene editing strategies for future CVI correction.

Oncogenic Mutations in the DNA-Binding Domain of FOXO1 that Disrupt Folding: Quantitative Insights from Experiments and Molecular Simulations

FOXO1, a member of the family of winged-helix motif Forkhead box (FOX) transcription factors, is the most abundantly expressed FOXO member in mature B cells. Sequencing of diffuse large B-cell lymphoma (DLBCL) tumors and cell lines identified specific mutations in the forkhead domain linked to loss of function. Differential scanning calorimetry and thermal shift assays were used to characterize how eight of these mutations affect the stability of the FOX domain. Mutations L183P and L183R were found to be particularly destabilizing. Electrophoresis mobility shift assays show these same mutations also disrupt FOXO1 binding to their canonical DNA sequences, suggesting that the loss of function is due to destabilization of the folded structure. Computational modeling of the effect of mutations on FOXO1 folding was performed using alchemical free energy perturbation (FEP), and a Markov model of the entire folding reaction was constructed from massively parallel molecular simulations, which predicts folding pathways involving the late folding of helix α3. Although FEP can qualitatively predict the destabilization from L183 mutations, we find that a simple hydrophobic transfer model, combined with estimates of unfolded-state solvent-accessible surface areas from molecular simulations, is able to more accurately predict changes in folding free energies due to mutations. These results suggest that the atomic detail provided by simulations is important for the accurate prediction of mutational effects on folding stability. Corresponding disease-associated mutations in other FOX family members support further experimental and computational studies of the folding mechanism of FOX domains.

Evolved DNA Duplex Readers for Strand-Asymmetrically Modified 5-Hydroxymethylcytosine/5-Methylcytosine CpG Dyads

5-Methylcytosine (mC) and 5-hydroxymethylcytosine (hmC), the two main epigenetic modifications of mammalian DNA, exist in symmetric and asymmetric combinations in the two strands of CpG dyads. However, revealing such combinations in single DNA duplexes is a significant challenge. Here, we evolve methyl-CpG-binding domains (MBDs) derived from MeCP2 by bacterial cell surface display, resulting in the first affinity probes for hmC/mC CpGs. One mutant has low nanomolar affinity for a single hmC/mC CpG, discriminates against all 14 other modified CpG dyads, and rivals the selectivity of wild-type MeCP2. Structural studies indicate that this protein has a conserved scaffold and recognizes hmC and mC with two dedicated sets of residues. The mutant allows us to selectively address and enrich hmC/mC-containing DNA fragments from genomic DNA backgrounds. We anticipate that this novel probe will be a versatile tool to unravel the function of hmC/mC marks in diverse aspects of chromatin biology.

Identification of Hotspot Residues in Binding of SARS-CoV-2 Spike and Human ACE2 Proteins

Coronaviruses are a large family of viruses that can cause respiratory infections with varying severity from common cold to severe diseases such as novel coronavirus disease (COVID-19). COVID-19 has been declared as a global pandemic by the World Health Organization on March 11, 2020 and with the development of vaccines it slowed down as of this date. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) uses its spike glycoprotein (Sgp) to bind human angiotensin-converting enzyme 2 (hACE2) receptor, and mediates membrane fusion and virus entry. The recognition of Sgp to human ACE2 and its high affinity for it has been of great importance since this provides the first step in viral entry to human cells. Therefore, it is crucial to identify key residues (hotspots) in this process. In this study, computational alanine scanning has been performed for Sgp and hACE2. The residues identified with significance in binding and other residues in close proximity were studied further through molecular mechanics-based protein binding free energy change prediction methods. Moreover, the interfacial residues in both proteins were investigated for their cooperative binding. Additionally, folding free energy changes upon mutation to Ala were calculated to assess their effect on stability of Sgp and hACE2. Our results taken together with findings from previous studies revealed the residues that are most significant and are relatively significant in binding of Sgp to human ACE2 protein.

A Methylation-Directed, Synthetic Pap Switch Based on Self-Complementary Regulatory DNA Reconstituted in an All E. coli Cell-Free Expression System

Pyelonephritis-associated pili (pap) enable migration of the uropathogenic Escherichia coli strain (UPEC) through the urinary tract. UPEC can switch between a stable 'ON phase' where the corresponding pap genes are expressed and a stable 'OFF phase' where their transcription is repressed. Hereditary DNA methylation of either one of two GATC motives within the regulatory region stabilizes the respective phase over many generations. The underlying molecular mechanism is only partly understood. Previous investigations suggest that in vivo phase-variation stability results from cooperative action of the transcriptional regulators Lrp and PapI. Here, we use an E. coli cell-free expression system to study molecular functions of the pap regulatory region based on a specially designed, synthetic construct flanked by two reporter genes encoding fluorescent proteins for simple readout. On the basis of our observations we suggest that besides Lrp, the conformation of the self-complementary regulatory DNA plays a strong role in the regulation of phase-variation. Our work not only contributes to better understand the phase variation mechanism, but it represents a successful start for mimicking stable, hereditary, and strong expression control based on methylation. The conformation of the regulatory DNA corresponds to a Holliday junction. Gene expression must be expected to respond if opposite arms of the junction are drawn outward.

MeCP2: The Genetic Driver of Rett Syndrome Epigenetics

Mutations in methyl CpG binding protein 2 (MeCP2) are the major cause of Rett syndrome (RTT), a rare neurodevelopmental disorder with a notable period of developmental regression following apparently normal initial development. Such MeCP2 alterations often result in changes to DNA binding and chromatin clustering ability, and in the stability of this protein. Among other functions, MeCP2 binds to methylated genomic DNA, which represents an important epigenetic mark with broad physiological implications, including neuronal development. In this review, we will summarize the genetic foundations behind RTT, and the variable degrees of protein stability exhibited by MeCP2 and its mutated versions. Also, past and emerging relationships that MeCP2 has with mRNA splicing, miRNA processing, and other non-coding RNAs (ncRNA) will be explored, and we suggest that these molecules could be missing links in understanding the epigenetic consequences incurred from genetic ablation of this important chromatin modifier. Importantly, although MeCP2 is highly expressed in the brain, where it has been most extensively studied, the role of this protein and its alterations in other tissues cannot be ignored and will also be discussed. Finally, the additional complexity to RTT pathology introduced by structural and functional implications of the two MeCP2 isoforms (MeCP2-E1 and MeCP2-E2) will be described. Epigenetic therapeutics are gaining clinical popularity, yet treatment for Rett syndrome is more complicated than would be anticipated for a purely epigenetic disorder, which should be taken into account in future clinical contexts.

Ab-initio binding of barnase–barstar with DelPhiForce steered Molecular Dynamics (DFMD) approach

Receptor–ligand interactions are involved in various biological processes, therefore understanding the binding mechanism and ability to predict the binding mode are essential for many biological investigations. While many computational methods exist to predict the 3D structure of the corresponding complex provided the knowledge of the monomers, here we use the newly developed DelPhiForce steered Molecular Dynamics (DFMD) approach to model the binding of barstar to barnase to demonstrate that first-principles methods are also capable of modeling the binding. Essential component of DFMD approach is enhancing the role of long-range electrostatic interactions to provide guiding force of the monomers toward their correct binding orientation and position. Thus, it is demonstrated that the DFMD can successfully dock barstar to barnase even if the initial positions and orientations of both are completely different from the correct ones. Thus, the electrostatics provides orientational guidance along with pulling force to deliver the ligand in close proximity to the receptor.

Computational and structural studies of MeCP2 and associated mutants

Rett Syndrome is a rare genetic disorder exclusively seen in girls. Approximately 95% of RTT cases is caused by mutations in the MeCP2 gene which codes for Methyl-CpG-binding protein 2 (MeCP2). In this review, first, a brief introductory review of Rett Syndrome, MeCP2 protein structure and function, mutation types and frequencies, and phenotype–genotype relationships were provided. After that, the current knowledge on the wild-type and mutant MeCP2 protein structure and dynamics as well as its binding to DNA is reviewed. The review particularly focuses on computational (such as molecular dynamics) and experimental (such as electrophoretic mobility shift assays) studies on the MeCP2 binding to different types of DNA as well as the computational and experimental (such as circular dichroism) studies on the stability changes upon mutations. In the end, a brief opinion on future outlook for further computational studies is provided.

SAAMBE-3D: Predicting Effect of Mutations on Protein-Protein Interactions

Maintaining wild type protein-protein interactions is essential for the normal function of cell and any mutation that alter their characteristics can cause disease. Therefore, the ability to correctly and quickly predict the effect of amino acid mutations is crucial for understanding disease effects and to be able to carry out genome-wide studies. Here, we report a new development of the SAAMBE method, SAAMBE-3D, which is a machine learning-based approach, resulting in accurate predictions and is extremely fast. It achieves the Pearson correlation coefficient ranging from 0.78 to 0.82 depending on the training protocol in benchmarking five-fold validation test against the SKEMPI v2.0 database and outperforms currently existing algorithms on various blind-tests. Furthermore, optimized and tested via five-fold cross-validation on the Cornell University dataset, the SAAMBE-3D achieves AUC of 1.0 and 0.96 on a homo and hereto-dimer test datasets. Another important feature of SAAMBE-3D is that it is very fast, it takes less than a fraction of a second to complete a prediction. SAAMBE-3D is available as a web server and as well as a stand-alone code, the last one being another important feature allowing other researchers to directly download the code and run it on their local computer. Combined all together, SAAMBE-3D is an accurate and fast software applicable for genome-wide studies to assess the effect of amino acid mutations on protein-protein interactions. The webserver and the stand-alone codes (SAAMBE-3D for predicting the change of binding free energy and SAAMBE-3D-DN for predicting if the mutation is disruptive or non-disruptive) are available.

The Molecular Basis of MeCP2 Function in the Brain

MeCP2 is a reader of the DNA methylome that occupies a large proportion of the genome due to its high abundance and the frequency of its target sites. It has been the subject of extensive study because of its link with 'MECP2-related disorders', of which Rett syndrome is the most prevalent. This review integrates evidence from patient mutation data with results of experimental studies using mouse models, cell lines and in vitro systems to critically evaluate our understanding of MeCP2 protein function. Recent evidence challenges the idea that MeCP2 is a multifunctional hub that integrates diverse processes to underpin neuronal function, suggesting instead that its primary role is to recruit the NCoR1/2 co-repressor complex to methylated sites in the genome, leading to dampening of gene expression.

Computational chemistry methods to investigate the effects caused by DNA variants linked with disease

Computational chemistry offers variety of tools to study properties of biological macromolecules. These tools vary in terms of levels of details from quantum mechanical treatment to numerous macroscopic approaches. Here, we provide a review of computational chemistry algorithms and tools for modeling the effects of genetic variations and their association with diseases. Particular emphasis is given on modeling the effects of missense mutations on stability, conformational dynamics, binding, hydrogen bond network, salt bridges, and pH-dependent properties of the corresponding macromolecules. It is outlined that the disease may be caused by alteration of one or several of above-mentioned biophysical characteristics, and a successful prediction of pathogenicity requires detailed analysis of how the alterations affect the function of involved macromolecules. The review provides a short list of most commonly used algorithms to predict the molecular effects of mutations as well.

Plasticity at the DNA recognition site of the MeCP2 mCG-binding domain

MeCP2 is an abundant protein, involved in transcriptional repression by binding to CG and non-CG methylated DNA. However, MeCP2 might also function as a transcription activator as MeCP2 is found bound to sparsely methylated promoters of actively expressed genes. Furthermore, Attachment Region Binding Protein (ARBP), the chicken ortholog of MeCP2, has been reported to bind to Matrix/scaffold attachment regions (MARs/SARs) DNA with an unmethylated 5'-CAC/GTG-3' consensus sequence. In our previous study, although we have systemically measured the binding abilities of MBDs to unmethylated CAC/GTG DNA and the complex structures reveal that the MBD2-MBD (MBD of MBD2) binds to the unmethylated CAC/GTG DNA by recognizing the complementary GTG trinucleotide, how the MeCP2-MBD (MBD of MeCP2) recognizes the unmethylated CAC/GTG DNA, especially the MARs DNA, is still unclear. In this study, we investigated the binding characteristics of MeCP2 in recognizing unmethylated 5'-CAC/GTG-3' motif containing DNA by binding and structural studies. We found that MeCP2-MBD binds to MARs DNA with a comparable binding affinity to mCG DNA, and the MeCP2-CAC/GTG complex structure revealed that MeCP2 residues R111 and R133 form base-specific interactions with the GTG motif. For comparison, we also determined crystal structures of the MeCP2-MBD bound to mCG and mCAC/GTG DNA, respectively. Together, these crystal structures illustrate the adaptability of the MeCP2-MBD toward the GTG motif as well as the mCG DNA, and also provide structural basis of a biological role of MeCP2 as a transcription activator and its disease implications in Rett syndrome.

Tyr120Asp mutation alters domain flexibility and dynamics of MeCP2 DNA binding domain leading to impaired DNA interaction: Atomistic characterization of a Rett syndrome causing mutation

Mutations in the X-linked MECP2 gene represent the main origin of Rett syndrome, causing a profound intellectual disability in females. MeCP2 is an epigenetic transcriptional regulator containing two main functional domains: a methyl-CpG binding domain (MBD) and a transcription repression domain (TRD). Over 600 pathogenic mutations were reported to affect the whole protein; almost half of missense mutations affect the MBD. Understanding the impact of these mutations on the MBD structure and interaction with DNA will foster the comprehension of their pathogenicity and possibly genotype/phenotype correlation studies. Herein, we use molecular dynamics simulations to obtain a detailed view of the dynamics of WT and mutated MBD in the presence and absence of DNA. The pathogenic mutation Y120D is used as paradigm for our studies. Further, since the Y120 residue was previously found to be a phosphorylation site, we characterize the dynamic profile of the MBD also in the presence of Y120 phosphorylation (pY120). We found that addition of a phosphate group to Y120 or mutation in aspartic acid affect domain mobility that samples an alternative conformational space with respect to the WT, leading to impaired ability to interact with DNA. Experimental assays showing a significant reduction in the binding affinity between the mutated MBD and the DNA confirmed our predictions.

Stem Cell Technology for (Epi)genetic Brain Disorders

Despite the enormous efforts of the scientific community over the years, effective therapeutics for many (epi)genetic brain disorders remain unidentified. The common and persistent failures to translate preclinical findings into clinical success are partially attributed to the limited efficiency of current disease models. Although animal and cellular models have substantially improved our knowledge of the pathological processes involved in these disorders, human brain research has generally been hampered by a lack of satisfactory humanized model systems. This, together with our incomplete knowledge of the multifactorial causes in the majority of these disorders, as well as a thorough understanding of associated (epi)genetic alterations, has been impeding progress in gaining more mechanistic insights from translational studies. Over the last years, however, stem cell technology has been offering an alternative approach to study and treat human brain disorders. Owing to this technology, we are now able to obtain a theoretically inexhaustible source of human neural cells and precursors in vitro that offer a platform for disease modeling and the establishment of therapeutic interventions. In addition to the potential to increase our general understanding of how (epi)genetic alterations contribute to the pathology of brain disorders, stem cells and derivatives allow for high-throughput drugs and toxicity testing, and provide a cell source for transplant therapies in regenerative medicine. In the current chapter, we will demonstrate the validity of human stem cell-based models and address the utility of other stem cell-based applications for several human brain disorders with multifactorial and (epi)genetic bases, including Parkinson's disease (PD), Alzheimer's disease (AD), fragile X syndrome (FXS), Angelman syndrome (AS), Prader-Willi syndrome (PWS), and Rett syndrome (RTT).

Binding Analysis of Methyl-CpG Binding Domain of MeCP2 and Rett Syndrome Mutations

Methyl-CpG binding protein 2 (MeCP2) binds to methylated cytosine in CpG island through its methyl-CpG binding domain (MBD). Here, the effects of the Rett syndrome-causing missense mutations on binding affinity of MBD to cytosine (C), methylcytosine (mC), hydroxymethylcytosine (hmC), formylcytosine (fC), and carboxylcytosine (caC) in CpG dinucleotide are investigated. MeCP2-MBD binds to mC-containing variants of double stranded CpG stronger than any other cytosine modified CpG with the strongest affinity to mC/mC. Thirteen MBD missense mutations show reduced binding affinity for mC/mC ranging with a 2-fold decrease for T158M to 88-fold for R111G. The binding affinities of these mutants to C/C are also reduced to various degrees except for T158M. Consistent with free energy perturbation analysis, correlation of binding affinity with protein unfolding allows for grouping mutations into three clusters. Correlation of the first cluster includes mutations that have a higher tendency to unfold and have lesser affinity to mC/mC and C/C. Mutations in the second cluster have similar structural stability but various affinities to mC/mC and C/C. R111G and A140V belong to the third cluster in which the loss of protein flexibility may underlie their reduction in binding affinity to mC/mC and C/C. Most notably, R111 emerges as the key structural element that modulates the specific contacts with mCpG. Implications of the results for the mCpG binding mechanism of MeCP2-MBD are discussed. These analyses provide new insights on the structure and function relationships in MeCP2-MBD and offer new clues to their roles in the pathology of Rett syndrome.

SAAFEC: Predicting the Effect of Single Point Mutations on Protein Folding Free Energy Using a Knowledge-Modified MM/PBSA Approach

Folding free energy is an important biophysical characteristic of proteins that reflects the overall stability of the 3D structure of macromolecules. Changes in the amino acid sequence, naturally occurring or made in vitro, may affect the stability of the corresponding protein and thus could be associated with disease. Several approaches that predict the changes of the folding free energy caused by mutations have been proposed, but there is no method that is clearly superior to the others. The optimal goal is not only to accurately predict the folding free energy changes, but also to characterize the structural changes induced by mutations and the physical nature of the predicted folding free energy changes. Here we report a new method to predict the Single Amino Acid Folding free Energy Changes (SAAFEC) based on a knowledge-modified Molecular Mechanics Poisson-Boltzmann (MM/PBSA) approach. The method is comprised of two main components: a MM/PBSA component and a set of knowledge based terms delivered from a statistical study of the biophysical characteristics of proteins. The predictor utilizes a multiple linear regression model with weighted coefficients of various terms optimized against a set of experimental data. The aforementioned approach yields a correlation coefficient of 0.65 when benchmarked against 983 cases from 42 proteins in the ProTherm database. Availability: the webserver can be accessed via http://compbio.clemson.edu/SAAFEC/.

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.