Integration of small-angle X-ray scattering data into structural modeling of proteins and their assemblies

A major challenge in structural biology is to determine the configuration of domains and proteins in multidomain proteins and assemblies, respectively. All available data should be considered to maximize the accuracy and precision of these models. Small-angle X-ray scattering (SAXS) efficiently provides low-resolution experimental data about the shapes of proteins and their assemblies. Thus, we integrated SAXS profiles into our software for modeling proteins and their assemblies by satisfaction of spatial restraints. Specifically, we modeled the quaternary structures of multidomain proteins with structurally defined rigid domains as well as quaternary structures of binary complexes of structurally defined rigid proteins. In addition to SAXS profiles and the component structures, we used stereochemical restraints and an atomic distance-dependent statistical potential. The scoring function is optimized by a biased Monte Carlo protocol, including quasi-Newton and simulated annealing schemes. The final prediction corresponds to the best scoring solution in the largest cluster of many independently calculated solutions. To quantify how well the quaternary structures are determined based on their SAXS profiles, we used a benchmark of 12 simulated examples as well as an experimental SAXS profile of the homotetramer d-xylose isomerase. Optimization of the SAXS-dependent scoring function generally results in accurate models if sufficiently precise approximations for the constituent rigid bodies are available; otherwise, the best scoring models can have significant errors. Thus, SAXS profiles can play a useful role in the structural characterization of proteins and assemblies if they are combined with additional data and used judiciously. Our integration of a SAXS profile into modeling by satisfaction of spatial restraints will facilitate further integration of different kinds of data for structure determination of proteins and their assemblies.     

Quaternary structure of KATP channel SUR2A nucleotide binding domains resolved by synchrotron radiation X-ray scattering

Heterodimeric nucleotide binding domains NBD1/NBD2 distinguish the ATP-binding cassette protein SUR2A, a recognized regulatory subunit of cardiac ATP-sensitive K⁺ (K_ATP) channels. The tandem function of these core domains ensures metabolism-dependent gating of the Kir6.2 channel pore, yet their structural arrangement has not been resolved. Here, purified monodisperse and interference-free recombinant particles were subjected to synchrotron radiation small-angle X-ray scattering (SAXS) in solution. Intensity function analysis of SAXS profiles resolved NBD1 and NBD2 as octamers. Implemented by ab initio simulated annealing, shape determination prioritized an oblong envelope wrapping NBD1 and NBD2 with respective dimensions of 168 × 80 × 37 Å³ and 175 × 81 × 37 Å³ based on symmetry constraints, validated by atomic force microscopy. Docking crystal structure homology models against SAXS data reconstructed the NBD ensemble surrounding an inner cleft suitable for Kir6.2 insertion. Human heart disease-associated mutations introduced in silico verified the criticality of the mapped protein–protein interface. The resolved quaternary structure delineates thereby a macromolecular arrangement of K_ATP channel SUR2A regulatory domains.     

iATTRACT: Simultaneous global and local interface optimization for protein–protein docking refinement

Protein‐protein interactions are abundant in the cell but to date structural data for a large number of complexes is lacking. Computational docking methods can complement experiments by providing structural models of complexes based on structures of the individual partners. A major caveat for docking success is accounting for protein flexibility. Especially, interface residues undergo significant conformational changes upon binding. This limits the performance of docking methods that keep partner structures rigid or allow limited flexibility. A new docking refinement approach, iATTRACT, has been developed which combines simultaneous full interface flexibility and rigid body optimizations during docking energy minimization. It employs an atomistic molecular mechanics force field for intermolecular interface interactions and a structure‐based force field for intramolecular contributions. The approach was systematically evaluated on a large protein‐protein docking benchmark, starting from an enriched decoy set of rigidly docked protein–protein complexes deviating by up to 15 Å from the native structure at the interface. Large improvements in sampling and slight but significant improvements in scoring/discrimination of near native docking solutions were observed. Complexes with initial deviations at the interface of up to 5.5 Å were refined to significantly better agreement with the native structure. Improvements in the fraction of native contacts were especially favorable, yielding increases of up to 70%. Proteins 2015; 83:248–258. © 2014 Wiley Periodicals, Inc.     

