Structure of the human MLH1 N-terminus: implications for predisposition to Lynch syndrome

Special Issue: "Molecular Dynamics Simulations and Structural Analysis of Protein Domains"

The 3D protein structure is the basis for all their biological functions [...].

Identification of a Catalytic Lysine Residue Conserved Among GHKL ATPases: MutL, GyrB, and MORC

DNA mismatch repair endonuclease MutL is a member of GHKL ATPase superfamily. Mutations of MutL homologs are causative of a hereditary cancer, Lynch syndrome. We characterized MutL homologs from human and a hyperthermophile, Aquifex aeolicus, (aqMutL) to reveal the catalytic mechanism for the ATPase activity. Although involvement of a basic residue had not been conceived in the catalytic mechanism, analysis of the pH dependence of the aqMutL ATPase activity revealed that the reaction is catalyzed by a residue with an alkaline pKa. Analyses of mutant aqMutLs showed that Lys79 is the catalytic residue, and the corresponding residues were confirmed to be critical for activities of human MutL homologs, on the basis of which a catalytic mechanism for MutL ATPase is proposed. These and other results described here would contribute to evaluating the pathogenicity of Lynch syndrome-associated missense mutations. Furthermore, it was confirmed that the catalytic lysine residue is conserved among DNA gyrases and microrchidia ATPases, other members of GHKL ATPases, indicating that the catalytic mechanism proposed here is applicable to these members of the superfamily.

Classification of MLH1 Missense VUS Using Protein Structure-Based Deep Learning-Ramachandran Plot-Molecular Dynamics Simulations Method

Pathogenic variation in DNA mismatch repair (MMR) gene MLH1 is associated with Lynch syndrome (LS), an autosomal dominant hereditary cancer. Of the 3798 MLH1 germline variants collected in the ClinVar database, 38.7% (1469) were missense variants, of which 81.6% (1199) were classified as Variants of Uncertain Significance (VUS) due to the lack of functional evidence. Further determination of the impact of VUS on MLH1 function is important for the VUS carriers to take preventive action. We recently developed a protein structure-based method named "Deep Learning-Ramachandran Plot-Molecular Dynamics Simulation (DL-RP-MDS)" to evaluate the deleteriousness of MLH1 missense VUS. The method extracts protein structural information by using the Ramachandran plot-molecular dynamics simulation (RP-MDS) method, then combines the variation data with an unsupervised learning model composed of auto-encoder and neural network classifier to identify the variants causing significant change in protein structure. In this report, we applied the method to classify 447 MLH1 missense VUS. We predicted 126/447 (28.2%) MLH1 missense VUS were deleterious. Our study demonstrates that DL-RP-MDS is able to classify the missense VUS based solely on their impact on protein structure.

Integration of deep learning with Ramachandran plot molecular dynamics simulation for genetic variant classification

Functional classification of genetic variants is a key for their clinical applications in patient care. However, abundant variant data generated by the next-generation DNA sequencing technologies limit the use of experimental methods for their classification. Here, we developed a protein structure and deep learning (DL)-based system for genetic variant classification, DL-RP-MDS, which comprises two principles: 1) Extracting protein structural and thermodynamics information using the Ramachandran plot-molecular dynamics simulation (RP-MDS) method, 2) combining those data with an unsupervised learning model of auto-encoder and a neural network classifier to identify the statistical significance patterns of the structural changes. We observed that DL-RP-MDS provided higher specificity than over 20 widely used in silico methods in classifying the variants of three DNA damage repair genes: TP53, MLH1, and MSH2. DL-RP-MDS offers a powerful platform for high-throughput genetic variant classification. The software and online application are available at https://genemutation.fhs.um.edu.mo/DL-RP-MDS/.

MutS functions as a clamp loader by positioning MutL on the DNA during mismatch repair

Highly conserved MutS and MutL homologs operate as protein dimers in mismatch repair (MMR). MutS recognizes mismatched nucleotides forming ATP-bound sliding clamps, which subsequently load MutL sliding clamps that coordinate MMR excision. Several MMR models envision static MutS-MutL complexes bound to mismatched DNA via a positively charged cleft (PCC) located on the MutL N-terminal domains (NTD). We show MutL-DNA binding is undetectable in physiological conditions. Instead, MutS sliding clamps exploit the PCC to position a MutL NTD on the DNA backbone, likely enabling diffusion-mediated wrapping of the remaining MutL domains around the DNA. The resulting MutL sliding clamp enhances MutH endonuclease and UvrD helicase activities on the DNA, which also engage the PCC during strand-specific incision/excision. These MutS clamp-loader progressions are significantly different from the replication clamp-loaders that attach the polymerase processivity factors β-clamp/PCNA to DNA, highlighting the breadth of mechanisms for stably linking crucial genome maintenance proteins onto DNA.

