A method for integrative structure determination of protein-protein complexes

Evidence for a compact σ70 conformation in vitro and in vivo

The initiation of transcription in Escherichia coli (E. coli) is facilitated by promoter specificity factors, also known as σ factors, which may bind a promoter only as part of a complex with RNA polymerase (RNAP). By performing in vitro cross-linking mass spectrometry (CL-MS) of apo-σ70, we reveal structural features suggesting a compact conformation compared to the known RNAP-bound extended conformation. Then, we validate the existence of the compact conformation using in vivo CL-MS by identifying cross-links similar to those found in vitro, which deviate from the extended conformation only during the stationary phase of bacterial growth. Conclusively, we provide information in support of a compact conformation of apo-σ70 that exists in live cells, which might represent a transcriptionally inactive form that can be activated upon binding to RNAP.

A deep learning model for predicting optimal distance range in crosslinking mass spectrometry data

Macromolecular assemblies play an important role in all cellular processes. While there has recently been significant progress in protein structure prediction based on deep learning, large protein complexes cannot be predicted with these approaches. The integrative structure modeling approach characterizes multi-subunit complexes by computational integration of data from fast and accessible experimental techniques. Crosslinking mass spectrometry is one such technique that provides spatial information about the proximity of crosslinked residues. One of the challenges in interpreting crosslinking datasets is designing a scoring function that, given a structure, can quantify how well it fits the data. Most approaches set an upper bound on the distance between Cα atoms of crosslinked residues and calculate a fraction of satisfied crosslinks. However, the distance spanned by the crosslinker greatly depends on the neighborhood of the crosslinked residues. Here, we design a deep learning model for predicting the optimal distance range for a crosslinked residue pair based on the structures of their neighborhoods. We find that our model can predict the distance range with the area under the receiver-operator curve of 0.86 and 0.7 for intra- and inter-protein crosslinks, respectively. Our deep scoring function can be used in a range of structure modeling applications.

Multi-state modeling of antibody-antigen complexes with SAXS profiles and deep-learning models

Antibodies are an established class of human therapeutics. Epitope characterization is an important part of therapeutic antibody discovery. However, structural characterization of antibody-antigen complexes remains challenging. On the one hand, X-ray crystallography or cryo-electron microscopy provide atomic resolution characterization of the epitope, but the data collection process is typically long and the success rate is low. On the other hand, computational methods for modeling antibody-antigen structures from the individual components frequently suffer from a high false positive rate, rarely resulting in a unique solution. Recent deep learning models for structure prediction are also successful in predicting protein-protein complexes. However, they do not perform well for antibody-antigen complexes. Small Angle X-ray Scattering (SAXS) is a reliable technique for rapid structural characterization of protein samples in solution albeit at low resolution. Here, we present an integrative approach for modeling antigen-antibody complexes using the antibody sequence, antigen structure, and experimentally determined SAXS profiles of the antibody, antigen, and the complex. The method models antibody structures using a novel deep-learning approach, NanoNet. The structures of the antibodies and antigens are represented using multiple 3D conformations to account for compositional and conformational heterogeneity of the protein samples that are used to collect the SAXS data. The complexes are predicted by integrating the SAXS profiles with scoring functions for protein-protein interfaces that are based on statistical potentials and antibody-specific deep-learning models. We validated the method via application to four Fab:EGFR and one Fab:PCSK9 antibody:antigen complexes with experimentally available SAXS datasets. The integrative approach returns accurate predictions (interface RMSD<4Å) in the top five predictions for four out of five complexes (respective interface RMSD values of 1.95, 2.18, 2.66 and 3.87Å), providing support for the utility of such a computational pipeline for epitope characterization during therapeutic antibody discovery.

Integrative structure determination reveals functional global flexibility for an ultra-multimodular arabinanase

AbnA is an extracellular GH43 α-L-arabinanase from Geobacillus stearothermophilus, a key bacterial enzyme in the degradation and utilization of arabinan. We present herein its full-length crystal structure, revealing the only ultra-multimodular architecture and the largest structure to be reported so far within the GH43 family. Additionally, the structure of AbnA appears to contain two domains belonging to new uncharacterized carbohydrate-binding module (CBM) families. Three crystallographic conformational states are determined for AbnA, and this conformational flexibility is thoroughly investigated further using the "integrative structure determination" approach, integrating molecular dynamics, metadynamics, normal mode analysis, small angle X-ray scattering, dynamic light scattering, cross-linking, and kinetic experiments to reveal large functional conformational changes for AbnA, involving up to ~100 Å movement in the relative positions of its domains. The integrative structure determination approach demonstrated here may apply also to the conformational study of other ultra-multimodular proteins of diverse functions and structures.

New Horizons in Structural Biology of Membrane Proteins: Experimental Evaluation of the Role of Conformational Dynamics and Intrinsic Flexibility

A plethora of membrane proteins are found along the cell surface and on the convoluted labyrinth of membranes surrounding organelles. Since the advent of various structural biology techniques, a sub-population of these proteins has become accessible to investigation at near-atomic resolutions. The predominant bona fide methods for structure solution, X-ray crystallography and cryo-EM, provide high resolution in three-dimensional space at the cost of neglecting protein motions through time. Though structures provide various rigid snapshots, only an amorphous mechanistic understanding can be inferred from interpolations between these different static states. In this review, we discuss various techniques that have been utilized in observing dynamic conformational intermediaries that remain elusive from rigid structures. More specifically we discuss the application of structural techniques such as NMR, cryo-EM and X-ray crystallography in studying protein dynamics along with complementation by conformational trapping by specific binders such as antibodies. We finally showcase the strength of various biophysical techniques including FRET, EPR and computational approaches using a multitude of succinct examples from GPCRs, transporters and ion channels.

