Local energetic frustration conservation in protein families and superfamilies.

Energetic local frustration offers a biophysical perspective to interpret the effects of sequence variability on protein families. Here we present a methodology to analyze local frustration patterns within protein families and superfamilies that allows us to uncover constraints related to stability and function, and identify differential frustration patterns in families with a common ancestry. We analyze these signals in very well studied protein families such as PDZ, SH3, ɑ and ß globins and RAS families. Recent advances in protein structure prediction make it possible to analyze a vast majority of the protein space. An automatic and unsupervised proteome-wide analysis on the SARS-CoV-2 virus demonstrates the potential of our approach to enhance our understanding of the natural phenotypic diversity of protein families beyond single protein instances. We apply our method to modify biophysical properties of natural proteins based on their family properties, as well as perform unsupervised analysis of large datasets to shed light on the physicochemical signatures of poorly characterized proteins such as the ones belonging to emergent pathogens.

Best Practices on QM/MM Simulations of Biological Systems.

During the second half of the 20th century, following structural biology hallmark works on DNA and proteins, biochemists shifted their questions from "what does this molecule look like?" to "how does this process work?". Prompted by the theoretical and practical developments in computational chemistry, this led to the emergence of biomolecular simulations and, along with the 2013 Nobel Prize in Chemistry, to the development of hybrid QM/MM methods. QM/MM methods are necessary whenever the problem we want to address involves chemical reactivity and/or a change in the system's electronic structure, with archetypal examples being the studies of an enzyme's reaction mechanism and a metalloprotein's active site. In the last decades QM/MM methods have seen an increasing adoption driven by their incorporation in widely used biomolecular simulation software. However, properly setting up a QM/MM simulation is not an easy task, and several issues need to be properly addressed to obtain meaningful results. In the present work, we describe both the theoretical concepts and practical issues that need to be considered when performing QM/MM simulations. We start with a brief historical perspective on the development of these methods and describe when and why QM/MM methods are mandatory. Then we show how to properly select and analyze the performance of the QM level of theory, the QM system size, and the position and type of the boundaries. We show the relevance of performing prior QM model system (or QM cluster) calculations in a vacuum and how to use the corresponding results to adequately calibrate those derived from QM/MM. We also discuss how to prepare the starting structure and how to select an adequate simulation strategy, including those based on geometry optimizations as well as free energy methods. In particular, we focus on the determination of free energy profiles using multiple steered molecular dynamics (MSMD) combined with Jarzynski's equation. Finally, we describe the results for two illustrative and complementary examples: the reaction performed by chorismate mutase and the study of ligand binding to hemoglobins. Overall, we provide many practical recommendations (or shortcuts) together with important conceptualizations that we hope will encourage more and more researchers to incorporate QM/MM studies into their research projects.

An allosteric switch ensures efficient unidirectional information transmission by the histidine kinase DesK from Bacillus subtilis.

Phosphorylation carries chemical information in biological systems. In two-component systems (TCSs), the sensor histidine kinase and the response regulator are connected through phosphoryl transfer reactions that may be uni- or bidirectional. Directionality enables the construction of complex regulatory networks that optimize signal propagation and ensure the forward flow of information. We combined x-ray crystallography, hybrid quantum mechanics/molecular mechanics (QM/MM) simulations, and systems-integrative kinetic modeling approaches to study phosphoryl flow through the Bacillus subtilis thermosensing TCS DesK-DesR. The allosteric regulation of the histidine kinase DesK was critical to avoid back transfer of phosphoryl groups and futile phosphorylation-dephosphorylation cycles by isolating phosphatase, autokinase, and phosphotransferase activities. Interactions between the kinase's ATP-binding domain and the regulator's receiver domain placed the regulator in two distinct positions in the phosphotransferase and phosphatase complexes, thereby determining whether a key glutamine residue in DesK was properly situated to assist in the dephosphorylation reaction. Moreover, an energetically unfavorable phosphotransferase conformation when DesK was not phosphorylated minimized reverse phosphoryl transfer. DesR dimerization and a dissociative phosphoryl transfer reaction also enforced the direction of phosphoryl flow. Shorter or longer distances between the phosphoryl acceptor and donor residues shifted the phosphoryl transfer equilibrium by modulating the stabilizing effect of the Mg2+ cofactor. These mechanisms control the directionality of signal transmission and show how structure-encoded allostery stores and transmits information in signaling systems.