Isolated MLH1 Loss by Immunohistochemistry Because of Benign Germline MLH1 Polymorphisms

Mismatch repair (MMR) immunohistochemistry (IHC) is frequently used to inform prognosis, select (immuno-)therapy, and identify patients for heritable cancer syndrome testing. However, false-negative and false-positive MMR IHC interpretations have been described.

A Rare MSH2 Variant as a Candidate Marker for Lynch Syndrome II Screening in Tunisia: A Case of Diffuse Gastric Carcinoma

Several syndromic forms of digestive cancers are known to predispose to early-onset gastric tumors such as Hereditary Diffuse Gastric Cancer (HDGC) and Lynch Syndrome (LS). LSII is an extracolonic cancer syndrome characterized by a tumor spectrum including gastric cancer (GC). In the current work, our main aim was to identify the mutational spectrum underlying the genetic predisposition to diffuse gastric tumors occurring in a Tunisian family suspected of both HDGC and LS II syndromes. We selected the index case “JI-021”, which was a woman diagnosed with a Diffuse Gastric Carcinoma and fulfilling the international guidelines for both HDGC and LSII syndromes. For DNA repair, a custom panel targeting 87 candidate genes recovering the four DNA repair pathways was used. Structural bioinformatics analysis was conducted to predict the effect of the revealed variants on the functional properties of the proteins. DNA repair genes panel screening identified two variants: a rare MSH2 c.728G>A classified as a variant with uncertain significance (VUS) and a novel FANCD2 variant c.1879G>T. The structural prediction model of the MSH2 variant and electrostatic potential calculation showed for the first time that MSH2 c.728G>A is likely pathogenic and is involved in the MSH2-MLH1 complex stability. It appears to affect the MSH2-MLH1 complex as well as DNA-complex stability. The c.1879G>T FANCD2 variant was predicted to destabilize the protein structure. Our results showed that the MSH2 p.R243Q variant is likely pathogenic and is involved in the MSH2-MLH1 complex stability, and molecular modeling analysis highlights a putative impact on the binding with MLH1 by disrupting the electrostatic potential, suggesting the revision of its status from VUS to likely pathogenic. This variant seems to be a shared variant in the Mediterranean region. These findings emphasize the importance of testing DNA repair genes for patients diagnosed with diffuse GC with suspicion of LSII and colorectal cancer allowing better clinical surveillance for more personalized medicine.

PMS2 variant results in loss of ATPase activity without compromising mismatch repair

Hereditary cancer syndromes account for approximately 5%-10% of all diagnosed cancer cases. Lynch syndrome (LS) is an autosomal dominant hereditary cancer condition that predisposes individuals to an elevated lifetime risk for developing colorectal, endometrial, and other cancers. LS results from a pathogenic mutation in one of four mismatch repair (MMR) genes (MSH2, MSH6, MLH1, and PMS2). The diagnosis of LS is often challenged by the identification of missense mutations, termed variants of uncertain significance, whose functional effect on the protein is not known. Of the eight PMS2 variants initially selected for this study, we identified a variant within the N-terminal domain where asparagine 335 is mutated to serine, p.Asn335Ser, which lacked ATPase activity, yet appears to be proficient in MMR. To expand our understanding of this functional dichotomy, we performed biophysical and structural studies, and noted that p.Asn335Ser binds to ATP but is unable to hydrolyze it to ADP. To examine the impact of p.Asn335Ser on MMR, we developed a novel in-cell fluorescent-based microsatellite instability reporter that revealed p.Asn335Ser maintained genomic stability. We conclude that in the absence of gross structural changes, PMS2 ATP hydrolysis is not necessary for proficient MMR and that the ATPase deficient p.Asn335Ser variant is likely benign.

Linker domain function predicts pathogenic MLH1 missense variants

The pathogenic consequences of 369 unique human HsMLH1 missense variants has been hampered by the lack of a detailed function in mismatch repair (MMR). Here single-molecule images show that HsMSH2-HsMSH6 provides a platform for HsMLH1-HsPMS2 to form a stable sliding clamp on mismatched DNA. The mechanics of sliding clamp progression solves a significant operational puzzle in MMR and provides explicit predictions for the distribution of clinically relevant HsMLH1 missense mutations.

DNA Mismatch Repair and its Role in Huntington's Disease

DNA mismatch repair (MMR) is a highly conserved genome stabilizing pathway that corrects DNA replication errors, limits chromosomal rearrangements, and mediates the cellular response to many types of DNA damage. Counterintuitively, MMR is also involved in the generation of mutations, as evidenced by its role in causing somatic triplet repeat expansion in Huntington’s disease (HD) and other neurodegenerative disorders. In this review, we discuss the current state of mechanistic knowledge of MMR and review the roles of key enzymes in this pathway. We also present the evidence for mutagenic function of MMR in CAG repeat expansion and consider mechanistic hypotheses that have been proposed. Understanding the role of MMR in CAG expansion may shed light on potential avenues for therapeutic intervention in HD.

