Atomistic simulations reveal impacts of missense mutations on the structure and function of SynGAP1.

De novo mutations in the synaptic GTPase activating protein (SynGAP) are associated with neurological disorders like intellectual disability, epilepsy, and autism. SynGAP is also implicated in Alzheimer's disease and cancer. Although pathogenic variants are highly penetrant in neurodevelopmental conditions, a substantial number of them are caused by missense mutations that are difficult to diagnose. Hence, in silico mutagenesis was performed for probing the missense effects within the N-terminal region of SynGAP structure. Through extensive molecular dynamics simulations, encompassing three 150-ns replicates for 211 variants, the impact of missense mutations on the protein fold was assessed. The effect of the mutations on the folding stability was also quantitatively assessed using free energy calculations. The mutations were categorized as potentially pathogenic or benign based on their structural impacts. Finally, the study introduces wild-type-SynGAP in complex with RasGTPase at the inner membrane, while considering the potential effects of mutations on these key interactions. This study provides structural perspective to the clinical assessment of SynGAP missense variants and lays the foundation for future structure-based drug discovery.

An outlook on structural biology after AlphaFold: tools, limits and perspectives.

AlphaFold and similar groundbreaking, AI-based tools, have revolutionized the field of structural bioinformatics, with their remarkable accuracy in ab-initio protein structure prediction. This success has catalyzed the development of new software and pipelines aimed at incorporating AlphaFold's predictions, often focusing on addressing the algorithm's remaining challenges. Here, we present the current landscape of structural bioinformatics shaped by AlphaFold, and discuss how the field is dynamically responding to this revolution, with new software, methods, and pipelines. While the excitement around AI-based tools led to their widespread application, it is essential to acknowledge that their practical success hinges on their integration into established protocols within structural bioinformatics, often neglected in the context of AI-driven advancements. Indeed, user-driven intervention is still as pivotal in the structure prediction process as in complementing state-of-the-art algorithms with functional and biological knowledge.

The fitness cost of spurious phosphorylation.

The fidelity of signal transduction requires the binding of regulatory molecules to their cognate targets. However, the crowded cell interior risks off-target interactions between proteins that are functionally unrelated. How such off-target interactions impact fitness is not generally known. Here, we use Saccharomyces cerevisiae to inducibly express tyrosine kinases. Because yeast lacks bona fide tyrosine kinases, the resulting tyrosine phosphorylation is biologically spurious. We engineered 44 yeast strains each expressing a tyrosine kinase, and quantitatively analysed their phosphoproteomes. This analysis resulted in ~30,000 phosphosites mapping to ~3500 proteins. The number of spurious pY sites generated correlates strongly with decreased growth, and we predict over 1000 pY events to be deleterious. However, we also find that many of the spurious pY sites have a negligible effect on fitness, possibly because of their low stoichiometry. This result is consistent with our evolutionary analyses demonstrating a lack of phosphotyrosine counter-selection in species with tyrosine kinases. Our results suggest that, alongside the risk for toxicity, the cell can tolerate a large degree of non-functional crosstalk as interaction networks evolve.

A Computational Approach in the Systematic Search of the Interaction Partners of Alternatively Spliced TREM2 Isoforms.

Alzheimer's disease is the most common form of dementia, characterized by the pathological accumulation of amyloid-beta (Aß) plaques and tau neurofibrillary tangles. Triggering receptor expressed on myeloid cells 2 (TREM2) is increasingly recognized as playing a central role in Aß clearance and microglia activation in AD. The TREM2 gene transcriptional product is alternatively spliced to produce three different protein isoforms. The canonical TREM2 isoform binds to DAP12 to activate downstream pathways. However, little is known about the function or interaction partners of the alternative TREM2 isoforms. The present study utilized a computational approach in a systematic search for new interaction partners of the TREM2 isoforms by integrating several state-of-the-art structural bioinformatics tools from initial large-scale screening to one-on-one corroborative modeling and eventual all-atom visualization. CD9, a cell surface glycoprotein involved in cell-cell adhesion and migration, was identified as a new interaction partner for two TREM2 isoforms, and CALM, a calcium-binding protein involved in calcium signaling, was identified as an interaction partner for a third TREM2 isoform, highlighting the potential role of cell adhesion and calcium regulation in AD.

