HADDOCK versus HADDOCK: new features and performance of HADDOCK2.0 on the CAPRI targets

Galaxy7TM: flexible GPCR-ligand docking by structure refinement

G-protein-coupled receptors (GPCRs) play important physiological roles related to signal transduction and form a major group of drug targets. Prediction of GPCR-ligand complex structures has therefore important implications to drug discovery. With previously available servers, it was only possible to first predict GPCR structures by homology modeling and then perform ligand docking on the model structures. However, model structures generated without explicit consideration of specific ligands of interest can be inaccurate because GPCR structures can be affected by ligand binding. The Galaxy7TM server, freely accessible at http://galaxy.seoklab.org/7TM, improves an input GPCR structure by simultaneous ligand docking and flexible structure refinement using GALAXY methods. The server shows better performance in both ligand docking and GPCR structure refinement than commonly used programs AutoDock Vina and Rosetta MPrelax, respectively.

Automated structure modeling of large protein assemblies using crosslinks as distance restraints

Crosslinking mass spectrometry is increasingly used for structural characterization of multisubunit protein complexes. Chemical crosslinking captures conformational heterogeneity, which typically results in conflicting crosslinks that cannot be satisfied in a single model, making detailed modeling a challenging task. Here we introduce an automated modeling method dedicated to large protein assemblies ('XL-MOD' software is available at http://aria.pasteur.fr/supplementary-data/x-links) that (i) uses a form of spatial restraints that realistically reflects the distribution of experimentally observed crosslinked distances; (ii) automatically deals with ambiguous and/or conflicting crosslinks and identifies alternative conformations within a Bayesian framework; and (iii) allows subunit structures to be flexible during conformational sampling. We demonstrate our method by testing it on known structures and available crosslinking data. We also crosslinked and modeled the 17-subunit yeast RNA polymerase III at atomic resolution; the resulting model agrees remarkably well with recently published cryoelectron microscopy structures and provides additional insights into the polymerase structure.

Structural modeling of G-protein coupled receptors: An overview on automatic web-servers

Despite the significant efforts and discoveries during the last few years in G protein-coupled receptor (GPCR) expression and crystallization, the receptors with known structures to date are limited only to a small fraction of human GPCRs. The lack of experimental three-dimensional structures of the receptors represents a strong limitation that hampers a deep understanding of their function. Computational techniques are thus a valid alternative strategy to model three-dimensional structures. Indeed, recent advances in the field, together with extraordinary developments in crystallography, in particular due to its ability to capture GPCRs in different activation states, have led to encouraging results in the generation of accurate models. This, prompted the community of modelers to render their methods publicly available through dedicated databases and web-servers. Here, we present an extensive overview on these services, focusing on their advantages, drawbacks and their role in successful applications. Future challenges in the field of GPCR modeling, such as the predictions of long loop regions and the modeling of receptor activation states are presented as well.

Viral Evolved Inhibition Mechanism of the RNA Dependent Protein Kinase PKR's Kinase Domain, a Structural Perspective

The protein kinase PKR activated by viral dsRNA, phosphorylates the eIF2α, which inhibit the mechanism of translation initiation. Viral evolved proteins mimicking the eIF2α block its phosphorylation and help in the viral replication. To decipher the molecular basis for the PKR’s substrate and inhibitor interaction mechanisms, we carried the molecular dynamics studies on the catalytic domain of PKR in complex with substrate eIF2α, and inhibitors TAT and K3L. The studies conducted show the altered domain movements of N lobe, which confers open and close state to the substrate-binding cavity. In addition, PKR exhibits variations in the secondary structural transition of the activation loop residues, and inter molecular contacts with the substrate and the inhibitors. Phosphorylation of the P+1 loop at the Thr-451 increases the affinity of the binding proteins exhibiting its role in the phosphorylation events. The implications of structural mechanisms uncovered will help to understand the basis of the evolution of the host-viral and the viral replication mechanisms.

Computational Study of Binding of μ-Conotoxin GIIIA to Bacterial Sodium Channels NaVAb and NaVRh

Structures of several voltage-gated sodium (NaV) channels from bacteria have been determined recently, but the same feat might not be achieved for the mammalian counterparts in the near future. Thus, at present, computational studies of the mammalian NaV channels have to be performed using homology models based on the bacterial crystal structures. A successful homology model for the mammalian NaV1.4 channel was recently constructed using the extensive mutation data for binding of μ-conotoxin GIIIA to NaV1.4, which was further validated through studies of binding of other μ-conotoxins and ion permeation. Understanding the similarities and differences between the bacterial and mammalian NaV channels is an important issue, and the NaV1.4-GIIIA system provides a good opportunity for such a comparison. To this end, we study the binding of GIIIA to the bacterial channels NaVAb and NaVRh. The complex structures are obtained using docking and molecular dynamics simulations, and the dissociation of GIIIA is studied through umbrella sampling simulations. The results are compared to those obtained from the NaV1.4-GIIIA system, and the differences in the binding modes arising from the changes in the selectivity filters are highlighted.

A Flexible Domain-Domain Hinge Promotes an Induced-fit Dominant Mechanism for the Loading of Guide-DNA into Argonaute Protein in Thermus thermophilus

Argonaute proteins (Ago) are core components of the RNA Induced Silencing Complex (RISC) that load and utilize small guide nucleic acids to silence mRNAs or cleave foreign DNAs. Despite the essential role of Ago in gene regulation and defense against virus, the molecular mechanism of guide-strand loading into Ago remains unclear. We explore such a mechanism in the bacterium Thermus thermophilus Ago (TtAgo), via a computational approach combining molecular dynamics, bias-exchange metadynamics, and protein-DNA docking. We show that apo TtAgo adopts multiple closed states that are unable to accommodate guide-DNA. Conformations able to accommodate the guide are beyond the reach of thermal fluctuations from the closed states. These results suggest an induced-fit dominant mechanism for guide-strand loading in TtAgo, drastically different from the two-step mechanism for human Ago 2 (hAgo2) identified in our previous study. Such a difference between TtAgo and hAgo2 is found to mainly originate from the distinct rigidity of their L1-PAZ hinge. Further comparison among known Ago structures from various species indicates that the L1-PAZ hinge may be flexible in general for prokaryotic Ago's but rigid for eukaryotic Ago's.

Molecular characterization of HIV-1 Nef and ACOT8 interaction: insights from in silico structural predictions and in vitro functional assays

HIV-1 Nef interacts with several cellular proteins, among which the human peroxisomal thioesterase 8 (ACOT8). This interaction may be involved in the endocytosis regulation of membrane proteins and might modulate lipid composition in membrane rafts. Nef regions involved in the interaction have been experimentally characterized, whereas structural details of the ACOT8 protein are unknown. The lack of structural information hampers the comprehension of the functional consequences of the complex formation during HIV-1 infection. We modelled, through in silico predictions, the ACOT8 structure and we observed a high charge complementarity between Nef and ACOT8 surfaces, which allowed the identification of the ACOT8 putative contact points involved in the interaction. The predictions were validated by in vitro assays through the development of ACOT8 deletion mutants. Coimmunoprecipitation and immunofluorescence analyses showed that ACOT8 Arg⁴⁵-Phe⁵⁵ and Arg⁸⁶-Pro⁹³ regions are involved in Nef association. In addition, K91S mutation abrogated the interaction with Nef, indicating that Lys⁹¹ plays a key role in the interaction. Finally, when associated with ACOT8, Nef may be preserved from degradation. These findings improve the comprehension of the association between HIV-1 Nef and ACOT8, helping elucidating the biological effect of their interaction.

A deep learning framework for modeling structural features of RNA-binding protein targets

RNA-binding proteins (RBPs) play important roles in the post-transcriptional control of RNAs. Identifying RBP binding sites and characterizing RBP binding preferences are key steps toward understanding the basic mechanisms of the post-transcriptional gene regulation. Though numerous computational methods have been developed for modeling RBP binding preferences, discovering a complete structural representation of the RBP targets by integrating their available structural features in all three dimensions is still a challenging task. In this paper, we develop a general and flexible deep learning framework for modeling structural binding preferences and predicting binding sites of RBPs, which takes (predicted) RNA tertiary structural information into account for the first time. Our framework constructs a unified representation that characterizes the structural specificities of RBP targets in all three dimensions, which can be further used to predict novel candidate binding sites and discover potential binding motifs. Through testing on the real CLIP-seq datasets, we have demonstrated that our deep learning framework can automatically extract effective hidden structural features from the encoded raw sequence and structural profiles, and predict accurate RBP binding sites. In addition, we have conducted the first study to show that integrating the additional RNA tertiary structural features can improve the model performance in predicting RBP binding sites, especially for the polypyrimidine tract-binding protein (PTB), which also provides a new evidence to support the view that RBPs may own specific tertiary structural binding preferences. In particular, the tests on the internal ribosome entry site (IRES) segments yield satisfiable results with experimental support from the literature and further demonstrate the necessity of incorporating RNA tertiary structural information into the prediction model. The source code of our approach can be found in https://github.com/thucombio/deepnet-rbp.

Structural basis for PHDVC5HCHNSD1-C2HRNizp1 interaction: implications for Sotos syndrome

Sotos syndrome is an overgrowth syndrome caused by mutations within the functional domains ofNSD1 gene coding for NSD1, a multidomain protein regulating chromatin structure and gene expression. In particular, PHDVC5HCHNSD1 tandem domain, composed by a classical (PHDV) and an atypical (C5HCH) plant homeo-domain (PHD) finger, is target of several pathological missense-mutations. PHDVC5HCHNSD1 is also crucial for NSD1-dependent transcriptional regulation and interacts with the C2HR domain of transcriptional repressor Nizp1 (C2HRNizp1)in vitro To get molecular insights into the mechanisms dictating the patho-physiological relevance of the PHD finger tandem domain, we solved its solution structure and provided a structural rationale for the effects of seven Sotos syndrome point-mutations. To investigate PHDVC5HCHNSD1 role as structural platform for multiple interactions, we characterized its binding to histone H3 peptides and to C2HRNizp1 by ITC and NMR. We observed only very weak electrostatic interactions with histone H3 N-terminal tails, conversely we proved specific binding to C2HRNizp1 We solved C2HRNizp1 solution structure and generated a 3D model of the complex, corroborated by site-directed mutagenesis. We suggest a mechanistic scenario where NSD1 interactions with cofactors such as Nizp1 are impaired by PHDVC5HCHNSD1 pathological mutations, thus impacting on the repression of growth-promoting genes, leading to overgrowth conditions.

Structure prediction and docking-based molecular insights of human YB-1 and nucleic acid interaction

Y-box-binding protein 1 (YB-1), a cold shock domain protein, is one of the most conserved nucleic acid-binding proteins. The multifunctional human YB-1 is a member of a large family of proteins with an evolutionary ancient cold shock domain. The presence of a cold shock domain is a specific feature of Y-box-binding proteins and allows attributing them to a wider group of proteins containing a cold shock domain. This protein is involved in a number of cellular processes including proliferation, differentiation and stress response. The YB-1 performs its function both in the cytoplasm and in the cell nucleus. In this study, we present the structure of full-length human YB-1 protein along with investigation of their nucleic acid-binding preferential. The study also focuses on biases for particular purine and pyrimidine bases. The overall goal of this study was to model and validate full-length YB-1 protein and to compare its nucleic acid-binding studies with previous reports.

Crystal Structure of a Ube2S-Ubiquitin Conjugate

Protein ubiquitination occurs through the sequential formation and reorganization of specific protein-protein interfaces. Ubiquitin-conjugating (E2) enzymes, such as Ube2S, catalyze the formation of an isopeptide linkage between the C-terminus of a "donor" ubiquitin and a primary amino group of an "acceptor" ubiquitin molecule. This reaction involves an intermediate, in which the C-terminus of the donor ubiquitin is thioester-bound to the active site cysteine of the E2 and a functionally important interface is formed between the two proteins. A docked model of a Ube2S-donor ubiquitin complex was generated previously, based on chemical shift mapping by NMR, and predicted contacts were validated in functional studies. We now present the crystal structure of a covalent Ube2S-ubiquitin complex. The structure contains an interface between Ube2S and ubiquitin in trans that resembles the earlier model in general terms, but differs in detail. The crystallographic interface is more hydrophobic than the earlier model and is stable in molecular dynamics (MD) simulations. Remarkably, the docked Ube2S-donor complex converges readily to the configuration seen in the crystal structure in 3 out of 8 MD trajectories. Since the crystallographic interface is fully consistent with mutational effects, this indicates that the structure provides an energetically favorable representation of the functionally critical Ube2S-donor interface.

Cryo-EM Data Are Superior to Contact and Interface Information in Integrative Modeling

Protein-protein interactions carry out a large variety of essential cellular processes. Cryo-electron microscopy (cryo-EM) is a powerful technique for the modeling of protein-protein interactions at a wide range of resolutions, and recent developments have caused a revolution in the field. At low resolution, cryo-EM maps can drive integrative modeling of the interaction, assembling existing structures into the map. Other experimental techniques can provide information on the interface or on the contacts between the monomers in the complex. This inevitably raises the question regarding which type of data is best suited to drive integrative modeling approaches. Systematic comparison of the prediction accuracy and specificity of the different integrative modeling paradigms is unavailable to date. Here, we compare EM-driven, interface-driven, and contact-driven integrative modeling paradigms. Models were generated for the protein docking benchmark using the ATTRACT docking engine and evaluated using the CAPRI two-star criterion. At 20 Å resolution, EM-driven modeling achieved a success rate of 100%, outperforming the other paradigms even with perfect interface and contact information. Therefore, even very low resolution cryo-EM data is superior in predicting heterodimeric and heterotrimeric protein assemblies. Our study demonstrates that a force field is not necessary, cryo-EM data alone is sufficient to accurately guide the monomers into place. The resulting rigid models successfully identify regions of conformational change, opening up perspectives for targeted flexible remodeling.

Structural characterization of ANGPTL8 (betatrophin) with its interacting partner lipoprotein lipase

Angiopoietin-like protein 8 (ANGPTL8) (also known as betatrophin) is a newly identified secretory protein with a potential role in autophagy, lipid metabolism and pancreatic beta-cell proliferation. Its structural characterization is required to enhance our current understanding of its mechanism of action which could help in identifying its receptor and/or other binding partners. Based on the physiological significance and necessity of exploring structural features of ANGPTL8, the present study is conducted with a specific aim to model the structure of ANGPTL8 and study its possible interactions with Lipoprotein Lipase (LPL). To the best of our knowledge, this is the first attempt to predict 3-dimensional (3D) structure of ANGPTL8. Three different approaches were used for modeling of ANGPTL8 including homology modeling, de-novo structure prediction and their amalgam which is then proceeded by structure verification using ERRATT, PROSA, Qmean and Ramachandran plot scores. The selected models of ANGPTL8 were further evaluated for protein-protein interaction (PPI) analysis with LPL using CPORT and HADDOCK server. Our results have shown that the crystal structure of iSH2 domain of Phosphatidylinositol 3-kinase (PI3K) p85β subunit (PDB entry: 3mtt) is a good candidate for homology modeling of ANGPTL8. Analysis of inter-molecular interactions between the structure of ANGPTL8 and LPL revealed existence of several non-covalent interactions. The residues of LPL involved in these interactions belong from its lid region, thrombospondin (TSP) region and heparin binding site which is suggestive of a possible role of ANGPTL8 in regulating the proteolysis, motility and localization of LPL. Besides, the conserved residues of SE1 region of ANGPTL8 formed interactions with the residues around the hinge region of LPL. Overall, our results support a model of inhibition of LPL by ANGPTL8 through the steric block of its catalytic site which will be further explored using wet lab studies in future.

A target-protection mechanism of antibiotic resistance at atomic resolution: insights into FusB-type fusidic acid resistance

Antibiotic resistance in clinically important bacteria can be mediated by proteins that physically associate with the drug target and act to protect it from the inhibitory effects of an antibiotic. We present here the first detailed structural characterization of such a target protection mechanism mediated through a protein-protein interaction, revealing the architecture of the complex formed between the FusB fusidic acid resistance protein and the drug target (EF-G) it acts to protect. Binding of FusB to EF-G induces conformational and dynamic changes in the latter, shedding light on the molecular mechanism of fusidic acid resistance.

Nuclear Magnetic Resonance Insight into the Multiple Glycosaminoglycan Binding Modes of the Link Module from Human TSG-6

Tumor necrosis factor-stimulated gene-6 (TSG-6) is a hyaluronan (HA)-binding protein that is essential for stabilizing and remodeling the extracellular matrix (ECM) during ovulation and inflammatory disease processes such as arthritis. The Link module, one of the domains of TSG-6, is responsible for binding hyaluronan and other glycosaminoglycans found in the ECM. In this study, we used a well-defined chondroitin sulfate (CS) hexasaccharide (ΔC444S) to determine the structure of the Link module, in solution, in its chondroitin sulfate-bound state. A variety of nuclear magnetic resonance techniques were employed, including chemical shift perturbation, residual dipolar couplings (RDCs), nuclear Overhauser effects, spin relaxation measurements, and paramagnetic relaxation enhancements from a spin-labeled analogue of ΔC444S. The binding site for ΔC444S on the Link module overlapped with that of HA. Surprisingly, ΔC444S binding induced dimerization of the Link module (as confirmed by analytical ultracentrifugation), and a second weak binding site that partially overlapped with a previously identified heparin site was detected. A dimer model was generated using chemical shift perturbations and RDCs as restraints in the docking program HADDOCK. We postulate that the molecular cross-linking enhanced by the multiple binding modes of the Link module might be critical for remodeling the ECM during inflammation/ovulation and might contribute to other functions of TSG-6.

Different Types of Molecular Docking Based on Variations of Interacting Molecules

Molecular docking plays an important role in drug discovery research by facilitating target identification, target validation, virtual screening for lead identification and lead optimization. Depending upon the nature of the disease of interest, targets can be either protein or DNA while drugs are mostly organic small molecules. Different types of molecular docking techniques like protein-protein or protein-DNA or protein-small molecule or DNA-small molecule are employed for achieving the above mentioned objectives. This chapter provides a clear idea of the position of molecular docking in drug discovery with detailed discussion on different types of molecular docking based on the varieties of interacting partners. Subsequently the authors provide a detailed list of tools that can be used for docking in drug discovery and discus some examples of molecular docking in drug discovery before concluding with a remark on future areas of improvement in molecular docking related to drug discovery.

An investigation on the interaction modes of a single-strand DNA aptamer and RBP4 protein: a molecular dynamic simulations approach

Monitoring and evaluating structural and functional alternations in RBP4 induced by its specific aptamer binding to design new aptamers for diagnostic and therapeutic purposes with reduced insulin resistance.

Production, Isolation, and Structural Analysis of Ligands and Receptors of the TGF-β Superfamily

The ability to understand the molecular mechanisms by which secreted signaling proteins of the TGF-β superfamily assemble their cell surface receptors into complexes to initiate downstream signaling is dependent upon the ability to determine atomic-resolution structures of the signaling proteins, the ectodomains of the receptors, and the complexes that they form. The structures determined to date have revealed major differences in the overall architecture of the signaling complexes formed by the TGF-βs and BMPs, which has provided insights as to how they have evolved to fulfill their distinct functions. Such studies, have however, only been applied to a few members of the TGF-β superfamily, which is largely due to the difficulty of obtaining milligram-scale quantities of highly homogenous preparations of the disulfide-rich signaling proteins and receptor ectodomains of the superfamily. Here we describe methods used to produce signaling proteins and receptor ectodomains of the TGF-β superfamily using bacterial and mammalian expression systems and procedures to purify them to homogeneity.

K, Subbarao N. Docking and molecular dynamics simulation study of EGFR1 with EGF-like peptides to understand molecular interactions

Epidermal growth factor receptors (EGFRs) are oncogenes, which regulate the expression of genes in various pathways, allowing cells to grow and divide.

Structural transition of ETS1 from an auto-inhibited to functional state upon association with the p16INK4anative and mutated promoter region

Detailed elucidation of structural changes invoked on transcriptional factors and their target genes upon their association is pivotal for understanding the genetic level regulations imposed in several diseases including ovarian cancer.

CXCL1/MGSA Is a Novel Glycosaminoglycan (GAG)-binding Chemokine: STRUCTURAL EVIDENCE FOR TWO DISTINCT NON-OVERLAPPING BINDING DOMAINS

In humans, the chemokine CXCL1/MGSA (hCXCL1) plays fundamental and diverse roles in pathophysiology, from microbial killing to cancer progression, by orchestrating the directed migration of immune and non-immune cells. Cellular trafficking is highly regulated and requires concentration gradients that are achieved by interactions with sulfated glycosaminoglycans (GAGs). However, very little is known regarding the structural basis underlying hCXCL1-GAG interactions. We addressed this by characterizing the binding of GAG heparin oligosaccharides to hCXCL1 using NMR spectroscopy. Binding experiments under conditions at which hCXCL1 exists as monomers and dimers indicate that the dimer is the high-affinity GAG ligand. NMR experiments and modeling studies indicate that lysine and arginine residues mediate binding and that they are located in two non-overlapping domains. One domain, consisting of N-loop and C-helical residues (defined as α-domain) has also been identified previously as the GAG-binding domain for the related chemokine CXCL8/IL-8. The second domain, consisting of residues from the N terminus, 40s turn, and third β-strand (defined as β-domain) is novel. Eliminating β-domain binding by mutagenesis does not perturb α-domain binding, indicating two independent GAG-binding sites. It is known that N-loop and N-terminal residues mediate receptor activation, and we show that these residues are also involved in extensive GAG interactions. We also show that the GAG-bound hCXCL1 completely occlude receptor binding. We conclude that hCXCL1-GAG interactions provide stringent control over regulating chemokine levels and receptor accessibility and activation, and that chemotactic gradients mediate cellular trafficking to the target site.

Solution NMR Spectroscopy Provides an Avenue for the Study of Functionally Dynamic Molecular Machines: The Example of Protein Disaggregation

Solution-based NMR spectroscopy has been an important tool for studying the structure and dynamics of relatively small proteins and protein complexes with aggregate molecular masses under approximately 50 kDa. The development of new experiments and labeling schemes, coupled with continued improvements in hardware, has significantly reduced this size limitation, enabling atomic-resolution studies of molecular machines in the 1 MDa range. In this Perspective, some of the important advances are highlighted in the context of studies of molecular chaperones involved in protein disaggregation. New insights into the structural biology of disaggregation obtained from NMR studies are described, focusing on the unique capabilities of the methodology for obtaining atomic-resolution descriptions of dynamic systems.

Interacting mechanism of ID3 HLH domain towards E2A/E12 transcription factor - An Insight through molecular dynamics and docking approach

Inhibitor of DNA binding protein 3 (ID3) has long been characterized as an oncogene that implicates its functional role through its Helix-Loop-Helix (HLH) domain upon protein-protein interaction. An insight into the dimerization brought by this domain helps in identifying the key residues that favor the mechanism behind it. Molecular dynamics (MD) simulations were performed for the HLH proteins ID3 and Transcription factor E2-alpha (E2A/E12) and their ensemble complex (ID3-E2A/E12) to gather information about the HLH domain region and its role in the interaction process. Further evaluation of the results by Principal Component Analysis (PCA) and Free Energy Landscape (FEL) helped in revealing residues of E2A/E12: Lys570, Ala595, Val598, and Ile599 and ID3: Glu53, Gln63, and Gln66 buried in their HLH motifs imparting key roles in dimerization process. Furthermore the T-pad analysis results helped in identifying the key fluctuations and conformational transitions using the intrinsic properties of the residues present in the domain region of the proteins thus specifying their crucial role towards molecular recognition. The study provides an insight into the interacting mechanism of the ID3-E2A/E12 complex and maps the structural transitions arising in the essential conformational space indicating the key structural changes within the helical regions of the motif. It thereby describes how the internal dynamics of the proteins might regulate their intrinsic structural features and its subsequent functionality.

Accurate refinement of docked protein complexes using evolutionary information and deep learning

One of the major challenges for protein docking methods is to accurately discriminate native-like structures from false positives. Docking methods are often inaccurate and the results have to be refined and re-ranked to obtain native-like complexes and remove outliers. In a previous work, we introduced AccuRefiner, a machine learning based tool for refining protein-protein complexes. Given a docked complex, the refinement tool produces a small set of refined versions of the input complex, with lower root-mean-square-deviation (RMSD) of atomic positions with respect to the native structure. The method employs a unique ranking tool that accurately predicts the RMSD of docked complexes with respect to the native structure. In this work, we use a deep learning network with a similar set of features and five layers. We show that a properly trained deep learning network can accurately predict the RMSD of a docked complex with 1.40 Å error margin on average, by approximating the complex relationship between a wide set of scoring function terms and the RMSD of a docked structure. The network was trained on 35000 unbound docking complexes generated by RosettaDock. We tested our method on 25 different putative docked complexes produced also by RosettaDock for five proteins that were not included in the training data. The results demonstrate that the high accuracy of the ranking tool enables AccuRefiner to consistently choose the refinement candidates with lower RMSD values compared to the coarsely docked input structures.

Transient collagen triple helix binding to a key metalloproteinase in invasion and development

Skeletal development and invasion by tumor cells depends on proteolysis of collagen by the pericellular metalloproteinase MT1-MMP. Its hemopexin-like (HPX) domain binds to collagen substrates to facilitate their digestion. Spin labeling and paramagnetic nuclear magnetic resonance (NMR) detection have revealed how the HPX domain docks to collagen I-derived triple helix. Mutations impairing triple-helical peptidase activity corroborate the interface. Saturation transfer difference NMR suggests rotational averaging around the longitudinal axis of the triple-helical peptide. Part of the interface emerges as unique and potentially targetable for selective inhibition. The triple helix crosses the junction of blades I and II at a 45° angle to the symmetry axis of the HPX domain, placing the scissile Gly∼Ile bond near the HPX domain and shifted ∼25 Å from MMP-1 complexes. This raises the question of the MT1-MMP catalytic domain folding over the triple helix during catalysis, a possibility accommodated by the flexibility between domains suggested by atomic force microscopy images.

Characterization of the interaction between Escherichia coli topoisomerase IV E subunit and an ATP competitive inhibitor

Bacterial topoisomerase IV (ParE) is essential for DNA replication and serves as an attractive target for antibacterial drug development. The X-ray structure of the N-terminal 24 kDa ParE, responsible for ATP binding has been solved. Due to the accessibility of structural information of ParE, many potent ParE inhibitors have been discovered. In this study, a pyridylurea lead molecule against ParE of Escherichia coli (eParE) was characterized with a series of biochemical and biophysical techniques. More importantly, solution NMR analysis of compound binding to eParE provides better understanding of the molecular interactions between the inhibitor and eParE.

The HADDOCK2.2 Web Server: User-Friendly Integrative Modeling of Biomolecular Complexes

The prediction of the quaternary structure of biomolecular macromolecules is of paramount importance for fundamental understanding of cellular processes and drug design. In the era of integrative structural biology, one way of increasing the accuracy of modeling methods used to predict the structure of biomolecular complexes is to include as much experimental or predictive information as possible in the process. This has been at the core of our information-driven docking approach HADDOCK. We present here the updated version 2.2 of the HADDOCK portal, which offers new features such as support for mixed molecule types, additional experimental restraints and improved protocols, all of this in a user-friendly interface. With well over 6000 registered users and 108,000 jobs served, an increasing fraction of which on grid resources, we hope that this timely upgrade will help the community to solve important biological questions and further advance the field. The HADDOCK2.2 Web server is freely accessible to non-profit users at http://haddock.science.uu.nl/services/HADDOCK2.2.

In situ structural analysis of the human nuclear pore complex

Nuclear pore complexes are fundamental components of all eukaryotic cells that mediate nucleocytoplasmic exchange. Determining their 110-megadalton structure imposes a formidable challenge and requires in situ structural biology approaches. Of approximately 30 nucleoporins (Nups), 15 are structured and form the Y and inner-ring complexes. These two major scaffolding modules assemble in multiple copies into an eight-fold rotationally symmetric structure that fuses the inner and outer nuclear membranes to form a central channel of ~60 nm in diameter. The scaffold is decorated with transport-channel Nups that often contain phenylalanine-repeat sequences and mediate the interaction with cargo complexes. Although the architectural arrangement of parts of the Y complex has been elucidated, it is unclear how exactly it oligomerizes in situ. Here we combine cryo-electron tomography with mass spectrometry, biochemical analysis, perturbation experiments and structural modelling to generate, to our knowledge, the most comprehensive architectural model of the human nuclear pore complex to date. Our data suggest previously unknown protein interfaces across Y complexes and to inner-ring complex members. We show that the transport-channel Nup358 (also known as Ranbp2) has a previously unanticipated role in Y-complex oligomerization. Our findings blur the established boundaries between scaffold and transport-channel Nups. We conclude that, similar to coated vesicles, several copies of the same structural building block--although compositionally identical--engage in different local sets of interactions and conformations.

Solution NMR characterization of chemokine CXCL8/IL-8 monomer and dimer binding to glycosaminoglycans: structural plasticity mediates differential binding interactions

Structural plasticity plays a major role in determining differential binding of CXCL8 monomer and dimer to glycosaminoglycans (GAGs) and that dimer is the high-affinity GAG ligand. We propose that these properties play important roles in orchestrating in vivo chemokine-mediated neutrophil function.

Chemokine CXCL8/interleukin-8 (IL-8) plays a crucial role in directing neutrophils and oligodendrocytes to combat infection/injury and tumour cells in metastasis development. CXCL8 exists as monomers and dimers and interaction of both forms with glycosaminoglycans (GAGs) mediate these diverse cellular processes. However, very little is known regarding the structural basis underlying CXCL8–GAG interactions. There are conflicting reports on the affinities, geometry and whether the monomer or dimer is the high-affinity GAG ligand. To resolve these issues, we characterized the binding of a series of heparin-derived oligosaccharides [heparin disaccharide (dp2), heparin tetrasaccharide (dp4), heparin octasaccharide (dp8) and heparin 14-mer (dp14)] to the wild-type (WT) dimer and a designed monomer using solution NMR spectroscopy. The pattern and extent of binding-induced chemical shift perturbation (CSP) varied between dimer and monomer and between longer and shorter oligosaccharides. NMR-based structural models show that different interaction modes coexist and that the nature of interactions varied between monomer and dimer and oligosaccharide length. MD simulations indicate that the binding interface is structurally plastic and provided residue-specific details of the dynamic nature of the binding interface. Binding studies carried out under conditions at which WT CXCL8 exists as monomers and dimers provide unambiguous evidence that the dimer is the high-affinity GAG ligand. Together, our data indicate that a set of core residues function as the major recognition/binding site, a set of peripheral residues define the various binding geometries and that the structural plasticity of the binding interface allows multiplicity of binding interactions. We conclude that structural plasticity most probably regulates in vivo CXCL8 monomer/dimer–GAG interactions and function.

Molecular view of ligands specificity for CAG repeats in anti-Huntington therapy

Huntington's disease is a fatal and devastating neurodegenerative genetic disorder for which there is currently no cure. It is characterized by Huntingtin protein's mRNA transcripts with 36 or more CAG repeats. Inhibiting the formation of pathological complexes between these expanded transcripts and target proteins may be a valuable strategy against the disease. Yet, the rational design of molecules specifically targeting the expanded CAG repeats is limited by the lack of structural information. Here, we use well-tempered metadynamics-based free energy calculations to investigate pose and affinity of two ligands targeting CAG repeats for which affinities have been previously measured. The first consists of two 4-guanidinophenyl rings linked by an ester group. It is the most potent ligand identified so far, with Kd = 60(30) nM. The second consists of a 4-phenyl dihydroimidazole and 4-1H-indole dihydroimidazole connected by a C-C bond (Kd = 700(80) nM). Our calculations reproduce the experimental affinities and uncover the recognition pattern between ligands' and their RNA target. They also provide a molecular basis for the markedly different affinity of the two ligands for CAG repeats as observed experimentally. These findings may pave the way for a structure-based hit-to-lead optimization to further improve ligand selectivity toward CAG repeat-containing mRNAs.

Solid-state NMR Study of the YadA Membrane-Anchor Domain in the Bacterial Outer Membrane

MAS-NMR was used to study the structure and dynamics at ambient temperatures of the membrane-anchor domain of YadA (YadA-M) in a pellet of the outer membrane of E. coli in which it was expressed. YadA is an adhesin from the pathogen Yersinia enterocolitica that is involved in interactions with the host cell, and it is a model protein for studying the autotransport process. Existing assignments were sucessfully transferred to a large part of the YadA-M protein in the E. coli lipid environment by using (13) C-(13) C DARR and PDSD spectra at different mixing times. The chemical shifts in most regions of YadA-M are unchanged relative to those in microcrystalline YadA-M preparations from which a structure has previously been solved, including the ASSA region that is proposed to be involved in transition-state hairpin formation for transport of the soluble domain. Comparisons of the dynamics between the microcrystalline and membrane-embedded samples indicate greater flexibility of the ASSA region in the outer-membrane preparation at physiological temperatures. This study will pave the way towards MAS-NMR structure determination of membrane proteins, and a better understanding of functionally important dynamic residues in native membrane environments.

Evolution of Macromolecular Docking Techniques: The Case Study of Nickel and Iron Metabolism in Pathogenic Bacteria

The interaction between macromolecules is a fundamental aspect of most biological processes. The computational techniques used to study protein-protein and protein-nucleic acid interactions have evolved in the last few years because of the development of new algorithms that allow the a priori incorporation, in the docking process, of experimentally derived information, together with the possibility of accounting for the flexibility of the interacting molecules. Here we review the results and the evolution of the techniques used to study the interaction between metallo-proteins and DNA operators, all involved in the nickel and iron metabolism of pathogenic bacteria, focusing in particular on Helicobacter pylori (Hp). In the first part of the article we discuss the methods used to calculate the structure of complexes of proteins involved in the activation of the nickel-dependent enzyme urease. In the second part of the article, we concentrate on two applications of protein-DNA docking conducted on the transcription factors HpFur (ferric uptake regulator) and HpNikR (nickel regulator). In both cases we discuss the technical expedients used to take into account the conformational variability of the multi-domain proteins involved in the calculations.

On the interaction of Helicobacter pylori NikR, a Ni(II)-responsive transcription factor, with the urease operator: in solution and in silico studies

Helicobacter pylori (Hp) is a carcinogen that relies on Ni(II) to survive in the extreme pH conditions of the human guts. The regulation of genes coding for Ni(II) enzymes and proteins is effected by the nickel-responsive transcription factor NikR, composed of a DNA-binding domain (DBD) and a metal-binding domain (MBD). The scope of this study is to obtain the molecular details of the HpNikR interaction with the urease operator OP ureA , in solution. The size of the full-length protein prevents the characterization of the HpNikR-OP ureA interaction using NMR. We thus investigated the two separate domains of HpNikR. The conservation of their oligomeric state was established by multiple-angle light scattering. Isothermal calorimetric titrations indicated that the thermodynamics of Ni(II) binding to the isolated MBD is independent of the presence of the adjacent DBDs. The NMR spectra of the isolated DBD support considerable conservation of its structural properties. The spectral perturbations induced on the DBD by OP ureA provided information useful to calculate a structural model of the HpNikR-OP ureA complex using a docking computational protocol. The NMR assignment of the residues involved in the protein-DNA interaction represents a starting point for the development of drugs potentially able to eradicate H. pylori infections. All evidences so far collected, in this and previous studies, consistently indicate that binding of Ni(II) to the MBD increases the HpNikR-DNA affinity by modulating the dynamic, and not the structural, properties of the protein, suggesting that the formation of a stable complex relies upon an induced fit mechanism.

Markov State Models Reveal a Two-Step Mechanism of miRNA Loading into the Human Argonaute Protein: Selective Binding followed by Structural Re-arrangement

Argonaute (Ago) proteins and microRNAs (miRNAs) are central components in RNA interference, which is a key cellular mechanism for sequence-specific gene silencing. Despite intensive studies, molecular mechanisms of how Ago recognizes miRNA remain largely elusive. In this study, we propose a two-step mechanism for this molecular recognition: selective binding followed by structural re-arrangement. Our model is based on the results of a combination of Markov State Models (MSMs), large-scale protein-RNA docking, and molecular dynamics (MD) simulations. Using MSMs, we identify an open state of apo human Ago-2 in fast equilibrium with partially open and closed states. Conformations in this open state are distinguished by their largely exposed binding grooves that can geometrically accommodate miRNA as indicated in our protein-RNA docking studies. miRNA may then selectively bind to these open conformations. Upon the initial binding, the complex may perform further structural re-arrangement as shown in our MD simulations and eventually reach the stable binary complex structure. Our results provide novel insights in Ago-miRNA recognition mechanisms and our methodology holds great potential to be widely applied in the studies of other important molecular recognition systems.

In RNA interference, Argonaute proteins and microRNAs together form the functional core that regulates the gene expression with high sequence specificity. Elucidating the detailed mechanism of molecular recognition between Argonaute proteins and microRNAs is thus important not only for the fundamental understanding of RNA interference, but also for the further development of microRNA-based therapeutic application. In this work, we propose a two-step model to understand the mechanism of microRNA loading into human Argonaute-2: selective binding followed by structural re-arrangement. Our model is based on the results from a combined approach of molecular dynamics simulations, Markov State Models and protein-RNA docking. In particular, we identify a metastable open state of apo hAgo2 in rapid equilibrium with other states. Some of conformations in this open state have largely exposed RNA binding groove that can accommodate microRNA. We further show that the initial Argonaute-microRNA binding complex undergoes structural re-arrangement to reach stable binary crystal structure. These results provide novel insights into the underlying mechanism of Argonaute-microRNA recognition. In addition, our method is readily applicable to the investigation of other complex molecular recognition events such as protein-protein interactions and protein-ligand binding.

Structural basis for cellobiose dehydrogenase action during oxidative cellulose degradation

A new paradigm for cellulose depolymerization by fungi focuses on an oxidative mechanism involving cellobiose dehydrogenases (CDH) and copper-dependent lytic polysaccharide monooxygenases (LPMO); however, mechanistic studies have been hampered by the lack of structural information regarding CDH. CDH contains a haem-binding cytochrome (CYT) connected via a flexible linker to a flavin-dependent dehydrogenase (DH). Electrons are generated from cellobiose oxidation catalysed by DH and shuttled via CYT to LPMO. Here we present structural analyses that provide a comprehensive picture of CDH conformers, which govern the electron transfer between redox centres. Using structure-based site-directed mutagenesis, rapid kinetics analysis and molecular docking, we demonstrate that flavin-to-haem interdomain electron transfer (IET) is enabled by a haem propionate group and that rapid IET requires a closed CDH state in which the propionate is tightly enfolded by DH. Following haem reduction, CYT reduces LPMO to initiate oxygen activation at the copper centre and subsequent cellulose depolymerization.

Computational approaches for designing potent and selective analogs of peptide toxins as novel therapeutics

Peptide toxins provide valuable therapeutic leads for many diseases. As they bind to their targets with high affinity, potency is usually ensured. However, toxins also bind to off-target receptors, causing potential side effects. Thus, a major challenge in generating drugs from peptide toxins is ensuring their specificity for their intended targets. Computational methods can play an important role in solving such design problems through construction of accurate models of receptor-toxin complexes and calculation of binding free energies. Here we review the computational methods used for this purpose and their application to toxins targeting ion channels. We describe ShK and HsTX1 toxins, high-affinity blockers of the voltage-gated potassium channel Kv1.3, which could be developed as therapeutic agents for autoimmune diseases.

Models for the Metal Transfer Complex of the N-Terminal Region of CusB and CusF

The tripartite CusCFBA pump in Escherichia coli is a very effective heavy metal extrusion system specific for Cu(I) and Ag(I). The N-terminal region of the membrane fusion protein CusB (CusB-NT) is highly disordered, and hence, experimentally characterizing its structure is challenging. In a previous study, this disorder was confirmed with molecular dynamics simulations, although some key structural elements were determined. It was experimentally shown that CusB-NT is fully functional in transferring the metal from the metallochaperone CusF. In this study, we docked these two entities together and formed two representative metal coordination modes, which consist of residues from both proteins. In this way, we created two potential CusB-NT/CusF complexes that share coordination of Cu(I) and thereby represent structural models for the metal transfer process. Each model complex was simulated for 4 μs. The previously observed structural disorder in CusB-NT disappeared upon complexation with CusF. The only differences between the two models occurred in the M21-M36 loop region of CusB-NT and the open flap of CusF: we observed the model with two CusB-NT methionine residues and a CusF methionine as the metal coordination site (termed "MMM") to be more stable than the model with a CusB-NT methionine, a CusF methionine, and a CusF histidine ligating the metal (termed "MMH"). The observed stability of the MMM model was probed for an additional 2 μs, yielding a total simulation time of 6 μs. We hypothesize that both MMM and MMH configurations might take part in the metal exchange process in which the MMH configuration would appear first and would be followed by the MMM configuration. Given the experimental finding of comparable binding affinities of CusB-NT and CusF, the increased stability of the MMM configuration might be a determinant for the transfer from CusF to CusB-NT. The metal would be transferred from the more CusF-dominated metal binding environment (MMH model) to a more CusB-dominated one (MMM model) in which the coordination environment is more stable. From the MMM model, the metal ion would ultimately be coordinated by the CusB methionines only, which would complete the Cu(I) transfer process.

Analysis and Ranking of Protein-Protein Docking Models Using Inter-Residue Contacts and Inter-Molecular Contact Maps

In view of the increasing interest both in inhibitors of protein-protein interactions and in protein drugs themselves, analysis of the three-dimensional structure of protein-protein complexes is assuming greater relevance in drug design. In the many cases where an experimental structure is not available, protein-protein docking becomes the method of choice for predicting the arrangement of the complex. However, reliably scoring protein-protein docking poses is still an unsolved problem. As a consequence, the screening of many docking models is usually required in the analysis step, to possibly single out the correct ones. Here, making use of exemplary cases, we review our recently introduced methods for the analysis of protein complex structures and for the scoring of protein docking poses, based on the use of inter-residue contacts and their visualization in inter-molecular contact maps. We also show that the ensemble of tools we developed can be used in the context of rational drug design targeting protein-protein interactions.

Heat Shock Protein 90 Associates with the Per-Arnt-Sim Domain of Heme-free Soluble Guanylate Cyclase: IMplications for Enzyme Maturation

Heat shock protein 90 (hsp90) drives heme insertion into the β1 subunit of soluble guanylate cyclase (sGC) β1, which enables it to associate with a partner sGCα1 subunit and mature into a nitric oxide (NO)-responsive active form. We utilized fluorescence polarization measurements and hydrogen-deuterium exchange mass spectrometry to define molecular interactions between the specific human isoforms hsp90β and apo-sGCβ1. hsp90β and its isolated M domain, but not its isolated N and C domains, bind with low micromolar affinity to a heme-free, truncated version of sGCβ1 (sGCβ1(1-359)-H105F). Surprisingly, hsp90β and its M domain bound to the Per-Arnt-Sim (PAS) domain of apo-sGC-β1(1-359), which lies adjacent to its heme-binding (H-NOX) domain. The interaction specifically involved solvent-exposed regions in the hsp90β M domain that are largely distinct from sites utilized by other hsp90 clients. The interaction strongly protected two regions of the sGCβ1 PAS domain and caused local structural relaxation in other regions, including a PAS dimerization interface and a segment in the H-NOX domain. Our results suggest a means by which the hsp90β interaction could prevent apo-sGCβ1 from associating with its partner sGCα1 subunit while enabling structural changes to assist heme insertion into the H-NOX domain. This mechanism would parallel that in other clients like the aryl hydrocarbon receptor and HIF1α, which also interact with hsp90 through their PAS domains to control protein partner and small ligand binding interactions.

Surfing the Protein-Protein Interaction Surface Using Docking Methods: Application to the Design of PPI Inhibitors

Blocking protein-protein interactions (PPI) using small molecules or peptides modulates biochemical pathways and has therapeutic significance. PPI inhibition for designing drug-like molecules is a new area that has been explored extensively during the last decade. Considering the number of available PPI inhibitor databases and the limited number of 3D structures available for proteins, docking and scoring methods play a major role in designing PPI inhibitors as well as stabilizers. Docking methods are used in the design of PPI inhibitors at several stages of finding a lead compound, including modeling the protein complex, screening for hot spots on the protein-protein interaction interface and screening small molecules or peptides that bind to the PPI interface. There are three major challenges to the use of docking on the relatively flat surfaces of PPI. In this review we will provide some examples of the use of docking in PPI inhibitor design as well as its limitations. The combination of experimental and docking methods with improved scoring function has thus far resulted in few success stories of PPI inhibitors for therapeutic purposes. Docking algorithms used for PPI are in the early stages, however, and as more data are available docking will become a highly promising area in the design of PPI inhibitors or stabilizers.

A Three-protein Charge Zipper Stabilizes a Complex Modulating Bacterial Gene Silencing

The Hha/YmoA nucleoid-associated proteins help selectively silence horizontally acquired genetic material, including pathogenicity and antibiotic resistance genes and their maintenance in the absence of selective pressure. Members of the Hha family contribute to gene silencing by binding to the N-terminal dimerization domain of H-NS and modifying its selectivity. Hha-like proteins and the H-NS N-terminal domain are unusually rich in charged residues, and their interaction is mostly electrostatic-driven but, nonetheless, highly selective. The NMR-based structural model of the complex between Hha/YmoA and the H-NS N-terminal dimerization domain reveals that the origin of the selectivity is the formation of a three-protein charge zipper with interdigitated complementary charged residues from Hha and the two units of the H-NS dimer. The free form of YmoA shows collective microsecond-millisecond dynamics that can by measured by NMR relaxation dispersion experiments and shows a linear dependence with the salt concentration. The number of residues sensing the collective dynamics and the population of the minor form increased in the presence of H-NS. Additionally, a single residue mutation in YmoA (D43N) abolished H-NS binding and the dynamics of the apo-form, suggesting the dynamics and binding are functionally related.

Structural basis of intramitochondrial phosphatidic acid transport mediated by Ups1-Mdm35 complex

Ups1 forms a complex with Mdm35 and is critical for the transport of phosphatidic acid (PA) from the mitochondrial outer membrane to the inner membrane. We report the crystal structure of the Ups1-Mdm35-PA complex and the functional characterization of Ups1-Mdm35 in PA binding and transfer. Ups1 features a barrel-like structure consisting of an antiparallel β-sheet and three α-helices. Mdm35 adopts a three-helical clamp-like structure to wrap around Ups1 to form a stable complex. The β-sheet and α-helices of Ups1 form a long tunnel-like pocket to accommodate the substrate PA, and a short helix α2 acts as a lid to cover the pocket. The hydrophobic residues lining the pocket and helix α2 are critical for PA binding and transfer. In addition, a hydrophilic patch on the surface of Ups1 near the PA phosphate-binding site also plays an important role in the function of Ups1-Mdm35. Our study reveals the molecular basis of the function of Ups1-Mdm35 and sheds new light on the mechanism of intramitochondrial phospholipid transport by the MSF1/PRELI family proteins.