POPSCOMP: an automated interaction analysis of biomolecular complexes

TNIK: A redox sensor in endothelial cell permeability

Dysregulation of endothelial barrier integrity can lead to vascular leak and potentially fatal oedema. TNF-α controls endothelial permeability during inflammation and requires the actin organizing Ezrin-Radixin-Moesin (ERM) proteins. We identified TRAF2 and NCK-interacting kinase (TNIK) as a kinase directly phosphorylating and activating ERM, specifically at the plasma membrane of primary human endothelial cells. TNIK mediates TNF-α-dependent cellular stiffness and paracellular gap formation in vitro and is essential in driving inflammatory oedema formation in vivo. Unlike its homologs, TNIK activity is negatively and reversibly regulated by H2O2-mediated oxidation of C202 within the kinase domain. TNIK oxidation results in intermolecular disulfide bond formation and loss of kinase activity. Pharmacologic inhibition of endogenous reactive oxygen species production in endothelial cells elevated TNIK-dependent ERM phosphorylation, endothelial cell contraction, and cell rounding. Together, we highlight an interplay between TNIK, ERM phosphorylation, and redox signalling in regulating TNF-induced endothelial cell permeability.

VCAb: a web-tool for structure-guided exploration of antibodies

Motivation

Effective responses against immune challenges require antibodies of different isotypes performing specific effector functions. Structural information on these isotypes is essential to engineer antibodies with desired physico-chemical features of their antigen-binding properties, and optimal developability as potential therapeutics. In silico mutational scanning profiles on antibody structures would further pinpoint candidate mutations for enhancing antibody stability and function. Current antibody structure databases lack consistent annotations of isotypes and structural coverage of 3D antibody structures, as well as computed deep mutation profiles.

Results

The V and C region bearing antibody (VCAb) web-tool is established to clarify these annotations and provides an accessible resource to facilitate antibody engineering and design. VCAb currently provides data on 7,166 experimentally determined antibody structures including both V and C regions from different species. Additionally, VCAb provides annotations of species and isotypes with numbering schemes applied. These information can be interactively queried or downloaded in batch.

Availability and implementation

VCAb is implemented as a R shiny application to enable interactive data interrogation. The online application is freely accessible https://fraternalilab.cs.ucl.ac.uk/VCAb/. The source code to generate the database and the online application is available open-source at https://github.com/Fraternalilab/VCAb.

Saturation genome editing of DDX3X clarifies pathogenicity of germline and somatic variation

Loss-of-function of DDX3X is a leading cause of neurodevelopmental disorders (NDD) in females. DDX3X is also a somatically mutated cancer driver gene proposed to have tumour promoting and suppressing effects. We perform saturation genome editing of DDX3X, testing in vitro the functional impact of 12,776 nucleotide variants. We identify 3432 functionally abnormal variants, in three distinct classes. We train a machine learning classifier to identify functionally abnormal variants of NDD-relevance. This classifier has at least 97% sensitivity and 99% specificity to detect variants pathogenic for NDD, substantially out-performing in silico predictors, and resolving up to 93% of variants of uncertain significance. Moreover, functionally-abnormal variants can account for almost all of the excess nonsynonymous DDX3X somatic mutations seen in DDX3X-driven cancers. Systematic maps of variant effects generated in experimentally tractable cell types have the potential to transform clinical interpretation of both germline and somatic disease-associated variation.

Structural basis of dimerization of chemokine receptors CCR5 and CXCR4

G protein-coupled receptors (GPCRs) are prominent drug targets responsible for extracellular-to-intracellular signal transduction. GPCRs can form functional dimers that have been poorly characterized so far. Here, we show the dimerization mechanism of the chemokine receptors CCR5 and CXCR4 by means of an advanced free-energy technique named coarse-grained metadynamics. Our results reproduce binding events between the GPCRs occurring in the minute timescale, revealing a symmetric and an asymmetric dimeric structure for each of the three investigated systems, CCR5/CCR5, CXCR4/CXCR4, and CCR5/CXCR4. The transmembrane helices TM4-TM5 and TM6-TM7 are the preferred binding interfaces for CCR5 and CXCR4, respectively. The identified dimeric states differ in the access to the binding sites of the ligand and G protein, indicating that dimerization may represent a fine allosteric mechanism to regulate receptor activity. Our study offers structural basis for the design of ligands able to modulate the formation of CCR5 and CXCR4 dimers and in turn their activity, with therapeutic potential against HIV, cancer, and immune-inflammatory diseases.

Systematic detection of tertiary structural modules in large RNAs and RNP interfaces by Tb-seq

Compact RNA structural motifs control many aspects of gene expression, but we lack methods for finding these structures in the vast expanse of multi-kilobase RNAs. To adopt specific 3-D shapes, many RNA modules must compress their RNA backbones together, bringing negatively charged phosphates into close proximity. This is often accomplished by recruiting multivalent cations (usually Mg2+), which stabilize these sites and neutralize regions of local negative charge. Coordinated lanthanide ions, such as terbium (III) (Tb3+), can also be recruited to these sites, where they induce efficient RNA cleavage, thereby revealing compact RNA 3-D modules. Until now, Tb3+ cleavage sites were monitored via low-throughput biochemical methods only applicable to small RNAs. Here we present Tb-seq, a high-throughput sequencing method for detecting compact tertiary structures in large RNAs. Tb-seq detects sharp backbone turns found in RNA tertiary structures and RNP interfaces, providing a way to scan transcriptomes for stable structural modules and potential riboregulatory motifs.

PDBe and PDBe-KB: Providing high-quality, up-to-date and integrated resources of macromolecular structures to support basic and applied research and education

The archiving and dissemination of protein and nucleic acid structures as well as their structural, functional and biophysical annotations is an essential task that enables the broader scientific community to conduct impactful research in multiple fields of the life sciences. The Protein Data Bank in Europe (PDBe; pdbe.org) team develops and maintains several databases and web services to address this fundamental need. From data archiving as a member of the Worldwide PDB consortium (wwPDB; wwpdb.org), to the PDBe Knowledge Base (PDBe-KB; pdbekb.org), we provide data, data-access mechanisms, and visualizations that facilitate basic and applied research and education across the life sciences. Here, we provide an overview of the structural data and annotations that we integrate and make freely available. We describe the web services and data visualization tools we offer, and provide information on how to effectively use or even further develop them. Finally, we discuss the direction of our data services, and how we aim to tackle new challenges that arise from the recent, unprecedented advances in the field of structure determination and protein structure modeling.