ModFOLD9: A Web Server for Independent Estimates of 3D Protein Model Quality.

Accurate models of protein tertiary structures are now available from numerous advanced prediction methods, although the accuracy of each method often varies depending on the specific protein target. Additionally, many models may still contain significant local errors. Therefore, reliable, independent model quality estimates are essential both for identifying errors and selecting the very best models for further biological investigations. ModFOLD9 is a leading independent server for detecting the local errors in models produced by any method, and it can accurately discriminate between high-quality models from multiple alternative approaches. ModFOLD9 incorporates several new scores from deep learning-based approaches, leading to greatly improved prediction accuracy compared with earlier versions of the server. ModFOLD9 is continuously independently benchmarked, and it is shown to be highly competitive with other public servers. ModFOLD9 is freely available at https://www.reading.ac.uk/bioinf/ModFOLD/.

Clinical and genetic spectrum of Ataxia Telangiectasia Tunisian patients: Bioinformatic analysis unveil mechanisms of ATM variants pathogenicity.

Ataxia Telangiectasia (AT) is a rare multisystemic neurodegenerative disease caused by biallelic mutations in the ATM gene. Few clinical studies on AT disease have been conducted in Tunisia, however, the mutational landscape is still undefined. Our aim is to determine the clinical and genetic spectrum of AT Tunisian patients and to explore the potential underlying mechanism of variant pathogenicity. Sanger sequencing was performed for nine AT patients. A comprehensive computational analysis was conducted to evaluate the possible pathogenic effect of ATM identified variants. Genetic screening of ATM gene has identified nine different variants from which six have not been previously reported. In silico analysis has predicted a pathogenic effect of identified mutations. This was corroborated by a structural bioinformatics study based on molecular modeling and docking for novel missense mutations. Our findings suggest a profound impact of identified mutations not only on the ATM protein stability, but also on the ATM-ligand interactions. Our study characterizes the mutational landscape of AT Tunisian patients which will allow to set up genetic counseling and prenatal diagnosis for families at risk and expand the spectrum of ATM variants worldwide. Furthermore, understanding the mechanism that underpin variant pathogenicity could provide further insights into disease pathogenesis.

Respiratory tract infections: an update on the complexity of bacterial diversity, therapeutic interventions and breakthroughs.

Respiratory tract infections (RTIs) have a significant impact on global health, especially among children and the elderly. The key bacterial pathogens Streptococcus pneumoniae, Haemophilus influenzae, Klebsiella pneumoniae, Staphylococcus aureus and non-fermenting Gram Negative bacteria such as Acinetobacter baumannii and Pseudomonas aeruginosa are most commonly associated with RTIs. These bacterial pathogens have evolved a diverse array of resistance mechanisms through horizontal gene transfer, often mediated by mobile genetic elements and environmental acquisition. Treatment failures are primarily due to antimicrobial resistance and inadequate bacterial engagement, which necessitates the development of alternative treatment strategies. To overcome this, our review mainly focuses on different virulence mechanisms and their resulting pathogenicity, highlighting different therapeutic interventions to combat resistance. To prevent the antimicrobial resistance crisis, we also focused on leveraging the application of artificial intelligence and machine learning to manage RTIs. Integrative approaches combining mechanistic insights are crucial for addressing the global challenge of antimicrobial resistance in respiratory infections.

Studying protein-protein interactions: Latest and most popular approaches.