A Lynch syndrome-associated mutation at a Bergerat ATP-binding fold destabilizes the structure of the DNA mismatch repair endonuclease MutL

In humans, mutations in genes encoding homologs of the DNA mismatch repair endonuclease MutL cause a hereditary cancer that is known as Lynch syndrome. Here, we determined the crystal structures of the N-terminal domain (NTD) of MutL from the thermophilic eubacterium Aquifex aeolicus (aqMutL) complexed with ATP analogs at 1.69-1.73 Å. The structures revealed significant structural similarities to those of a human MutL homolog, postmeiotic segregation increased 2 (PMS2). We introduced five Lynch syndrome-associated mutations clinically found in human PMS2 into the aqMutL NTD and investigated the protein stability, ATPase activity, and DNA-binding ability of these protein variants. Among the mutations studied, the most unexpected results were obtained for the residue Ser34. Ser34 (Ser46 in PMS2) is located at a previously identified Bergerat ATP-binding fold. We found that the S34I aqMutL NTD retains ATPase and DNA-binding activities. Interestingly, CD spectrometry and trypsin-limited proteolysis indicated the disruption of a secondary structure element of the S34I NTD, destabilizing the overall structure of the aqMutL NTD. In agreement with this, the recombinant human PMS2 S46I NTD was easily digested in the host Escherichia coli cells. Moreover, other mutations resulted in reduced DNA-binding or ATPase activity. In summary, using the thermostable aqMutL protein as a model molecule, we have experimentally determined the effects of the mutations on MutL endonuclease; we discuss the pathological effects of the corresponding mutations in human PMS2.

Probing the DNA-binding center of the MutL protein from the Escherichia coli mismatch repair system via crosslinking and Förster resonance energy transfer

As no crystal structure of full-size MutL bound to DNA has been obtained up to date, in the present work we used crosslinking and Förster resonance energy transfer (FRET) assays for probing the putative DNA-binding center of MutL from Escherichia coli. Several single-cysteine MutL variants (scMutL) were used for site-specific crosslinking or fluorophore modification. The crosslinking efficiency between scMutL proteins and mismatched DNA modified with thiol-reactive probes correlated with the distances from the Cys residues to the DNA calculated from a model of MutS-MutL-DNA complex. FRET-based investigation of DNA binding with different scMutL variants clearly showed that the highest signals were detected for the variants MutL(T218C) and MutL(A251C) indicating closeness of the positions 218 and 251 to DNA in the MutL-DNA complex. Indeed, the Cys218 and Cys251 of scMutL were crosslinked to the reactive DNA with the highest yield demonstrating their proximity to DNA in the MutL-DNA complex. The presence of MutS increased the yield of conjugate formation between the MutL variants and the modified DNA due to tighter MutL-DNA interactions caused by MutS binding to MutL.

Computational and cellular studies reveal structural destabilization and degradation of MLH1 variants in Lynch syndrome

Defective mismatch repair leads to increased mutation rates, and germline loss-of-function variants in the repair component MLH1 cause the hereditary cancer predisposition disorder known as Lynch syndrome. Early diagnosis is important, but complicated by many variants being of unknown significance. Here we show that a majority of the disease-linked MLH1 variants we studied are present at reduced cellular levels. We show that destabilized MLH1 variants are targeted for chaperone-assisted proteasomal degradation, resulting also in degradation of co-factors PMS1 and PMS2. In silico saturation mutagenesis and computational predictions of thermodynamic stability of MLH1 missense variants revealed a correlation between structural destabilization, reduced steady-state levels and loss-of-function. Thus, we suggest that loss of stability and cellular degradation is an important mechanism underlying many MLH1 variants in Lynch syndrome. Combined with analyses of conservation, the thermodynamic stability predictions separate disease-linked from benign MLH1 variants, and therefore hold potential for Lynch syndrome diagnostics.

Stochastic Processes and Component Plasticity Governing DNA Mismatch Repair

DNA mismatch repair (MMR) is a DNA excision-resynthesis process that principally enhances replication fidelity. Highly conserved MutS (MSH) and MutL (MLH/PMS) homologs initiate MMR and in higher eukaryotes act as DNA damage sensors that can trigger apoptosis. MSH proteins recognize mismatched nucleotides, whereas the MLH/PMS proteins mediate multiple interactions associated with downstream MMR events including strand discrimination and strand-specific excision that are initiated at a significant distance from the mismatch. Remarkably, the biophysical functions of the MLH/PMS proteins have been elusive for decades. Here we consider recent observations that have helped to define the mechanics of MLH/PMS proteins and their role in choreographing MMR. We highlight the stochastic nature of DNA interactions that have been visualized by single-molecule analysis and the plasticity of protein complexes that employ thermal diffusion to complete the progressions of MMR.

