Exome Sequencing has a high diagnostic rate in sporadic congenital hypopituitarism and reveals novel candidate genes.

CONTEXT: The pituitary gland is key for childhood growth, puberty, and metabolism. Pituitary dysfunction is associated with a spectrum of phenotypes, from mild to severe. Congenital Hypopituitarism (CH) is the most commonly reported pediatric endocrine dysfunction with an incidence of 1:4000, yet low rates of genetic diagnosis have been reported. OBJECTIVE: We aimed to unveil the genetic etiology of CH in a large cohort of patients from Argentina. METHODS: We performed whole exome sequencing of 137 unrelated cases of CH, the largest cohort examined with this method to date. RESULTS: Of the 137 cases, 19.1% and 16% carried pathogenic or likely pathogenic variants in known and new genes, respectively, while 28.2% carried variants of uncertain significance. This high yield was achieved through the integration of broad gene panels (genes described in animal models and/or other disorders), an unbiased candidate gene screen with a new bioinformatics pipeline (including genes high loss of function intolerance), and analysis of copy number variants. Three novel findings emerged. First, the most prevalent affected gene encodes the cell adhesion factor ROBO1. Affected children had a spectrum of phenotypes, consistent with a role beyond pituitary stalk interruption syndrome. Second, we found that CHD7 mutations also produce a phenotypic spectrum, not always associated with full CHARGE syndrome. Third, we add new evidence of pathogenicity in the genes PIBF1 and TBC1D32, and report 13 novel candidate genes associated with CH (e.g. PTPN6, ARID5B). CONCLUSION: Overall, these results provide an unprecedented insight into the diverse genetic etiology of hypopituitarism.

Bacterial cytochrome P450s: a bioinformatics odyssey of substrate discovery.

Bacterial P450 cytochromes (BacCYPs) are versatile heme-containing proteins responsible for oxidation reactions on a wide range of substrates, contributing to the production of valuable natural products with limitless biotechnological potential. While the sequencing of microbial genomes has provided a wealth of BacCYP sequences, functional characterization lags behind, hindering our understanding of their roles. This study employs a comprehensive approach to predict BacCYP substrate specificity, bridging the gap between sequence and function. We employed an integrated approach combining sequence and functional data analysis, genomic context exploration, 3D structural modeling with molecular docking, and phylogenetic clustering. The research begins with an in-depth analysis of BacCYP sequence diversity and structural characteristics, revealing conserved motifs and recurrent residues in the active site. Phylogenetic analysis identifies distinct groups within the BacCYP family based on sequence similarity. However, our study reveals that sequence alone does not consistently predict substrate specificity, necessitating additional perspectives. The study delves into the genetic context of BacCYPs, utilizing neighboring gene information to infer potential substrates, a method proven very effective in many cases. Molecular docking is employed to assess BacCYP-substrate interactions, confirming potential substrates and providing insights into selectivity. Finally, a comprehensive strategy is proposed for predicting BacCYP substrates, involving all the evaluated approaches. The effectiveness of this strategy is demonstrated with two case studies, highlighting its potential for substrate discovery.

Globin phylogeny, evolution and function, the newest update.

Our globin census update allows us to refine our vision of globin origin, evolution, and structure to function relationship in the context of the currently accepted tree of life. The modern globin domain originates as a single domain, three-over-three α-helical folded structure before the diversification of the kingdoms of life (Bacteria, Archaea, Eukarya). Together with the diversification of prokaryotes, three monophyletic globin families (M, S, and T) emerged, most likely in Proteobacteria and Actinobacteria, displaying specific sequence and structural features, and spread by vertical and horizontal gene transfer, most probably already present in the last universal common ancestor (LUCA). Non-globin domains were added, and eventually lost again, creating multi-domain structures in key branches of M- (FHb and Adgb) and the vast majority of S globins, which with their coevolved multi-domain architectures, have predominantly "sensor" functions. Single domain T-family globins diverged into four major groups and most likely display functions related to reactive nitrogen and oxygen species (RNOS) chemistry, as well as oxygen storage/transport which drives the evolution of its major branches with their characteristic key distal residues (B10, E11, E7, and G8). M-family evolution also lead to distinctive major types (FHb and Fgb, Ngb, Adgb, GbX vertebrate Gbs), and shows the shift from high oxygen affinity controlled by TyrB10-Gln/AsnE11 likely related to RNOS chemistry in microorganisms, to a moderate oxygen affinity storage/transport function controlled by hydrophobic B10/E11-HisE7 in multicellular animals.

