Accurate SAXS Profile Computation and its Assessment by Contrast Variation Experiments

AbstractA major challenge in structural biology is to characterize structures of proteins and their assemblies in solution. At low resolution, such a characterization may be achieved by small angle x-ray scattering (SAXS). Because SAXS analyses often require comparing profiles calculated from many atomic models against those determined by experiment, rapid and accurate profile computation from molecular structures is needed. We developed fast open-source x-ray scattering (FoXS) for profile computation. To match the experimental profile within the experimental noise, FoXS explicitly computes all interatomic distances and implicitly models the first hydration layer of the molecule. For assessing the accuracy of the modeled hydration layer, we performed contrast variation experiments for glucose isomerase and lysozyme, and found that FoXS can accurately represent density changes of this layer. The hydration layer model was also compared with a SAXS profile calculated for the explicit water molecules in the high-resolution structures of glucose isomerase and lysozyme. We tested FoXS on eleven protein, one DNA, and two RNA structures, revealing superior accuracy and speed versus CRYSOL, AquaSAXS, the Zernike polynomials-based method, and Fast-SAXS-pro. In addition, we demonstrated a significant correlation of the SAXS score with the accuracy of a structural model. Moreover, FoXS utility for analyzing heterogeneous samples was demonstrated for intrinsically flexible XLF-XRCC4 filaments and Ligase III-DNA complex. FoXS is extensively used as a standalone web server as a component of integrative structure determination by programs IMP, Chimera, and BILBOMD, as well as in other applications that require rapidly and accurately calculated SAXS profiles

PatchDock and SymmDock: servers for rigid and symmetric docking

Here, we describe two freely available web servers for molecular docking. The PatchDock method performs structure prediction of protein–protein and protein–small molecule complexes. The SymmDock method predicts the structure of a homomultimer with cyclic symmetry given the structure of the monomeric unit. The inputs to the servers are either protein PDB codes or uploaded protein structures. The services are available at . The methods behind the servers are very efficient, allowing large-scale docking experiments

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/

The X-ray structure of human calbindin-D28K: an improved model

Calbindin-D28K is a widely expressed calcium-buffering cytoplasmic protein that is involved in many physiological processes. It has been shown to interact with other proteins, suggesting a role as a calcium sensor. Many of the targets of calbindin-D28K are of therapeutic interest: for example, inositol monophosphatase, the putative target of lithium therapy in bipolar disorder. Presented here is the first crystal structure of human calbindin-D28K. There are significant deviations in the tertiary structure when compared with the NMR structure of rat calbindin-D28K (PDB entry 2g9b), despite 98% sequence identity. Smallangle X-ray scattering (SAXS) indicates that the crystal structure better predicts the properties of calbindin-D28K in solution compared with the NMR structure. Here, the first direct visualization of the calcium-binding properties of calbindinD28K is presented. Four of the six EF-hands that make up the secondary structure of the protein contain a calcium-binding site. Two distinct conformations of the N-terminal EF-hand calcium-binding site were identified using long-wavelength calcium single-wavelength anomalous dispersion (SAD). This flexible region has previously been recognized as a protein–protein interaction interface. SAXS data collected in both the presence and absence of calcium indicate that there are no large structural differences in the globular structure of calbindin-D28K between the calcium-loaded and unloaded proteins

A real-time proximity querying algorithm for haptic-based molecular docking

Intermolecular binding underlies every metabolic and regulatory processes of the cell, and the therapeutic and pharmacological properties of drugs. Molecular docking systems model and simulate these interactions in silico and allow us to study the binding process. Haptic-based docking provides an immersive virtual docking environment where the user can interact with and guide the molecules to their binding pose. Moreover, it allows human perception, intuition and knowledge to assist and accelerate the docking process, and reduces incorrect binding poses. Crucial for interactive docking is the real-time calculation of interaction forces. For smooth and accurate haptic exploration and manipulation, force-feedback cues have to be updated at a rate of 1 kHz. Hence, force calculations must be performed within 1ms. To achieve this, modern haptic-based docking approaches often utilize pre-computed force grids and linear interpolation. However, such grids are time-consuming to pre-compute (especially for large molecules), memory hungry, can induce rough force transitions at cell boundaries and cannot be applied to flexible docking. Here we propose an efficient proximity querying method for computing intermolecular forces in real time. Our motivation is the eventual development of a haptic-based docking solution that can model molecular flexibility. Uniquely in a haptics application we use octrees to decompose the 3D search space in order to identify the set of interacting atoms within a cut-off distance. Force calculations are then performed on this set in real time. The implementation constructs the trees dynamically, and computes the interaction forces of large molecular structures (i.e. consisting of thousands of atoms) within haptic refresh rates. We have implemented this method in an immersive, haptic-based, rigid-body, molecular docking application called Haptimol_RD. The user can use the haptic device to orientate the molecules in space, sense the interaction forces on the device, and guide the molecules to their binding pose. Haptimol_RD is designed to run on consumer level hardware, i.e. there is no need for specialized/proprietary hardware

Identification of the N-terminal Peptide Binding Site of Glucose-regulated Protein 94

Because the stress protein GRP94 can augment presentation of peptides to T cells, it is important to define how it, as well as all other HSP90 family members, binds peptides. Having previously shown that the N-terminal half of GRP94 can account for the peptide binding activity of the full-length protein, we now locate this binding site by testing predictions of a molecular docking model. The best predicted site was on the opposite face of the β sheet from the pan-HSP90 radicicol-binding pocket, in close proximity to a deep hydrophobic pocket. The peptide and radicicol-binding sites are distinct, as shown by the ability of a radicicol-refractive mutant to bind peptide. When the fluorophore acrylodan is attached to Cys(117)within the hydrophobic pocket, its fluorescence is reduced upon peptide binding, consistent with proximity of the two ligands. Substitution of His(125), which contacts the bound peptide, compromises peptide-binding activity. We conclude that peptide binds to the concave face of the β sheet of the N-terminal domain, where binding is regulated during the action cycle of the chaperone

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/