Alessandro Pandini

Ligand binding and activation of the Ah receptor

The Ah receptor (AhR) is a ligand-dependent transcription factor that can be activated by structurally diverse synthetic and naturally-occurring chemicals. Although a significant amount of information is available with respect to the planar aromatic hydrocarbon AhR ligands, the actual spectrum of chemicals that can bind to and activate the AhR is only now being elucidated. In addition, the lack of information regarding the actual three-dimensional structure of the AhR ligand binding domain (LBD) has hindered detailed analysis of the molecular mechanisms by which these ligands bind to and active AhR signal transduction. In this review we describe the current state of knowledge with respect to naturally occurring AhR ligands and present and discuss the first theoretical model of the AhR LBD based on crystal structures of homologous PAS family members.

Detection of allosteric signal transmission by information-theoretic analysis of protein dynamics

Allostery offers a highly specific way to modulate protein function. Therefore, understanding this mechanism is of increasing interest for protein science and drug discovery. However, allosteric signal transmission is difficult to detect experimentally and to model because it is often mediated by local structural changes propagating along multiple pathways. To address this, we developed a method to identify communication pathways by an information‐theoretical analysis of molecular dynamics simulations. Signal propagation was described as information exchange through a network of correlated local motions, modeled as transitions between canonical states of protein fragments. The method was used to describe allostery in two‐component regulatory systems. In particular, the transmission from the allosteric site to the signaling surface of the receiver domain NtrC was shown to be mediated by a layer of hub residues. The location of hubs preferentially connected to the allosteric site was found in close agreement with key residues experimentally identified as involved in the signal transmission. The comparison with the networks of the homologues CheY and FixJ highlighted similarities in their dynamics. In particular, we showed that a preorganized network of fragment connections between the allosteric and functional sites exists already in the inactive state of all three proteins.—Pandini, A., Fornili, A., Fraternali, F., Kleinjung, J. Detection of allosteric signal transmission by information‐theoretic analysis of protein dynamics. FASEB J. 26, 868–881 (2012). www.fasebj.org

Predicting the accuracy of protein-ligand docking on homology models

Ligand–protein docking is increasingly used in Drug Discovery. The initial limitations imposed by a reduced availability of target protein structures have been overcome by the use of theoretical models, especially those derived by homology modeling techniques. While this greatly extended the use of docking simulations, it also introduced the need for general and robust criteria to estimate the reliability of docking results given the model quality. To this end, a large‐scale experiment was performed on a diverse set including experimental structures and homology models for a group of representative ligand–protein complexes. A wide spectrum of model quality was sampled using templates at different evolutionary distances and different strategies for target–template alignment and modeling. The obtained models were scored by a selection of the most used model quality indices. The binding geometries were generated using AutoDock, one of the most common docking programs. An important result of this study is that indeed quantitative and robust correlations exist between the accuracy of docking results and the model quality, especially in the binding site. Moreover, state‐of‐the‐art indices for model quality assessment are already an effective tool for an a priori prediction of the accuracy of docking experiments in the context of groups of proteins with conserved structural characteristics.

Comparative analysis of homology models of the AH receptor ligand binding domain: verification of structure-function predictions by site-directed mutagenesis of a nonfunctional receptor.

The aryl hydrocarbon receptor (AHR) is a ligand-dependent transcription factor that mediates the biological and toxic effects of a wide variety of structurally diverse chemicals, including the toxic environmental contaminant 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). While significant interspecies differences in AHR ligand binding specificity, selectivity, and response have been observed, the structural determinants responsible for those differences have not been determined, and homology models of the AHR ligand-binding domain (LBD) are available for only a few species. Here we describe the development and comparative analysis of homology models of the LBD of 16 AHRs from 12 mammalian and nonmammalian species and identify the specific residues contained within their ligand binding cavities. The ligand-binding cavity of the fish AHR exhibits differences from those of mammalian and avian AHRs, suggesting a slightly different TCDD binding mode. Comparison of the internal cavity in the LBD model of zebrafish (zf) AHR2, which binds TCDD with high affinity, to that of zfAHR1a, which does not bind TCDD, revealed that the latter has a dramatically shortened binding cavity due to the side chains of three residues (Tyr296, Thr386, and His388) that reduce the amount of internal space available to TCDD. Mutagenesis of two of these residues in zfAHR1a to those present in zfAHR2 (Y296H and T386A) restored the ability of zfAHR1a to bind TCDD and to exhibit TCDD-dependent binding to DNA. These results demonstrate the importance of these two amino acids and highlight the predictive potential of comparative analysis of homology models from diverse species. The availability of these AHR LBD homology models will facilitate in-depth comparative studies of AHR ligand binding and ligand-dependent AHR activation and provide a novel avenue for examining species-specific differences in AHR responsiveness.