Gp29 LysA of mycobacteriophage TM4 can hydrolyze peptidoglycan through an N-acetyl-muramoyl-L-alanine amidase activity.

Bacteriophage endolysins are crucial for progeny release at the end of the lytic cycle. Mycobacteriophage's genomes carry a lysin A essential gene, whose product cleaves the peptidoglycan (PG) layer and a lysin B, coding for an esterase, that cleaves the linkage between the mycolic acids and the arabinogalactan-PG complex. Lysin A mycobacteriophage proteins are highly modular and in gp29 (LysA) of phage TM4 three distinctive domains were identified. By bioinformatics analysis the central module was previously found to be similar to an amidase-2 domain family with an N-acetylmuramoyl -L-alanine amidase activity. We demonstrated experimentally that purified LysA is able to lyse a suspension of Micrococcus lysodeikticus and can promote cell lysis when expressed in E. coli and Mycobacterium smegmatis. After incubation of LysA with MDP (Muramyl dipeptide, N-acetyl-muramyl-L-alanyl-D-isoglutamine) we detected the presence of N-acetylmuramic acid (NAcMur) and L-Ala- D- isoGlutamine (L-Ala-D-isoGln) corroborating the proposed muramidase activity of this enzyme. This protein was stabilized at acidic pH in the presence of Zn consistent with the increase of the enzymatic activity under these conditions. By homology modeling, we predicted that the Zn ion is coordinated by His 226, His 335, and Asp 347 and we also identified the amino acid Glu 290 as the catalytic residue. LysA activity was completely abolished in derived mutants on these key residues, suggesting that the PG hydrolysis solely relies on the central domain of the protein.

A Novel SLC5A5 Variant Reveals the Crucial Role of Kinesin Light Chain 2 in Thyroid Hormonogenesis.

CONTEXT: Iodide transport defect (ITD) (Online Mendelian Inheritance in Man No. 274400) is an uncommon cause of dyshormonogenic congenital hypothyroidism due to loss-of-function variants in the SLC5A5 gene, which encodes the sodium/iodide symporter (NIS), causing deficient iodide accumulation in thyroid follicular cells. OBJECTIVE: This work aims to determine the molecular basis of a patient's ITD clinical phenotype. METHODS: The propositus was diagnosed with dyshormonogenic congenital hypothyroidism with minimal 99mTc-pertechnetate accumulation in a eutopic thyroid gland. The propositus SLC5A5 gene was sequenced. Functional in vitro characterization of the novel NIS variant was performed. RESULTS: Sanger sequencing revealed a novel homozygous missense p.G561E NIS variant. Mechanistically, the G561E substitution reduces iodide uptake, because targeting of G561E NIS to the plasma membrane is reduced. Biochemical analyses revealed that G561E impairs the recognition of an adjacent tryptophan-acidic motif by the kinesin-1 subunit kinesin light chain 2 (KLC2), interfering with NIS maturation beyond the endoplasmic reticulum, and reducing iodide accumulation. Structural bioinformatic analysis suggests that G561E shifts the equilibrium of the unstructured tryptophan-acidic motif toward a more structured conformation unrecognizable to KLC2. Consistently, knockdown of Klc2 causes defective NIS maturation and consequently decreases iodide accumulation in rat thyroid cells. Morpholino knockdown of klc2 reduces thyroid hormone synthesis in zebrafish larvae leading to a hypothyroid state as revealed by expression profiling of key genes related to the hypothalamic-pituitary-thyroid axis. CONCLUSION: We report a novel NIS pathogenic variant associated with dyshormonogenic congenital hypothyroidism. Detailed molecular characterization of G561E NIS uncovered the significance of KLC2 in thyroid physiology.

Comprehensive Identification of Pathogenic Gene Variants in Patients With Neuroendocrine Disorders.

PURPOSE: Congenital hypopituitarism (CH) can present in isolation or with other birth defects. Mutations in multiple genes can cause CH, and the use of a genetic screening panel could establish the prevalence of mutations in known and candidate genes for this disorder. It could also increase the proportion of patients that receive a genetic diagnosis. METHODS: We conducted target panel genetic screening using single-molecule molecular inversion probes sequencing to assess the frequency of mutations in known hypopituitarism genes and new candidates in Argentina. We captured genomic deoxyribonucleic acid from 170 pediatric patients with CH, either alone or with other abnormalities. We performed promoter activation assays to test the functional effects of patient variants in LHX3 and LHX4. RESULTS: We found variants classified as pathogenic, likely pathogenic, or with uncertain significance in 15.3% of cases. These variants were identified in known CH causative genes (LHX3, LHX4, GLI2, OTX2, HESX1), in less frequently reported genes (FOXA2, BMP4, FGFR1, PROKR2, PNPLA6) and in new candidate genes (BMP2, HMGA2, HNF1A, NKX2-1). CONCLUSION: In this work, we report the prevalence of mutations in known CH genes in Argentina and provide evidence for new candidate genes. We show that CH is a genetically heterogeneous disease with high phenotypic variation and incomplete penetrance, and our results support the need for further gene discovery for CH. Identifying population-specific pathogenic variants will improve the capacity of genetic data to predict eventual clinical outcomes.