Antihyperlipidemic effect of selected pyrimidine derivatives mediated through multiple pathways

Hyperlipidemia is worth-mentioning risk factor in quickly expanding atherosclerosis, myocardial infarction, and stroke. This study attempted to determine effectiveness of selected pyrimidine derivatives: 5-(3-Hydroxybenzylidene)-2, 4, 6(1H, 3H, 5H)-pyrimidinetrione (SR-5), 5-(4-Hydroxybenzylidene)-2, 4, 6(1H, 3H, 5H)-pyrimidinetrione (SR-8), 5-(3-Chlorobenzylidene)-2, 4, 6(1H, 3H, 5H)-pyrimidinetrione (SR-9), and 5-(4-Chlorobenzylidene)-2, 4, 6(1H, 3H, 5H)-pyrimidinetrione (SR-10) against hyperlipidemia. In silico results revealed that SR-5, SR-8, SR-9, and SR-10 exhibited high affinity with 3-hydroxy-3-methylglutaryl coenzyme A (HMGCoA) possessing binding energy values of -8.2, -8.4, -8.6, and -9.5 Kcal/mol, respectively, and moderate (<-8 Kcal/mol) against other selected targets. In vivo findings showed that test drugs (25 and 50 mg/Kg) significantly decreased HFD rat total cholesterol, triglycerides, low-density lipoprotein, very-low-density lipoprotein, atherogenic index, coronary risk index, alkaline phosphatase, aspartate transaminase, alanine transaminase, and bilirubin and increased high-density lipoprotein (p < 0.05, p < 0.01, p < 0.001 vs HFD group). In animal liver tissues, SR-5, SR-8, SR-9, and SR-10 inhibited HMGCoA reductase enzyme, enhanced glutathione-s-transferase, reduced glutathione, catalase levels, improved cellular architecture in histopathological examination, and decreased expression of inflammatory markers: cyclo-oxygenase 2, tumor necrosis factor alpha, phosphorylated c-Jun N-terminal kinase, and phosphorylated-nuclear factor kappa B, evidenced in immunohistochemistry and enzyme-linked immunosorbent assay molecular investigations. This study indicates that SR-5, SR-8, SR-9, and SR-10 exhibit antihyperlipidemic action, mediated possibly through HMGCoA inhibition, hepatoprotection, antioxidant, and anti-inflammatory pathways.

Integrative proteomics identifies thousands of distinct, multi-epitope, and high-affinity nanobodies

The antibody immune response is essential for the survival of mammals. However, we still lack a systematic understanding of the antibody repertoire. Here, we developed a proteomic strategy to survey, at an unprecedented scale, the landscape of antigen-engaged, circulating camelid heavy-chain antibodies, whose minimal binding fragments are called VHH antibodies or nanobodies. The sensitivity and robustness of this approach were validated with three antigens spanning orders of magnitude in immune responses; thousands of distinct, high-affinity nanobody families were reliably identified and quantified. Using high-throughput structural modeling, cross-linking mass spectrometry, mutagenesis, and deep learning, we mapped and analyzed the epitopes of >100,000 antigen-nanobody complexes. Our results revealed a surprising diversity of ultrahigh-affinity camelid nanobodies for specific antigen binding on various dominant epitope clusters. Nanobodies utilize both shape and charge complementarity to enable highly selective antigen binding. Interestingly, we found that nanobody-antigen binding can mimic conserved intracellular protein-protein interactions. A record of this paper's Transparent Peer Review process is included in the Supplemental information.

From integrative structural biology to cell biology

Integrative modeling is an increasingly important tool in structural biology, providing structures by combining data from varied experimental methods and prior information. As a result, molecular architectures of large, heterogeneous, and dynamic systems, such as the ∼52-MDa Nuclear Pore Complex, can be mapped with useful accuracy, precision, and completeness. Key challenges in improving integrative modeling include expanding model representations, increasing the variety of input data and prior information, quantifying a match between input information and a model in a Bayesian fashion, inventing more efficient structural sampling, as well as developing better model validation, analysis, and visualization. In addition, two community-level challenges in integrative modeling are being addressed under the auspices of the Worldwide Protein Data Bank (wwPDB). First, the impact of integrative structures is maximized by PDB-Development, a prototype wwPDB repository for archiving, validating, visualizing, and disseminating integrative structures. Second, the scope of structural biology is expanded by linking the wwPDB resource for integrative structures with archives of data that have not been generally used for structure determination but are increasingly important for computing integrative structures, such as data from various types of mass spectrometry, spectroscopy, optical microscopy, proteomics, and genetics. To address the largest of modeling problems, a type of integrative modeling called metamodeling is being developed; metamodeling combines different types of input models as opposed to different types of data to compute an output model. Collectively, these developments will facilitate the structural biology mindset in cell biology and underpin spatiotemporal mapping of the entire cell.