PPIs, or protein-protein interactions, are essential for many biological processes. According to the findings, abnormal PPIs have been linked to several diseases, such as cancer and infectious and neurological disorders. Consequently, focusing on PPIs is a path toward disease treatment and a crucial tool for producing novel medications. Many methods exist to investigate PPIs, including low- and high-throughput studies. Since many PPIs have been discovered using in vitro and in vivo experimental approaches, the use of computational methods to predict PPIs has grown due to the expanding scale of PPI data and the intrinsic complexity of interacting mechanisms. Recognizing PPI networks offers a systematic means of predicting protein functions, and pathways that are included. These investigations can help uncover the underlying molecular mechanisms of complex phenotypes and clarify the biological processes related to health and diseases. Therefore, our goal in this study is to provide an overview of the latest and most popular approaches for investigating PPIs. We also overview some important clinical approaches based on the PPIs and how these interactions can be targeted.

Geometric descriptors for beta turns.

Beta turns, in which the protein backbone abruptly changes direction over four amino acid residues, are the most common type of protein secondary structure after alpha helices and beta sheets and play key structural and functional roles. Previous work has produced classification systems for turn geometry at multiple levels of precision, but these operate in backbone dihedral-angle (Ramachandran) space, and the absence of a local Euclidean-space coordinate system and structural alignment for turns, or of any systematic Euclidean-space characterization of turn backbone shape, presents challenges for the visualization, comparison and analysis of the wide range of turn conformations and the design of turns and the structures that incorporate them. This work derives a turn-local coordinate system that implicitly aligns turns, together with a set of geometric descriptors that characterize the bulk BB shapes of turns and describe modes of structural variation not explicitly captured by existing systems. These modes are shown to be meaningful by the demonstration of clear relationships between descriptor values and the electrostatic energy of the beta-turn H-bond, the overrepresentations of key side-chain motifs, and the structural contexts of turns. Geometric turn descriptors complement Ramachandran-space classifications, and they can be used to select turn structures for compatibility with particular side-chain interactions or contexts. Potential applications include protein design and other tasks in which an enhanced Euclidean-space characterization of turns may improve understanding or performance. The web-based tools ExploreTurns, MapTurns, and ProfileTurn, available at www.betaturn.com, incorporate turn-local coordinates and turn descriptors and demonstrate their utility.

miRVim: Three-dimensional miRNA Structure Database.

INTRODUCTION: MicroRNAs (miRNAs), a distinct category of non-coding RNAs, exert multifaceted regulatory functions in a variety of organisms, including humans, animals, and plants. The inventory of identified miRNAs stands at approximately 60,000 among all species, and 1,926 in Homo sapiens manifest miRNA expression. METHOD: Their theranostic role has been explored by researchers over the last few decades, positioning them as prominent therapeutic targets as our understanding of RNA targeting advances. However, the limited availability of experimentally determined miRNA structures has constrained drug discovery efforts relying on virtual screening or computational methods, including machine learning and artificial intelligence. RESULTS: To address this lacuna, miRVim has been developed, providing a repository of human miRNA structures derived from both two-dimensional (MXFold2, CentroidFold, and RNAFold) and three-dimensional (RNAComposer and 3dRNA) structure prediction algorithms, in addition to experimentally available structures from the RCSB PDB repository. miRVim contains 13,971 predicted secondary structures and 17,045 predicted three-dimensional structures, filling the gap of unavailability of miRNA structure data bank. This database aims to facilitate computational data analysis for drug discovery, opening new avenues for advancing technologies, such as machine learning-based predictions in the field of RNA biology. CONCLUSION: The publicly accessible structures provided by miRVim, available at https://mirna.in/miRVim, offer a valuable resource for the research community, advancing the field of miRNA-related computational analysis and drug discovery.

Molecular Origins of the Mendelian Rare Diseases Reviewed by Orpha.net: A Structural Bioinformatics Investigation.

The study of rare diseases is important not only for the individuals affected but also for the advancement of medical knowledge and a deeper understanding of human biology and genetics. The wide repertoire of structural information now available from reliable and accurate prediction methods provides the opportunity to investigate the molecular origins of most of the rare diseases reviewed in the Orpha.net database. Thus, it has been possible to analyze the topology of the pathogenic missense variants found in the 2515 proteins involved in Mendelian rare diseases (MRDs), which form the database for our structural bioinformatics study. The amino acid substitutions responsible for MRDs showed different mutation site distributions at different three-dimensional protein depths. We then highlighted the depth-dependent effects of pathogenic variants for the 20,061 pathogenic variants that are present in our database. The results of this structural bioinformatics investigation are relevant, as they provide additional clues to mitigate the damage caused by MRD.

