Match_Motif: A rapid computational tool to assist in protein–protein interaction design

In order to generate protein assemblies with a desired function, the rational design of protein–protein binding interfaces is of significant interest. Approaches based on random mutagenesis or directed evolution may involve complex experimental selection procedures. Also, molecular modeling approaches to design entirely new proteins and interactions with partner molecules can involve large computational efforts and screening steps. In order to simplify at least the initial effort for designing a putative binding interface between two proteins the Match_Motif approach has been developed. It employs the large collection of known protein–protein complex structures to suggest interface modifications that may lead to improved binding for a desired input interaction geometry. The approach extracts interaction motifs based on the backbone structure of short (four residues) segments and the relative arrangement with respect to short segments on the partner protein. The interaction geometry is used to search through a database of such motifs in known stable bound complexes. All matches are rapidly identified (within a few seconds) and collected and can be used to guide changes in the interface that may lead to improved binding. In the output, an alternative interface structure is also proposed based on the frequency of occurrence of side chains at a given interface position in all matches and based on sterical considerations. Applications of the procedure to known complex structures and alternative arrangements are presented and discussed. The program, data files, and example applications can be downloaded from https://www.groups.ph.tum.de/t38/downloads/.     

Two distinct classes of cochaperones compete for the EEVD motif in heat shock protein 70 to tune its chaperone activities

Chaperones of the heat shock protein 70 (Hsp70) family engage in protein–protein interactions with many cochaperones. One “hotspot” for cochaperone binding is the EEVD motif, found at the extreme C terminus of cytoplasmic Hsp70s. This motif is known to bind tetratricopeptide repeat domain cochaperones, such as the E3 ubiquitin ligase CHIP. In addition, the EEVD motif also interacts with a structurally distinct domain that is present in class B J-domain proteins, such as DnaJB4. These observations suggest that CHIP and DnaJB4 might compete for binding to Hsp70’s EEVD motif; however, the molecular determinants of such competition are not clear. Using a collection of EEVD-derived peptides, including mutations and truncations, we explored which residues are critical for binding to both CHIP and DnaJB4. These results revealed that some features, such as the C-terminal carboxylate, are important for both interactions. However, CHIP and DnaJB4 also had unique preferences, especially at the isoleucine position immediately adjacent to the EEVD. Finally, we show that competition between these cochaperones is important in vitro, as DnaJB4 limits the ubiquitination activity of the Hsp70–CHIP complex, whereas CHIP suppresses the client refolding activity of the Hsp70–DnaJB4 complex. Together, these data suggest that the EEVD motif has evolved to support diverse protein–protein interactions, such that competition between cochaperones may help guide whether Hsp70-bound proteins are folded or degraded.     

Letting go: Deep computational modeling insights into pH-dependent calcium affinity

Calcium and other cofactors can feature as key additions to a molecular interface, to the extent that the cofactor is completely buried in the bound state. How can such an interaction be regulated then? The answer: By facilitating a switch through an allosteric network. Although a number of unbinding mechanisms are being characterized, an extensive computational study by Joswig et al. reveals a detailed model for the pattern recognition receptor langerin.     

Frequent Assembly of Chimeric Complexes in the Protein Interaction Network of an Interspecies Yeast Hybrid

Hybrids between species often show extreme phenotypes, including some that take place at the molecular level. In this study, we investigated the phenotypes of an interspecies diploid hybrid in terms of protein–protein interactions inferred from protein correlation profiling. We used two yeast species, Saccharomyces cerevisiae and Saccharomyces uvarum, which are interfertile, but yet have proteins diverged enough to be differentiated using mass spectrometry. Most of the protein–protein interactions are similar between hybrid and parents, and are consistent with the assembly of chimeric complexes, which we validated using an orthogonal approach for the prefoldin complex. We also identified instances of altered protein–protein interactions in the hybrid, for instance, in complexes related to proteostasis and in mitochondrial protein complexes. Overall, this study uncovers the likely frequent occurrence of chimeric protein complexes with few exceptions, which may result from incompatibilities or imbalances between the parental proteomes.     