ClusPro FMFT-SAXS: Ultra-fast Filtering Using Small-Angle X-ray Scattering Data in Protein Docking

We have recently demonstrated that incorporation of small-angle X-ray scattering (SAXS)-based filtering in our heavily used docking server ClusPro improves docking results. However, the filtering step is time consuming, since ≈ 10⁵ conformations have to be sequentially processed. At the same time, we have demonstrated the possibility of ultra-fast systematic energy evaluation for all rigid body orientations of two proteins, by sampling using Fast Manifold Fourier Transform (FMFT), if energies are represented as a combination of convolution-like expressions. Here we present a novel FMFT-based algorithm FMFT-SAXS for massive SAXS computation on multiple conformations of a protein complex. This algorithm exploits the convolutional form of SAXS calculation function. FMFT-SAXS allows computation of SAXS profiles for millions of conformations in a matter of minutes, providing an opportunity to explore the whole conformational space of two interacting proteins. We demonstrate the application of the new FMFT-SAXS approach to significantly speed up SAXS filtering step in our current docking protocol (1 to 2 orders of magnitude faster, running in several minutes on a modern 16-core CPU) without loss of accuracy. This is demonstrated on the benchmark set as well as on the experimental data. The new approach is available as a part of ClusPro server (https://beta.cluspro.org) and as an open source C library (https://bitbucket.org/abc-group/libfmftsaxs).     

Improving ranking of models for protein complexes with side chain modeling and atomic potentials

An atomically detailed potential for docking pairs of proteins is derived using mathematical programming. A refinement algorithm that builds atomically detailed models of the complex and combines coarse grained and atomic scoring is introduced. The refinement step consists of remodeling the interface side chains of the top scoring decoys from rigid docking followed by a short energy minimization. The refined models are then re‐ranked using a combination of coarse grained and atomic potentials. The docking algorithm including the refinement and re‐ranking, is compared favorably to other leading docking packages like ZDOCK, Cluspro, and PATCHDOCK, on the ZLAB 3.0 Benchmark and a test set of 30 novel complexes. A detailed analysis shows that coarse grained potentials perform better than atomic potentials for realistic unbound docking (where the exact structures of the individual bound proteins are unknown), probably because atomic potentials are more sensitive to local errors. Nevertheless, the atomic potential captures a different signal from the residue potential and as a result a combination of the two scores provides a significantly better prediction than each of the approaches alone. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

The utility of geometrical and chemical restraint information extracted from predicted ligand-binding sites in protein structure refinement

Exhaustive exploration of molecular interactions at the level of complete proteomes requires efficient and reliable computational approaches to protein function inference. Ligand docking and ranking techniques show considerable promise in their ability to quantify the interactions between proteins and small molecules. Despite the advances in the development of docking approaches and scoring functions, the genome-wide application of many ligand docking/screening algorithms is limited by the quality of the binding sites in theoretical receptor models constructed by protein structure prediction. In this study, we describe a new template-based method for the local refinement of ligand-binding regions in protein models using remotely related templates identified by threading. We designed a Support Vector Regression (SVR) model that selects correct binding site geometries in a large ensemble of multiple receptor conformations. The SVR model employs several scoring functions that impose geometrical restraints on the Cα positions, account for the specific chemical environment within a binding site and optimize the interactions with putative ligands. The SVR score is well correlated with the RMSD from the native structure; in 47% (70%) of the cases, the Pearson’s correlation coefficient is >0.5 (>0.3). When applied to weakly homologous models, the average heavy atom, local RMSD from the native structure of the top-ranked (best of top five) binding site geometries is 3.1 Å (2.9 Å) for roughly half of the targets; this represents a 0.1 (0.3) Å average improvement over the original predicted structure. Focusing on the subset of strongly conserved residues, the average heavy atom RMSD is 2.6 Å (2.3 Å). Furthermore, we estimate the upper bound of template-based binding site refinement using only weakly related proteins to be ～2.6 Å RMSD. This value also corresponds to the plasticity of the ligand-binding regions in distant homologues. The Binding Site Refinement (BSR) approach is available to the scientific community as a web server that can be accessed at http://cssb.biology.gatech.edu/bsr/.     

In silico and in vitro methods to identify ebola virus VP35-dsRNA inhibitors

Ebola virus continues to be problematic as sporadic outbreaks in Africa continue to arise, and as terrorist organizations have considered the virus for bioterrorism use. Several proteins within the virus have been targeted for antiviral chemotherapy, including VP35, a dsRNA binding protein that promotes viral replication, protects dsRNA from degradation, and prevents detection of the viral genome by immune complexes. To augment the scope of our antiviral research, we have now employed molecular modeling techniques to enrich the population of compounds for further testing in vitro. In the initial docking of a static VP35 structure with an 80,000 compound library, 40 compounds were selected, of which four compounds inhibited VP35 with IC₅₀ <200 μM, with the best compounds having an IC₅₀ of 20 μM. By superimposing 26 VP35 structures, we determined four aspartic acid residues were highly flexible and the docking was repeated under flexible parameters. Of 14 compounds chosen for testing, five compounds inhibited VP35 with IC₅₀ <200 μM and one compound with an IC₅₀ of 4 μM. These studies demonstrate the value of docking in silico for enriching compounds for testing in vitro, and specifically using multiple structures as a guide for detecting flexibility and provide a foundation for further development of small molecule inhibitors directed towards VP35.     

Pre-analytical factors in clinical proteomics investigations: Impact of ex vivo protein modifications for multiple sclerosis biomarker discovery

Serum proteome investigations have raised an incredible interest in the research of novel molecular biomarker, nevertheless few of the proposed evidences have been translated to the clinical practice. One of the limiting factors has been the lack of generally accepted guidelines for clinical proteomics studies and the lack of a robust analytical and pre-analytical ground for the proposed classification models. Pre-analytical issues may results in a deep impact for biomarker discovery campaign. In this study we present a systematic evaluation of sample storage and sampling conditions for clinical proteomics investigations.We have developed and validated a linear MALDI-TOF-MS protein profiling method to explore the low protein molecular weight region (5–20 kDa) of serum samples. Data normalization and processing was performed using optimise peak detection routine (LIMPIC) able to describe each group under investigation. Data were acquired either from healthy volunteers and from multiple sclerosis patients in order to highlight ex vivo protein profile alteration related to different physio-pathological conditions. Our data showed critical conditions for serum protein profiles depending on storage times and temperatures: 23 °C, 4 °C, ? 20 °C and ? 80 °C. We demonstrated that upon a ? 20 °C short term storage, characteristic degradation profiles are associated with different clinical groups. Protein signals were further identified after preparative HPLC separation by peptide sequencing on a nanoLC-Q-TOF TANDEM mass spectrometer. Apolipoprotein A-IV and complement C3 protein fragments, transthyretin and the oxidized isoforms in different apolipoprotein species represent the major molecular features of such a degradation pattern.     

3D-QSAR (CoMFA,CoMSIA), molecular docking and molecular dynamics simulations study of 6-aryl-5-cyano-pyrimidine derivatives to explore the structure requirements of LSD1 inhibitors

Recently, Histone Lysine Specific Demethylase 1 (LSD1) was regarded as a promising anticancer target for the novel drug discovery. And several small molecules as LSD1 inhibitors in different structures have been reported. In this work, we carried out a molecular modeling study on the 6-aryl-5-cyano-pyrimidine fragment LSD1 inhibitors using three-dimensional quantitative structure–activity relationship (3D-QSAR), molecular docking and molecular dynamics simulations. Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were used to generate 3D-QSAR models. The results show that the best CoMFA model has q² = 0.802, r²_ncv = 0.979, and the best CoMSIA model has q² = 0.799, r²_ncv = 0.982. The electrostatic, hydrophobic and H-bond donor fields play important roles in the models. Molecular docking studies predict the binding mode and the interactions between the ligand and the receptor protein. Molecular dynamics simulations results reveal that the complex of the ligand and the receptor protein are stable at 300 K. All the results can provide us more useful information for our further drug design.     

Calorimetric and spectroscopic investigations of the thermal denaturation of wild type nitrite reductase

Nitrite reductase (NiR) is a multicopper protein, with a trimeric structure containing two types of copper site: type 1 is present in each subunit whereas type 2 is localized at the subunits interface. The paper reports on the thermal behaviour of wild type NiR from Alcaligenes faecalis S-6. The temperature-induced changes of the copper centres are characterized by optical spectroscopy and electron paramagnetic resonance spectroscopy, and by establishing the thermal stability by differential scanning calorimetry. The calorimetric profile of the enzyme shows a single endothermic peak with maximum heat absorption at T_m ≈ 100 °C, revealing an exceptional thermal stability. The thermal transition is irreversible and the scan rate dependence of the calorimetric trace indicates that the denaturation of NiR is kinetically controlled. The divergence of the activation energy values determined by different methods is used as a criterion for the inapplicability of the one-step irreversible model. The best fit of the DSC profiles is obtained when the classical Lumry–Eyring model, N ? U ? F, is considered. The simulation results indicate that the irreversible step prevails on the reversible one. Moreover, it is found that the conformational changes within the type-1 copper environments precede the denaturation of the whole protein. No evidence of protein dissociation within the temperature range investigated was observed.     

IKKβ inhibitors identification part I: Homology model assisted structure based virtual screening

Control of NF-κB release through the inhibition of IKKβ has been identified as a potential target for the treatment of inflammatory and autoimmune diseases. We have employed structure based virtual screening scheme to identify lead like molecule from ChemDiv database. Homology models of IKKβ enzyme were developed based on the crystal structures of four kinases. The efficiency of the homology model has been validated at different levels. Docking of known inhibitors library revealed the possible binding mode of inhibitors. Besides, the docking sequence analyses results indicate the responsibility of Glu172 in selectivity. Structure based virtual screening of ChemDiv database has yielded 277 hits. Top scoring 75 compounds were selected and purchased for the IKKβ enzyme inhibition test. From the combined approach of virtual screening followed by biological screening, we have identified six novel compounds that can work against IKKβ, in which 1 compound had highest inhibition rate 82.09% at 10 μM and IC₅₀ 1.76 μM and 5 compounds had 25.35–48.80% inhibition.     

Biased coarse-grained molecular dynamics simulation approach for flexible fitting of X-ray structure into cryo electron microscopy maps

Several approaches have been introduced to interpret, in terms of high-resolution structure, low-resolution structural data as obtained from cryo-EM. As conformational changes are often observed in biological molecules, these techniques need to take into account the flexibility of proteins. Flexibility has been described in terms of movement between rigid domains and between rigid secondary structure elements, which present some limitations for studying dynamical properties. Normal mode analysis has also been used, but is limited to medium resolution data. All-atom molecular dynamics fitting techniques are more appropriate to fit structures into higher-resolution data as full protein flexibility is considered, but are cumbersome in terms of computational time. Here, we introduce a coarse-grained approach; a Go-model was used to represent biological molecules, combined with biased molecular dynamics to reproduce accurately conformational transitions. Illustrative examples on simulated data are shown. Accurate fittings can be obtained for resolution ranging from 5 to 20 Å. The approach was also tested on experimental data of Elongation Factor G and Escherichia coli RNA polymerase, where its validity is compared to previous models obtained from different techniques. This comparison demonstrates that quantitative flexible techniques, as opposed to manual docking, need to be considered to interpret low-resolution data.     

BCL::SAXS: GPU accelerated Debye method for computation of small angle X‐ray scattering profiles

Small angle X‐ray scattering (SAXS) is an experimental technique used for structural characterization of macromolecules in solution. Here, we introduce BCL::SAXS—an algorithm designed to replicate SAXS profiles from rigid protein models at different levels of detail. We first show our derivation of BCL::SAXS and compare our results with the experimental scattering profile of hen egg white lysozyme. Using this protein we show how to generate SAXS profiles representing: (1) complete models, (2) models with approximated side chain coordinates, and (3) models with approximated side chain and loop region coordinates. We evaluated the ability of SAXS profiles to identify a correct protein topology from a non‐redundant benchmark set of proteins. We find that complete SAXS profiles can be used to identify the correct protein by receiver operating characteristic (ROC) analysis with an area under the curve (AUC) > 99%. We show how our approximation of loop coordinates between secondary structure elements improves protein recognition by SAχS for protein models without loop regions and side chains. Agreement with SAXS data is a necessary but not sufficient condition for structure determination. We conclude that experimental SAXS data can be used as a filter to exclude protein models with large structural differences from the native. Proteins 2015; 83:1500–1512. © 2015 Wiley Periodicals, Inc.     

Improving protein template recognition by using small-angle x-ray scattering profiles

Small-angle x-ray scattering (SAXS) is able to extract low-resolution protein shape information without requiring a specific crystal formation. However, it has found little use in atomic-level protein structure determination due to the uncertainty of residue-level structural assignment. We developed a new algorithm, SAXSTER, to couple the raw SAXS data with protein-fold-recognition algorithms and thus improve template-based protein-structure predictions. We designed nine different matching scoring functions of template and experimental SAXS profiles. The logarithm of the integrated correlation score showed the best template recognition ability and had the highest correlation with the true template modeling (TM)-score of the target structures. We tested the method in large-scale protein-fold-recognition experiments and achieved significant improvements in prioritizing the best template structures. When SAXSTER was applied to the proteins of asymmetric SAXS profile distributions, the average TM-score of the top-ranking templates increased by 18% after homologous templates were excluded, which corresponds to a p-value < 10⁻⁹ in Student's t-test. These data demonstrate a promising use of SAXS data to facilitate computational protein structure modeling, which is expected to work most efficiently for proteins of irregular global shape and/or multiple-domain protein complexes.     

Finding rigid bodies in protein structures: Application to flexible fitting into cryoEM maps

We present RIBFIND, a method for detecting flexibility in protein structures via the clustering of secondary structural elements (SSEs) into rigid bodies. To test the usefulness of the method in refining atomic structures within cryoEM density we incorporated it into our flexible fitting protocol (Flex-EM). Our benchmark includes 13 pairs of protein structures in two conformations each, one of which is represented by a corresponding cryoEM map. Refining the structures in simulated and experimental maps at the 5–15 Å resolution range using rigid bodies identified by RIBFIND shows a significant improvement over using individual SSEs as rigid bodies. For the 15 Å resolution simulated maps, using RIBFIND-based rigid bodies improves the initial fits by 40.64% on average, as compared to 26.52% when using individual SSEs. Furthermore, for some test cases we show that at the sub-nanometer resolution range the fits can be further improved by applying a two-stage refinement protocol (using RIBFIND-based refinement followed by an SSE-based refinement). The method is stand-alone and could serve as a general interactive tool for guiding flexible fitting into EM maps.     

Structures of xyloglucans in primary cell walls of gymnosperms,monilophytes (ferns sensu lato) and lycophytes

Little is known about the structures of the xyloglucans in the primary cell walls of vascular plants (tracheophytes) other than angiosperms. Xyloglucan structures were examined in 13 species of gymnosperms, 13 species of monilophytes (ferns sensu lato), and two species of lycophytes. Wall preparations were obtained, extracted with 6 M sodium hydroxide, and the extracts treated with a xyloglucan-specific endo-(1 → 4)-β-glucanase preparation. The oligosaccharides released were analysed by matrix-assisted laser-desorption ionisation time-of-flight mass spectrometry and by high-performance anion-exchange chromatography. The xyloglucan oligosaccharide profiles from the gymnosperm walls were similar to those from the walls of most eudicotyledons and non-commelinid monocotyledons, indicating that the xyloglucans were fucogalactoxyloglucans, containing the fucosylated units XXFG and XLFG. The xyloglucan oligosaccharide profiles for six of the monilophyte species were similar to those of the gymnosperms, indicating they were also fucogalactoxyloglucans. Phylogenetically, these monilophyte species were from both basal and more derived orders. However, the profiles for the other monilophyte species showed various significant differences, including additional oligosaccharides. In three of the species, these additional oligosaccharides contained arabinosyl residues which were most abundant in the profile of Equisetum hyemale. The two species of lycophytes examined, Selaginella kraussiana and Lycopodium cernuum, had quite different xyloglucan oligosaccharide profiles, but neither were fucogalactoxyloglucans. The S. kraussiana profile had abundant oligosaccharides containing arabinosyl residues. The L. cernuum profile indicated the xyloglucan had a very complex structure.     

Origin of defects in assembled supramolecular structures of sweet potato starches with different amylopectin chain-length distribution

A combined approach of fluorophore-assisted capillary electrophoresis (FACEL), high-sensitivity differential scanning calorimetry (DSC), wide-angle X-ray scattering (WAXS), small-angle X-ray scattering (SAXS), and light (LM) and scanning electron microscopy (SEM) was applied to study the effects of changes in amylopectin chain-length distribution on the assembly structures of sweet potato starches with similar amylose levels. It was shown that unlike ordinary sweet potato starch, starch extracted from Quick Sweet cultivar of sweet potato had anomalous high level of amylopectin chains with a degree of polymerization (DP) 6–12. Joint analysis of the obtained data revealed that amylopectin chains with DP 10–24 are, apparently, the dominant material for the formation of supramolecular structures in starch granules. In contrast, amylopectin chains with DP < 10 facilitated the formation of defects within crystalline lamellae. An increase in relative content of amylopectin chains with DP < 10 is accompanied by the correlated structural alterations manifested at all levels of starch granule organization (crystalline lamellae, amylopectin clusters, semi-crystalline growth rings, and granule morphology). Thus, the short amylopectin chains with DP < 10 were considered as an origin of the defectiveness in starch supramolecular structures.     

Interaction of prodigiosin with HSA and β-Lg: Spectroscopic and molecular docking studies

Human serum albumin (HSA) and bovine β-lactoglobulin (β-Lg) are both introduced as blood and oral carrier scaffolds with high affinity for a wide range of pharmaceutical compounds. Prodigiosin, a natural three pyrrolic compound produced by Serratia marcescens, exhibits many pharmaceutical properties associated with health benefits. In the present study, the interaction of prodigiosin with HSA and β-Lg was investigated using fluorescence spectroscopy, circular dichroism (CD) and computational docking. Prodigiosin interacts with the Sudlow’s site I of HSA and the calyx of β-Lg with association constant of 4.41 × 10⁴ and 1.99 × 10⁴ M⁻¹ to form 1:1 and 2:3 complexes at 300 K, respectively. The results indicated that binding of prodigiosin to HSA and β-Lg caused strong fluorescence quenching of both proteins through static quenching mechanism. Electrostatic and hydrophobic interactions are the major forces in the stability of PG–HSA complex with enthalpy- and entropy-driving mode, although the formation of prodigiosin–β-Lg complex is entropy-driven hydrophobic associations. CD spectra showed slight conformational changes in both proteins due to the binding of prodigiosin. Moreover, the ligand displacement assay, pH-dependent interaction and protein–ligand docking study confirmed that the prodigiosin binds to residues located in the subdomain IIA and IIIA of HSA and central calyx of β-Lg.     

New coumarin derivatives: Design,synthesis and use as inhibitors of hMAO

A series new 2H-chromene-3-carboxamides (4a–4i) and S-2H-chromene-3-carbothioates (5j–5t) were synthesized and evaluated as monoamine oxidase A and B inhibitors. Among them, compound 5k (IC₅₀ = 0.21 μM, IC_{50 iproniazid} = 7.65 μM) showed the most activity and higher MAO-B selectivity (189.2-fold vs 1-fold) with respect to the MAO-A isoform. The need to clarify at a 3D level some important molecular aspects of discussed SAR, we undertaked a number of docking simulations to better assess. The steric effect was analyzed interms of both posing and scoring by investigating the nature of the binding interactions. The docking results of active compound 5k with hMAO-B complex indicated that conserved residue ILE 199 was important for ligand binding via Sigma–Pi interaction.