FireDock: a web server for fast interaction refinement in molecular docking†

Structural details of protein–protein interactions are invaluable for understanding and deciphering biological mechanisms. Computational docking methods aim to predict the structure of a protein–protein complex given the structures of its single components. Protein flexibility and the absence of robust scoring functions pose a great challenge in the docking field. Due to these difficulties most of the docking methods involve a two-tier approach: coarse global search for feasible orientations that treats proteins as rigid bodies, followed by an accurate refinement stage that aims to introduce flexibility into the process. The FireDock web server, presented here, is the first web server for flexible refinement and scoring of protein–protein docking solutions. It includes optimization of side-chain conformations and rigid-body orientation and allows a high-throughput refinement. The server provides a user-friendly interface and a 3D visualization of the results. A docking protocol consisting of a global search by PatchDock and a refinement by FireDock was extensively tested. The protocol was successful in refining and scoring docking solution candidates for cases taken from docking benchmarks. We provide an option for using this protocol by automatic redirection of PatchDock candidate solutions to the FireDock web server for refinement. The FireDock web server is available at http://bioinfo3d.cs.tau.ac.il/FireDock/

Rigidity and flexibility of biological networks

The network approach became a widely used tool to understand the behaviour of complex systems in the last decade. We start from a short description of structural rigidity theory. A detailed account on the combinatorial rigidity analysis of protein structures, as well as local flexibility measures of proteins and their applications in explaining allostery and thermostability is given. We also briefly discuss the network aspects of cytoskeletal tensegrity. Finally, we show the importance of the balance between functional flexibility and rigidity in protein-protein interaction, metabolic, gene regulatory and neuronal networks. Our summary raises the possibility that the concepts of flexibility and rigidity can be generalized to all networks.Comment: 21 pages, 4 figures, 1 tabl

Digital monitoring and compensation of MDL based on higher-order Poincaré spheres

We present a digital technique able to monitor and compensate for the mode-dependent losses (MDL) in space-division multiplexing (SDM) transmission systems. The working principle of the technique is based on the analysis of the received signal samples in the higher-order Poincaré spheres (HoPs). When an arbitrary pair of tributaries is represented in the respective HoPs, the effect of the MDL can be modeled as a simple shift of the constellation points in a such sphere. Therefore, the MDL can be estimated by computing those shifts over all the HoPs and the induced signal distortions can be compensated by re-centering all the constellations in the respective HoPs. It should be highlighted that the proposed technique is scalable with an arbitrary number of spatial channels, modulation format agonistic and free of training sequences. The HoPs-based MDL monitoring (compensation) technique allows the MDL estimation (compensation) up to values of ≈ 6 dB. The proposed technique can partially compensate the MDL distortion, making a MDL sensitive algorithm in an insensitive one. When the proposed technique assists a HoPs-based space-demultiplexing algorithm, it provides signal-to-noise ratio (SNR) enhancements of 2, 4 and 8 dB for PM-QPSK, PM-16QAM and PM-64QAM signals, respectively, for the particular case of a SDM-based transmission system with a spatial diversity of 2 and 2 dB of MDL.publishe

The Solution Structures of Two Human IgG1 Antibodies Show Conformational Stability and Accommodate Their C1q and FcγR Ligands.

The human IgG1 antibody subclass shows distinct properties compared with the IgG2, IgG3, and IgG4 subclasses and is the most exploited subclass in therapeutic antibodies. It is the most abundant subclass, has a half-life as long as that of IgG2 and IgG4, binds the FcγR receptor, and activates complement. There is limited structural information on full-length human IgG1 because of the challenges of crystallization. To rectify this, we have studied the solution structures of two human IgG1 6a and 19a monoclonal antibodies in different buffers at different temperatures. Analytical ultracentrifugation showed that both antibodies were predominantly monomeric, with sedimentation coefficients s20,w (0) of 6.3-6.4 S. Only a minor dimer peak was observed, and the amount was not dependent on buffer conditions. Solution scattering showed that the x-ray radius of gyration Rg increased with salt concentration, whereas the neutron Rg values remained unchanged with temperature. The x-ray and neutron distance distribution curves P(r) revealed two peaks, M1 and M2, whose positions were unchanged in different buffers to indicate conformational stability. Constrained atomistic scattering modeling revealed predominantly asymmetric solution structures for both antibodies with extended hinge structures. Both structures were similar to the only known crystal structure of full-length human IgG1. The Fab conformations in both structures were suitably positioned to permit the Fc region to bind readily to its FcγR and C1q ligands without steric clashes, unlike human IgG4. Our molecular models for human IgG1 explain its immune activities, and we discuss its stability and function for therapeutic applications

GPU.proton.DOCK: Genuine Protein Ultrafast proton equilibria consistent DOCKing