Mechanistic insights into the inhibition of the CRISPR-Cas surveillance complex by anti-CRISPR protein AcrIF13

Clustered regularly interspaced short palindromic repeats (CRISPRs) and CRISPR-associated (Cas) proteins provide prokaryotes with nucleic acid–based adaptive immunity against infections of mobile genetic elements, including phages. To counteract this immune process, phages have evolved various anti-CRISPR (Acr) proteins which deactivate CRISPR-Cas–based immunity. However, the mechanisms of many of these Acr-mediated inhibitions are not clear. Here, we report the crystal structure of AcrIF13 and explore its inhibition mechanism. The structure of AcrIF13 is unique and displays a negatively charged surface. Additionally, biochemical studies identified that AcrIF13 interacts with the type I-F CRISPR-Cas surveillance complex (Csy complex) to block target DNA recognition and that the Cas5f-8f tail and Cas7.6f subunit of the Csy complex are specific binding targets of AcrIF13. Further mutational studies demonstrated that several negatively charged residues of AcrIF13 and positively charged residues of Cas8f and Cas7f of the Csy complex are involved in AcrIF13–Csy binding. Together, our findings provide mechanistic insights into the inhibition mechanism of AcrIF13 and further suggest the prevalence of the function of Acr proteins as DNA mimics.     

Structural and mechanistic insights into the inhibition of type I-F CRISPR-Cas system by anti-CRISPR protein AcrIF23

Prokaryotes evolved clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated (Cas) proteins as a kind of adaptive immune defense against mobile genetic elements including harmful phages. To counteract this defense, many mobile genetic elements in turn encode anti-CRISPR proteins (Acrs) to inactivate the CRISPR-Cas system. While multiple mechanisms of Acrs have been uncovered, it remains unknown whether other mechanisms are utilized by uncharacterized Acrs. Here, we report a novel mechanism adopted by recently identified AcrIF23. We show that AcrIF23 interacts with the Cas2/3 helicase-nuclease in the type I-F CRISPR-Cas system, similar to AcrIF3. The structure of AcrIF23 demonstrated a novel fold and structure-based mutagenesis identified a surface region of AcrIF23 involved in both Cas2/3-binding and its inhibition capacity. Unlike AcrIF3, however, we found AcrIF23 only potently inhibits the DNA cleavage activity of Cas2/3 but does not hinder the recruitment of Cas2/3 to the CRISPR RNA-guided surveillance complex (the Csy complex). Also, in contrast to AcrIF3 which hinders substrate DNA recognition by Cas2/3, we show AcrIF23 promotes DNA binding to Cas2/3. Taken together, our study identifies a novel anti-CRISPR mechanism used by AcrIF23 and highlights the diverse mechanisms adopted by Acrs.     

Efficient molecular docking of NMR structures: application to HIV-1 protease

Docking ligands into an ensemble of NMR conformers is essential to structure-based drug discovery if only NMR structures are available for the target. However, sequentially docking ligands into each NMR conformer through standard single-receptor-structure docking, referred to as sequential docking, is computationally expensive for large-scale database screening because of the large number of NMR conformers involved. Recently, we developed an efficient ensemble docking algorithm to consider protein structural variations in ligand binding. The algorithm simultaneously docks ligands into an ensemble of protein structures and achieves comparable performance to sequential docking without significant increase in computational time over single-structure docking. Here, we applied this algorithm to docking with NMR structures. The HIV-1 protease was used for validation in terms of docking accuracy and virtual screening. Ensemble docking of the NMR structures identified 91% of the known inhibitors under the criterion of RMSD < 2.0 A for the best-scored conformation, higher than the average success rate of single docking of individual crystal structures (66%). In the virtual screening test, on average, ensemble docking of the NMR structures obtained higher enrichments than single-structure docking of the crystal structures. In contrast, docking of either the NMR minimized average structure or a single NMR conformer performed less satisfactorily on both binding mode prediction and virtual screening, indicating that a single NMR structure may not be suitable for docking calculations. The success of ensemble docking of the NMR structures suggests an efficient alternative method for standard single docking of crystal structures and for considering protein flexibility.     