The interplay between serine proteases and caspase-1 regulates the autophagy-mediated secretion of Interleukin-1 beta in human neutrophils.

Neutrophils play major roles against bacteria and fungi infections not only due to their microbicide properties but also because they release mediators like Interleukin-1 beta (IL-1ß) that contribute to orchestrate the inflammatory response. This cytokine is a leaderless protein synthesized in the cytoplasm as a precursor (pro-IL-1ß) that is proteolytically processed to its active isoform and released from human neutrophils by secretory autophagy. In most myeloid cells, pro-IL-1ß is processed by caspase-1 upon inflammasome activation. Here we employed neutrophils from both healthy donors and patients with a gain-of-function (GOF) NLRP3-mutation to dissect IL-1ß processing in these cells. We found that although caspase-1 is required for IL-1ß secretion, it undergoes rapid inactivation, and instead, neutrophil serine proteases play a key role in pro-IL-1ß processing. Our findings bring to light distinctive features of the regulation of caspase-1 activity in human neutrophils and reveal new molecular mechanisms that control human neutrophil IL-1ß secretion.

Solvent Sites Improve Docking Performance of Protein-Protein Complexes and Protein-Protein Interface-Targeted Drugs.

Protein-protein interactions (PPIs) are essential, and modulating their function through PPI-targeted drugs is an important research field. PPI sites are shallow protein surfaces readily accessible to the solvent, thus lacking a proper pocket to fit a drug, while their lack of endogenous ligands prevents drug design by chemical similarity. The development of PPI-blocking compounds is, therefore, a tough challenge. Mixed solvent molecular dynamics has been shown to reveal protein-ligand interaction hot spots in protein active sites by identifying solvent sites (SSs). Furthermore, our group has shown that SSs significantly improve protein-ligand docking. In the present work, we extend our analysis to PPI sites. In particular, we analyzed water, ethanol, and phenol-derived sites in terms of their capacity to predict protein-drug and protein-protein interactions. Subsequently, we show how this information can be incorporated to improve both protein-ligand and protein-protein docking. Finally, we highlight the presence of aromatic clusters as key elements of the corresponding interactions.

Delayed Puberty Due to a WDR11 Truncation at Its N-Terminal Domain Leading to a Mild Form of Ciliopathy Presenting With Dissociated Central Hypogonadism: Case Report.

Pubertal delay in males is frequently due to constitutional delay of growth and puberty, but pathologic hypogonadism should be considered. After general illnesses and primary testicular failure are ruled out, the main differential diagnosis is central (or hypogonadotropic) hypogonadism, resulting from a defective function of the gonadotropin-releasing hormone (GnRH)/gonadotropin axis. Ciliopathies arising from defects in non-motile cilia are responsible for developmental disorders affecting the sense organs and the reproductive system. WDR11-mediated signaling in non-motile cilia is critical for fetal development of GnRH neurons. Only missense variants of WDR11 have been reported to date in patients with central hypogonadism, suggesting that nonsense variants could lead to more complex phenotypes. We report the case of a male patient presenting with delayed puberty due to Kallmann syndrome (central hypogonadism associated with hyposmia) in whom the next-generation sequencing analysis identified a novel heterozygous base duplication, leading to a frameshift and a stop codon in the N-terminal region of WDR11. The variant was predicted to undergo nonsense-mediated decay and classified as probably pathogenic following the American College of Medical Genetics and Genomics (ACMG) criteria. This is the first report of a variant in the WDR11 N-terminal region predicted to lead to complete expression loss that, contrary to expectations, led to a mild form of ciliopathy resulting in isolated Kallmann syndrome.

Cosolvent Sites-Based Discovery of Mycobacterium Tuberculosis Protein Kinase G Inhibitors.

