DIMA 3.0: Domain Interaction Map

Domain Interaction MAp (DIMA, available at http://webclu.bio.wzw.tum.de/dima) is a database of predicted and known interactions between protein domains. It integrates 5807 structurally known interactions imported from the iPfam and 3did databases and 46 900 domain interactions predicted by four computational methods: domain phylogenetic profiling, domain pair exclusion algorithm correlated mutations and domain interaction prediction in a discriminative way. Additionally predictions are filtered to exclude those domain pairs that are reported as non-interacting by the Negatome database. The DIMA Web site allows to calculate domain interaction networks either for a domain of interest or for entire organisms, and to explore them interactively using the Flash-based Cytoscape Web software

Structural Characterization of the Integrin Αiibβ3 Transmembrane and Cytosolic Domains

Integrin’s are the principal cell surface receptors that link the cytoskeleton to the extracellular matrix. They exist in active conformations that can bind extracellular ligands and resting conformations that cannot. The platelet integrin αIIbβ3 is a prototypical regulated integrin that is resting on a circulating platelet and becomes activated to adhere the platelet to the vascular endothelium or subendothelial matrix. The integrin is composed of α and β subunits and each subunit contains a single transmembrane helix that form an α/β heterodimer in the resting state. Additionally, each subunit contains a cytosolic domain that binds signaling proteins that affect the resting-active equilibrium. Activation signals are transduced across the membrane by separating the transmembrane heterodimer. The structure of the resting integrin αIIbβ3’s transmembrane and cytosolic domains was characterized by molecular modeling and NMR spectroscopy. First, software was developed to model transmembrane helix dimers using experimental mutagenesis results as a modeling restraint. Next, the αIIb/β3 transmembrane heterodimer was modeled and the model was compared to published experimental data and other published models. The model correlated well with experimental findings and converged on the same structure as other top performing models, suggesting this conformation approximates the native interface. The model’s interface includes αIIb residue Met987 and β3 residue Leu712. These residues were mutated to cysteine to crosslink peptides corresponding to the αIIb and β3 cytosolic tails, and the disulfide-linked construct was probed by NMR spectroscopy. NMR revealed that the αIIb and β3 cytosolic tails have a dynamic interface. The αIIb subunit is natively unstructured and the β3 subunit consists of a hydrophobic helix followed by two amphiphilic helices. The amphiphilic portions of β3 include domains that interact with cytosolic proteins, but the membrane embedding of its hydrophobic faces sequesters some of the interacting residues. This result suggests that the integrin’s resting-active equilibrium is coupled to an equilibrium between membrane embedded and solvent exposed conformations of the β3 cytosolic tail, providing new insight into integrin activation

Deconvolving mutational patterns of poliovirus outbreaks reveals its intrinsic fitness landscape.

Vaccination has essentially eradicated poliovirus. Yet, its mutation rate is higher than that of viruses like HIV, for which no effective vaccine exists. To investigate this, we infer a fitness model for the poliovirus viral protein 1 (vp1), which successfully predicts in vitro fitness measurements. This is achieved by first developing a probabilistic model for the prevalence of vp1 sequences that enables us to isolate and remove data that are subject to strong vaccine-derived biases. The intrinsic fitness constraints derived for vp1, a capsid protein subject to antibody responses, are compared with those of analogous HIV proteins. We find that vp1 evolution is subject to tighter constraints, limiting its ability to evade vaccine-induced immune responses. Our analysis also indicates that circulating poliovirus strains in unimmunized populations serve as a reservoir that can seed outbreaks in spatio-temporally localized sub-optimally immunized populations

Mutantelec: AnIn Silicomutation simulation platform for comparative electrostatic potential profiling of proteins

The electrostatic potential plays a key role in many biological processes like determining the affinity of a ligand to a given protein target, and they are responsible for the catalytic activity of many enzymes. Understanding the effect that amino acid mutations will have on the electrostatic potential of a protein, will allow a thorough understanding of which residues are the most important in a protein. MutantElec, is a friendly web application for in silico generation of site-directed mutagenesis of proteins and the comparison of electrostatic potential between the wild type protein and the mutant(s), based on the three-dimensional structure of the protein. The effect of the mutation is evaluated using different approach to the traditional surface map. MutantElec provides a graphical display of the results that allows the visualization of changes occurring at close distance from the mutation and thus uncovers the local and global impact of a specific chang

Cyclin D-Cdk4,6 Drives Cell-Cycle Progression via the Retinoblastoma Protein's C-Terminal Helix.

The cyclin-dependent kinases Cdk4 and Cdk6 form complexes with D-type cyclins to drive cell proliferation. A well-known target of cyclin D-Cdk4,6 is the retinoblastoma protein Rb, which inhibits cell-cycle progression until its inactivation by phosphorylation. However, the role of Rb phosphorylation by cyclin D-Cdk4,6 in cell-cycle progression is unclear because Rb can be phosphorylated by other cyclin-Cdks, and cyclin D-Cdk4,6 has other targets involved in cell division. Here, we show that cyclin D-Cdk4,6 docks one side of an alpha-helix in the Rb C terminus, which is not recognized by cyclins E, A, and B. This helix-based docking mechanism is shared by the p107 and p130 Rb-family members across metazoans. Mutation of the Rb C-terminal helix prevents its phosphorylation, promotes G1 arrest, and enhances Rb's tumor suppressive function. Our work conclusively demonstrates that the cyclin D-Rb interaction drives cell division and expands the diversity of known cyclin-based protein docking mechanisms

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. © 2013 American Chemical Society

Selection of sequence motifs and generative Hopfield-Potts models for protein familiesilies

Statistical models for families of evolutionary related proteins have recently gained interest: in particular pairwise Potts models, as those inferred by the Direct-Coupling Analysis, have been able to extract information about the three-dimensional structure of folded proteins, and about the effect of amino-acid substitutions in proteins. These models are typically requested to reproduce the one- and two-point statistics of the amino-acid usage in a protein family, {\em i.e.}~to capture the so-called residue conservation and covariation statistics of proteins of common evolutionary origin. Pairwise Potts models are the maximum-entropy models achieving this. While being successful, these models depend on huge numbers of {\em ad hoc} introduced parameters, which have to be estimated from finite amount of data and whose biophysical interpretation remains unclear. Here we propose an approach to parameter reduction, which is based on selecting collective sequence motifs. It naturally leads to the formulation of statistical sequence models in terms of Hopfield-Potts models. These models can be accurately inferred using a mapping to restricted Boltzmann machines and persistent contrastive divergence. We show that, when applied to protein data, even 20-40 patterns are sufficient to obtain statistically close-to-generative models. The Hopfield patterns form interpretable sequence motifs and may be used to clusterize amino-acid sequences into functional sub-families. However, the distributed collective nature of these motifs intrinsically limits the ability of Hopfield-Potts models in predicting contact maps, showing the necessity of developing models going beyond the Hopfield-Potts models discussed here.Comment: 26 pages, 16 figures, to app. in PR

First-Step Mutations for Adaptation at Elevated Temperature Increase Capsid Stability in a Virus

The relationship between mutation, protein stability and protein function plays a central role in molecular evolution. Mutations tend to be destabilizing, including those that would confer novel functions such as host-switching or antibiotic resistance. Elevated temperature may play an important role in preadapting a protein for such novel functions by selecting for stabilizing mutations. In this study, we test the stability change conferred by single mutations that arise in a G4-like bacteriophage adapting to elevated temperature. The vast majority of these mutations map to interfaces between viral coat proteins, suggesting they affect protein-protein interactions. We assess their effects by estimating thermodynamic stability using molecular dynamic simulations and measuring kinetic stability using experimental decay assays. The results indicate that most, though not all, of the observed mutations are stabilizing