Resilience of death: intrinsic disorder in proteins involved in the programmed cell death

It is recognized now that intrinsically disordered proteins (IDPs), which do not have unique 3D structures as a whole or in noticeable parts, constitute a significant fraction of any given proteome. IDPs are characterized by an astonishing structural and functional diversity that defines their ability to be universal regulators of various cellular pathways. Programmed cell death (PCD) is one of the most intricate cellular processes where the cell uses specialized cellular machinery and intracellular programs to kill itself. This cell-suicide mechanism enables metazoans to control cell numbers and to eliminate cells that threaten the animal''s survival. PCD includes several specific modules, such as apoptosis, autophagy, and programmed necrosis (necroptosis). These modules are not only tightly regulated but also intimately interconnected and are jointly controlled via a complex set of protein–protein interactions. To understand the role of the intrinsic disorder in controlling and regulating the PCD, several large sets of PCD-related proteins across 28 species were analyzed using a wide array of modern bioinformatics tools. This study indicates that the intrinsic disorder phenomenon has to be taken into consideration to generate a complete picture of the interconnected processes, pathways, and modules that determine the essence of the PCD. We demonstrate that proteins involved in regulation and execution of PCD possess substantial amount of intrinsic disorder. We annotate functional roles of disorder across and within apoptosis, autophagy, and necroptosis processes. Disordered regions are shown to be implemented in a number of crucial functions, such as protein–protein interactions, interactions with other partners including nucleic acids and other ligands, are enriched in post-translational modification sites, and are characterized by specific evolutionary patterns. We mapped the disorder into an integrated network of PCD pathways and into the interactomes of selected proteins that are involved in the p53-mediated apoptotic signaling pathway.     

3D cryo-EM imaging of bacterial flagella: Novel structural and mechanistic insights into cell motility

Bacterial flagella are nanomachines that enable cells to move at high speeds. Comprising 25 and more different types of proteins, the flagellum is a large supramolecular assembly organized into three widely conserved substructures: a basal body including the rotary motor, a connecting hook, and a long filament. The whole flagellum from Escherichia coli weighs ∼20 MDa, without considering its filament portion, which is by itself a ∼1.6 GDa structure arranged as a multimer of ∼30,000 flagellin protomers. Breakthroughs regarding flagellar structure and function have been achieved in the last few years, mainly because of the revolutionary improvements in 3D cryo-EM methods. This review discusses novel structures and mechanistic insights derived from such high-resolution studies, advancing our understanding of each one of the three major flagellar segments. The rotation mechanism of the motor has been unveiled with unprecedented detail, showing a two-cogwheel machine propelled by a Brownian ratchet device. In addition, by imaging the flagellin-like protomers that make up the hook in its native bent configuration, their unexpected conformational plasticity challenges the paradigm of a two-state conformational rearrangement mechanism for flagellin-fold proteins. Finally, imaging of the filaments of periplasmic flagella, which endow Spirochete bacteria with their singular motility style, uncovered a strikingly asymmetric protein sheath that coats the flagellin core, challenging the view of filaments as simple homopolymeric structures that work as freely whirling whips. Further research will shed more light on the functional details of this amazing nanomachine, but our current understanding has definitely come a long way.     

Chemical cross-linking and mass spectrometry to determine the subunit interaction network in a recombinant human SAGA HAT subcomplex