Computer-aided drug discovery methods play a major role in the development of therapeutically important small molecules, but their performance needs to be improved. Molecular dynamics simulations in mixed solvents are useful in understanding protein-ligand recognition and improving molecular docking predictions. In this work, we used ethanol as a cosolvent to find relevant interactions for ligands toward protein kinase G, an essential protein of Mycobacterium tuberculosis (Mtb). We validated the hot spots by screening a database of fragment-like compounds and another one of known kinase inhibitors. Next, we performed a pharmacophore-guided docking simulation and found three low micromolar inhibitors, including one with a novel chemical scaffold that we expanded to four derivative compounds. Binding affinities were characterized by intrinsic fluorescence quenching assays, isothermal titration calorimetry, and the analysis of melting curves. The predicted binding mode was confirmed by X-ray crystallography. Finally, the compounds significantly inhibited the viability of Mtb in infected THP-1 macrophages.

Type 1 Insulin-Like Growth Factor Receptor Nuclear Localization in High-Grade Glioma Cells Enhances Motility, Metabolism, and In Vivo Tumorigenesis.

Gliomas are the most frequent solid tumors in children. Among these, high-grade gliomas are less common in children than in adults, though they are similar in their aggressive clinical behavior. In adults, glioblastoma is the most lethal tumor of the central nervous system. Insulin-like growth factor 1 receptor (IGF1R) plays an important role in cancer biology, and its nuclear localization has been described as an adverse prognostic factor in different tumors. Previously, we have demonstrated that, in pediatric gliomas, IGF1R nuclear localization is significantly associated with high-grade tumors, worst clinical outcome, and increased risk of death. Herein we explore the role of IGF1R intracellular localization by comparing two glioblastoma cell lines that differ only in their IGF1R capacity to translocate to the nucleus. In vitro, IGF1R nuclear localization enhances glioblastoma cell motility and metabolism without affecting their proliferation. In vivo, IGF1R has the capacity to translocate to the nucleus and allows not only a higher proliferation rate and the earlier development of tumors but also renders the cells sensitive to OSI906 therapy. With this work, we provide evidence supporting the implications of the presence of IGF1R in the nucleus of glioma cells and a potential therapeutic opportunity for patients harboring gliomas with IGF1R nuclear localization.

Mycobacterium tuberculosis DosS binds H2S through its Fe3+ heme iron to regulate the DosR dormancy regulon.

Mycobacterium tuberculosis (Mtb) senses and responds to host-derived gasotransmitters NO and CO via heme-containing sensor kinases DosS and DosT and the response regulator DosR. Hydrogen sulfide (H2S) is an important signaling molecule in mammals, but its role in Mtb physiology is unclear. We have previously shown that exogenous H2S can modulate expression of genes in the Dos dormancy regulon via an unknown mechanism(s). Here, we test the hypothesis that Mtb senses and responds to H2S via the DosS/T/R system. Using UV-Vis and EPR spectroscopy, we show that H2S binds directly to the ferric (Fe3+) heme of DosS (KDapp = 5.30 µM) but not the ferrous (Fe2+) form. No interaction with DosT(Fe2+-O2) was detected. We found that the binding of sulfide can slowly reduce the DosS heme iron to the ferrous form. Steered Molecular Dynamics simulations show that H2S, and not the charged HS- species, can enter the DosS heme pocket. We also show that H2S increases DosS autokinase activity and subsequent phosphorylation of DosR, and H2S-mediated increases in Dos regulon gene expression is lost in Mtb lacking DosS. Finally, we demonstrate that physiological levels of H2S in macrophages can induce DosR regulon genes via DosS. Overall, these data reveal a novel mechanism whereby Mtb senses and responds to a third host gasotransmitter, H2S, via DosS(Fe3+). These findings highlight the remarkable plasticity of DosS and establish a new paradigm for how bacteria can sense multiple gasotransmitters through a single heme sensor kinase.

Specificity and Reactivity of Mycobacterium tuberculosis Serine/Threonine Kinases PknG and PknB.

