Unbiased Protein Interface Prediction Based on Ligand Diversity Quantification

Proteins interact with each other to perform essential functions in cells. Consequently, identification of their binding interfaces can provide key information for drug design. Here, we introduce Weighted Protein Interface Prediction (WePIP), an original framework which predicts protein interfaces from homologous complexes. WePIP takes advantage of a novel weighted score which is not only based on structural neighbours' information but, unlike current state-of-the-art methods, also takes into consideration the nature of their interaction partners. Experimental validation demonstrates that our weighted schema significantly improves prediction performance. In particular, we have established a major contribution to ligand diversity quantification. Moreover, application of our framework on a standard dataset shows WePIP performance compares favourably with other state of the art methods

Virtual screening for inhibitors of the human TSLP:TSLPR interaction

The pro-inflammatory cytokine thymic stromal lymphopoietin (TSLP) plays a pivotal role in the pathophysiology of various allergy disorders that are mediated by type 2 helper T cell (Th2) responses, such as asthma and atopic dermatitis. TSLP forms a ternary complex with the TSLP receptor (TSLPR) and the interleukin-7-receptor subunit alpha (IL-7Ra), thereby activating a signaling cascade that culminates in the release of pro-inflammatory mediators. In this study, we conducted an in silico characterization of the TSLP: TSLPR complex to investigate the drugability of this complex. Two commercially available fragment libraries were screened computationally for possible inhibitors and a selection of fragments was subsequently tested in vitro. The screening setup consisted of two orthogonal assays measuring TSLP binding to TSLPR: a BLI-based assay and a biochemical assay based on a TSLP: alkaline phosphatase fusion protein. Four fragments pertaining to diverse chemical classes were identified to reduce TSLP: TSLPR complex formation to less than 75% in millimolar concentrations. We have used unbiased molecular dynamics simulations to develop a Markov state model that characterized the binding pathway of the most interesting compound. This work provides a proof-ofprinciple for use of fragments in the inhibition of TSLP: TSLPR complexation

A versatile Halo- and SNAP-tagged BMP/TGFβ receptor library for quantification of cell surface ligand binding

TGFβs, BMPs and Activins regulate numerous developmental and homeostatic processes and signal through hetero-tetrameric receptor complexes composed of two types of serine/threonine kinase receptors. Each of the 33 different ligands possesses unique affinities towards specific receptor types. However, the lack of specific tools hampered simultaneous testing of ligand binding towards all BMP/TGFβ receptors. Here we present a N-terminally Halo- and SNAP-tagged TGFβ/BMP receptor library to visualize receptor complexes in dual color. In combination with fluorescently labeled ligands, we established a Ligand Surface Binding Assay (LSBA) for optical quantification of receptor-dependent ligand binding in a cellular context. We highlight that LSBA is generally applicable to test (i) binding of different ligands such as Activin A, TGFβ1 and BMP9, (ii) for mutant screens and (iii) evolutionary comparisons. This experimental set-up opens opportunities for visualizing ligand-receptor binding dynamics, essential to determine signaling specificity and is easily adaptable for other receptor signaling pathways

Scoring docking conformations using predicted protein interfaces

BACKGROUND: Since proteins function by interacting with other molecules, analysis of protein-protein interactions is essential for comprehending biological processes. Whereas understanding of atomic interactions within a complex is especially useful for drug design, limitations of experimental techniques have restricted their practical use. Despite progress in docking predictions, there is still room for improvement. In this study, we contribute to this topic by proposing T-PioDock, a framework for detection of a native-like docked complex 3D structure. T-PioDock supports the identification of near-native conformations from 3D models that docking software produced by scoring those models using binding interfaces predicted by the interface predictor, Template based Protein Interface Prediction (T-PIP). RESULTS: First, exhaustive evaluation of interface predictors demonstrates that T-PIP, whose predictions are customised to target complexity, is a state-of-the-art method. Second, comparative study between T-PioDock and other state-of-the-art scoring methods establishes T-PioDock as the best performing approach. Moreover, there is good correlation between T-PioDock performance and quality of docking models, which suggests that progress in docking will lead to even better results at recognising near-native conformations. CONCLUSION: Accurate identification of near-native conformations remains a challenging task. Although availability of 3D complexes will benefit from template-based methods such as T-PioDock, we have identified specific limitations which need to be addressed. First, docking software are still not able to produce native like models for every target. Second, current interface predictors do not explicitly consider pairwise residue interactions between proteins and their interacting partners which leaves ambiguity when assessing quality of complex conformations

Using single-cell technologies to map the human immune system - implications for nephrology.

Advances in single-cell technologies are transforming our understanding of cellular identity. For instance, the application of single-cell RNA sequencing and mass cytometry technologies to the study of immune cell populations in blood, secondary lymphoid organs and the renal tract is helping researchers to map the complex immune landscape within the kidney, define cell ontogeny and understand the relationship of kidney-resident immune cells with their circulating counterparts. These studies also provide insights into the interactions of immune cell populations with neighbouring epithelial and endothelial cells in health, and across a range of kidney diseases and cancer. These data have translational potential and will aid the identification of drug targets and enable better prediction of off-target effects. The application of single-cell technologies to clinical renal biopsy samples, or even cells within urine, will improve diagnostic accuracy and assist with personalized prognostication for patients with various kidney diseases. A comparison of immune cell types in peripheral blood and secondary lymphoid organs in healthy individuals and in patients with systemic autoimmune diseases that affect the kidney will also help to unravel the mechanisms that underpin the breakdown in self-tolerance and propagation of autoimmune responses. Together, these immune cell atlases have the potential to transform nephrology

