Pathogenic variants damage cell composition and single cell transcription in cardiomyopathies.

Pathogenic variants in genes that cause dilated cardiomyopathy (DCM) and arrhythmogenic cardiomyopathy (ACM) convey high risks for the development of heart failure through unknown mechanisms. Using single-nucleus RNA sequencing, we characterized the transcriptome of 880,000 nuclei from 18 control and 61 failing, nonischemic human hearts with pathogenic variants in DCM and ACM genes or idiopathic disease. We performed genotype-stratified analyses of the ventricular cell lineages and transcriptional states. The resultant DCM and ACM ventricular cell atlas demonstrated distinct right and left ventricular responses, highlighting genotype-associated pathways, intercellular interactions, and differential gene expression at single-cell resolution. Together, these data illuminate both shared and distinct cellular and molecular architectures of human heart failure and suggest candidate therapeutic targets.

Cells of the adult human heart.

Cardiovascular disease is the leading cause of death worldwide. Advanced insights into disease mechanisms and therapeutic strategies require a deeper understanding of the molecular processes involved in the healthy heart. Knowledge of the full repertoire of cardiac cells and their gene expression profiles is a fundamental first step in this endeavour. Here, using state-of-the-art analyses of large-scale single-cell and single-nucleus transcriptomes, we characterize six anatomical adult heart regions. Our results highlight the cellular heterogeneity of cardiomyocytes, pericytes and fibroblasts, and reveal distinct atrial and ventricular subsets of cells with diverse developmental origins and specialized properties. We define the complexity of the cardiac vasculature and its changes along the arterio-venous axis. In the immune compartment, we identify cardiac-resident macrophages with inflammatory and protective transcriptional signatures. Furthermore, analyses of cell-to-cell interactions highlight different networks of macrophages, fibroblasts and cardiomyocytes between atria and ventricles that are distinct from those of skeletal muscle. Our human cardiac cell atlas improves our understanding of the human heart and provides a valuable reference for future studies.

Structural-Functional Features of the Thyrotropin Receptor: A Class A G-Protein-Coupled Receptor at Work.

The thyroid-stimulating hormone receptor (TSHR) is a member of the glycoprotein hormone receptors, a sub-group of class A G-protein-coupled receptors (GPCRs). TSHR and its endogenous ligand thyrotropin (TSH) are of essential importance for growth and function of the thyroid gland and proper function of the TSH/TSHR system is pivotal for production and release of thyroid hormones. This receptor is also important with respect to pathophysiology, such as autoimmune (including ophthalmopathy) or non-autoimmune thyroid dysfunctions and cancer development. Pharmacological interventions directly targeting the TSHR should provide benefits to disease treatment compared to currently available therapies of dysfunctions associated with the TSHR or the thyroid gland. Upon TSHR activation, the molecular events conveying conformational changes from the extra- to the intracellular side of the cell across the membrane comprise reception, conversion, and amplification of the signal. These steps are highly dependent on structural features of this receptor and its intermolecular interaction partners, e.g., TSH, antibodies, small molecules, G-proteins, or arrestin. For better understanding of signal transduction, pathogenic mechanisms such as autoantibody action and mutational modifications or for developing new pharmacological strategies, it is essential to combine available structural data with functional information to generate homology models of the entire receptor. Although so far these insights are fragmental, in the past few decades essential contributions have been made to investigate in-depth the involved determinants, such as by structure determination via X-ray crystallography. This review summarizes available knowledge (as of December 2016) concerning the TSHR protein structure, associated functional aspects, and based on these insights we suggest several receptor complex models. Moreover, distinct TSHR properties will be highlighted in comparison to other class A GPCRs to understand the molecular activation mechanisms of this receptor comprehensively. Finally, limitations of current knowledge and lack of information are discussed highlighting the need for intensified efforts toward TSHR structure elucidation.

Molecular dissection of colorectal cancer in pre-clinical models identifies biomarkers predicting sensitivity to EGFR inhibitors.

Colorectal carcinoma represents a heterogeneous entity, with only a fraction of the tumours responding to available therapies, requiring a better molecular understanding of the disease in precision oncology. To address this challenge, the OncoTrack consortium recruited 106 CRC patients (stages I-IV) and developed a pre-clinical platform generating a compendium of drug sensitivity data totalling >4,000 assays testing 16 clinical drugs on patient-derived in vivo and in vitro models. This large biobank of 106 tumours, 35 organoids and 59 xenografts, with extensive omics data comparing donor tumours and derived models provides a resource for advancing our understanding of CRC. Models recapitulate many of the genetic and transcriptomic features of the donors, but defined less complex molecular sub-groups because of the loss of human stroma. Linking molecular profiles with drug sensitivity patterns identifies novel biomarkers, including a signature outperforming RAS/RAF mutations in predicting sensitivity to the EGFR inhibitor cetuximab.