Large protein complex interfaces have evolved to promote cotranslational assembly

Assembly pathways of protein complexes should be precise and efficient to minimise misfolding and unwanted interactions with other proteins in the cell. One way to achieve this efficiency is by seeding assembly pathways during translation via the cotranslational assembly of subunits. While recent evidence suggests that such cotranslational assembly is widespread, little is known about the properties of protein complexes associated with the phenomenon. Here, using a combination of proteome-specific protein complex structures and publicly available ribosome profiling data, we show that cotranslational assembly is particularly common between subunits that form large intermolecular interfaces. To test whether large interfaces have evolved to promote cotranslational assembly, as opposed to cotranslational assembly being a non-adaptive consequence of large interfaces, we compared the sizes of first and last translated interfaces of heteromeric subunits in bacterial, yeast, and human complexes. When considering all together, we observe the N-terminal interface to be larger than the C-terminal interface 54% of the time, increasing to 64% when we exclude subunits with only small interfaces, which are unlikely to cotranslationally assemble. This strongly suggests that large interfaces have evolved as a means to maximise the chance of successful cotranslational subunit binding.

Energy Landscapes for Base-Flipping in a Model DNA Duplex

We explore the process of base-flipping for four central bases, adenine, guanine, cytosine, and thymine, in a deoxyribonucleic acid (DNA) duplex using the energy landscape perspective. NMR imino-proton exchange and fluorescence correlation spectroscopy studies have been used in previous experiments to obtain lifetimes for bases in paired and extrahelical states. However, the difference of almost 4 orders of magnitude in the base-flipping rates obtained by the two methods implies that they are exploring different pathways and possibly different open states. Our results support the previous suggestion that minor groove opening may be favored by distortions in the DNA backbone and reveal links between sequence effects and the direction of opening, i.e., whether the base flips toward the major or the minor groove side. In particular, base flipping along the minor groove pathway was found to align toward the 5' side of the backbone. We find that bases align toward the 3' side of the backbone when flipping along the major groove pathway. However, in some cases for cytosine and thymine, the base flipping along the major groove pathway also aligns toward the 5' side. The sequence effect may be caused by the polar interactions between the flipping-base and its neighboring bases on either of the strands. For guanine flipping toward the minor groove side, we find that the equilibrium constant for opening is large compared to flipping via the major groove. We find that the estimated rates of base opening, and hence the lifetimes of the closed state, obtained for thymine flipping through small and large angles along the major groove differ by 6 orders of magnitude, whereas for thymine flipping through small angles along the minor groove and large angles along the major groove, the rates differ by 3 orders of magnitude.

PDBe-KB: collaboratively defining the biological context of structural data