Versatile and multivalent nanobodies efficiently neutralize SARS-CoV-2

Cost-effective, efficacious therapeutics are urgently needed to combat the COVID-19 pandemic. In this study, we used camelid immunization and proteomics to identify a large repertoire of highly potent neutralizing nanobodies (Nbs) to the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein receptor binding domain (RBD). We discovered Nbs with picomolar to femtomolar affinities that inhibit viral infection at concentrations below the nanograms-per-milliliter level, and we determined a structure of one of the most potent Nbs in complex with the RBD. Structural proteomics and integrative modeling revealed multiple distinct and nonoverlapping epitopes and indicated an array of potential neutralization mechanisms. We bioengineered multivalent Nb constructs that achieved ultrahigh neutralization potency (half-maximal inhibitory concentration as low as 0.058 ng/ml) and may prevent mutational escape. These thermostable Nbs can be rapidly produced in bulk from microbes and resist lyophilization and aerosolization.

Genetic interaction mapping informs integrative structure determination of protein complexes

Determining structures of protein complexes is crucial for understanding cellular functions. Here, we describe an integrative structure determination approach that relies on in vivo measurements of genetic interactions. We construct phenotypic profiles for point mutations crossed against gene deletions or exposed to environmental perturbations, followed by converting similarities between two profiles into an upper bound on the distance between the mutated residues. We determine the structure of the yeast histone H3-H4 complex based on ~500,000 genetic interactions of 350 mutants. We then apply the method to subunits Rpb1-Rpb2 of yeast RNA polymerase II and subunits RpoB-RpoC of bacterial RNA polymerase. The accuracy is comparable to that based on chemical cross-links; using restraints from both genetic interactions and cross-links further improves model accuracy and precision. The approach provides an efficient means to augment integrative structure determination with in vivo observations.

Combining Experimental Data and Computational Methods for the Non-Computer Specialist

Experimental methods are indispensable for the study of the function of biological macromolecules, not just as static structures, but as dynamic systems that change conformation, bind partners, perform reactions, and respond to different stimulus. However, providing a detailed structural interpretation of the results is often a very challenging task. While experimental and computational methods are often considered as two different and separate approaches, the power and utility of combining both is undeniable. The integration of the experimental data with computational techniques can assist and enrich the interpretation, providing new detailed molecular understanding of the systems. Here, we briefly describe the basic principles of how experimental data can be combined with computational methods to obtain insights into the molecular mechanism and expand the interpretation through the generation of detailed models.

Versatile, Multivalent Nanobody Cocktails Efficiently Neutralize SARS-CoV-2

The outbreak of COVID-19 has severely impacted global health and the economy. Cost-effective, highly efficacious therapeutics are urgently needed. Here, we used camelid immunization and proteomics to identify a large repertoire of highly potent neutralizing nanobodies (Nbs) to the SARS-CoV-2 spike (S) protein receptor-binding domain (RBD). We discovered multiple elite Nbs with picomolar to femtomolar affinities that inhibit viral infection at sub-ng/ml concentration, more potent than some of the best human neutralizing antibodies. We determined a crystal structure of such an elite neutralizing Nb in complex with RBD. Structural proteomics and integrative modeling revealed multiple distinct and non-overlapping epitopes and indicated an array of potential neutralization mechanisms. Structural characterization facilitated the bioengineering of novel multivalent Nb constructs into multi-epitope cocktails that achieved ultrahigh neutralization potency (IC50s as low as 0.058 ng/ml) and may prevent mutational escape. These thermostable Nbs can be rapidly produced in bulk from microbes and resist lyophilization, and aerosolization. These promising agents are readily translated into efficient, cost-effective, and convenient therapeutics to help end this once-in-a-century health crisis.

Structural insights on binding mechanism of CAD complexes (CPSase, ATCase and DHOase)

Pyrimidine biosynthetic pathway enzymes constitute an important target for the development of antitumor drugs. To understand the role of binding mechanisms underlying the inborn errors of pyrimidine biosynthetic pathway, structure and function of enzymes have been analyzed. Pyrimidine biosynthetic pathway is initiated by CAD enzymes that harbor the first three enzymatic activities facilitated by Carbamoyl Phosphate Synthetase (CPSase), Aspartate Transcarbamoylase (ATCase) and Dihydroorotase (DHOase). While being an attractive therapeutic target, the lack of data driven us to study the CPSase (CarA and CarB) and its mode of binding to ATCase and DHOase which are the major limitation for its structural optimization. Understanding the binding mode of CPSase, ATCase and DHOase could help to identify the potential interface hotspot residues that favor the mechanism behind it. The mechanistic insight into the CAD complexes were achieved through Molecular modeling, Protein-Protein docking, Alanine scanning and Molecular dynamics (MD) Studies. The hotspot residues present in the CarB region of carboxy phosphate and carbamoyl phosphate synthetic domains are responsible for the assembly of CAD (CPSase-ATCase-DHOase) complexes. Overall analysis suggests that the identified hotspot residues were confirmed by alanine scanning and important for the regulation of pyrimidine biosynthesis. MD simulations analysis provided the prolonged stability of the interacting complexes. The present study reveals the novel hotspot residues such as Glu134, Glu147, Glu154, Asp266, Lys269, Glu274, Asp333, Trp459, Asp526, Asp528, Glu533, Glu544, Glu546, Glu800, Val855, Asp877, Tyr884 and Gln919 which could be targeted for structure-based inhibitor design to potentiate the CAD mediated regulation of aggressive tumors.Communicated by Ramaswamy H. Sarma.