Genomics and drug profiling of fatal TCF3-HLF-positive acute lymphoblastic leukemia identifies recurrent mutation patterns and therapeutic options.

TCF3-HLF-positive acute lymphoblastic leukemia (ALL) is currently incurable. Using an integrated approach, we uncovered distinct mutation, gene expression and drug response profiles in TCF3-HLF-positive and treatment-responsive TCF3-PBX1-positive ALL. We identified recurrent intragenic deletions of PAX5 or VPREB1 in constellation with the fusion of TCF3 and HLF. Moreover somatic mutations in the non-translocated allele of TCF3 and a reduction of PAX5 gene dosage in TCF3-HLF ALL suggest cooperation within a restricted genetic context. The enrichment for stem cell and myeloid features in the TCF3-HLF signature may reflect reprogramming by TCF3-HLF of a lymphoid-committed cell of origin toward a hybrid, drug-resistant hematopoietic state. Drug response profiling of matched patient-derived xenografts revealed a distinct profile for TCF3-HLF ALL with resistance to conventional chemotherapeutics but sensitivity to glucocorticoids, anthracyclines and agents in clinical development. Striking on-target sensitivity was achieved with the BCL2-specific inhibitor venetoclax (ABT-199). This integrated approach thus provides alternative treatment options for this deadly disease.

Recurrent somatic alterations of FGFR1 and NTRK2 in pilocytic astrocytoma.

Pilocytic astrocytoma, the most common childhood brain tumor, is typically associated with mitogen-activated protein kinase (MAPK) pathway alterations. Surgically inaccessible midline tumors are therapeutically challenging, showing sustained tendency for progression and often becoming a chronic disease with substantial morbidities. Here we describe whole-genome sequencing of 96 pilocytic astrocytomas, with matched RNA sequencing (n = 73), conducted by the International Cancer Genome Consortium (ICGC) PedBrain Tumor Project. We identified recurrent activating mutations in FGFR1 and PTPN11 and new NTRK2 fusion genes in non-cerebellar tumors. New BRAF-activating changes were also observed. MAPK pathway alterations affected all tumors analyzed, with no other significant mutations identified, indicating that pilocytic astrocytoma is predominantly a single-pathway disease. Notably, we identified the same FGFR1 mutations in a subset of H3F3A-mutated pediatric glioblastoma with additional alterations in the NF1 gene. Our findings thus identify new potential therapeutic targets in distinct subsets of pilocytic astrocytoma and childhood glioblastoma.

Development of immune-specific interaction potentials and their application in the multi-agent-system VaccImm.

Peptide vaccination in cancer therapy is a promising alternative to conventional methods. However, the parameters for this personalized treatment are difficult to access experimentally. In this respect, in silico models can help to narrow down the parameter space or to explain certain phenomena at a systems level. Herein, we develop two empirical interaction potentials specific to B-cell and T-cell receptor complexes and validate their applicability in comparison to a more general potential. The interaction potentials are applied to the model VaccImm which simulates the immune response against solid tumors under peptide vaccination therapy. This multi-agent system is derived from another immune system simulator (C-ImmSim) and now includes a module that enables the amino acid sequence of immune receptors and their ligands to be taken into account. The multi-agent approach is combined with approved methods for prediction of major histocompatibility complex (MHC)-binding peptides and the newly developed interaction potentials. In the analysis, we critically assess the impact of the different modules on the simulation with VaccImm and how they influence each other. In addition, we explore the reasons for failures in inducing an immune response by examining the activation states of the immune cell populations in detail.In summary, the present work introduces immune-specific interaction potentials and their application to the agent-based model VaccImm which simulates peptide vaccination in cancer therapy.

GPCR-SSFE: a comprehensive database of G-protein-coupled receptor template predictions and homology models.