The Protein Data Bank in Europe - Knowledge Base (PDBe-KB, https://pdbe-kb.org) is an open collaboration between world-leading specialist data resources contributing functional and biophysical annotations derived from or relevant to the Protein Data Bank (PDB). The goal of PDBe-KB is to place macromolecular structure data in their biological context by developing standardised data exchange formats and integrating functional annotations from the contributing partner resources into a knowledge graph that can provide valuable biological insights. Since we described PDBe-KB in 2019, there have been significant improvements in the variety of available annotation data sets and user functionality. Here, we provide an overview of the consortium, highlighting the addition of annotations such as predicted covalent binders, phosphorylation sites, effects of mutations on the protein structure and energetic local frustration. In addition, we describe a library of reusable web-based visualisation components and introduce new features such as a bulk download data service and a novel superposition service that generates clusters of superposed protein chains weekly for the whole PDB archive.

Pathogenic missense protein variants affect different functional pathways and proteomic features than healthy population variants

Missense variants are present amongst the healthy population, but some of them are causative of human diseases. A classification of variants associated with “healthy” or “diseased” states is therefore not always straightforward. A deeper understanding of the nature of missense variants in health and disease, the cellular processes they may affect, and the general molecular principles which underlie these differences is essential to offer mechanistic explanations of the true impact of pathogenic variants. Here, we have formalised a statistical framework which enables robust probabilistic quantification of variant enrichment across full-length proteins, their domains, and 3D structure-defined regions. Using this framework, we validate and extend previously reported trends of variant enrichment in different protein structural regions (surface/core/interface). By examining the association of variant enrichment with available functional pathways and transcriptomic and proteomic (protein half-life, thermal stability, abundance) data, we have mined a rich set of molecular features which distinguish between pathogenic and population variants: Pathogenic variants mainly affect proteins involved in cell proliferation and nucleotide processing and are enriched in more abundant proteins. Additionally, rare population variants display features closer to common than pathogenic variants. We validate the association between these molecular features and variant pathogenicity by comparing against existing in silico variant impact annotations. This study provides molecular details into how different proteins exhibit resilience and/or sensitivity towards missense variants and provides the rationale to prioritise variant-enriched proteins and protein domains for therapeutic targeting and development. The ZoomVar database, which we created for this study, is available at fraternalilab.kcl.ac.uk/ZoomVar. It allows users to programmatically annotate missense variants with protein structural information and to calculate variant enrichment in different protein structural regions.

How do can one improve the classification of genetic variants as harmful or harmless? This study uses a robust statistical analysis to exploit the interplay between protein structure, proteomic measurements and functional pathways to enable better discrimination between missense variants in health and disease.

Understanding the structural details of APOBEC3-DNA interactions using graph-based representations

Human APOBEC3 (A3; apolipoprotein B mRNA editing catalytic polypeptide-like 3) is a family of seven enzymes involved in generating mutations in nascent reverse transcripts of many retroviruses, as well as the human genome in a range of cancer types. The structural details of the interaction between A3 proteins and DNA molecules are only available for a few family members. Here we use homology modelling techniques to address the difference in structural coverage of human A3 enzymes interacting with different DNA substrates. A3-DNA interfaces are represented as residue networks ("graphs"), based on which features at these interfaces are compared and quantified. We demonstrate that graph-based representations are effective in highlighting structural features of A3-DNA interfaces. By large-scale in silico mutagenesis of the bound DNA chain, we predicted the preference of substrate DNA sequence for multiple A3 domains. These data suggested that computational modelling approaches could contribute in the exploration of the structural basis for sequence specificity in A3 substrate selection, and demonstrated the utility of graph-based approaches in evaluating a large number of structural models generated in silico.

•

APOBEC3(A3)-DNA structures have been resolved with modified deaminase domains.

•

Structural modelling of interaction between wild-type A3 domains and DNA substrates.

•

Graph-based representations reveal structural differences across A3-DNA interfaces.

•

Using in silico mutagenesis we compared substrate preference of multiple A3 domains.

•

Graph-based approaches can efficiently compare a large number of structural models.

Mining the PDB for Tractable Cases Where X-ray Crystallography Combined with Fragment Screens Can Be Used to Systematically Design Protein-Protein Inhibitors: Two Test Cases Illustrated by IL1β-IL1R and p38α-TAB1 Complexes

Nowadays, it is possible to combine X-ray crystallography and fragment screening in a medium throughput fashion to chemically probe the surfaces used by proteins to interact and use the outcome of the screens to systematically design protein-protein inhibitors. To prove it, we first performed a bioinformatics analysis of the Protein Data Bank protein complexes, which revealed over 400 cases where the crystal lattice of the target in the free form is such that large portions of the interacting surfaces are free from lattice contacts and therefore accessible to fragments during soaks. Among the tractable complexes identified, we then performed single fragment crystal screens on two particular interesting cases: the Il1β-ILR and p38α-TAB1 complexes. The result of the screens showed that fragments tend to bind in clusters, highlighting the small-molecule hotspots on the surface of the target protein. In most of the cases, the hotspots overlapped with the binding sites of the interacting proteins.

PDBe-KB: a community-driven resource for structural and functional annotations

The Protein Data Bank in Europe-Knowledge Base (PDBe-KB, https://pdbe-kb.org) is a community-driven, collaborative resource for literature-derived, manually curated and computationally predicted structural and functional annotations of macromolecular structure data, contained in the Protein Data Bank (PDB). The goal of PDBe-KB is two-fold: (i) to increase the visibility and reduce the fragmentation of annotations contributed by specialist data resources, and to make these data more findable, accessible, interoperable and reusable (FAIR) and (ii) to place macromolecular structure data in their biological context, thus facilitating their use by the broader scientific community in fundamental and applied research. Here, we describe the guidelines of this collaborative effort, the current status of contributed data, and the PDBe-KB infrastructure, which includes the data exchange format, the deposition system for added value annotations, the distributable database containing the assembled data, and programmatic access endpoints. We also describe a series of novel web-pages-the PDBe-KB aggregated views of structure data-which combine information on macromolecular structures from many PDB entries. We have recently released the first set of pages in this series, which provide an overview of available structural and functional information for a protein of interest, referenced by a UniProtKB accession.

Robotic assisted generation of 2'-deoxy-2'-fluoro-modifed RNA aptamers - High performance enabling strategies in aptamer selection

Aptamer selection is a laborious procedure, requiring expertise and significant resources. These characteristics limit the accessibility of researchers to these molecular tools. We describe a selection procedure, making use of a robotic system that allows the fully automated selection of RNA and 2'deoxy-2'-fluoro pyrimidine RNA aptamers. The platform offers a rapid access to aptamers for basic research and development, therefore opening the path to aptamer-based systemic analysis of proteomes in biological settings.

Inhibitor-induced HER2-HER3 heterodimerisation promotes proliferation through a novel dimer interface

While targeted therapy against HER2 is an effective first-line treatment in HER2+ breast cancer, acquired resistance remains a clinical challenge. The pseudokinase HER3, heterodimerisation partner of HER2, is widely implicated in the resistance to HER2-mediated therapy. Here, we show that lapatinib, an ATP-competitive inhibitor of HER2, is able to induce proliferation cooperatively with the HER3 ligand neuregulin. This counterintuitive synergy between inhibitor and growth factor depends on their ability to promote atypical HER2-HER3 heterodimerisation. By stabilising a particular HER2 conformer, lapatinib drives HER2-HER3 kinase domain heterocomplex formation. This dimer exists in a head-to-head orientation distinct from the canonical asymmetric active dimer. The associated clustering observed for these dimers predisposes to neuregulin responses, affording a proliferative outcome. Our findings provide mechanistic insights into the liabilities involved in targeting kinases with ATP-competitive inhibitors and highlight the complex role of protein conformation in acquired resistance.

Conserved asymmetry underpins homodimerization of Dicer-associated double-stranded RNA-binding proteins

Double-stranded RNA-binding domains (dsRBDs) are commonly found in modular proteins that interact with RNA. Two varieties of dsRBD exist: canonical Type A dsRBDs interact with dsRNA, while non-canonical Type B dsRBDs lack RNA-binding residues and instead interact with other proteins. In higher eukaryotes, the microRNA biogenesis enzyme Dicer forms a 1:1 association with a dsRNA-binding protein (dsRBP). Human Dicer associates with HIV TAR RNA-binding protein (TRBP) or protein activator of PKR (PACT), while Drosophila Dicer-1 associates with Loquacious (Loqs). In each case, the interaction involves a region of the protein that contains a Type B dsRBD. All three dsRBPs are reported to homodimerize, with the Dicer-binding region implicated in self-association. We report that these dsRBD homodimers display structural asymmetry and that this unusual self-association mechanism is conserved from flies to humans. We show that the core dsRBD is sufficient for homodimerization and that mutation of a conserved leucine residue abolishes self-association. We attribute differences in the self-association properties of Loqs, TRBP and PACT to divergence of the composition of the homodimerization interface. Modifications that make TRBP more like PACT enhance self-association. These data are examined in the context of miRNA biogenesis and the protein/protein interaction properties of Type B dsRBDs.

Investigating Ebola virus pathogenicity using molecular dynamics

BACKGROUND

Ebolaviruses have been known to cause deadly disease in humans for 40 years and have recently been demonstrated in West Africa to be able to cause large outbreaks. Four Ebolavirus species cause severe disease associated with high mortality in humans. Reston viruses are the only Ebolaviruses that do not cause disease in humans. Conserved amino acid changes in the Reston virus protein VP24 compared to VP24 of other Ebolaviruses have been suggested to alter VP24 binding to host cell karyopherins resulting in impaired inhibition of interferon signalling, which may explain the difference in human pathogenicity. Here we used protein structural analysis and molecular dynamics to further elucidate the interaction between VP24 and KPNA5.

RESULTS

As a control experiment, we compared the interaction of wild-type and R137A-mutant (known to affect KPNA5 binding) Ebola virus VP24 with KPNA5. Results confirmed that the R137A mutation weakens direct VP24-KPNA5 binding and enables water molecules to penetrate at the interface. Similarly, Reston virus VP24 displayed a weaker interaction with KPNA5 than Ebola virus VP24, which is likely to reduce the ability of Reston virus VP24 to prevent host cell interferon signalling.

CONCLUSION

Our results provide novel molecular detail on the interaction of Reston virus VP24 and Ebola virus VP24 with human KPNA5. The results indicate a weaker interaction of Reston virus VP24 with KPNA5 than Ebola virus VP24, which is probably associated with a decreased ability to interfere with the host cell interferon response. Hence, our study provides further evidence that VP24 is a key player in determining Ebolavirus pathogenicity.

Structural modeling of the AhR:ARNT complex in the bHLH-PASA-PASB region elucidates the key determinants of dimerization

Elucidation of the dimerization process of the aryl hydrocarbon receptor (AhR) with the AhR nuclear translocator (ARNT) is crucial for understanding the mechanisms underlying the functional activity of AhR, including mediation of the toxicity of environmental contaminants. In this work, for the first time a structural model of the AhR:ARNT dimer encompassing the entire bHLH-PASA-PASB domain region is proposed. It is developed by using a template-based modeling approach, relying on the recently available crystallographic structures of two dimers of homologous systems in the bHLH-PAS family of proteins: the CLOCK:BMAL1 and the HIF2α:ARNT heterodimers. The structural and energetic characteristics of the modeled AhR:ARNT protein-protein interface are determined by evaluating the variations in solvent accessible surface area, the total binding free energy and the per-residue free energy contributions obtained by the MM-GBSA method and the Energy Decomposition Analysis. The analyses of the intricate network of inter-domain interactions at the dimerization interfaces provide insights into the key determinants of dimerization. These are confirmed by comparison of the computational findings with the available experimental mutagenesis and functional analysis data. The results presented here on the AhR:ARNT dimer structure and interactions provide a framework to start analyzing the mechanism of AhR transformation into its functional DNA binding form.

On the role of residue phosphorylation in 14-3-3 partners: AANAT as a case study

Twenty years ago, a novel concept in protein structural biology was discovered: the intrinsically disordered regions (IDRs). These regions remain largely unstructured under native conditions and the more are studied, more properties are attributed to them. Possibly, one of the most important is their ability to conform a new type of protein-protein interaction. Besides the classical domain-to-domain interactions, IDRs follow a 'fly-casting' model including 'induced folding'. Unfortunately, it is only possible to experimentally explore initial and final states. However, the complete movie of conformational changes of protein regions and their characterization can be addressed by in silico experiments. Here, we simulate the binding of two proteins to describe how the phosphorylation of a single residue modulates the entire process. 14-3-3 protein family is considered a master regulator of phosphorylated proteins and from a modern point-of-view, protein phosphorylation is a three component system, with writers (kinases), erasers (phosphatases) and readers. This later biological role is attributed to the 14-3-3 protein family. Our molecular dynamics results show that phosphorylation of the key residue Thr31 in a partner of 14-3-3, the aralkylamine N-acetyltransferase, releases the fly-casting mechanism during binding. On the other hand, the non-phosphorylation of the same residue traps the proteins, systematically and repeatedly driving the simulations into wrong protein-protein conformations.

Radical-induced purine lesion formation is dependent on DNA helical topology

Herein we report the quantification of purine lesions arising from gamma-radiation sourced hydroxyl radicals (HO^•) on tertiary dsDNA helical forms of supercoiled (SC), open circular (OC), and linear (L) conformation, along with single-stranded folded and non-folded sequences of guanine-rich DNA in selected G-quadruplex structures. We identify that DNA helical topology and folding plays major, and unexpected, roles in the formation of 8-oxo-7,8-dihydro-2'-deoxyguanosine (8-oxo-dG) and 8-oxo-7,8-dihydro-2'-deoxyadenosine (8-oxo-dA), along with tandem-type purine lesions 5',8-cyclo-2'-deoxyguanosine (5',8-cdG) and 5',8-cyclo-2'-deoxyadenosine (5',8-cdA). SC, OC, and L dsDNA conformers together with folded and non-folded G-quadruplexes d[TGGGGT]₄ (TG4T), d[AGGG(TTAGGG)₃] (Tel22), and the mutated tel24 d[TTGGG(TTAGGG)₃A] (mutTel24) were exposed to HO^• radicals and purine lesions were then quantified via stable isotope dilution LC-MS/MS analysis. Purine oxidation in dsDNA follows L > OC ≫ SC indicating greater damage towards the extended B-DNA topology. Conversely, G-quadruplex sequences were significantly more resistant toward purine oxidation in their unfolded states as compared with G-tetrad folded topologies; this effect is confirmed upon comparative analysis of Tel22 (∼50% solution folded) and mutTel24 (∼90% solution folded). In an effort to identify the accessibly of hydroxyl radicals to quadruplex purine nucleobases, G-quadruplex solvent cavities were then modeled at 1.33 Å with evidence suggesting that folded G-tetrads may act as potential oxidant traps to protect against chromosomal DNA damage.

Homology-Based Modeling of Universal Stress Protein from Listeria innocua Up-Regulated under Acid Stress Conditions

An Universal Stress Protein (USP) expressed under acid stress condition by Listeria innocua ATCC 33090 was investigated. The USP was up-regulated not only in the stationary phase but also during the exponential growth phase. The three dimensional (3D) structure of USP was predicted using a combined proteomic and bioinformatics approach. Phylogenetic analysis showed that the USP from Listeria detected in our study was distant from the USPs of other bacteria (such as Pseudomonas spp., Escherichia coli, Salmonella spp.) and clustered in a separate and heterogeneous class including several USPs from Listeria spp. and Lactobacillus spp. An important information on the studied USP was obtained from the 3D-structure established through the homology modeling procedure. In detail, the Model_USP-691 suggested that the investigated USP had a homo-tetrameric quaternary structure. Each monomer presented an architecture analogous to the Rossmann-like α/β-fold with five parallel β-strands, and four α-helices. The analysis of monomer-monomer interfaces and quality of the structure alignments confirmed the model reliability. In fact, the structurally and sequentially conserved hydrophobic residues of the β-strand 5 (in particular the residues V¹⁴⁶ and V¹⁴⁸) were involved in the inter-chains contact. Moreover, the highly conserved residues I¹³⁹ and H¹⁴¹ in the region α4 were involved in the dimer association and functioned as hot spots into monomer–monomer interface assembly. The hypothetical assembly of dimers was also supported by the large interface area and by the negative value of solvation free energy gain upon interface interaction. Finally, the structurally conserved ATP-binding motif G-2X-G-9X-G(S/T-N) suggested for a putative role of ATP in stabilizing the tetrameric assembly of the USP. Therefore, the results obtained from a multiple approach, consisting in the application of kinetic, proteomic, phylogenetic and modeling analyses, suggest that Listeria USP could be considered a new type of ATP-binding USP involved in the response to acid stress condition during the exponential growth phase.

Evidence that the Rhodopsin Kinase (GRK1) N-Terminus and the Transducin Gα C-Terminus Interact with the Same "Hydrophobic Patch" on Rhodopsin TM5

Phosphorylation of G protein-coupled receptors (GPCRs) terminates their ability to couple with and activate G proteins by increasing their affinity for arrestins. Unfortunately, detailed information regarding how GPCRs interact with the kinases responsible for their phosphorylation is still limited. Here, we purified fully functional GPCR kinase 1 (GRK1) using a rapid method and used it to gain insights into how this important kinase interacts with the GPCR rhodopsin. Specifically, we find that GRK1 uses the same site on rhodopsin as the transducin (Gt) Gtα C-terminal tail and the arrestin "finger loop", a cleft formed in the cytoplasmic face of the receptor upon activation. Our studies also show GRK1 requires two conserved residues located in this cleft (L226 and V230) that have been shown to be required for Gt activation due to their direct interactions with hydrophobic residues on the Gα C-terminal tail. Our data and modeling studies are consistent with the idea that all three proteins (Gt, GRK1, and visual arrestin) bind, at least in part, in the same site on rhodopsin and interact with the receptor through a similar hydrophobic contact-driven mechanism.

PinSnps: structural and functional analysis of SNPs in the context of protein interaction networks

Summary: We present a practical computational pipeline to readily perform data analyses of protein–protein interaction networks by using genetic and functional information mapped onto protein structures. We provide a 3D representation of the available protein structure and its regions (surface, interface, core and disordered) for the selected genetic variants and/or SNPs, and a prediction of the mutants’ impact on the protein as measured by a range of methods. We have mapped in total 2587 genetic disorder-related SNPs from OMIM, 587 873 cancer-related variants from COSMIC, and 1 484 045 SNPs from dbSNP. All result data can be downloaded by the user together with an R-script to compute the enrichment of SNPs/variants in selected structural regions.

Availability and Implementation: PinSnps is available as open-access service at http://fraternalilab.kcl.ac.uk/PinSnps/

Contact:franca.fraternali@kcl.ac.uk

Supplementary information:Supplementary data are available at Bioinformatics online.

Assembly of Influenza Hemagglutinin Fusion Peptides in a Phospholipid Bilayer by Coarse-grained Computer Simulations

Membrane fusion is critical to eukaryotic cellular function and crucial to the entry of enveloped viruses such as influenza and human immunodeficiency virus. Influenza viral entry in the host cell is mediated by a 20–23 amino acid long sequence, called the fusion peptide (FP). Recently, possible structures for the fusion peptide (ranging from an inverted V shaped α-helical structure to an α-helical hairpin, or to a complete α-helix) and their implication in the membrane fusion initiation have been proposed. Despite the large number of studies devoted to the structure of the FP, the mechanism of action of this peptide remains unclear with several mechanisms having been suggested, including the induction of local disorder, promoting membrane curvature, and/or altering local membrane composition. In recent years, several research groups have employed atomistic and/or coarse-grained molecular dynamics (MD) simulations to investigate the matter. In all previous works, the behavior of a single FP monomer was studied, while in this manuscript, we use a simplified model of a tripeptide (TP) monomer of FP (TFP) instead of a single FP monomer because each Influenza Hemagglutinin contains three FP molecules in the biological system. In this manuscript we report findings targeted at understanding the fusogenic properties and the collective behavior of these trimers of FP peptides on a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine model membrane. Here we show how the TFP monomers self-assemble into differently sized oligomers in the presence of the membrane. We measure the perturbation to the structure of the phospholipid membrane caused by the presence of these TFP oligomers. Our work (i) shows how self-assembly of TFP in the presence of the membrane induces non negligible deformation to the membrane and (ii) could be a useful starting point to stimulate discussion and further work targeted to fusion pore formation.

Protein-protein interaction networks studies and importance of 3D structure knowledge

Protein-protein interaction networks (PPINs) are a powerful tool to study biological processes in living cells. In this review, we present the progress of PPIN studies from abstract to more detailed representations. We will focus on 3D interactome networks, which offer detailed information at the atomic level. This information can be exploited in understanding not only the underlying cellular mechanisms, but also how human variants and disease-causing mutations affect protein functions and complexes' stability. Recent studies have used structural information on PPINs to also understand the molecular mechanisms of binding partner selection. We will address the challenges in generating 3D PPINs due to the restricted number of solved protein structures. Finally, some of the current use of 3D PPINs will be discussed, highlighting their contribution to the studies in genotype-phenotype relationships and in the optimization of targeted studies to design novel chemical compounds for medical treatments.

Design and application of implicit solvent models in biomolecular simulations

•

Implicit solvent replaces explicit water by a potential of mean force.

•

Popular models are SASA, VOL and Generalized Born.

•

Implicit solvent is used in MD, protein modelling, folding, design, prediction and drug screening.

•

Large-scale simulations allow for parametrisation via force matching.

•

Application to nucleic acids and membranes is challenging.

We review implicit solvent models and their parametrisation by introducing the concepts and recent devlopments of the most popular models with a focus on parametrisation via force matching. An overview of recent applications of the solvation energy term in protein dynamics, modelling, design and prediction is given to illustrate the usability and versatility of implicit solvation in reproducing the physical behaviour of biomolecular systems. Limitations of implicit modes are discussed through the example of more challenging systems like nucleic acids and membranes.

Structure and stability insights into tumour suppressor p53 evolutionary related proteins

The p53 family of genes and their protein products, namely, p53, p63 and p73, have over one billion years of evolutionary history. Advances in computational biology and genomics are enabling studies of the complexities of the molecular evolution of p53 protein family to decipher the underpinnings of key biological conditions spanning from cancer through to various metabolic and developmental disorders and facilitate the design of personalised medicines. However, a complete understanding of the inherent nature of the thermodynamic and structural stability of the p53 protein family is still lacking. This is due, to a degree, to the lack of comprehensive structural information for a large number of homologous proteins and to an incomplete knowledge of the intrinsic factors responsible for their stability and how these might influence function. Here we investigate the thermal stability, secondary structure and folding properties of the DNA-binding domains (DBDs) of a range of proteins from the p53 family using biophysical methods. While the N- and the C-terminal domains of the p53 family show sequence diversity and are normally targets for post-translational modifications and alternative splicing, the central DBD is highly conserved. Together with data obtained from Molecular Dynamics simulations in solution and with structure based homology modelling, our results provide further insights into the molecular properties of evolutionary related p53 proteins. We identify some marked structural differences within the p53 family, which could account for the divergence in biological functions as well as the subtleties manifested in the oligomerization properties of this family.

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.

Replica exchange improves sampling in low-resolution docking stage of RosettaDock

Many protein-protein docking protocols are based on a shotgun approach, in which thousands of independent random-start trajectories minimize the rigid-body degrees of freedom. Another strategy is enumerative sampling as used in ZDOCK. Here, we introduce an alternative strategy, ReplicaDock, using a small number of long trajectories of temperature replica exchange. We compare replica exchange sampling as low-resolution stage of RosettaDock with RosettaDock's original shotgun sampling as well as with ZDOCK. A benchmark of 30 complexes starting from structures of the unbound binding partners shows improved performance for ReplicaDock and ZDOCK when compared to shotgun sampling at equal or less computational expense. ReplicaDock and ZDOCK consistently reach lower energies and generate significantly more near-native conformations than shotgun sampling. Accordingly, they both improve typical metrics of prediction quality of complex structures after refinement. Additionally, the refined ReplicaDock ensembles reach significantly lower interface energies and many previously hidden features of the docking energy landscape become visible when ReplicaDock is applied.

CONS-COCOMAPS: a novel tool to measure and visualize the conservation of inter-residue contacts in multiple docking solutions

Background

The development of accurate protein-protein docking programs is making this kind of simulations an effective tool to predict the 3D structure and the surface of interaction between the molecular partners in macromolecular complexes. However, correctly scoring multiple docking solutions is still an open problem. As a consequence, the accurate and tedious screening of many docking models is usually required in the analysis step.

Methods

All the programs under CONS-COCOMAPS have been written in python, taking advantage of python libraries such as SciPy and Matplotlib. CONS-COCOMAPS is freely available as a web tool at the URL:

http://www.molnac.unisa.it/BioTools/conscocomaps/.

Results

Here we presented CONS-COCOMAPS, a novel tool to easily measure and visualize the consensus in multiple docking solutions. CONS-COCOMAPS uses the conservation of inter-residue contacts as an estimate of the similarity between different docking solutions. To visualize the conservation, CONS-COCOMAPS uses intermolecular contact maps.

Conclusions

The application of CONS-COCOMAPS to test-cases taken from recent CAPRI rounds has shown that it is very efficient in highlighting even a very weak consensus that often is biologically meaningful.

Protein-water interactions in MD simulations: POPS/POPSCOMP solvent accessibility analysis, solvation forces and hydration sites

The effects of solvation on molecular recognition are investigated from different perspectives, ranging from methods to analyse explicit solvent dynamical behaviour at the protein surface to methods for the implicit treatment of solvent effects associated with the conformational behaviour of biomolecules. The here presented implicit solvation method is based on an analytical approximation of the Solvent Accessible Surface Area (SASA) of solute molecules, which is computationally efficient and easy to parametrise. The parametrised SASA solvation method is discussed in the light of protein design and ligand binding studies. The POPS program for the SASA computation on single molecules and complex interfaces is described in detail. Explicit solvent behaviour is described here in the form of solvent density maps at the protein surface. We highlight the usefulness of that approach in defining the organisation of specific water molecules at functional sites and in determining hydrophobicity scores for the identification of potential interaction patches.

COCOMAPS: a web application to analyze and visualize contacts at the interface of biomolecular complexes

SUMMARY

Herein we present COCOMAPS, a novel tool for analyzing, visualizing and comparing the interface in protein-protein and protein-nucleic acids complexes. COCOMAPS combines traditional analyses and 3D visualization of the interface with the effectiveness of intermolecular contact maps.

AVAILABILITY

COCOMAPS is accessible as a public web tool at http://www.molnac.unisa.it/BioTools/cocomaps

CONTACT

romina.oliva@uniparthenope.it; lcavallo@unisa.it.

Comparing interfacial dynamics in protein-protein complexes: an elastic network approach

Background

The transient, or permanent, association of proteins to form organized complexes is one of the most common mechanisms of regulation of biological processes. Systematic physico-chemical studies of the binding interfaces have previously shown that a key mechanism for the formation/stabilization of dimers is the steric and chemical complementarity of the two semi-interfaces. The role of the fluctuation dynamics at the interface of the interacting subunits, although expectedly important, proved more elusive to characterize. The aim of the present computational study is to gain insight into salient dynamics-based aspects of protein-protein interfaces.

Results

The interface dynamics was characterized by means of an elastic network model for 22 representative dimers covering three main interface types. The three groups gather dimers sharing the same interface but with good (type I) or poor (type II) similarity of the overall fold, or dimers sharing only one of the semi-interfaces (type III). The set comprises obligate dimers, which are complexes for which no structural representative of the free form(s) is available. Considerations were accordingly limited to bound and unbound forms of the monomeric subunits of the dimers. We proceeded by first computing the mobility of amino acids at the interface of the bound forms and compare it with the mobility of (i) other surface amino acids (ii) interface amino acids in the unbound forms. In both cases different dynamic patterns were observed across interface types and depending on whether the interface belongs to an obligate or non-obligate complex.

Conclusions

The comparative investigation indicated that the mobility of amino acids at the dimeric interface is generally lower than for other amino acids at the protein surface. The change in interfacial mobility upon removing "in silico" the partner monomer (unbound form) was next found to be correlated with the interface type, size and obligate nature of the complex. In particular, going from the unbound to the bound forms, the interfacial mobility is noticeably reduced for dimers with type I interfaces, while it is largely unchanged for type II ones. The results suggest that these structurally- and biologically-different types of interfaces are stabilized by different balancing mechanisms between enthalpy and conformational entropy.

Two modules in the BRC repeats of BRCA2 mediate structural and functional interactions with the RAD51 recombinase

The breast and ovarian cancer suppressor protein BRCA2 controls the RAD51 recombinase in reactions that lead to homologous DNA recombination (HDR). BRCA2 binds RAD51 via eight conserved BRC repeat motifs of approximately 35 amino acids, each with a varying capacity to bind RAD51. BRC repeats both promote and inhibit RAD51 assembly on different DNA substrates to regulate HDR, but the structural basis for these functions is unclear. Here, we demarcate two tetrameric clusters of hydrophobic residues in the BRC repeats, interacting with distinct pockets in RAD51, and show that the co-location of both modules within a single BRC repeat is necessary for BRC–RAD51 binding and function. The two modules comprise the sequence FxxA, known to inhibit RAD51 assembly by blocking the oligomerization interface, and a previously unrecognized tetramer with the consensus sequence LFDE, which binds to a RAD51 pocket distinct from this interface. The LFDE motif is essential in BRC repeats for modes of RAD51 binding both permissive and inhibitory to RAD51 oligomerization. Targeted insertion of point mutations in RAD51 that disrupt the LFDE-binding pocket impair its assembly at DNA damage sites in living cells. Our findings suggest a model for the modular architecture of BRC repeats that provides fresh insight into the mechanisms regulating homologous DNA recombination.

Structural basis of nucleotide exchange and client binding by the Hsp70 cochaperone Bag2

Cochaperones are essential for Hsp70- and Hsc70-mediated folding of proteins and include nucleotide-exchange factors (NEFs) that assist protein folding by accelerating ADP-ATP exchange on Hsp70. The cochaperone Bag2 binds misfolded Hsp70 clients and also acts as an NEF, but the molecular basis for its function is unclear. We show that, rather than being a member of the Bag domain family, Bag2 contains a new type of Hsp70 NEF domain, which we call the 'brand new bag' (BNB) domain. Free and Hsc70-bound crystal structures of Bag2-BNB show its dimeric structure, in which a flanking linker helix and loop bind to Hsc70 to promote nucleotide exchange. NMR analysis demonstrates that the client binding sites and Hsc70-interaction sites of the Bag2-BNB overlap, and that Hsc70 can displace clients from Bag2-BNB, indicating a distinct mechanism for the regulation of Hsp70-mediated protein folding by Bag2.

Crystal structures of the pro-inflammatory cytokine interleukin-23 and its complex with a high-affinity neutralizing antibody

Interleukin (IL)-23 is a pro-inflammatory cytokine playing a key role in the pathogenesis of several autoimmune and inflammatory diseases. We have determined the crystal structures of the heterodimeric p19-p40 IL-23 and its complex with the Fab (antigen-binding fragment) of a neutralizing antibody at 2.9 and 1.9 A, respectively. The IL-23 structure closely resembles that of IL-12. They share the common p40 subunit, and IL-23 p19 overlaps well with IL-12 p35. Along the hydrophilic heterodimeric interface, fewer charged residues are involved for IL-23 compared with IL-12. The binding site of the Fab is located exclusively on the p19 subunit, and comparison with published cytokine-receptor structures suggests that it overlaps with the IL-23 receptor binding site.

Molecular interactions of ASPP1 and ASPP2 with the p53 protein family and the apoptotic promoters PUMA and Bax

The apoptosis stimulating p53 proteins, ASPP1 and ASPP2, are the first two common activators of the p53 protein family that selectively enable the latter to regulate specific apoptotic target genes, which facilitates yes yet unknown mechanisms for discrimination between cell cycle arrest and apoptosis. To better understand the interplay between ASPP- and p53-family of proteins we investigated the molecular interactions between them using biochemical methods and structure-based homology modelling. The data demonstrate that: (i) the binding of ASPP1 and ASPP2 to p53, p63 and p73 is direct; (ii) the C-termini of ASPP1 and ASPP2 interact with the DNA-binding domains of p53 protein family with dissociation constants, K_d, in the lower micro-molar range; (iii) the stoichiometry of binding is 1:1; (iv) the DNA-binding domains of p53 family members are sufficient for these protein–protein interactions; (v) EMSA titrations revealed that while tri-complex formation between ASPPs, p53 family of proteins and PUMA/Bax is mutually exclusive, ASPP2 (but not ASPP1) formed a complex with PUMA (but not Bax) and displaced p53 and p73. The structure-based homology modelling revealed subtle differences between ASPP2 and ASPP1 and together with the experimental data provide novel mechanistic insights.

Structure of a tyrosyl-tRNA synthetase splicing factor bound to a group I intron RNA

The 'RNA world' hypothesis holds that during evolution the structural and enzymatic functions initially served by RNA were assumed by proteins, leading to the latter's domination of biological catalysis. This progression can still be seen in modern biology, where ribozymes, such as the ribosome and RNase P, have evolved into protein-dependent RNA catalysts ('RNPzymes'). Similarly, group I introns use RNA-catalysed splicing reactions, but many function as RNPzymes bound to proteins that stabilize their catalytically active RNA structure. One such protein, the Neurospora crassa mitochondrial tyrosyl-tRNA synthetase (TyrRS; CYT-18), is bifunctional and both aminoacylates mitochondrial tRNA(Tyr) and promotes the splicing of mitochondrial group I introns. Here we determine a 4.5-A co-crystal structure of the Twort orf142-I2 group I intron ribozyme bound to splicing-active, carboxy-terminally truncated CYT-18. The structure shows that the group I intron binds across the two subunits of the homodimeric protein with a newly evolved RNA-binding surface distinct from that which binds tRNA(Tyr). This RNA binding surface provides an extended scaffold for the phosphodiester backbone of the conserved catalytic core of the intron RNA, allowing the protein to promote the splicing of a wide variety of group I introns. The group I intron-binding surface includes three small insertions and additional structural adaptations relative to non-splicing bacterial TyrRSs, indicating a multistep adaptation for splicing function. The co-crystal structure provides insight into how CYT-18 promotes group I intron splicing, how it evolved to have this function, and how proteins could have incrementally replaced RNA structures during the transition from an RNA world to an RNP world.

Conformational stability and DNA binding energetics of the rat thyroid transcription factor 1 homeodomain

The conformational stability of the rat thyroid transcription factor 1 homeodomain, TTF-1HD, has been investigated by means of circular dichroism (CD) and differential scanning calorimetry (DSC) measurements at pH 5.0 as a function of KCl concentration. Thermal unfolding of TTF-1HD is a reversible two-state transition. The protein is not stable against temperature, showing a denaturation temperature of 32 degrees C in the absence of salt and 50 degrees C at 75 mM KCl. The binding energetics of TTF-1HD to its target DNA sequence has been characterized by means of isothermal titration calorimetry (ITC) measurements, complemented with CD data. At 25 degrees C, pH 5.0 and 75 mM KCl, the binding constant amounts to 1.5 x 10(8)M(-1) and the binding enthalpy change amounts to -41 kJ mol(-1). The process is enthalpy driven, but also the entropy change is favorable to complex formation. To gain a molecular level understanding of the interactions determining the association of TTF-1HD to the target DNA sequence structural information would be requested, but it is not yet available. Therefore, structural models of two complexes, TTF-1HD with the target DNA sequence and TTF-1HD with a modified DNA sequence, have been constructed by using as a template the NMR structure of the complex between NK-2 HD and its target DNA, and by performing molecular dynamics simulations 3.5 ns long. Analysis of these models allows one to shed light on the origin of the DNA binding specificity characteristic of TTF-1HD.

Interaction of malaria parasite-inhibitory antibodies with the merozoite surface protein MSP1(19) by computational docking

Merozoite surface protein 1 (MSP1) of the malaria parasite Plasmodium falciparum is an important vaccine candidate antigen. Antibodies specific for the C-terminal maturation product, MSP1(19), have been shown to inhibit erythrocyte invasion and parasite growth. Specific monoclonal antibodies react with conformational epitopes contained within the two EGF-like domains that constitute the antigen MSP1(19). To gain greater insight into the inhibitory process, the authors selected two strongly inhibitory antibodies (designated 12.8 and 12.10) and modeled their structures by homology. Computational docking was used to generate antigen-antibody complexes and a selection filter based on NMR data was applied to obtain plausible models. Molecular Dynamics simulations of the selected complexes were performed to evaluate the role of specific side chains in the binding. Favorable complexes were obtained that complement the NMR data in defining specific binding sites. These models can provide valuable guidelines for future experimental work that is devoted to the understanding of the action mechanism of invasion-inhibitory antibodies.

Probing the structure of the Caulobacter crescentus ribosome with chemical labeling and mass spectrometry

The ribosomal proteins of Caulobacter crescentus were amidinated before and after disassembly of the organelle and the results analyzed by mass spectrometry. Comparison with structural information from previous X-ray crystal studies of other bacterial ribosomes provides insight about the C. crescentus ribosome. In total, 47 of the 54 proteins present in the ribosome of C. crescentus were detected after labeling. The extent of derivatization for each protein is strongly dependent on the solvent accessibility of its target residues. Proteins of the ribosome stalk, which are known to be largely solvent-accessible, were labeled quite extensively. In striking contrast, other proteins that are known to be highly shielded in their subunits were labeled at very few of their potential sites. Furthermore, evidence that protein L12 binds to the ribosome via its N-terminal domain is consistent with previous findings.

Interactions of the C2 domain of human factor V with a model membrane

Activated coagulation Factor V is an important cofactor of the coagulation cascade that catalyzes the formation of the prothrombinase complex on the surface of membranes rich in phosphatidyl-L-serine (PS). Here we report molecular dynamics simulations of the two crystallographic structures (the open and closed conformations) of domain C2 of coagulation Factor V (FaVC2). The calculations were performed in water (1.5 ns for each conformation) and in the presence of a neutral phospholipid bilayer model (POPE; 10 ns for each conformation) in order to describe the dynamics of the free (plasma circulating) and membrane bound forms of FaVC2. Water simulations confirmed the hypothesis that the plasma circulating form is in the closed conformation. In contrast, the membrane simulations showed that both conformations are energetically compatible with membrane binding. We have investigated the mechanism, the dynamics, and the energetics of the binding process. Our data are consistent with published estimates of the immersion depth of the aromatic residues (W26 and W27), and with mutagenesis studies involving specific residues located on the spikes at the bottom of the FaVC2 structure. Electrostatic interactions between the phospholipid head groups and hydrophilic residues at the bottom of the structure play a key role in the binding process by creating a large number of hydrogen bonds that anchor the protein to the membrane. The simulations identified a stable phospholipid binding pocket reminiscent of a previously suggested PS interaction site. Our structural data could contribute to the design of potential inhibitors able to disrupt membrane association.