Carbohydrate-Binding Capability and Functional Conformational Changes of AbnE, an Arabino-oligosaccharide Binding Protein

ABC importers are membrane proteins responsible for the transport of nutrients into the cells of prokaryotes. Although the structures of ABC importers vary, all contain four conserved domains: two nucleotide-binding domains (NBDs), which bind and hydrolyze ATP, and two transmembrane domains (TMDs), which help translocate the substrate. ABC importers are also dependent on an additional protein component, a high-affinity substrate-binding protein (SBP) that specifically binds the target ligand for delivery to the appropriate ABC transporter. AbnE is a SBP belonging to the ABC importer for arabino-oligosaccharides in the Gram-positive thermophilic bacterium Geobacillus stearothermophilus. Using isothermal titration calorimetry (ITC), purified AbnE was shown to bind medium-sized arabino-oligosaccharides, in the range of arabino-triose (A3) to arabino-octaose (A8), all with K_d values in the nanomolar range. We describe herein the 3D structure of AbnE in its closed conformation in complex with a wide range of arabino-oligosaccharide substrates (A2-A8). These structures provide the basis for the detailed structural analysis of the AbnE-sugar complexes, and together with complementary quantum chemical calculations, site-specific mutagenesis, and isothermal titration calorimetry (ITC) experiments, provide detailed insights into the AbnE-substrate interactions involved. Small-angle X-ray scattering (SAXS) experiments and normal mode analysis (NMA) are used to study the conformational changes of AbnE, and these data, taken together, suggest clues regarding its binding mode to the full ABC importer.

Neighborhood Preference of Amino Acids in Protein Structures and its Applications in Protein Structure Assessment

Amino acids form protein 3D structures in unique manners such that the folded structure is stable and functional under physiological conditions. Non-specific and non-covalent interactions between amino acids exhibit neighborhood preferences. Based on structural information from the protein data bank, a statistical energy function was derived to quantify amino acid neighborhood preferences. The neighborhood of one amino acid is defined by its contacting residues, and the energy function is determined by the neighboring residue types and relative positions. The neighborhood preference of amino acids was exploited to facilitate structural quality assessment, which was implemented in the neighborhood preference program NEPRE. The source codes are available via https://github.com/LiuLab-CSRC/NePre.

Collision Cross Section Calculations Using HPCCS

A technical overview of the High Performance Collision Cross Section (HPCCS) software for accurate and efficient calculations of collision cross sections for molecular ions ranging from small organic molecules to large protein complexes is presented. The program uses helium or nitrogen as buffer gas with considerable gains in computer time compared to publicly available codes under the Trajectory Method approximation. HPCCS is freely available under the Academic Use License at https://github.com/cepid-cces/hpccs .

Proper evaluation of chemical cross-linking-based spatial restraints improves the precision of modeling homo-oligomeric protein complexes

Background

The function of oligomeric proteins is inherently linked to their quaternary structure. In the absence of high-resolution data, low-resolution information in the form of spatial restraints can significantly contribute to the precision and accuracy of structural models obtained using computational approaches. To obtain such restraints, chemical cross-linking coupled with mass spectrometry (XL-MS) is commonly used. However, the use of XL-MS in the modeling of protein complexes comprised of identical subunits (homo-oligomers) is often hindered by the inherent ambiguity of intra- and inter-subunit connection assignment.

Results

We present a comprehensive evaluation of (1) different methods for inter-residue distance calculations, and (2) different approaches for the scoring of spatial restraints. Our results show that using Solvent Accessible Surface distances (SASDs) instead of Euclidean distances (EUCs) greatly reduces the assignation ambiguity and delivers better modeling precision. Furthermore, ambiguous connections should be considered as inter-subunit only when the intra-subunit alternative exceeds the distance threshold. Modeling performance can also be improved if symmetry, characteristic for most homo-oligomers, is explicitly defined in the scoring function.

Conclusions

Our findings provide guidelines for proper evaluation of chemical cross-linking-based spatial restraints in modeling homo-oligomeric protein complexes, which could facilitate structural characterization of this important group of proteins.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-3032-x) contains supplementary material, which is available to authorized users.

Archiving and disseminating integrative structure models