Mycobacterium tuberculosis (Mtb), the causative agent of Tuberculosis, has 11 eukaryotic-like serine/threonine protein kinases, which play essential roles in cell growth, signal transduction, and pathogenesis. Protein kinase G (PknG) regulates the carbon and nitrogen metabolism by phosphorylation of the glycogen accumulation regulator (GarA) protein at Thr21. Protein kinase B (PknB) is involved in cell wall synthesis and cell shape, as well as phosphorylates GarA but at Thr22. While PknG seems to be constitutively activated and recognition of GarA requires phosphorylation in its unstructured tail, PknB activation is triggered by phosphorylation of its activation loop, which allows binding of the forkhead-associated domain of GarA. In the present work, we used molecular dynamics and quantum-mechanics/molecular mechanics simulations of the catalytically competent complex and kinase activity assays to understand PknG/PknB specificity and reactivity toward GarA. Two hydrophobic residues in GarA, Val24 and Phe25, seem essential for PknG binding and allow specificity for Thr21 phosphorylation. On the other hand, phosphorylated residues in PknB bind Arg26 in GarA and regulate its specificity for Thr22. We also provide a detailed analysis of the free energy profile for the phospho-transfer reaction and show why PknG has a constitutively active conformation not requiring priming phosphorylation in contrast to PknB. Our results provide new insights into these two key enzymes relevant for Mtb and the mechanisms of serine/threonine phosphorylation in bacteria.

MotSASi: Functional short linear motifs (SLiMs) prediction based on genomic single nucleotide variants and structural data.

Short linear motifs (SLiMs) are key to cell physiology mediating reversible protein-protein interactions. Precise identification of SLiMs remains a challenge, being the main drawback of most bioinformatic prediction tools, their low specificity (high number of false positives). An important, usually overlooked, aspect is the relation between SLiMs mutations and disease. The presence of variants in each residue position can be used to assess the relevance of the corresponding residue(s) for protein function, and its (in)tolerance to change. In the present work, we combined sequence variant information and structural analysis of the energetic impact of single amino acid substitution (SAS) in SLiM-Receptor complex structure, and showed that it improves prediction of true functional SLiMs. Our strategy is based on building a SAS tolerance matrix that shows, for each position, whether one of the possible 19 SAS is tolerated or not. Herein we present the MotSASi strategy and analyze in detail 3 SLiMs involved in intracellular protein trafficking (phospho-independent tyrosine-based motif (NPx[Y/F]), type 1 PDZ-binding motif ([S/T]x[V/I/L]COOH) and tryptophan-acidic motif ([L/M]xW[D/E])). Our results show that inclusion of variant and structure information improves both prediction of true SLiMs and rejection of false positives, while also allowing better classification of variants inside SLiMs, a result with a direct impact in clinical genomics.

A Remote Secondary Binding Pocket Promotes Heteromultivalent Targeting of DC-SIGN.

Dendritic cells (DC) are antigen-presenting cells coordinating the interplay of the innate and the adaptive immune response. The endocytic C-type lectin receptors DC-SIGN and Langerin display expression profiles restricted to distinct DC subtypes and have emerged as prime targets for next-generation immunotherapies and anti-infectives. Using heteromultivalent liposomes copresenting mannosides bearing aromatic aglycones with natural glycan ligands, we serendipitously discovered striking cooperativity effects for DC-SIGN+ but not for Langerin+ cell lines. Mechanistic investigations combining NMR spectroscopy with molecular docking and molecular dynamics simulations led to the identification of a secondary binding pocket for the glycomimetics. This pocket, located remotely of DC-SIGN's carbohydrate bindings site, can be leveraged by heteromultivalent avidity enhancement. We further present preliminary evidence that the aglycone allosterically activates glycan recognition and thereby contributes to DC-SIGN-specific cell targeting. Our findings have important implications for both translational and basic glycoscience, showcasing heteromultivalent targeting of DCs to improve specificity and supporting potential allosteric regulation of DC-SIGN and CLRs in general.

High-throughput splicing assays identify missense and silent splice-disruptive POU1F1 variants underlying pituitary hormone deficiency.

Pituitary hormone deficiency occurs in ∼1:4,000 live births. Approximately 3% of the cases are due to mutations in the alpha isoform of POU1F1, a pituitary-specific transcriptional activator. We found four separate heterozygous missense variants in unrelated individuals with hypopituitarism that were predicted to affect a minor isoform, POU1F1 beta, which can act as a transcriptional repressor. These variants retain repressor activity, but they shift splicing to favor the expression of the beta isoform, resulting in dominant-negative loss of function. Using a high-throughput splicing reporter assay, we tested 1,070 single-nucleotide variants in POU1F1. We identified 96 splice-disruptive variants, including 14 synonymous variants. In separate cohorts, we found two additional synonymous variants nominated by this screen that co-segregate with hypopituitarism. This study underlines the importance of evaluating the impact of variants on splicing and provides a catalog for interpretation of variants of unknown significance in POU1F1.