Understanding the way how proteins interact with each other to form transient or stable protein complexes is a key aspect in structural biology. In this study, we combined chemical cross-linking with mass spectrometry to determine the binding stoichiometry and map the protein–protein interaction network of a human SAGA HAT subcomplex. MALDI-MS equipped with high mass detection was used to follow the cross-linking reaction using bis[sulfosuccinimidyl] suberate (BS3) and confirm the heterotetrameric stoichiometry of the specific stabilized subcomplex. Cross-linking with isotopically labeled BS3 d0-d4 followed by trypsin digestion allowed the identification of intra- and intercross-linked peptides using two dedicated search engines: pLink and xQuest. The identified interlinked peptides suggest a strong network of interaction between GCN5, ADA2B and ADA3 subunits; SGF29 is interacting with GCN5 and ADA3 but not with ADA2B. These restraint data were combined to molecular modeling and a low-resolution interacting model for the human SAGA HAT subcomplex could be proposed, illustrating the potential of an integrative strategy using cross-linking and mass spectrometry for addressing the structural architecture of multiprotein complexes.     

Computational design of small molecular modulators of protein–protein interactions with a novel thermodynamic cycle: Allosteric inhibitors of HIV‐1 integrase

Targeting protein–protein interactions for therapeutic development involves designing small molecules to either disrupt or enhance a known PPI. For this purpose, it is necessary to compute reliably the effect of chemical modifications of small molecules on the protein–protein association free energy. Here we present results obtained using a novel thermodynamic free energy cycle, for the rational design of allosteric inhibitors of HIV‐1 integrase (ALLINI) that act specifically in the early stage of the infection cycle. The new compounds can serve as molecular probes to dissect the multifunctional mechanisms of ALLINIs to inform the discovery of new allosteric inhibitors. The free energy protocol developed here can be more broadly applied to study quantitatively the effects of small molecules on modulating the strengths of protein–protein interactions.     

Deamidation drives molecular aging of the SARS-CoV-2 spike protein receptor-binding motif

The spike protein is the main protein component of the SARS-CoV-2 virion surface. The spike receptor-binding motif mediates recognition of the human angiotensin-converting enzyme 2 receptor, a critical step in infection, and is the preferential target for spike-neutralizing antibodies. Posttranslational modifications of the spike receptor-binding motif have been shown to modulate viral infectivity and host immune response, but these modifications are still being explored. Here we studied asparagine deamidation of the spike protein, a spontaneous event that leads to the appearance of aspartic and isoaspartic residues, which affect both the protein backbone and its charge. We used computational prediction and biochemical experiments to identify five deamidation hotspots in the SARS-CoV-2 spike protein. Asparagine residues 481 and 501 in the receptor-binding motif deamidate with a half-life of 16.5 and 123 days at 37 °C, respectively. Deamidation is significantly slowed at 4 °C, indicating a strong dependence of spike protein molecular aging on environmental conditions. Deamidation of the spike receptor-binding motif decreases the equilibrium constant for binding to the human angiotensin-converting enzyme 2 receptor more than 3.5-fold, yet its high conservation pattern suggests some positive effect on viral fitness. We propose a model for deamidation of the full SARS-CoV-2 virion illustrating how deamidation of the spike receptor-binding motif could lead to the accumulation on the virion surface of a nonnegligible chemically diverse spike population in a timescale of days. Our findings provide a potential mechanism for molecular aging of the spike protein with significant consequences for understanding virus infectivity and vaccine development.     

The 14-3-3γ isoform binds to and regulates the localization of endoplasmic reticulum (ER) membrane protein TMCC3 for the reticular network of the ER

The reticular network of the endoplasmic reticulum (ER) is formed by connecting ER tubules through three-way junctions and undergoes constant remodeling through formation and loss of the three-way junctions. Transmembrane and coiled-coil domain family 3 (TMCC3), an ER membrane protein localizing at three-way junctions, has been shown to positively regulate formation of the reticular ER network. However, elements that negatively regulate TMCC3 localization have not been characterized. In this study, we report that 14-3-3γ, a phospho-serine/phospho-threonine-binding protein involved in various signal transduction pathways, is a negative regulator of TMCC3. We demonstrate that overexpression of 14-3-3γ reduced localization of TMCC3 to three-way junctions and decreased the number of three-way junctions. TMCC3 bound to 14-3-3γ through the N terminus and had deduced 14-3-3 binding motifs. Additionally, we determined that a TMCC3 mutant substituting alanine for serine to be phosphorylated in the binding motif reduced binding to 14-3-3γ. The TMCC3 mutant was more prone than wildtype TMCC3 to localize at three-way junctions in the cells overexpressing 14-3-3γ. Furthermore, the TMCC3 mutant rescued the ER sheet expansion caused by TMCC3 knockdown less than wild-type TMCC3. Taken together, these results indicate that 14-3-3γ binding negatively regulates localization of TMCC3 to the three-way junctions for the proper reticular ER network, implying that the negative regulation of TMCC3 by 14-3-3γ would underlie remodeling of the reticular network of the ER.     