BACKGROUND: G protein-coupled receptors (GPCRs) transduce a wide variety of extracellular signals to within the cell and therefore have a key role in regulating cell activity and physiological function. GPCR malfunction is responsible for a wide range of diseases including cancer, diabetes and hyperthyroidism and a large proportion of drugs on the market target these receptors. The three dimensional structure of GPCRs is important for elucidating the molecular mechanisms underlying these diseases and for performing structure-based drug design. Although structural data are restricted to only a handful of GPCRs, homology models can be used as a proxy for those receptors not having crystal structures. However, many researchers working on GPCRs are not experienced homology modellers and are therefore unable to benefit from the information that can be gleaned from such three-dimensional models. Here, we present a comprehensive database called the GPCR-SSFE, which provides initial homology models of the transmembrane helices for a large variety of family A GPCRs. DESCRIPTION: Extending on our previous theoretical work, we have developed an automated pipeline for GPCR homology modelling and applied it to a large set of family A GPCR sequences. Our pipeline is a fragment-based approach that exploits available family A crystal structures. The GPCR-SSFE database stores the template predictions, sequence alignments, identified sequence and structure motifs and homology models for 5025 family A GPCRs. Users are able to browse the GPCR dataset according to their pharmacological classification or search for results using a UniProt entry name. It is also possible for a user to submit a GPCR sequence that is not contained in the database for analysis and homology model building. The models can be viewed using a Jmol applet and are also available for download along with the alignments. CONCLUSIONS: The data provided by GPCR-SSFE are useful for investigating general and detailed sequence-structure-function relationships of GPCRs, performing structure-based drug design and for better understanding the molecular mechanisms underlying disease-associated mutations in GPCRs. The effectiveness of our multiple template and fragment approach is demonstrated by the accuracy of our predicted homology models compared to recently published crystal structures.

From molecular details of the interplay between transmembrane helices of the thyrotropin receptor to general aspects of signal transduction in family a G-protein-coupled receptors (GPCRs).

Transmembrane helices (TMHs) 5 and 6 are known to be important for signal transduction by G-protein-coupled receptors (GPCRs). Our aim was to characterize the interface between TMH5 and TMH6 of the thyrotropin receptor (TSHR) to gain molecular insights into aspects of signal transduction and regulation. A proline at TMH5 position 5.50 is highly conserved in family A GPCRs and causes a twist in the helix structure. Mutation of the TSHR-specific alanine (Ala-5935·5°) at this position to proline resulted in a 20-fold reduction of cell surface expression. This indicates that TMH5 in the TSHR might have a conformation different from most other family A GPCRs by forming a regular α-helix. Furthermore, linking our own and previous data from directed mutagenesis with structural information led to suggestions of distinct pairs of interacting residues between TMH5 and TMH6 that are responsible for stabilizing either the basal or the active state. Our insights suggest that the inactive state conformation is constrained by a core set of polar interactions among TMHs 2, 3, 6, and 7 and in contrast that the active state conformation is stabilized mainly by non-polar interactions between TMHs 5 and 6. Our findings might be relevant for all family A GPCRs as supported by a statistical analysis of residue properties between the TMHs of a vast number of GPCR sequences.

Meet me halfway: when genomics meets structural bioinformatics.

The DNA sequencing technology developed by Frederick Sanger in the 1970s established genomics as the basis of comparative genetics. The recent invention of next-generation sequencing (NGS) platform has added a new dimension to genome research by generating ultra-fast and high-throughput sequencing data in an unprecedented manner. The advent of NGS technology also provides the opportunity to study genetic diseases where sequence variants or mutations are sought to establish a causal relationship with disease phenotypes. However, it is not a trivial task to seek genetic variants responsible for genetic diseases and even harder for complex diseases such as diabetes and cancers. In such polygenic diseases, multiple genes and alleles, which can exist in healthy individuals, come together to contribute to common disease phenotypes in a complex manner. Hence, it is desirable to have an approach that integrates omics data with both knowledge of protein structure and function and an understanding of networks/pathways, i.e. functional genomics and systems biology; in this way, genotype-phenotype relationships can be better understood. In this review, we bring this 'bottom-up' approach alongside the current NGS-driven genetic study of genetic variations and disease aetiology. We describe experimental and computational techniques for assessing genetic variants and their deleterious effects on protein structure and function.

Principles and determinants of G-protein coupling by the rhodopsin-like thyrotropin receptor.