From Genome to Drugs: New Approaches in Antimicrobial Discovery.

Decades of successful use of antibiotics is currently challenged by the emergence of increasingly resistant bacterial strains. Novel drugs are urgently required but, in a scenario where private investment in the development of new antimicrobials is declining, efforts to combat drug-resistant infections become a worldwide public health problem. Reasons behind unsuccessful new antimicrobial development projects range from inadequate selection of the molecular targets to a lack of innovation. In this context, increasingly available omics data for multiple pathogens has created new drug discovery and development opportunities to fight infectious diseases. Identification of an appropriate molecular target is currently accepted as a critical step of the drug discovery process. Here, we review how diverse layers of multi-omics data in conjunction with structural/functional analysis and systems biology can be used to prioritize the best candidate proteins. Once the target is selected, virtual screening can be used as a robust methodology to explore molecular scaffolds that could act as inhibitors, guiding the development of new drug lead compounds. This review focuses on how the advent of omics and the development and application of bioinformatics strategies conduct a "big-data era" that improves target selection and lead compound identification in a cost-effective and shortened timeline.

Structural modeling of a novel membrane-bound globin-coupled sensor in Geobacter sulfurreducens.

Globin-coupled sensors (GCS) usually consist of three domains: a sensor/globin, a linker, and a transmitter domain. The globin domain (GD), activated by ligand binding and/or redox change, induces an intramolecular signal transduction resulting in a response of the transmitter domain. Depending on the nature of the transmitter domain, GCSs can have different activities and functions, including adenylate and di-guanylate cyclase, histidine kinase activity, aerotaxis and/or oxygen sensing function. The gram-negative delta-proteobacterium Geobacter sulfurreducens expresses a protein with a GD covalently linked to a four transmembrane domain, classified, by sequence similarity, as GCS (GsGCS). While its GD is fully characterized, not so its transmembrane domain, which is rarely found in the globin superfamily. In the present work, GsGCS was characterized spectroscopically and by native ion mobility-mass spectrometry in combination with cryo-electron microscopy. Although lacking high resolution, the oligomeric state and the electron density map were valuable for further rational modeling of the full-length GsGCS structure. This model demonstrates that GsGCS forms a transmembrane domain-driven tetramer with minimal contact between the GDs and with the heme groups oriented outward. This organization makes an intramolecular signal transduction less likely. Our results, including the auto-oxidation rate and redox potential, suggest a potential role for GsGCS as redox sensor or in a membrane-bound e-/H+ transfer. As such, GsGCS might act as a player in connecting energy production to the oxidation of organic compounds and metal reduction. Database searches indicate that GDs linked to a four or seven helices transmembrane domain occur more frequently than expected.

Biased Docking for Protein-Ligand Pose Prediction.

The interaction between a protein and its ligands is one of the basic and most important processes in biological chemistry. Docking methods aim to predict the molecular 3D structure of protein-ligand complexes starting from coordinates of the protein and the ligand separately. They are widely used in both industry and academia, especially in the context of drug development projects. AutoDock4 is one of the most popular docking tools and, as for any docking method, its performance is highly system dependent. Knowledge about specific protein-ligand interactions on a particular target can be used to successfully overcome this limitation. Here, we describe how to apply the AutoDock Bias protocol, a simple and elegant strategy that allows users to incorporate target-specific information through a modified scoring function that biases the ligand structure towards those poses (or conformations) that establish selected interactions. We discuss two examples using different bias sources. In the first, we show how to steer dockings towards interactions derived from crystal structures of the receptor with different ligands; in the second example, we define and apply hydrophobic biases derived from Molecular Dynamics simulations in mixed solvents. Finally, we discuss general concepts of biased docking, its performance in pose prediction, and virtual screening campaigns as well as other potential applications.

Oculocutaneous albinism type 1B associated with a functionally significant tyrosinase gene polymorphism detected with Whole Exome Sequencing.