Drug to genome to drug: a computational large-scale chemogenomics screening for novel drug candidates against sporotrichosis.

Sporotrichosis is recognized as the predominant subcutaneous mycosis in South America, attributed to pathogenic species within the Sporothrix genus. Notably, in Brazil, Sporothrix brasiliensis emerges as the principal species, exhibiting significant sapronotic, zoonotic and enzootic epidemic potential. Consequently, the discovery of novel therapeutic agents for the treatment of sporotrichosis is imperative. The present study is dedicated to the repositioning of pharmaceuticals for sporotrichosis therapy. To achieve this goal, we designed a pipeline with the following steps: (a) compilation and preparation of Sporothrix genome data; (b) identification of orthologous proteins among the species; (c) identification of homologous proteins in publicly available drug-target databases; (d) selection of Sporothrix essential targets using validated genes from Saccharomyces cerevisiae; (e) molecular modeling studies; and (f) experimental validation of selected candidates. Based on this approach, we were able to prioritize eight drugs for in vitro experimental validation. Among the evaluated compounds, everolimus and bifonazole demonstrated minimum inhibitory concentration (MIC) values of 0.5 µg/mL and 4.0 µg/mL, respectively. Subsequently, molecular docking studies suggest that bifonazole and everolimus may target specific proteins within S. brasiliensis- namely, sterol 14-α-demethylase and serine/threonine-protein kinase TOR, respectively. These findings shed light on the potential binding affinities and binding modes of bifonazole and everolimus with their probable targets, providing a preliminary understanding of the antifungal mechanism of action of these compounds. In conclusion, our research advances the understanding of the therapeutic potential of bifonazole and everolimus, supporting their further investigation as antifungal agents for sporotrichosis in prospective hit-to-lead and preclinical investigations.

EGG: Accuracy Estimation of Individual Multimeric Protein Models Using Deep Energy-Based Models and Graph Neural Networks.

Reliable and accurate methods of estimating the accuracy of predicted protein models are vital to understanding their respective utility. Discerning how the quaternary structure conforms can significantly improve our collective understanding of cell biology, systems biology, disease formation, and disease treatment. Accurately determining the quality of multimeric protein models is still computationally challenging, as the space of possible conformations is significantly larger when proteins form in complex with one another. Here, we present EGG (energy and graph-based architectures) to assess the accuracy of predicted multimeric protein models. We implemented message-passing and transformer layers to infer the overall fold and interface accuracy scores of predicted multimeric protein models. When evaluated with CASP15 targets, our methods achieved promising results against single model predictors: fourth and third place for determining the highest-quality model when estimating overall fold accuracy and overall interface accuracy, respectively, and first place for determining the top three highest quality models when estimating both overall fold accuracy and overall interface accuracy.

Harnessing the 3D-Beacons Network: A Comprehensive Guide to Accessing and Displaying Protein Structure Data.

Recent advancements in protein structure determination and especially in protein structure prediction techniques have led to the availability of vast amounts of macromolecular structures. However, the accessibility and integration of these structures into scientific workflows are hindered by the lack of standardization among publicly available data resources. To address this issue, we introduced the 3D-Beacons Network, a unified platform that aims to establish a standardized framework for accessing and displaying protein structure data. In this article, we highlight the importance of standardized approaches for accessing protein structure data and showcase the capabilities of 3D-Beacons. We describe four protocols for finding and accessing macromolecular structures from various specialist data resources via 3D-Beacons. First, we describe three scenarios for programmatically accessing and retrieving data using the 3D-Beacons API. Next, we show how to perform sequence-based searches to find structures from model providers. Then, we demonstrate how to search for structures and fetch them directly into a workflow using JalView. Finally, we outline the process of facilitating access to data from providers interested in contributing their structures to the 3D-Beacons Network. © 2024 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Programmatic access to the 3D-Beacons API Basic Protocol 2: Sequence-based search using the 3D-Beacons API Basic Protocol 3: Accessing macromolecules from 3D-Beacons with JalView Basic Protocol 4: Enhancing data accessibility through 3D-Beacons.