Structural insights into the role of the WW2 domain on tandem WW–PPxY motif interactions of oxidoreductase WWOX

Class I WW domains are present in many proteins of various functions and mediate protein interactions by binding to short linear PPxY motifs. Tandem WW domains often bind peptides with multiple PPxY motifs, but the interplay of WW–peptide interactions is not always intuitive. The WW domain–containing oxidoreductase (WWOX) harbors two WW domains: an unstable WW1 capable of PPxY binding and stable WW2 that cannot bind PPxY. The WW2 domain has been suggested to act as a WW1 domain chaperone, but the underlying mechanism of its chaperone activity remains to be revealed. Here, we combined NMR, isothermal calorimetry, and structural modeling to elucidate the roles of both WW domains in WWOX binding to its PPxY-containing substrate ErbB4. Using NMR, we identified an interaction surface between these two domains that supports a WWOX conformation compatible with peptide substrate binding. Isothermal calorimetry and NMR measurements also indicated that while binding affinity to a single PPxY motif is marginally increased in the presence of WW2, affinity to a dual-motif peptide increases 10-fold. Furthermore, we found WW2 can directly bind double-motif peptides using its canonical binding site. Finally, differential binding of peptides in mutagenesis experiments was consistent with a parallel N- to C-terminal PPxY tandem motif orientation in binding to the WW1–WW2 tandem domain, validating structural models of the interaction. Taken together, our results reveal the complex nature of tandem WW-domain organization and substrate binding, highlighting the contribution of WWOX WW2 to both protein stability and target binding.     

Integration of New Genes into Cellular Networks,and Their Structural Maturation

It has been recently discovered that new genes can originate de novo from noncoding DNA, and several biological traits including expression or sequence composition form a continuum from noncoding sequences to conserved genes. In this article, using yeast genes I test whether the integration of new genes into cellular networks and their structural maturation shows such a continuum by analyzing their changes with gene age. I show that 1) The number of regulatory, protein–protein, and genetic interactions increases continuously with gene age, although with very different rates. New regulatory interactions emerge rapidly within a few million years, while the number of protein–protein and genetic interactions increases slowly, with a rate of 2–2.25 × 10⁻⁸/year and 4.8 × 10⁻⁸/year, respectively. 2) Gene essentiality evolves relatively quickly: the youngest essential genes appear in proto-genes ∼14 MY old. 3) In contrast to interactions, the secondary structure of proteins and their robustness to mutations indicate that new genes face a bottleneck in their evolution: proto-genes are characterized by high β-strand content, high aggregation propensity, and low robustness against mutations, while conserved genes are characterized by lower strand content and higher stability, most likely due to the higher probability of gene loss among young genes and accumulation of neutral mutations.     

Targeting the cyclin-binding groove site to inhibit the catalytic activity of CDK2/cyclin A complex using p27KIP1-derived peptidomimetic inhibitors