GPU.proton.DOCK (Genuine Protein Ultrafast proton equilibria consistent DOCKing) is a state of the art service for in silico prediction of protein–protein interactions via rigorous and ultrafast docking code. It is unique in providing stringent account of electrostatic interactions self-consistency and proton equilibria mutual effects of docking partners. GPU.proton.DOCK is the first server offering such a crucial supplement to protein docking algorithms—a step toward more reliable and high accuracy docking results. The code (especially the Fast Fourier Transform bottleneck and electrostatic fields computation) is parallelized to run on a GPU supercomputer. The high performance will be of use for large-scale structural bioinformatics and systems biology projects, thus bridging physics of the interactions with analysis of molecular networks. We propose workflows for exploring in silico charge mutagenesis effects. Special emphasis is given to the interface-intuitive and user-friendly. The input is comprised of the atomic coordinate files in PDB format. The advanced user is provided with a special input section for addition of non-polypeptide charges, extra ionogenic groups with intrinsic pKa values or fixed ions. The output is comprised of docked complexes in PDB format as well as interactive visualization in a molecular viewer. GPU.proton.DOCK server can be accessed at http://gpudock.orgchm.bas.bg/

SwarmDock and the Use of Normal Modes in Protein-Protein Docking

Here is presented an investigation of the use of normal modes in protein-protein docking, both in theory and in practice. Upper limits of the ability of normal modes to capture the unbound to bound conformational change are calculated on a large test set, with particular focus on the binding interface, the subset of residues from which the binding energy is calculated. Further, the SwarmDock algorithm is presented, to demonstrate that the modelling of conformational change as a linear combination of normal modes is an effective method of modelling flexibility in protein-protein docking

Surfactant protein D inhibits HIV-1 infection of target cells via interference with gp120-CD4 interaction and modulates pro-inflammatory cytokine production

© 2014 Pandit et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.Surfactant Protein SP-D, a member of the collectin family, is a pattern recognition protein, secreted by mucosal epithelial cells and has an important role in innate immunity against various pathogens. In this study, we confirm that native human SP-D and a recombinant fragment of human SP-D (rhSP-D) bind to gp120 of HIV-1 and significantly inhibit viral replication in vitro in a calcium and dose-dependent manner. We show, for the first time, that SP-D and rhSP-D act as potent inhibitors of HIV-1 entry in to target cells and block the interaction between CD4 and gp120 in a dose-dependent manner. The rhSP-D-mediated inhibition of viral replication was examined using three clinical isolates of HIV-1 and three target cells: Jurkat T cells, U937 monocytic cells and PBMCs. HIV-1 induced cytokine storm in the three target cells was significantly suppressed by rhSP-D. Phosphorylation of key kinases p38, Erk1/2 and AKT, which contribute to HIV-1 induced immune activation, was significantly reduced in vitro in the presence of rhSP-D. Notably, anti-HIV-1 activity of rhSP-D was retained in the presence of biological fluids such as cervico-vaginal lavage and seminal plasma. Our study illustrates the multi-faceted role of human SPD against HIV-1 and potential of rhSP-D for immunotherapy to inhibit viral entry and immune activation in acute HIV infection. © 2014 Pandit et al.The work (Project no. 2011-16850) was supported by Medical Innovation Fund of Indian Council of Medical Research, New Delhi, India (www.icmr.nic.in/)

Protein–protein HADDocking using exclusively pseudocontact shifts

In order to enhance the structure determination process of macromolecular assemblies by NMR, we have implemented long-range pseudocontact shift (PCS) restraints into the data-driven protein docking package HADDOCK. We demonstrate the efficiency of the method on a synthetic, yet realistic case based on the lanthanide-labeled N-terminal ε domain of the E. coli DNA polymerase III (ε186) in complex with the HOT domain. Docking from the bound form of the two partners is swiftly executed (interface RMSDs < 1 Å) even with addition of very large amount of noise, while the conformational changes of the free form still present some challenges (interface RMSDs in a 3.1–3.9 Å range for the ten lowest energy complexes). Finally, using exclusively PCS as experimental information, we determine the structure of ε186 in complex with the HOT-homologue θ subunit of the E. coli DNA polymerase III

Developing a Comparative Docking Protocol for the Prediction of Peptide Selectivity Profiles: Investigation of Potassium Channel Toxins

During the development of selective peptides against highly homologous targets, a reliable tool is sought that can predict information on both mechanisms of binding and relative affinities. These tools must first be tested on known profiles before application on novel therapeutic candidates. We therefore present a comparative docking protocol in HADDOCK using critical motifs, and use it to “predict” the various selectivity profiles of several major αKTX scorpion toxin families versus Kv1.1, Kv1.2 and Kv1.3. By correlating results across toxins of similar profiles, a comprehensive set of functional residues can be identified. Reasonable models of channel-toxin interactions can be then drawn that are consistent with known affinity and mutagenesis. Without biological information on the interaction, HADDOCK reproduces mechanisms underlying the universal binding of αKTX-2 toxins, and Kv1.3 selectivity of αKTX-3 toxins. The addition of constraints encouraging the critical lysine insertion confirms these findings, and gives analogous explanations for other families, including models of partial pore-block in αKTX-6. While qualitatively informative, the HADDOCK scoring function is not yet sufficient for accurate affinity-ranking. False minima in low-affinity complexes often resemble true binding in high-affinity complexes, despite steric/conformational penalties apparent from visual inspection. This contamination significantly complicates energetic analysis, although it is usually possible to obtain correct ranking via careful interpretation of binding-well characteristics and elimination of false positives. Aside from adaptations to the broader potassium channel family, we suggest that this strategy of comparative docking can be extended to other channels of interest with known structure, especially in cases where a critical motif exists to improve docking effectiveness