The MLH1 ATPase domain is needed for suppressing aberrant formation of interstitial telomeric sequences

Genome instability gives rise to cancer. MLH1, commonly known for its important role in mismatch repair (MMR), DNA damage signaling and double-strand break (DSB) repair, safeguards genome stability. Recently we have reported a novel role of MLH1 in preventing aberrant formation of interstitial telomeric sequences (ITSs) at intra-chromosomal regions. Deficiency in MLH1, in particular its N-terminus, leads to an increase of ITSs. Here, we identify that the ATPase activity in the MLH1 N-terminal domain is important for suppressing the formation of ITSs. The ATPase activity is also needed for recruiting MLH1 to DSBs. Moreover, defective ATPase activity of MLH1 causes an increase in micronuclei formation. Our results highlight the crucial role of MLH1's ATPase domain in preventing the aberrant formation of telomeric sequences at the intra-chromosomal regions and preserving genome stability.

Human MLH1 suppresses the insertion of telomeric sequences at intra-chromosomal sites in telomerase-expressing cells

Aberrant formation of interstitial telomeric sequences (ITSs) promotes genome instabilities. However, it is unclear how aberrant ITS formation is suppressed in human cells. Here, we report that MLH1, a key protein involved in mismatch repair (MMR), suppresses telomeric sequence insertion (TSI) at intra-chromosomal regions. The frequency of TSI can be elevated by double-strand break (DSB) inducer and abolished by ATM/ATR inhibition. Suppression of TSI requires MLH1 recruitment to DSBs, indicating that MLH1's role in DSB response/repair is important for suppressing TSI. Moreover, TSI requires telomerase activity but is independent of the functional status of p53 and Rb. Lastly, we show that TSI is associated with chromosome instabilities including chromosome loss, micronuclei formation and chromosome breakage that are further elevated by replication stress. Our studies uncover a novel link between MLH1, telomerase, telomere and genome stability.

Crystal structure and DNA-binding property of the ATPase domain of bacterial mismatch repair endonuclease MutL from Aquifex aeolicus

DNA mismatch repair (MMR) system corrects mismatched bases that are generated mainly by DNA replication errors. The repair system excises the error-containing single-stranded region and enables the re-synthesis of the strand. In the early reactions of MMR, MutL endonuclease incises the newly-synthesized/error-containing strand of the duplex to initiate the downstream excision reaction. MutL endonuclease consists of the N-terminal ATPase and C-terminal endonuclease domains. In this study, we report the crystal structure of the ATPase domain of MutL endonuclease from Aquifex aeolicus. The overall structure of the domain was similar to those of human MutL homologs and Escherichia coli MutL, although E. coli MutL has no endonuclease activity. The ATPase domain was comprised of two subdomains: the N-terminal ATP-binding subdomain and the C-terminal α-β sandwich subdomain. Site-directed mutagenesis experiment identified DNA-interacting eight basic amino acid residues, which were distributed across both the two subdomains and formed a DNA-binding cleft. Docking simulation between the structures of the ATPase and endonuclease domains generated a reliable model structure for the full-length A. aeolicus MutL, which satisfies our previous result of small-angle X-ray scattering analysis. On the basis of the model structure and further experimental results, we concluded that the two separate DNA-binding sites in the full-length A. aeolicus MutL simultaneously bind a dsDNA molecule.

Functional analysis of rare variants in mismatch repair proteins augments results from computation-based predictive methods

The cancer-predisposing Lynch Syndrome (LS) arises from germline mutations in DNA mismatch repair (MMR) genes, predominantly MLH1, MSH2, MSH6, and PMS2. A major challenge for clinical diagnosis of LS is the frequent identification of variants of uncertain significance (VUS) in these genes, as it is often difficult to determine variant pathogenicity, particularly for missense variants. Generic programs such as SIFT and PolyPhen-2, and MMR gene-specific programs such as PON-MMR and MAPP-MMR, are often used to predict deleterious or neutral effects of VUS in MMR genes. We evaluated the performance of multiple predictive programs in the context of functional biologic data for 15 VUS in MLH1, MSH2, and PMS2. Using cell line models, we characterized VUS predicted to range from neutral to pathogenic on mRNA and protein expression, basal cellular viability, viability following treatment with a panel of DNA-damaging agents, and functionality in DNA damage response (DDR) signaling, benchmarking to wild-type MMR proteins. Our results suggest that the MMR gene-specific classifiers do not always align with the experimental phenotypes related to DDR. Our study highlights the importance of complementary experimental and computational assessment to develop future predictors for the assessment of VUS.