The NRGTEN Python package: an extensible toolkit for coarse-grained normal mode analysis of proteins, nucleic acids, small molecules and their complexes

Summary: Coarse-grained normal mode analysis (NMA) is a fast computational technique to study the dynamics of biomolecules. Here we present the Najmanovich Research Group Toolkit for Elastic Networks (NRGTEN). NRGTEN is a Python toolkit that implements four different NMA models in addition to popular and novel metrics to benchmark and measure properties from these models. Furthermore, the toolkit is available as a public Python package and is easily extensible for the development or implementation of additional NMA models. The inclusion of the ENCoM model (Elastic Network Contact Model) developed in our group within NRGTEN is noteworthy, owing to its account for the specific chemical nature of atomic interactions. This makes possible some unique predictions of the effect of mutations, such as on stability (via changes in vibrational entropy differences), on the transition probability between different conformational states or on the flexibility profile of the whole macromolecule/complex (to study allostery and signalling). In addition, all NMA models can be used to generate conformational ensembles from a starting structure to aid in protein-protein, protein-ligand or other docking studies among applications. NRGTEN is freely available via a public Python package which can be easily installed on any modern machine and includes a detailed user guide hosted online. Availability and implementation: https://github.com/gregorpatof/nrgten_package/ Contact: [email protected]

IsoCleft Finder – a web-based tool for the detection and analysis of protein binding-site geometric and chemical similarities

IsoCleft Finder is a web-based tool for the detection of local geometric and chemical similarities between potential small-molecule binding cavities and a non-redundant dataset of ligand-bound known small-molecule binding-sites. The non-redundant dataset developed as part of this study is composed of 7339 entries representing unique Pfam/PDB-ligand (hetero group code) combinations with known levels of cognate ligand similarity. The query cavity can be uploaded by the user or detected automatically by the system using existing PDB entries as well as user-provided structures in PDB format. In all cases, the user can refine the definition of the cavity interactively via a browser-based Jmol 3D molecular visualization interface. Furthermore, users can restrict the search to a subset of the dataset using a cognate-similarity threshold. Local structural similarities are detected using the IsoCleft software and ranked according to two criteria (number of atoms in common and Tanimoto score of local structural similarity) and the associated Z-score and p-value measures of statistical significance. The results, including predicted ligands, target proteins, similarity scores, number of atoms in common, etc., are shown in a powerful interactive graphical interface. This interface permits the visualization of target ligands superimposed on the query cavity and additionally provides a table of pairwise ligand topological similarities. Similarities between top scoring ligands serve as an additional tool to judge the quality of the results obtained. We present several examples where IsoCleft Finder provides useful functional information. IsoCleft Finder results are complementary to existing approaches for the prediction of protein function from structure, rational drug design and x-ray crystallography. IsoCleft Finder can be found at: http://bcb.med.usherbrooke.ca/isocleftfinder

SPEAR: Systematic ProtEin AnnotatoR

Summary We present SPEAR, a lightweight and rapid SARS-CoV-2 variant annotation and scoring tool, for identifying mutations contributing to potential immune escape and transmissibility (ACE2 binding) at point of sequencing. SPEAR can be used in the field to evaluate genomic surveillance results in real-time and features a powerful interactive data visualisation report. Availability and implementation SPEAR and documentation are freely available on GitHub: https://github.com/m-crown/SPEAR and is implemented in Python and installable via Conda environment. Supplemental Supplementary data are available at Bioinformatics online

Toll-like receptor stimulation differentially regulates vasoactive intestinal peptide type 2 receptor in macrophages