Functionally activated cyclin-dependent kinase 2 (CDK2)/cyclin A complex has been validated as an interesting therapeutic target to develop the efficient antineoplastic drug based on the cell cycle arrest. Cyclin A binds to CDK2 and activates the kinases as well as recruits the substrate and inhibitors using a hydrophobic cyclin-binding groove (CBG). Blocking the cyclin substrate recruitment on CBG is an alternative approach to override the specificity hurdle of the currently available ATP site targeting CDK2 inhibitors. Greater understanding of the interaction of CDK2/cyclin A complex with p27 (negative regulator) reveals that the Leu-Phe-Gly (LFG) motif region of p27 binds with the CBG site of cyclin A to arrest the malignant cell proliferation that induces apoptosis. In the present study, Replacement with Partial Ligand Alternatives through Computational Enrichment (REPLACE) drug design strategies have been applied to acquire LFG peptide-derived peptidomimetics library. The peptidomimetics function is equivalent with respect to substrate p27 protein fashion but does not act as an ATP antagonist. The combined approach of molecular docking, molecular dynamics (MD), and molecular electrostatic potential and ADME/T prediction were carried out to evaluate the peptidomimetics. Resultant interaction and electrostatic potential maps suggested that smaller substituent is desirable at the position of phenyl ring to interact with Trp217, Arg250, and Gln254 residues in the active site. The best docked poses were refined by the MD simulations which resulted in conformational changes. After equilibration, the structure of the peptidomimetic and receptor complex was stable. The results revealed that the various substrate protein-derived peptidomimetics could serve as perfect leads against CDK2 protein.

Electronic supplementary material

The online version of this article (doi:10.1007/s12154-014-0124-y) contains supplementary material, which is available to authorized users.     

Seeing beyond the limit: A guide to choosing the right super-resolution microscopy technique

Super-resolution microscopy has become an increasingly popular and robust tool across the life sciences to study minute cellular structures and processes. However, with the increasing number of available super-resolution techniques has come an increased complexity and burden of choice in planning imaging experiments. Choosing the right super-resolution technique to answer a given biological question is vital for understanding and interpreting biological relevance. This is an often-neglected and complex task that should take into account well-defined criteria (e.g., sample type, structure size, imaging requirements). Trade-offs in different imaging capabilities are inevitable; thus, many researchers still find it challenging to select the most suitable technique that will best answer their biological question. This review aims to provide an overview and clarify the concepts underlying the most commonly available super-resolution techniques as well as guide researchers through all aspects that should be considered before opting for a given technique.     

Cellular retinoid-binding proteins transfer retinoids to human cytochrome P450 27C1 for desaturation

Cytochrome P450 27C1 (P450 27C1) is a retinoid desaturase expressed in the skin that catalyzes the formation of 3,4-dehydroretinoids from all-trans retinoids. Within the skin, retinoids are important regulators of proliferation and differentiation. In vivo, retinoids are bound to cellular retinol-binding proteins (CRBPs) and cellular retinoic acid–binding proteins (CRABPs). Interaction with these binding proteins is a defining characteristic of physiologically relevant enzymes in retinoid metabolism. Previous studies that characterized the catalytic activity of human P450 27C1 utilized a reconstituted in vitro system with free retinoids. However, it was unknown whether P450 27C1 could directly interact with holo-retinoid-binding proteins to receive all-trans retinoid substrates. To assess this, steady-state kinetic assays were conducted with free all-trans retinoids and holo-CRBP-1, holo-CRABP-1, and holo-CRABP-2. For holo-CRBP-1 and holo-CRABP-2, the k_cat/K_m values either decreased 5-fold or were equal to the respective free retinoid values. The k_cat/K_m value for holo-CRABP-1, however, decreased ∼65-fold in comparison with reactions with free all-trans retinoic acid. These results suggest that P450 27C1 directly accepts all-trans retinol and retinaldehyde from CRBP-1 and all-trans retinoic acid from CRABP-2, but not from CRABP-1. A difference in substrate channeling between CRABP-1 and CRABP-2 was also supported by isotope dilution experiments. Analysis of retinoid transfer from holo-CRABPs to P450 27C1 suggests that the decrease in k_cat observed in steady-state kinetic assays is due to retinoid transfer becoming rate-limiting in the P450 27C1 catalytic cycle. Overall, these results illustrate that, like the CYP26 enzymes involved in retinoic acid metabolism, P450 27C1 interacts with cellular retinoid-binding proteins.