The Transcription Factor NF-κB Mediates Thyrotropin-Stimulated Expression of Thyroid Differentiation Markers.

Background: The nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB) transcription factor is a key regulator of cell survival, proliferation, and gene expression. Although activation of NF-κB signaling in thyroid follicular cells after thyrotropin (TSH) receptor (TSHR) engagement has been reported, the downstream signaling leading to NF-κB activation remains unexplored. Here, we sought to elucidate the mechanisms that regulate NF-κB signaling activation in response to TSH stimulation. Methods: Fisher rat-derived thyroid cell lines and primary cultures of NF-κB essential modulator (NEMO)-deficient mice thyrocytes were used as models. Signaling pathways leading to the activation of NF-κB were investigated by using chemical inhibitors and phospho-specific antibodies. Luciferase reporter gene assays and site-directed mutagenesis were used to monitor NF-κB-dependent gene transcriptional activity and the expression of thyroid differentiation markers was assessed by reverse transcription quantitative polymerase chain reaction and Western blot, respectively. Chromatin immunoprecipitation (ChIP) was carried out to investigate NF-κB subunit p65 DNA binding, and small interfering RNA (siRNA)-mediated gene knockdown approaches were used for studying gene function. Results: Using thyroid cell lines, we observed that TSH treatment leads to protein kinase C (PKC)-mediated canonical NF-κB p65 subunit nuclear expression. Moreover, TSH stimulation phosphorylated the kinase TAK-1, and its knockdown abolished TSH-induced NF-κB transcriptional activity. TSH induced the transcriptional activity of the NF-κB subunit p65 in a protein kinase A (PKA)-dependent phosphorylation at Ser-276. In addition, p65 phosphorylation at Ser-276 induced acetyl transferase p300 recruitment, leading to its acetylation on Lys-310 and thereby enhancing its transcriptional activity. Evaluation of the role played by NF-κB in thyroid physiology demonstrated that the canonical NF-κB inhibitor BAY 11-7082 reduced TSH-induced expression of thyroid differentiation markers. The involvement of NF-κB signaling in thyroid physiology was confirmed by assessing the TSH-induced gene expression in primary cultures of NEMO-deficient mice thyrocytes. ChIP and the knockdown experiments revealed that p65 is a nuclear effector of TSH actions, inducing the transcripcional expression of thyroid differentiation markers. Conclusions: Taken together, our results point to NF-κB being a pivotal mediator in the TSH-induced thyroid follicular cell differentiation, a relevant finding with potential physiological and pathophysiological implications.

A Carboxy-Terminal Monoleucine-Based Motif Participates in the Basolateral Targeting of the Na+/I- Symporter.

The Na+/iodide (I-) symporter (NIS), a glycoprotein expressed at the basolateral plasma membrane of thyroid follicular cells, mediates I- accumulation for thyroid hormonogenesis and radioiodide therapy for differentiated thyroid carcinoma. However, differentiated thyroid tumors often exhibit lower I- transport than normal thyroid tissue (or even undetectable I- transport). Paradoxically, the majority of differentiated thyroid cancers show intracellular NIS expression, suggesting abnormal targeting to the plasma membrane. Therefore, a thorough understanding of the mechanisms that regulate NIS plasma membrane transport would have multiple implications for radioiodide therapy. In this study, we show that the intracellularly facing carboxy-terminus of NIS is required for the transport of the protein to the plasma membrane. Moreover, the carboxy-terminus contains dominant basolateral information. Using internal deletions and site-directed mutagenesis at the carboxy-terminus, we identified a highly conserved monoleucine-based sorting motif that determines NIS basolateral expression. Furthermore, in clathrin adaptor protein (AP)-1B-deficient cells, NIS sorting to the basolateral plasma membrane is compromised, causing the protein to also be expressed at the apical plasma membrane. Computer simulations suggest that the AP-1B subunit σ1 recognizes the monoleucine-based sorting motif in NIS carboxy-terminus. Although the mechanisms by which NIS is intracellularly retained in thyroid cancer remain elusive, our findings may open up avenues for identifying molecular targets that can be used to treat radioiodide-refractory thyroid tumors that express NIS intracellularly.