Vasoactive intestinal peptide (VIP) was originally isolated as a vasodilator intestinal peptide, then as a neuropeptide. In the immune system, VIP is described as an endogenous macrophage-deactivating factor. VIP exerts its immunological actions in a paracrine and/or autocrine manner, through specific receptors. However, very little is known about the molecular regulation of VIP type 2 receptor (VPAC2) in the immune system. We now report that different toll-like receptor (TLR) ligands selectively regulate the VPAC2 receptor gene and show a gene repression system controlled by key protein kinase signalling cascades in macrophages. VPAC2 gene expression is regulated by gram-positive (TLR2 ligands) and gram-negative bacteria wall constituents (TLR4 ligands). Moreover, VPAC2 is tightly regulated: TLR2- or TLR2/6- but not TLR2/1-mediated mechanisms are responsible for the induction of VPAC2. TLR stimulation by viral or bacterial nucleic acids did not modify the VPAC2 mRNA levels. Remarkably, imiquimod – a synthetic TLR7 ligand – led to a potent up-regulation of VPAC2 gene expression. TLR5 stimulation by flagellin present in gram-positive and gram-negative bacteria did not affect VPAC2 mRNA. The p38 mitogen-activated protein kinase (MAPK) activity accounted for the TLR4-mediated induction of VPAC2 gene expression. Surprisingly, our data strongly suggest for the first time a tightly repressed control of VPAC2 mRNA induction by elements downstream of MAPK kinase 1/2, PI3K/Akt, and particularly Jun-NH2-terminal kinase signalling pathways

Altruism can proliferate through group/kin selection despite high random gene flow

The ways in which natural selection can allow the proliferation of cooperative behavior have long been seen as a central problem in evolutionary biology. Most of the literature has focused on interactions between pairs of individuals and on linear public goods games. This emphasis led to the conclusion that even modest levels of migration would pose a serious problem to the spread of altruism in group structured populations. Here we challenge this conclusion, by analyzing evolution in a framework which allows for complex group interactions and random migration among groups. We conclude that contingent forms of strong altruism can spread when rare under realistic group sizes and levels of migration. Our analysis combines group-centric and gene-centric perspectives, allows for arbitrary strength of selection, and leads to extensions of Hamilton's rule for the spread of altruistic alleles, applicable under broad conditions.Comment: 5 pages, 2 figures. Supplementary material with 50 pages and 26 figure

Structural and Chemical Profiling of the Human Cytosolic Sulfotransferases

The human cytosolic sulfotransfases (hSULTs) comprise a family of 12 phase II enzymes involved in the metabolism of drugs and hormones, the bioactivation of carcinogens, and the detoxification of xenobiotics. Knowledge of the structural and mechanistic basis of substrate specificity and activity is crucial for understanding steroid and hormone metabolism, drug sensitivity, pharmacogenomics, and response to environmental toxins. We have determined the crystal structures of five hSULTs for which structural information was lacking, and screened nine of the 12 hSULTs for binding and activity toward a panel of potential substrates and inhibitors, revealing unique “chemical fingerprints” for each protein. The family-wide analysis of the screening and structural data provides a comprehensive, high-level view of the determinants of substrate binding, the mechanisms of inhibition by substrates and environmental toxins, and the functions of the orphan family members SULT1C3 and SULT4A1. Evidence is provided for structural “priming” of the enzyme active site by cofactor binding, which influences the spectrum of small molecules that can bind to each enzyme. The data help explain substrate promiscuity in this family and, at the same time, reveal new similarities between hSULT family members that were previously unrecognized by sequence or structure comparison alone

The structure of the first representative of Pfam family PF06475 reveals a new fold with possible involvement in glycolipid metabolism

PA1994, a Pfam PF06475 (DUF1089) family homolog from P. aeruginosa, reveals remote similarities to lipoprotein localization factors and a conserved putative glycolipid-binding site

To hit or not to hit, that is the question -genome-wide structure-based druggability predictions for <i>pseudomonas aeruginosa </i>proteins

Pseudomonas aeruginosa is a Gram-negative bacterium known to cause opportunistic infections in immune-compromised or immunosuppressed individuals that often prove fatal. New drugs to combat this organism are therefore sought after. To this end, we subjected the gene products of predicted perturbative genes to structure-based druggability predictions using DrugPred. Making this approach suitable for large-scale predictions required the introduction of new methods for calculation of descriptors, development of a workflow to identify suitable pockets in homologous proteins and establishment of criteria to obtain valid druggability predictions based on homologs. We were able to identify 29 perturbative proteins of P. aeruginosa that may contain druggable pockets, including some of them with no or no drug-like inhibitors deposited in ChEMBL. These proteins form promising novel targets for drug discovery against P. aeruginosa