Deep learning in structural bioinformatics: current applications and future perspectives.

In this review article, we explore the transformative impact of deep learning (DL) on structural bioinformatics, emphasizing its pivotal role in a scientific revolution driven by extensive data, accessible toolkits and robust computing resources. As big data continue to advance, DL is poised to become an integral component in healthcare and biology, revolutionizing analytical processes. Our comprehensive review provides detailed insights into DL, featuring specific demonstrations of its notable applications in bioinformatics. We address challenges tailored for DL, spotlight recent successes in structural bioinformatics and present a clear exposition of DL-from basic shallow neural networks to advanced models such as convolution, recurrent, artificial and transformer neural networks. This paper discusses the emerging use of DL for understanding biomolecular structures, anticipating ongoing developments and applications in the realm of structural bioinformatics.

Enhancing MD simulations: ASGARD's automated analysis for GROMACS.

Molecular Dynamics (MD) simulations are essential in analyzing the physical movement of molecules, with GROMACS being a widely recognized open-source package for this purpose. However, conducting analyses individually in GROMACS can take excessive time and effort. Addressing this challenge, we introduce ASGARD, an innovative workflow designed to streamline and automate the analysis of MD simulation of protein or protein-ligand complex. Unlike the traditional, manual approach, ASGARD enables researchers to generate comprehensive analyses with a single command line, significantly accelerating the research process and avoiding the laborious task of manual report generation. This tool automatically performs a range of analyses post-simulation, including system stability and flexibility assessments through RMSD Fluctuation and Distribution calculations. It further provides dynamic analysis using SASA, DSSP method graphs, and various interaction analyses. A key feature of ASGARD is its user-friendly design; it requires no additional installations or dependencies, making it highly accessible for researchers. In conclusion, ASGARD simplifies the MD simulation analysis process and substantially enhances efficiency and productivity in molecular research by providing an integrated, one-command analysis solution.Communicated by Ramaswamy H. Sarma.

Native dynamics and allosteric responses in PTP1B probed by high-resolution HDX-MS.

Protein tyrosine phosphatase 1B (PTP1B) is a validated therapeutic target for obesity, diabetes, and certain types of cancer. In particular, allosteric inhibitors hold potential for therapeutic use, but an incomplete understanding of conformational dynamics and allostery in this protein has hindered their development. Here, we interrogate solution dynamics and allosteric responses in PTP1B using high-resolution hydrogen-deuterium exchange mass spectrometry (HDX-MS), an emerging and powerful biophysical technique. Using HDX-MS, we obtain a detailed map of backbone amide exchange that serves as a proxy for the solution dynamics of apo PTP1B, revealing several flexible loops interspersed among more constrained and rigid regions within the protein structure, as well as local regions that exchange faster than expected from their secondary structure and solvent accessibility. We demonstrate that our HDX rate data obtained in solution adds value to estimates of conformational heterogeneity derived from a pseudo-ensemble constructed from ~200 crystal structures of PTP1B. Furthermore, we report HDX-MS maps for PTP1B with active-site versus allosteric small-molecule inhibitors. These maps suggest distinct and widespread effects on protein dynamics relative to the apo form, including changes in locations distal (>35 Å) from the respective ligand binding sites. These results illuminate that allosteric inhibitors of PTP1B can induce unexpected changes in dynamics that extend beyond the previously understood allosteric network. Together, our data suggest a model of BB3 allostery in PTP1B that combines conformational restriction of active-site residues with compensatory liberation of distal residues that aid in entropic balancing. Overall, our work showcases the potential of HDX-MS for elucidating aspects of protein conformational dynamics and allosteric effects of small-molecule ligands and highlights the potential of integrating HDX-MS alongside other complementary methods, such as room-temperature X-ray crystallography, NMR spectroscopy, and molecular dynamics simulations, to guide the development of new therapeutics.