Structural Insights into the HWE Histidine Kinase Family: The Brucella Blue Light-Activated Histidine Kinase Domain.

In response to light, as part of a two-component system, the Brucella blue light-activated histidine kinase (LOV-HK) increases its autophosphorylation, modulating the virulence of this microorganism. The Brucella histidine kinase (HK) domain belongs to the HWE family, for which there is no structural information. The HWE family is exclusively present in proteobacteria and usually coupled to a wide diversity of light sensor domains. This work reports the crystal structure of the Brucella HK domain, which presents two different dimeric assemblies in the asymmetric unit: one similar to the already described canonical parallel homodimers (C) and the other, an antiparallel non-canonical (NC) dimer, each with distinct relative subdomain orientations and dimerization interfaces. Contrary to these crystallographic structures and unlike other HKs, in solution, the Brucella HK domain is monomeric and still active, showing an astonishing instability of the dimeric interface. Despite this instability, using cross-linking experiments, we show that the C dimer is the functionally relevant species. Mutational analysis demonstrates that the autophosphorylation activity occurs in cis. The different relative subdomain orientations observed for the NC and C states highlight the large conformational flexibility of the HK domain. Through the analysis of these alternative conformations by means of molecular dynamics simulations, we also propose a catalytic mechanism for Brucella LOV-HK.

Efficient Calculation of Enzyme Reaction Free Energy Profiles Using a Hybrid Differential Relaxation Algorithm: Application to Mycobacterial Zinc Hydrolases.

Determination of the free energy profile for an enzyme reaction mechanism is of primordial relevance, paving the way for our understanding of the enzyme's catalytic power at the molecular level. Although hybrid, mostly DFT-based, QM/MM methods have been extensively applied to this type of studies, achieving accurate and statistically converged results at a moderate computational cost is still an open challenge. Recently, we have shown that accurate results can be achieved in less computational time, combining Jarzynski's relationship with a hybrid differential relaxation algorithm (HyDRA), which allows partial relaxation of the solvent during the nonequilibrium steering of the reaction. In this work, we have applied this strategy to study two mycobacterial zinc hydrolases. Mycobacterium tuberculosis infections are still a worldwide problem and thus characterization and validation of new drug targets is an intense field of research. Among possible drug targets, recently two essential zinc hydrolases, MshB (Rv1170) and MA-amidase (Rv3717), have been proposed and structurally characterized. Although possible mechanisms have been proposed by analogy to the widely studied human Zn hydrolases, several key issues, particularly those related to Zn coordination sphere and its role in catalysis, remained unanswered. Our results show that mycobacterial Zn hydrolases share a basic two-step mechanism. First, the attacking water becomes deprotonated by the conserved base and establishes the new C-O bond leading to a tetrahedral intermediate. The intermediate requires moderate reorganization to allow for proton transfer to the amide N and C-N bond breaking to occur in the second step. Zn ion plays a key role in stabilizing the tetrahedral intermediate and balancing the negative charge of the substrate during hydroxide ion attack. Finally, comparative analysis of other Zn hydrolases points to a convergent mechanistic evolution.

The NtrY/X two-component system of Brucella spp. acts as a redox sensor and regulates the expression of nitrogen respiration enzymes.

Brucella spp. are facultative intracellular bacteria pathogenic for many mammalian species including humans, causing a disease called brucellosis. Learning how Brucella adapts to its intracellular niche is crucial for understanding its pathogenesis mechanism, allowing for the development of new and more effective vaccines and treatments against brucellosis. Brucella pathogenesis resides mostly in its ability to adapt to the harsh environmental conditions encountered during host infection such as the oxygen depletion. The mechanism by which Brucella senses the oxygen tension and triggers its environmental adaptation is unknown. In this work we show that the Brucella abortus NtrY/NtrX two-component system is involved in oxygen sensing through a haem group contained in a Per-ARNT-SIM (PAS) domain of the NtrY histidine kinase. The NtrY haem iron can be reduced to the ferrous form and is rapidly oxidized to the ferric form in presence of oxygen. Importantly, we show that the oxidation state of the haem iron modulates the autokinase activity, being the anoxygenic reduced ferrous form the signalling state of NtrY. Also, we show that ntrY gene expression increases under low oxygen tension and that NtrY transfers its signal to its cognate response regulator NtrX, regulating in this way the expression of nitrogen respiration enzymes. Based on these findings, we postulate that NtrY acts as a redox sensor in Brucella spp.