In this study we wanted to gain insights into selectivity mechanisms between G-protein-coupled receptors (GPCR) and different subtypes of G-proteins. The thyrotropin receptor (TSHR) binds G-proteins promiscuously and activates both Gs (cAMP) and Gq (IP). Our goal was to dissect selectivity patterns for both pathways in the intracellular region of this receptor. We were particularly interested in the participation of poorly investigated receptor parts.We systematically investigated the amino acids of intracellular loop (ICL) 1 and helix 8 using site-directed mutagenesis alongside characterization of cAMP and IP accumulation. This approach was guided by a homology model of activated TSHR in complex with heterotrimeric Gq, using the X-ray structure of opsin with a bound G-protein peptide as a structural template.We provide evidence that ICL1 is significantly involved in G-protein activation and our model suggests potential interactions with subunits G alpha as well as G betagamma. Several amino acid substitutions impaired both IP and cAMP accumulation. Moreover, we found a few residues in ICL1 (L440, T441, H443) and helix 8 (R687) that are sensitive for Gq but not for Gs activation. Conversely, not even one residue was found that selectively affects cAMP accumulation only. Together with our previous mutagenesis data on ICL2 and ICL3 we provide here the first systematically completed map of potential interfaces between TSHR and heterotrimeric G-protein. The TSHR/Gq-heterotrimer complex is characterized by more selective interactions than the TSHR/Gs complex. In fact the receptor interface for binding Gs is a subset of that for Gq and we postulate that this may be true for other GPCRs coupling these G-proteins. Our findings support that G-protein coupling and preference is dominated by specific structural features at the intracellular region of the activated GPCR but is completed by additional complementary recognition patterns between receptor and G-protein subtypes.

Signaling-sensitive amino acids surround the allosteric ligand binding site of the thyrotropin receptor.

The thyrotropin receptor [thyroid-stimulating hormone receptor (TSHR)], a G-protein-coupled receptor (GPCR), is endogenously activated by thyrotropin, which binds to the extracellular region of the receptor. We previously identified a low-molecular-weight (LMW) agonist of the TSHR and predicted its allosteric binding pocket within the receptor's transmembrane domain. Because binding of the LMW agonist probably disrupts interactions or leads to formation of new interactions among amino acid residues surrounding the pocket, we tested whether mutation of residues at these positions would lead to constitutive signaling activity. Guided by molecular modeling, we performed site-directed mutagenesis of 24 amino acids in this spatial region, followed by functional characterization of the mutant receptors in terms of expression and signaling, measured as cAMP accumulation. We found that mutations V421I, Y466A, T501A, L587V, M637C, M637W, S641A, Y643F, L645V, and Y667A located in several helices exhibit constitutive activity. Of note is mutation M637W at position 6.48 in transmembrane helix 6, which has a significant effect on the interaction of the receptor with the LMW agonist. In summary, we found that a high proportion of residues in several helices surrounding the allosteric binding site of LMW ligands in the TSHR when mutated lead to constitutively active receptors. Our findings of signaling-sensitive residues in this region of the transmembrane bundle may be of general importance as this domain appears to be evolutionarily retained among GPCRs.

Structural bioinformatics mutation analysis reveals genotype-phenotype correlations in von Hippel-Lindau disease and suggests molecular mechanisms of tumorigenesis.

Mutations in the VHL gene lead to von Hippel-Lindau (VHL) disease, a clinically heterogeneous cancer syndrome. Here, we use software and database tools to understand and predict the phenotypes associated with missense mutations in the VHL gene product, pVHL. The protein product pVHL is known to interact with elongin B, elongin C, and the HIF substrate. By analyzing known and predicted interaction sites and predictions of thermodynamic stability change upon mutation, we generate new hypotheses regarding the molecular etiology of renal cell carcinoma (RCC) and pheochromocytoma (PCC) in VHL disease. We find that the molecular causes of RCC and PCC appear to be decoupled. RCC may arise through two distinct mechanisms: disruption of HIF interactions or binding at the elongin B interface. PCC is triggered by mutations which disrupt interactions at the elongin C binding site. These findings have important implications for VHL disease and for nonfamilial RCC, because most cases of clear cell RCC are linked with VHL inactivation. Additionally, predicting effects of genetic variation will be critical as genetic sequencing accelerates; the analytical strategy presented here may elucidate other systems as further data on genetic variation become available.

AIRE's CARD revealed, a new structure for central tolerance provokes transcriptional plasticity.

Developing T cells encounter peripheral self-antigens in the thymus in order to delete autoreactive clones. It is now known that the autoimmune regulator protein (AIRE), which is expressed in thymic medullary epithelial cells, plays a key role in regulating the thymic transcription of these peripheral tissue-specific antigens. Mutations in the AIRE gene are associated with a severe multiorgan autoimmune syndrome (APECED), and autoimmune reactivities are manifest in AIRE-deficient mice. Functional AIRE protein is expressed as distinct nuclear puncta, although no structural basis existed to explain their relevance to disease. In addressing the cell biologic basis for APECED, we made the unexpected discovery that an AIRE mutation hot spot lies in a caspase recruitment domain. Combined homology modeling and in vitro data now show how APECED mutations influence the activity of this transcriptional regulator. We also provide novel in vivo evidence for AIRE's association with a global transcription cofactor, which may underlie AIRE's focal, genome-wide, alteration of the transcriptome.