Databases of ligand-binding pockets and protein-ligand interactions.

Many research groups and institutions have created a variety of databases curating experimental and predicted data related to protein-ligand binding. The landscape of available databases is dynamic, with new databases emerging and established databases becoming defunct. Here, we review the current state of databases that contain binding pockets and protein-ligand binding interactions. We have compiled a list of such databases, fifty-three of which are currently available for use. We discuss variation in how binding pockets are defined and summarize pocket-finding methods. We organize the fifty-three databases into subgroups based on goals and contents, and describe standard use cases. We also illustrate that pockets within the same protein are characterized differently across different databases. Finally, we assess critical issues of sustainability, accessibility and redundancy.

Identification of type VI secretion system effector-immunity pairs using structural bioinformatics.

The type VI secretion system (T6SS) is an important mediator of microbe-microbe and microbe-host interactions. Gram-negative bacteria use the T6SS to inject T6SS effectors (T6Es), which are usually proteins with toxic activity, into neighboring cells. Antibacterial effectors have cognate immunity proteins that neutralize self-intoxication. Here, we applied novel structural bioinformatic tools to perform systematic discovery and functional annotation of T6Es and their cognate immunity proteins from a dataset of 17,920 T6SS-encoding bacterial genomes. Using structural clustering, we identified 517 putative T6E families, outperforming sequence-based clustering. We developed a logistic regression model to reliably quantify protein-protein interaction of new T6E-immunity pairs, yielding candidate immunity proteins for 231 out of the 517 T6E families. We used sensitive structure-based annotation which yielded functional annotations for 51% of the T6E families, again outperforming sequence-based annotation. Next, we validated four novel T6E-immunity pairs using basic experiments in E. coli. In particular, we showed that the Pfam domain DUF3289 is a homolog of Colicin M and that DUF943 acts as its cognate immunity protein. Furthermore, we discovered a novel T6E that is a structural homolog of SleB, a lytic transglycosylase, and identified a specific glutamate that acts as its putative catalytic residue. Overall, this study applies novel structural bioinformatic tools to T6E-immunity pair discovery, and provides an extensive database of annotated T6E-immunity pairs.

Uncovering the secrets of resistance: An introduction to computational methods in infectious disease research.

Antimicrobial resistance (AMR) is a growing global concern with significant implications for infectious disease control and therapeutics development. This chapter presents a comprehensive overview of computational methods in the study of AMR. We explore the prevalence and statistics of AMR, underscoring its alarming impact on public health. The role of AMR in infectious disease outbreaks and its impact on therapeutics development are discussed, emphasizing the need for novel strategies. Resistance mutations are pivotal in AMR, enabling pathogens to evade antimicrobial treatments. We delve into their importance and contribution to the spread of AMR. Experimental methods for quantitatively evaluating resistance mutations are described, along with their limitations. To address these challenges, computational methods provide promising solutions. We highlight the advantages of computational approaches, including rapid analysis of large datasets and prediction of resistance profiles. A comprehensive overview of computational methods for studying AMR is presented, encompassing genomics, proteomics, structural bioinformatics, network analysis, and machine learning algorithms. The strengths and limitations of each method are briefly outlined. Additionally, we introduce ResScan-design, our own computational method, which employs a protein (re)design protocol to identify potential resistance mutations and adaptation signatures in pathogens. Case studies are discussed to showcase the application of ResScan in elucidating hotspot residues, understanding underlying mechanisms, and guiding the design of effective therapies. In conclusion, we emphasize the value of computational methods in understanding and combating AMR. Integration of experimental and computational approaches can expedite the discovery of innovative antimicrobial treatments and mitigate the threat posed by AMR.