GSATools: analysis of allosteric communication and functional local motions using a structural alphabet

Motivation: GSATools is a free software package to analyze conformational ensembles and to detect functional motions in proteins by means of a structural alphabet. The software integrates with the widely used GROMACS simulation package and can generate a range of graphical outputs. Three applications can be supported: (i) investigation of the conformational variability of local structures; (ii) detection of allosteric communication; and (iii) identification of local regions that are critical for global functional motions. These analyses provide insights into the dynamics of proteins and allow for targeted design of functional mutants in theoretical and experimental studies. Availability: The C source code of the GSATools, along with a set of pre-compiled binaries, is freely available under GNU General Public License from http://mathbio.nimr.mrc.ac.uk/wiki/GSATools. Contact: alessandro.pandini@kcl.ac.uk or jkleinj@nimr.mrc.ac.uk Supplementary information: Supplementary data are available at Bioinformatics online.

Specialized dynamical properties of promiscuous residues revealed by simulated conformational ensembles

The ability to interact with different partners is one of the most important features in proteins. Proteins that bind a large number of partners (hubs) have been often associated with intrinsic disorder. However, many examples exist of hubs with an ordered structure, and evidence of a general mechanism promoting promiscuity in ordered proteins is still elusive. An intriguing hypothesis is that promiscuous binding sites have specific dynamical properties, distinct from the rest of the interface and pre-existing in the protein isolated state. Here, we present the first comprehensive study of the intrinsic dynamics of promiscuous residues in a large protein data set. Different computational methods, from coarse-grained elastic models to geometry-based sampling methods and to full-atom Molecular Dynamics simulations, were used to generate conformational ensembles for the isolated proteins. The flexibility and dynamic correlations of interface residues with a different degree of binding promiscuity were calculated and compared considering side chain and backbone motions, the latter both on a local and on a global scale. The study revealed that (a) promiscuous residues tend to be more flexible than nonpromiscuous ones, (b) this additional flexibility has a higher degree of organization, and (c) evolutionary conservation and binding promiscuity have opposite effects on intrinsic dynamics. Findings on simulated ensembles were also validated on ensembles of experimental structures extracted from the Protein Data Bank (PDB). Additionally, the low occurrence of single nucleotide polymorphisms observed for promiscuous residues indicated a tendency to preserve binding diversity at these positions. A case study on two ubiquitin-like proteins exemplifies how binding promiscuity in evolutionary related proteins can be modulated by the fine-tuning of the interface dynamics. The interplay between promiscuity and flexibility highlighted here can inspire new directions in protein–protein interaction prediction and design methods.

Conformational and functional analysis of molecular dynamics trajectories by self-organising maps.

BackgroundMolecular dynamics (MD) simulations are powerful tools to investigate the conformational dynamics of proteins that is often a critical element of their function. Identification of functionally relevant conformations is generally done clustering the large ensemble of structures that are generated. Recently, Self-Organising Maps (SOMs) were reported performing more accurately and providing more consistent results than traditional clustering algorithms in various data mining problems. We present a novel strategy to analyse and compare conformational ensembles of protein domains using a two-level approach that combines SOMs and hierarchical clustering.ResultsThe conformational dynamics of the α-spectrin SH3 protein domain and six single mutants were analysed by MD simulations. The Cαs Cartesian coordinates of conformations sampled in the essential space were used as input data vectors for SOM training, then complete linkage clustering was performed on the SOM prototype vectors. A specific protocol to optimize a SOM for structural ensembles was proposed: the optimal SOM was selected by means of a Taguchi experimental design plan applied to different data sets, and the optimal sampling rate of the MD trajectory was selected. The proposed two-level approach was applied to single trajectories of the SH3 domain independently as well as to groups of them at the same time. The results demonstrated the potential of this approach in the analysis of large ensembles of molecular structures: the possibility of producing a topological mapping of the conformational space in a simple 2D visualisation, as well as of effectively highlighting differences in the conformational dynamics directly related to biological functions.ConclusionsThe use of a two-level approach combining SOMs and hierarchical clustering for conformational analysis of structural ensembles of proteins was proposed. It can easily be extended to other study cases and to conformational ensembles from other sources.