Thirty years of molecular dynamics simulations on posttranslational modifications of proteins

Posttranslational modifications (PTMs) are an integral component to how cells respond to perturbation. While experimental advances have enabled improved PTM identification capabilities, the same throughput for characterizing how structural changes caused by PTMs equate to altered physiological function has not been maintained. In this Perspective, we cover the history of computational modeling and molecular dynamics simulations which have characterized the structural implications of PTMs. We distinguish results from different molecular dynamics studies based upon the timescales simulated and analysis approaches used for PTM characterization. Lastly, we offer insights into how opportunities for modern research efforts on in silico PTM characterization may proceed given current state-of-the-art computing capabilities and methodological advancements.Comment: 64 pages, 11 figure

Solvents to Fragments to Drugs: MD Applications in Drug Design

Simulations of molecular dynamics (MD) are playing an increasingly important role in structure-based drug discovery (SBDD). Here we review the use of MD for proteins in aqueous solvation, organic/aqueous mixed solvents (MDmix) and with small ligands, to the classic SBDD problems: Binding mode and binding free energy predictions. The simulation of proteins in their condensed state reveals solvent structures and preferential interaction sites (hot spots) on the protein surface. The information provided by water and its cosolvents can be used very effectively to understand protein ligand recognition and to improve the predictive capability of well-established methods such as molecular docking. The application of MD simulations to the study of the association of proteins with drug-like compounds is currently only possible for specific cases, as it remains computationally very expensive and labor intensive. MDmix simulations on the other hand, can be used systematically to address some of the common tasks in SBDD. With the advent of new tools and faster computers we expect to see an increase in the application of mixed solvent MD simulations to a plethora of protein targets to identify new drug candidates

Solvents to fragments to drugs: MD applications in drug design

Simulations of molecular dynamics (MD) are playing an increasingly important role in structure-based drug discovery (SBDD). Here we review the use of MD for proteins in aqueous solvation, organic/aqueous mixed solvents (MDmix) and with small ligands, to the classic SBDD problems: Binding mode and binding free energy predictions. The simulation of proteins in their condensed state reveals solvent structures and preferential interaction sites (hot spots) on the protein surface. The information provided by water and its cosolvents can be used very effectively to understand protein ligand recognition and to improve the predictive capability of well-established methods such as molecular docking. The application of MD simulations to the study of the association of proteins with drug-like compounds is currently only possible for specific cases, as it remains computationally very expensive and labor intensive. MDmix simulations on the other hand, can be used systematically to address some of the common tasks in SBDD. With the advent of new tools and faster computers we expect to see an increase in the application of mixed solvent MD simulations to a plethora of protein targets to identify new drug candidates.Fil: Defelipe, Lucas Alfredo. Universidad de Buenos Aires; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales; ArgentinaFil: Arcon, Juan Pablo. Universidad de Buenos Aires; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales; ArgentinaFil: Modenutti, Carlos Pablo. Universidad de Buenos Aires; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales; ArgentinaFil: Marti, Marcelo Adrian. Universidad de Buenos Aires; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales; ArgentinaFil: Turjanski, Adrian. Universidad de Buenos Aires; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales; ArgentinaFil: Barril, Xavier. Institucio Catalana de Recerca I Estudis Avancats

Deciphering Chemomechanical Couplings in Proteins Using Microsecond-level Molecular Dynamics Simulations

All-atom molecular dynamics (MD) simulations combine the high temporal resolution of experimental methods like smFRET and spatial resolution of methods like x-ray crystallography, to provide a detailed dynamic picture of biomolecular processes. Here, microsecond-level atomistic MD simulations have been used to characterize chemomechanical couplings in human fibroblast growth factor 1 (hFGF1) and the spike proteins of SARS CoV-1 and SARS-CoV-2. hFGF1 is a globular signaling protein that is involved in several physiological processes ranging from cell proliferation to wound healing. Experimental studies have previously described the low proteolytic and thermal stability of hFGF1, in addition to the stabilizing role of heparin. Here, a conformational change in the hFGF1 heparin-binding pocket that occurs only when heparin is absent, is described for the first time. Comparisons with experimental data indicate that this conformational transition is implicated in the low thermal stability of hFGF1. Unique electrostatic interactions that contribute to heparin-mediated stabilization are also described. This work also describes a novel binding affinity estimation approach involving restrained umbrella sampling simulations. The absolute binding affinity for the hFGF1-heparin interaction determined using this approach is in very good agreement with data from isothermal titration calorimetry (ITC) experiments. This binding affinity study revealed that restraining ligand orientation is essential for effective sampling along a protein-ligand distance collective variable.The differential dynamic behavior of the SARS-CoV-1 and CoV-2 spike proteins is also described in this work. Spike protein activation is the first step in the “effective binding” process leading to interaction with the human ACE2 receptor. This study shows that the active form of the CoV-1 spike protein is less stable than that of the CoV-2 spike protein and that the energy barriers associated with activation and inactivation are higher in CoV-2. A “pseudo-inactive” state of the CoV-1 spike protein is described for the first time, wherein the N-terminal domain (NTD) interacts with the receptor-binding domain (RBD). This highlights the potential role of the NTD in spike protein inactivation. The relatively slower kinetics of spike protein activation and inactivation in CoV-2 indicate that it might spend more time bound to the ACE2 receptor than CoV-1, which in turn might provide an explanation for the higher transmissibility of CoV-2