Whole Genome Sequencing Prioritizes CHEK2, EWSR1, and TIAM1 as Possible Predisposition Genes for Familial Non-Medullary Thyroid Cancer.

Familial inheritance in non-medullary thyroid cancer (NMTC) is an area that has yet to be adequately explored. Despite evidence suggesting strong familial clustering of non-syndromic NMTC, known variants still account for a very small percentage of the genetic burden. In a recent whole genome sequencing (WGS) study of five families with several NMTCs, we shortlisted promising variants with the help of our in-house developed Familial Cancer Variant Prioritization Pipeline (FCVPPv2). Here, we report potentially disease-causing variants in checkpoint kinase 2 (CHEK2), Ewing sarcoma breakpoint region 1 (EWSR1) and T-lymphoma invasion and metastasis-inducing protein 1 (TIAM1) in one family. Performing WGS on three cases, one probable case and one healthy individual in a family with familial NMTC left us with 112254 variants with a minor allele frequency of less than 0.1%, which was reduced by pedigree-based filtering to 6368. Application of the pipeline led to the prioritization of seven coding and nine non-coding variants from this family. The variant identified in CHEK2, a known tumor suppressor gene involved in DNA damage-induced DNA repair, cell cycle arrest, and apoptosis, has been previously identified as a germline variant in breast and prostate cancer and has been functionally validated by Roeb et al. in a yeast-based assay to have an intermediate effect on protein function. We thus hypothesized that this family may harbor additional disease-causing variants in other functionally related genes. We evaluated two further variants in EWSR1 and TIAM1 with promising in silico results and reported interaction in the DNA-damage repair pathway. Hence, we propose a polygenic mode of inheritance in this family. As familial NMTC is considered to be more aggressive than its sporadic counterpart, it is important to identify such susceptibility genes and their associated pathways. In this way, the advancement of personalized medicine in NMTC patients can be fostered. We also wish to reopen the discussion on monogenic vs polygenic inheritance in NMTC and instigate further development in this area of research.

Conformational Dynamics and Cooperativity Drive the Specificity of a Protein-Ligand Interaction.

Molecular recognition is critical for the fidelity of signal transduction in biology. Conversely, the disruption of protein-protein interactions can lead to disease. Thus, comprehension of the molecular determinants of specificity is essential for understanding normal biological signaling processes and for the development of precise therapeutics. Although high-resolution structures have provided atomic details of molecular interactions, much less is known about the influence of cooperativity and conformational dynamics. Here, we used the Tiam2 PSD-95/Dlg/ZO-1 (PDZ) domain and a quadruple mutant (QM), engineered by swapping the identity of four residues important for specificity in the Tiam1 PDZ into the Tiam2 PDZ domain, as a model system to investigate the role of cooperativity and dynamics in PDZ ligand specificity. Surprisingly, equilibrium binding experiments found that the ligand specificity of the Tiam2 QM was switched to that of the Tiam1 PDZ. NMR-based studies indicated that Tiam2 QM PDZ, but not other mutants, had extensive microsecond to millisecond motions distributed throughout the entire domain suggesting structural cooperativity between the mutated residues. Thermodynamic analyses revealed energetic cooperativity between residues in distinct specificity subpockets that was dependent upon the identity of the ligand, indicating a context-dependent binding mechanism. Finally, isothermal titration calorimetry experiments showed distinct entropic signatures along the mutational trajectory from the Tiam2 wild-type to the QM PDZ domain. Collectively, our studies provide unique insights into how structure, conformational dynamics, and thermodynamics combine to modulate ligand-binding specificity and have implications for the evolution, regulation, and design of protein-ligand interactions.

Adaptive landscape flattening in amino acid sequence space for the computational design of protein:peptide binding.

For the high throughput design of protein:peptide binding, one must explore a vast space of amino acid sequences in search of low binding free energies. This complex problem is usually addressed with either simple heuristic scoring or expensive sequence enumeration schemes. Far more efficient than enumeration is a recent Monte Carlo approach that adaptively flattens the energy landscape in sequence space of the unbound peptide and provides formally exact binding free energy differences. The method allows the binding free energy to be used directly as the design criterion. We propose several improvements that allow still more efficient sampling and can address larger design problems. They include the use of Replica Exchange Monte Carlo and landscape flattening for both the unbound and bound peptides. We used the method to design peptides that bind to the PDZ domain of the Tiam1 signaling protein and could serve as inhibitors of its activity. Four peptide positions were allowed to mutate freely. Almost 75 000 peptide variants were processed in two simulations of 109 steps each that used 1 CPU hour on a desktop machine. 96% of the theoretical sequence space was sampled. The relative binding free energies agreed qualitatively with values from experiment. The sampled sequences agreed qualitatively with an experimental library of Tiam1-binding peptides. The main assumption limiting accuracy is the fixed backbone approximation, which could be alleviated in future work by using increased computational resources and multi-backbone designs.

The phosphomimetic mutation of syndecan-4 binds and inhibits Tiam1 modulating Rac1 activity in PDZ interaction-dependent manner.

The small GTPases of the Rho family comprising RhoA, Rac1 and Cdc42 function as molecular switches controlling several essential biochemical pathways in eukaryotic cells. Their activity is cycling between an active GTP-bound and an inactive GDP-bound conformation. The exchange of GDP to GTP is catalyzed by guanine nucleotide exchange factors (GEFs). Here we report a novel regulatory mechanism of Rac1 activity, which is controlled by a phosphomimetic (Ser179Glu) mutant of syndecan-4 (SDC4). SDC4 is a ubiquitously expressed transmembrane, heparan sulfate proteoglycan. In this study we show that the Ser179Glu mutant binds strongly Tiam1, a Rac1-GEF reducing Rac1-GTP by 3-fold in MCF-7 breast adenocarcinoma cells. Mutational analysis unravels the PDZ interaction between SDC4 and Tiam1 is indispensable for the suppression of the Rac1 activity. Neither of the SDC4 interactions is effective alone to block the Rac1 activity, on the contrary, lack of either of interactions can increase the activity of Rac1, therefore the Rac1 activity is the resultant of the inhibitory and stimulatory effects. In addition, SDC4 can bind and tether RhoGDI1 (GDP-dissociation inhibitor 1) to the membrane. Expression of the phosphomimetic SDC4 results in the accumulation of the Rac1-RhoGDI1 complex. Co-immunoprecipitation assays (co-IP-s) reveal that SDC4 can form complexes with RhoGDI1. Together, the regulation of the basal activity of Rac1 is fine tuned and SDC4 is implicated in multiple ways.

Computational Design of the Tiam1 PDZ Domain and Its Ligand Binding.

PDZ domains direct protein-protein interactions and serve as models for protein design. Here, we optimized a protein design energy function for the Tiam1 and Cask PDZ domains that combines a molecular mechanics energy, Generalized Born solvent, and an empirical unfolded state model. Designed sequences were recognized as PDZ domains by the Superfamily fold recognition tool and had similarity scores comparable to natural PDZ sequences. The optimized model was used to redesign the two PDZ domains, by gradually varying the chemical potential of hydrophobic amino acids; the tendency of each position to lose or gain a hydrophobic character represents a novel hydrophobicity index. We also redesigned four positions in the Tiam1 PDZ domain involved in peptide binding specificity. The calculated affinity differences between designed variants reproduced experimental data and suggest substitutions with altered specificities.