Background: Oculocutaneous albinism (OCA) is a Mendelian disorder characterized by hypopigmentation of the skin, hair, and eyes, hypoplastic fovea, and low vision, known to be caused by mutations in the Tyrosinase (TYR) gene. Among the known TYR variants, some reduce but do not completely eliminate tyrosinase activity, allowing residual production of melanin and resulting in a contradictory assignment as either pathogenic or benign, preventing a precise clinical diagnostic.Materials and Methods: In the present work, we performed Whole Exome Sequencing and subsequent Sanger sequencing in a young male clinically diagnosed with OCA.Results: Whole-exome sequencing analysis revealed the identification of two variants in trans in TYR. The first, corresponds to a known pathogenic variant G47D, while the second S192Y, was considered a polymorphism due to its relatively high frequency in the European population.Conclusion: The lack of other pathogenic variants in TYR, the reported reduced enzymatic activity (ca. 40% respect to wt) for S192Y, together with the structural in-silico analysis strongly suggest that both reported variants are jointly disease-causing and that S192Y should be considered as likely pathogenic, especially when it is found in trans with a null variant.

Clamping, bending, and twisting inter-domain motions in the misfold-recognizing portion of UDP-glucose: Glycoprotein glucosyltransferase.

UDP-glucose:glycoprotein glucosyltransferase (UGGT) flags misfolded glycoproteins for ER retention. We report crystal structures of full-length Chaetomium thermophilum UGGT (CtUGGT), two CtUGGT double-cysteine mutants, and its TRXL2 domain truncation (CtUGGT-ΔTRXL2). CtUGGT molecular dynamics (MD) simulations capture extended conformations and reveal clamping, bending, and twisting inter-domain movements. We name "Parodi limit" the maximum distance on the same glycoprotein between a site of misfolding and an N-linked glycan that can be reglucosylated by monomeric UGGT in vitro, in response to recognition of misfold at that site. Based on the MD simulations, we estimate the Parodi limit as around 70-80 Å. Frequency distributions of distances between glycoprotein residues and their closest N-linked glycosylation sites in glycoprotein crystal structures suggests relevance of the Parodi limit to UGGT activity in vivo. Our data support a "one-size-fits-all adjustable spanner" UGGT substrate recognition model, with an essential role for the UGGT TRXL2 domain.

Inflammatory cutaneous lesions and pulmonary manifestations in a new patient with autosomal recessive ISG15 deficiency case report.

Interferon-stimulated gene 15 (ISG15) was the first ubiquitin-like modifier protein identified that acts by protein conjugation (ISGylation) and is thought to modulate IFN-induced inflammation. Here, we report a new patient from a non-consanguineous Argentinian family, who was followed for recurrent ulcerative skin lesions, cerebral calcifications and lung disease. Whole Exome Sequencing (WES) revealed two novel compound heterozygous variants (c.285del and c.299_312del, NM_005101.4 GRCh37(hg19), both classified as pathogenic according to ACMG criteria) in the ISG15 gene, resulting in a complete deficiency due to disruption of the second ubiquitin domain of the corresponding protein. The clinical phenotype of this patient is unique given the presence of recurrent pulmonary manifestations and the absence of mycobacterial infections, thus resulting in a phenotype distinct from that previously described in patients with biallelic loss-of-function (LOF) ISG15 variants. This case highlights the role of ISG15 as an immunomodulating factor whose LOF variants result in heterogeneous clinical presentations.

Systemic Type I IFN Inflammation in Human ISG15 Deficiency Leads to Necrotizing Skin Lesions.

Most monogenic disorders have a primary clinical presentation. Inherited ISG15 deficiency, however, has manifested with two distinct presentations to date: susceptibility to mycobacterial disease and intracranial calcifications from hypomorphic interferon-II (IFN-II) production and excessive IFN-I response, respectively. Accordingly, these patients were managed for their infectious and neurologic complications. Herein, we describe five new patients with six novel ISG15 mutations presenting with skin lesions who were managed for dermatologic disease. Cellularly, we denote striking specificity to the IFN-I response, which was previously assumed to be universal. In peripheral blood, myeloid cells display the most robust IFN-I signatures. In the affected skin, IFN-I signaling is observed in the keratinocytes of the epidermis, endothelia, and the monocytes and macrophages of the dermis. These findings define the specific cells causing circulating and dermatologic inflammation and expand the clinical spectrum of ISG15 deficiency to dermatologic presentations as a third phenotype co-dominant to the infectious and neurologic manifestations.