MinSet: a general approach to derive maximally representative database subsets by using fragment dictionaries and its application to the SCOP database

MOTIVATION The size of current protein databases is a challenge for many Bioinformatics applications, both in terms of processing speed and information redundancy. It may be therefore desirable to efficiently reduce the database of interest to a maximally representative subset. RESULTS The MinSet method employs a combination of a Suffix Tree and a Genetic Algorithm for the generation, selection and assessment of database subsets. The approach is generally applicable to any type of string-encoded data, allowing for a drastic reduction of the database size whilst retaining most of the information contained in the original set. We demonstrate the performance of the method on a database of protein domain structures encoded as strings. We used the SCOP40 domain database by translating protein structures into character strings by means of a structural alphabet and by extracting optimized subsets according to an entropy score that is based on a constant-length fragment dictionary. Therefore, optimized subsets are maximally representative for the distribution and range of local structures. Subsets containing only 10% of the SCOP structure classes show a coverage of >90% for fragments of length 1-4. AVAILABILITY http://mathbio.nimr.mrc.ac.uk/~jkleinj/MinSet. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

The Gearbox of the bacterial flagellar motor switch

Summary Switching of flagellar motor rotation sense dictates bacterial chemotaxis. Multi-subunit FliM-FliG rotor rings couple signal protein binding in FliM with reversal of a distant FliG C-terminal (FliGC) helix involved in stator contacts. Subunit dynamics were examined in conformer ensembles generated by molecular simulations from the X-ray structures. Principal component analysis extracted collective motions. Interfacial loop immobilization by complex formation coupled elastic fluctuations of the FliM middle (FliMM) and FliG middle (FliGM) domains. Coevolved mutations captured interfacial dynamics as well as contacts. FliGM rotation was amplified via two central hinges to the FliGC helix. Intrinsic flexibility, reported by the FliGMC ensembles, reconciled conformers with opposite FliGC helix orientations. FliG domain stacking deformed the inter-domain linker and reduced flexibility; but conformational changes were not triggered by engineered linker deletions that cause a rotation-locked phenotype. These facts suggest that binary rotation states arise from conformational selection by stacking interactions.

BCR-ABL residues interacting with ponatinib are critical to preserve the tumorigenic potential of the oncoprotein

Patients with chronic myeloid leukemia in whom tyrosine kinase inhibitors (TKIs) fail often present mutations in the BCR‐ABL catalytic domain. We noticed a lack of substitutions involving 4 amino acids (E286, M318, I360, and D381) that form hydrogen bonds with ponatinib. We therefore introduced mutations in each of these residues, either preserving or altering their physicochemical properties. We found that E286, M318, I360, and D381 are dispensable for ABL and BCR‐ABL protein stability but are critical for preserving catalytic activity. Indeed, only a “conservative” I360T substitution retained kinase proficiency and transforming potential. Molecular dynamics simulations of BCR‐ABLI360T revealed differences in both helix αC dynamics and protein‐correlated motions, consistent with a modified ATP‐binding pocket. Nevertheless, this mutant remained sensitive to ponatinib, imatinib, and dasatinib. These results suggest that changes in the 4 BCR‐ABL residues described here would be selected against by a lack of kinase activity or by maintained responsiveness to TKIs. Notably, amino acids equivalent to those identified in BCR‐ABL are conserved in 51% of human tyrosine kinases. Hence, these residues may represent an appealing target for the design of pharmacological compounds that would inhibit additional oncogenic tyrosine kinases while avoiding the emergence of resistance due to point mutations.—Buffa, P., Romano, C., Pandini, A., Massimino, M., Tirrò, E., Di Raimondo, F., Manzella, L., Fraternali, F., Vigneri, P. G. BCR‐ABL residues interacting with ponatinib are critical to preserve the tumorigenic potential of the oncoprotein. FASEB J. 28, 1221–1236 (2014). www.fasebj.org