Limitations in the applicability, accuracy, and precision of individual structure characterization methods can sometimes be overcome via an integrative modeling approach that relies on information from all available sources, including all available experimental data and prior models. The open-source Integrative Modeling Platform (IMP) is one piece of software that implements all computational aspects of integrative modeling. To maximize the impact of integrative structures, the coordinates should be made publicly available, as is already the case for structures based on X-ray crystallography, NMR spectroscopy, and electron microscopy. Moreover, the associated experimental data and modeling protocols should also be archived, such that the original results can easily be reproduced. Finally, it is essential that the integrative structures are validated as part of their publication and deposition. A number of research groups have already developed software to implement integrative modeling and have generated a number of structures, prompting the formation of an Integrative/Hybrid Methods Task Force. Following the recommendations of this task force, the existing PDBx/mmCIF data representation used for atomic PDB structures has been extended to address the requirements for archiving integrative structural models. This IHM-dictionary adds a flexible model representation, including coarse graining, models in multiple states and/or related by time or other order, and multiple input experimental information sources. A prototype archiving system called PDB-Dev ( https://pdb-dev.wwpdb.org ) has also been created to archive integrative structural models, together with a Python library to facilitate handling of integrative models in PDBx/mmCIF format.

Single-Molecule 3D Images of "Hole-Hole" IgG1 Homodimers by Individual-Particle Electron Tomography

The engineering of immunoglobulin-G molecules (IgGs) is of wide interest for improving therapeutics, for example by modulating the activity or multiplexing the specificity of IgGs to recognize more than one antigen. Optimization of engineered IgG requires knowledge of three-dimensional (3D) structure of synthetic IgG. However, due to flexible nature of the molecules, their structural characterization is challenging. Here, we use our reported individual-particle electron tomography (IPET) method with optimized negative-staining (OpNS) for direct 3D reconstruction of individual IgG hole-hole homodimer molecules. The hole-hole homodimer is an undesired variant generated during the production of a bispecific antibody using the knob-into-hole heterodimer technology. A total of 64 IPET 3D density maps at ~15 Å resolutions were reconstructed from 64 individual molecules, revealing 64 unique conformations. In addition to the known Y-shaped conformation, we also observed an unusual X-shaped conformation. The 3D structure of the X-shaped conformation contributes to our understanding of the structural details of the interaction between two heavy chains in the Fc domain. The IPET approach, as an orthogonal technique to characterize the 3D structure of therapeutic antibodies, provides insight into the 3D structural variety and dynamics of heterogeneous IgG molecules.

Integrative Structure Modeling: Overview and Assessment

Integrative structure modeling computationally combines data from multiple sources of information with the aim of obtaining structural insights that are not revealed by any single approach alone. In the first part of this review, we survey the commonly used sources of structural information and the computational aspects of model building. Throughout the past decade, integrative modeling was applied to various biological systems, with a focus on large protein complexes. Recent progress in the field of cryo-electron microscopy (cryo-EM) has resolved many of these complexes to near-atomic resolution. In the second part of this review, we compare a range of published integrative models with their higher-resolution counterparts with the aim of critically assessing their accuracy. This comparison gives a favorable view of integrative modeling and demonstrates its ability to yield accurate and informative results. We discuss possible roles of integrative modeling in the new era of cryo-EM and highlight future challenges and directions.

Integrating NMR, SAXS, and Atomistic Simulations: Structure and Dynamics of a Two-Domain Protein

Multidomain proteins with two or more independently folded functional domains are prevalent in nature. Whereas most multidomain proteins are linked linearly in sequence, roughly one-tenth possess domain insertions where a guest domain is implanted into a loop of a host domain, such that the two domains are connected by a pair of interdomain linkers. Here, we characterized the influence of the interdomain linkers on the structure and dynamics of a domain-insertion protein in which the guest LysM domain is inserted into a central loop of the host CVNH domain. Expanding upon our previous crystallographic and NMR studies, we applied SAXS in combination with NMR paramagnetic relaxation enhancement to construct a structural model of the overall two-domain system. Although the two domains have no fixed relative orientation, certain orientations were found to be preferred over others. We also assessed the accuracies of molecular mechanics force fields in modeling the structure and dynamics of tethered multidomain proteins by integrating our experimental results with microsecond-scale atomistic molecular dynamics simulations. In particular, our evaluation of two different combinations of the latest force fields and water models revealed that both combinations accurately reproduce certain structural and dynamical properties, but are inaccurate for others. Overall, our study illustrates the value of integrating experimental NMR and SAXS studies with long timescale atomistic simulations for characterizing structural ensembles of flexibly linked multidomain systems.

Epitope Mapping Using Yeast Display and Next Generation Sequencing

Monoclonal antibodies are the largest class of therapeutic proteins due in part to their ability to bind an antigen with a high degree of affinity and specificity. A precise determination of their epitope is important for gaining insights into their therapeutic mechanism of action and to help differentiate antibodies that bind the same antigen. Here, we describe a method to precisely and efficiently map the epitopes of multiple antibodies in parallel over the course of just several weeks. This approach is based on a combination of rational library design, yeast surface display, and next generation DNA sequencing and provides quantitative insights into the epitope residues most critical for the antibody-antigen interaction. As an example, we will use this method to map the epitopes of several antibodies that neutralize alpha toxin from Staphylococcus aureus.

LightDock: a new multi-scale approach to protein-protein docking

Motivation

Computational prediction of protein-protein complex structure by docking can provide structural and mechanistic insights for protein interactions of biomedical interest. However, current methods struggle with difficult cases, such as those involving flexible proteins, low-affinity complexes or transient interactions. A major challenge is how to efficiently sample the structural and energetic landscape of the association at different resolution levels, given that each scoring function is often highly coupled to a specific type of search method. Thus, new methodologies capable of accommodating multi-scale conformational flexibility and scoring are strongly needed.

Results

We describe here a new multi-scale protein-protein docking methodology, LightDock, capable of accommodating conformational flexibility and a variety of scoring functions at different resolution levels. Implicit use of normal modes during the search and atomic/coarse-grained combined scoring functions yielded improved predictive results with respect to state-of-the-art rigid-body docking, especially in flexible cases.

Availability and implementation

The source code of the software and installation instructions are available for download at https://life.bsc.es/pid/lightdock/.

Contact

juanf@bsc.es.

Supplementary information

Supplementary data are available at Bioinformatics online.

Chemical cross-linking in the structural analysis of protein assemblies

For decades, chemical cross-linking of proteins has been an established method to study protein interaction partners. The chemical cross-linking approach has recently been revived by mass spectrometric analysis of the cross-linking reaction products. Chemical cross-linking and mass spectrometric analysis (CXMS) enables the identification of residues that are close in three-dimensional (3D) space but not necessarily close in primary sequence. Therefore, this approach provides medium resolution information to guide de novo structure prediction, protein interface mapping and protein complex model building. The robustness and compatibility of the CXMS approach with multiple biochemical methods have made it especially appealing for challenging systems with multiple biochemical compositions and conformation states. This review provides an overview of the CXMS approach, describing general procedures in sample processing, data acquisition and analysis. Selection of proper chemical cross-linking reagents, strategies for cross-linked peptide identification, and successful application of CXMS in structural characterization of proteins and protein complexes are discussed.

Integrating Cross-Linking Experiments with Ab Initio Protein-Protein Docking

Ab initio protein-protein docking algorithms often rely on experimental data to identify the most likely complex structure. We integrated protein-protein docking with the experimental data of chemical cross-linking followed by mass spectrometry. We tested our approach using 19 cases that resulted from an exhaustive search of the Protein Data Bank for protein complexes with cross-links identified in our experiments. We implemented cross-links as constraints based on Euclidean distance or void-volume distance. For most test cases, the rank of the top-scoring near-native prediction was improved by at least twofold compared with docking without the cross-link information, and the success rate for the top 5 predictions nearly tripled. Our results demonstrate the delicate balance between retaining correct predictions and eliminating false positives. Several test cases had multiple components with distinct interfaces, and we present an approach for assigning cross-links to the interfaces. Employing the symmetry information for these cases further improved the performance of complex structure prediction.

Modeling Structure and Dynamics of Protein Complexes with SAXS Profiles

Small-angle X-ray scattering (SAXS) is an increasingly common and useful technique for structural characterization of molecules in solution. A SAXS experiment determines the scattering intensity of a molecule as a function of spatial frequency, termed SAXS profile. SAXS profiles can be utilized in a variety of molecular modeling applications, such as comparing solution and crystal structures, structural characterization of flexible proteins, assembly of multi-protein complexes, and modeling of missing regions in the high-resolution structure. Here, we describe protocols for modeling atomic structures based on SAXS profiles. The first protocol is for comparing solution and crystal structures including modeling of missing regions and determination of the oligomeric state. The second protocol performs multi-state modeling by finding a set of conformations and their weights that fit the SAXS profile starting from a single-input structure. The third protocol is for protein-protein docking based on the SAXS profile of the complex. We describe the underlying software, followed by demonstrating their application on interleukin 33 (IL33) with its primary receptor ST2 and DNA ligase IV-XRCC4 complex.

IRaPPA: information retrieval based integration of biophysical models for protein assembly selection

MOTIVATION

In order to function, proteins frequently bind to one another and form 3D assemblies. Knowledge of the atomic details of these structures helps our understanding of how proteins work together, how mutations can lead to disease, and facilitates the designing of drugs which prevent or mimic the interaction.

RESULTS

Atomic modeling of protein-protein interactions requires the selection of near-native structures from a set of docked poses based on their calculable properties. By considering this as an information retrieval problem, we have adapted methods developed for Internet search ranking and electoral voting into IRaPPA, a pipeline integrating biophysical properties. The approach enhances the identification of near-native structures when applied to four docking methods, resulting in a near-native appearing in the top 10 solutions for up to 50% of complexes benchmarked, and up to 70% in the top 100.

AVAILABILITY AND IMPLEMENTATION

IRaPPA has been implemented in the SwarmDock server ( http://bmm.crick.ac.uk/∼SwarmDock/ ), pyDock server ( http://life.bsc.es/pid/pydockrescoring/ ) and ZDOCK server ( http://zdock.umassmed.edu/ ), with code available on request.

CONTACT

moal@ebi.ac.uk.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Improved metrics for comparing structures of macromolecular assemblies determined by 3D electron-microscopy

Recent developments in 3-dimensional electron microcopy (3D-EM) techniques and a concomitant drive to look at complex molecular structures, have led to a rapid increase in the amount of volume data available for biomolecules. This creates a demand for better methods to analyse the data, including improved scores for comparison, classification and integration of data at different resolutions. To this end, we developed and evaluated a set of scoring functions that compare 3D-EM volumes. To test our scores we used a benchmark set of volume alignments derived from the Electron Microscopy Data Bank. We find that the performance of different scores vary with the map-type, resolution and the extent of overlap between volumes. Importantly, adding the overlap information to the local scoring functions can significantly improve their precision and accuracy in a range of resolutions. A combined score involving the local mutual information and overlap (LMI_OV) performs best overall, irrespective of the map category, resolution or the extent of overlap, and we recommend this score for general use. The local mutual information score itself is found to be more discriminatory than cross-correlation coefficient for intermediate-to-low resolution maps or when the map size and density distribution differ significantly. For comparing map surfaces, we implemented two filters to detect the surface points, including one based on the 'extent of surface exposure'. We show that scores that compare surfaces are useful at low resolutions and for maps with evident surface features. All the scores discussed are implemented in TEMPy (http://tempy.ismb.lon.ac.uk/).

Performance of ZDOCK and IRAD in CAPRI rounds 28-34

We report the performance of our protein-protein docking pipeline, including the ZDOCK rigid-body docking algorithm, on 19 targets in CAPRI rounds 28-34. Following the docking step, we reranked the ZDOCK predictions using the IRAD scoring function, pruned redundant predictions, performed energy landscape analysis, and utilized our interface prediction approach RCF. In addition, we applied constraints to the search space based on biological information that we culled from the literature, which increased the chance of making a correct prediction. For all but two targets we were able to find and apply biological information and we found the information to be highly accurate, indicating that effective incorporation of biological information is an important component for protein-protein docking. Proteins 2017; 85:408-416. © 2016 Wiley Periodicals, Inc.

Insights into HIV-1 proviral transcription from integrative structure and dynamics of the Tat:AFF4:P-TEFb:TAR complex

HIV-1 Tat hijacks the human superelongation complex (SEC) to promote proviral transcription. Here we report the 5.9 Å structure of HIV-1 TAR in complex with HIV-1 Tat and human AFF4, CDK9, and CycT1. The TAR central loop contacts the CycT1 Tat-TAR recognition motif (TRM) and the second Tat Zn²⁺-binding loop. Hydrogen-deuterium exchange (HDX) shows that AFF4 helix 2 is stabilized in the TAR complex despite not touching the RNA, explaining how it enhances TAR binding to the SEC 50-fold. RNA SHAPE and SAXS data were used to help model the extended (Tat Arginine-Rich Motif) ARM, which enters the TAR major groove between the bulge and the central loop. The structure and functional assays collectively support an integrative structure and a bipartite binding model, wherein the TAR central loop engages the CycT1 TRM and compact core of Tat, while the TAR major groove interacts with the extended Tat ARM.

DOI:http://dx.doi.org/10.7554/eLife.15910.001

Diffraction Techniques in Structural Biology

A detailed understanding of chemical and biological function and the mechanisms underlying the molecular activities ultimately requires atomic-resolution structural data. Diffraction-based techniques such as single-crystal X-ray crystallography, electron microscopy, and neutron diffraction are well established and they have paved the road to the stunning successes of modern-day structural biology. The major advances achieved in the last twenty years in all aspects of structural research, including sample preparation, crystallization, the construction of synchrotron and spallation sources, phasing approaches, and high-speed computing and visualization, now provide specialists and nonspecialists alike with a steady flow of molecular images of unprecedented detail. The present unit combines a general overview of diffraction methods with a detailed description of the process of a single-crystal X-ray structure determination experiment, from chemical synthesis or expression to phasing and refinement, analysis, and quality control. For novices it may serve as a stepping-stone to more in-depth treatises of the individual topics. Readers relying on structural information for interpreting functional data may find it a useful consumer guide. © 2016 by John Wiley & Sons, Inc.

Prion Protein-Antibody Complexes Characterized by Chromatography-Coupled Small-Angle X-Ray Scattering

Aberrant self-assembly, induced by structural misfolding of the prion proteins, leads to a number of neurodegenerative disorders. In particular, misfolding of the mostly α-helical cellular prion protein (PrP(C)) into a β-sheet-rich disease-causing isoform (PrP(Sc)) is the key molecular event in the formation of PrP(Sc) aggregates. The molecular mechanisms underlying the PrP(C)-to-PrP(Sc) conversion and subsequent aggregation remain to be elucidated. However, in persistently prion-infected cell-culture models, it was shown that treatment with monoclonal antibodies against defined regions of the prion protein (PrP) led to the clearing of PrP(Sc) in cultured cells. To gain more insight into this process, we characterized PrP-antibody complexes in solution using a fast protein liquid chromatography coupled with small-angle x-ray scattering (FPLC-SAXS) procedure. High-quality SAXS data were collected for full-length recombinant mouse PrP [denoted recPrP(23-230)] and N-terminally truncated recPrP(89-230), as well as their complexes with each of two Fab fragments (HuM-P and HuM-R1), which recognize N- and C-terminal epitopes of PrP, respectively. In-line measurements by fast protein liquid chromatography coupled with SAXS minimized data artifacts caused by a non-monodispersed sample, allowing structural analysis of PrP alone and in complex with Fab antibodies. The resulting structural models suggest two mechanisms for how these Fabs may prevent the conversion of PrP(C) into PrP(Sc).

FoXS, FoXSDock and MultiFoXS: Single-state and multi-state structural modeling of proteins and their complexes based on SAXS profiles

Small Angle X-ray Scattering (SAXS) is an increasingly common and useful technique for structural characterization of molecules in solution. A SAXS experiment determines the scattering intensity of a molecule as a function of spatial frequency, termed SAXS profile. Here, we describe three web servers for modeling atomic structures based on SAXS profiles. FoXS (Fast X-Ray Scattering) rapidly computes a SAXS profile of a given atomistic model and fits it to an experimental profile. FoXSDock docks two rigid protein structures based on a SAXS profile of their complex. MultiFoXS computes a population-weighted ensemble starting from a single input structure by fitting to a SAXS profile of the protein in solution. We describe the interfaces and capabilities of the servers (salilab.org/foxs), followed by demonstrating their application on Interleukin-33 (IL-33) and its primary receptor ST2.

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.

DksA regulates RNA polymerase in Escherichia coli through a network of interactions in the secondary channel that includes Sequence Insertion 1

Sensing and responding to nutritional status is a major challenge for microbial life. In Escherichia coli, the global response to amino acid starvation is orchestrated by guanosine-3',5'-bisdiphosphate and the transcription factor DksA. DksA alters transcription by binding to RNA polymerase and allosterically modulating its activity. Using genetic analysis, photo-cross-linking, and structural modeling, we show that DksA binds and acts upon RNA polymerase through prominent features of both the nucleotide-access secondary channel and the active-site region. This work is, to our knowledge, the first demonstration of a molecular function for Sequence Insertion 1 in the β subunit of RNA polymerase and significantly advances our understanding of how DksA binds to RNA polymerase and alters transcription.

NPPD: A Protein-Protein Docking Scoring Function Based on Dyadic Differences in Networks of Hydrophobic and Hydrophilic Amino Acid Residues

Protein-protein docking (PPD) predictions usually rely on the use of a scoring function to rank docking models generated by exhaustive sampling. To rank good models higher than bad ones, a large number of scoring functions have been developed and evaluated, but the methods used for the computation of PPD predictions remain largely unsatisfactory. Here, we report a network-based PPD scoring function, the NPPD, in which the network consists of two types of network nodes, one for hydrophobic and the other for hydrophilic amino acid residues, and the nodes are connected when the residues they represent are within a certain contact distance. We showed that network parameters that compute dyadic interactions and those that compute heterophilic interactions of the amino acid networks thus constructed allowed NPPD to perform well in a benchmark evaluation of 115 PPD scoring functions, most of which, unlike NPPD, are based on some sort of protein-protein interaction energy. We also showed that NPPD was highly complementary to these energy-based scoring functions, suggesting that the combined use of conventional scoring functions and NPPD might significantly improve the accuracy of current PPD predictions.

Emerging applications of small angle solution scattering in structural biology

Small angle solution X-ray and neutron scattering recently resurfaced as powerful tools to address an array of biological problems including folding, intrinsic disorder, conformational transitions, macromolecular crowding, and self or hetero-assembling of biomacromolecules. In addition, small angle solution scattering complements crystallography, nuclear magnetic resonance spectroscopy, and other structural methods to aid in the structure determinations of multidomain or multicomponent proteins or nucleoprotein assemblies. Neutron scattering with hydrogen/deuterium contrast variation, or X-ray scattering with sucrose contrast variation to a certain extent, is a convenient tool for characterizing the organizations of two-component systems such as a nucleoprotein or a lipid-protein assembly. Time-resolved small and wide-angle solution scattering to study biological processes in real time, and the use of localized heavy-atom labeling and anomalous solution scattering for applications as FRET-like molecular rulers, are amongst promising newer developments. Despite the challenges in data analysis and interpretation, these X-ray/neutron solution scattering based approaches hold great promise for understanding a wide variety of complex processes prevalent in the biological milieu.

Inferring the microscopic surface energy of protein-protein interfaces from mutation data

Mutations at protein-protein recognition sites alter binding strength by altering the chemical nature of the interacting surfaces. We present a simple surface energy model, parameterized with empirical ΔΔG values, yielding mean energies of -48 cal mol(-1) Å(-2) for interactions between hydrophobic surfaces, -51 to -80 cal mol(-1) Å(-2) for surfaces of complementary charge, and 66-83 cal mol(-1) Å(-2) for electrostatically repelling surfaces, relative to the aqueous phase. This places the mean energy of hydrophobic surface burial at -24 cal mol(-1) Å(-2) . Despite neglecting configurational entropy and intramolecular changes, the model correlates with empirical binding free energies of a functionally diverse set of rigid-body interactions (r = 0.66). When used to rerank docking poses, it can place near-native solutions in the top 10 for 37% of the complexes evaluated, and 82% in the top 100. The method shows that hydrophobic burial is the driving force for protein association, accounting for 50-95% of the cohesive energy. The model is available open-source from http://life.bsc.es/pid/web/surface_energy/ and via the CCharpPPI web server http://life.bsc.es/pid/ccharppi/.

Sequence co-evolution gives 3D contacts and structures of protein complexes

Protein-protein interactions are fundamental to many biological processes. Experimental screens have identified tens of thousands of interactions, and structural biology has provided detailed functional insight for select 3D protein complexes. An alternative rich source of information about protein interactions is the evolutionary sequence record. Building on earlier work, we show that analysis of correlated evolutionary sequence changes across proteins identifies residues that are close in space with sufficient accuracy to determine the three-dimensional structure of the protein complexes. We evaluate prediction performance in blinded tests on 76 complexes of known 3D structure, predict protein-protein contacts in 32 complexes of unknown structure, and demonstrate how evolutionary couplings can be used to distinguish between interacting and non-interacting protein pairs in a large complex. With the current growth of sequences, we expect that the method can be generalized to genome-wide elucidation of protein-protein interaction networks and used for interaction predictions at residue resolution.