Combinatorial docking approach for structure prediction of large proteins and multi-molecular assemblies

Protein folding and protein binding are similar processes. In both, structural units combinatorially associate with each other. In the case of folding, we mostly handle relatively small units, building blocks or domains, that are covalently linked. In the case of multi-molecular binding, the subunits are relatively large and are associated only by non-covalent bonds. Experimentally, the difficulty in the determination of the structures of such large assemblies increases with the complex size and the number of components it contains. Computationally, the prediction of the structures of multi-molecular complexes has largely not been addressed, probably owing to the magnitude of the combinatorial complexity of the problem. Current docking algorithms mostly target prediction of pairwise interactions. Here our goal is to predict the structures of multi-unit associations, whether these are chain-connected as in protein folding, or separate disjoint molecules in the assemblies. We assume that the structures of the single units are known, either through experimental determination or modeling. Our aim is to combinatorially assemble these units to predict their structure. To address this problem we have developed CombDock. CombDock is a combinatorial docking algorithm for the structural units assembly problem. Below, we briefly describe the algorithm and present examples of its various applications to folding and to multi-molecular assemblies. To test the robustness of the algorithm, we use inaccurate models of the structural units, derived either from crystal structures of unbound molecules or from modeling of the target sequences. The algorithm has been able to predict near-native arrangements of the input structural units in almost all of the cases, suggesting that a combinatorial approach can overcome the imperfect shape complementarity caused by the inaccuracy of the models. In addition, we further show that through a combinatorial docking strategy it is possible to enhance the predictions of pairwise interactions involved in a multi-molecular assembly.     

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.     

Accuracy of bound peptide structures determined by exchange transferred nuclear Overhauser data: A simulation study

The exchange-transferred NOE method to determine the three-dimensional structure of peptides bound to proteins, or other macromolecular systems, is becoming increasingly important in drug design efforts and for large or multicomponent assemblies, such as membrane receptors, where structural analysis of the full system is intractable. The exchange-transferred nuclear Overhauser effect spectroscopy (etNOESY) method allows the determination of the bound-state conformation of the peptide from the intra-molecular NOE interactions between ligand protons. Because only ligand–ligand NOEs are generally observable, the etNOESY method is restricted to fewer NOEs per residue than direct protein structure determination. In addition, the averaging of relaxation rates between free and bound states affects the measured cross-peak intensities, and possibly the accuracy of distance estimates. Accordingly, the study reported here was conducted to examine the conditions required to define a reliable structure. The program CORONA was used to simulate etNOE data using a rate-matrix including magnetic relaxation and exchange rates for two peptide–protein complexes derived from the reference complex of cAMP-dependent protein kinase ligated with a 24-residue inhibitor peptide. The results indicate that reasonably accurate peptide structures can be determined with relatively few NOE interactions when the interactions occur between non-neighboring residues. The reliability of the structural result is suggested from the pattern of NOE interactions. A structure with an accuracy of approximately 1.3 Å rms difference for the main-chain atoms can be obtained when etNOE interactions between non-neighboring residues occur over the length of the peptide. The global precision is higher (approximately 0.9 Å rms difference) but is not correlated to global accuracy. A local definition of precision along the backbone appears to be a good indicator of the local accuracy.     

Protein-nucleic acid complexes and the role of mass spectrometry in their structure determination

Mass spectrometry is now established as a powerful tool for the study of the stoichiometry, interactions, dynamics, and subunit architecture of large protein assemblies and their subcomplexes. Recent evidence has suggested that the 3D structure of protein complexes can be maintained intact in the gas phase, highlighting the potential of ion mobility to contribute to structural biology. A key challenge is to integrate the compositional and structural information from ion mobility mass spectrometry with molecular modelling approaches to produce 3D models of intact protein complexes. In this review, we focus on the mass spectrometry of protein-nucleic acid assemblies with particular attention to the application of ion mobility, an emerging technique in structural studies. We also discuss the challenges that lie ahead for the full integration of ion mobility mass spectrometry with structural biology.     

Protein-nucleic acid complexes and the role of mass spectrometry in their structure determination

Mass spectrometry is now established as a powerful tool for the study of the stoichiometry, interactions, dynamics, and subunit architecture of large protein assemblies and their subcomplexes. Recent evidence has suggested that the 3D structure of protein complexes can be maintained intact in the gas phase, highlighting the potential of ion mobility to contribute to structural biology. A key challenge is to integrate the compositional and structural information from ion mobility mass spectrometry with molecular modelling approaches to produce 3D models of intact protein complexes. In this review, we focus on the mass spectrometry of protein-nucleic acid assemblies with particular attention to the application of ion mobility, an emerging technique in structural studies. We also discuss the challenges that lie ahead for the full integration of ion mobility mass spectrometry with structural biology.     

Ab initio modeling of the herpesvirus VP26 core domain assessed by CryoEM density

Efforts in structural biology have targeted the systematic determination of all protein structures through experimental determination or modeling. In recent years, 3-D electron cryomicroscopy (cryoEM) has assumed an increasingly important role in determining the structures of these large macromolecular assemblies to intermediate resolutions (6–10 Å). While these structures provide a snapshot of the assembly and its components in well-defined functional states, the resolution limits the ability to build accurate structural models. In contrast, sequence-based modeling techniques are capable of producing relatively robust structural models for isolated proteins or domains. In this work, we developed and applied a hybrid modeling approach, utilizing cryoEM density and ab initio modeling to produce a structural model for the core domain of a herpesvirus structural protein, VP26. Specifically, this method, first tested on simulated data, utilizes the cryoEM density map as a geometrical constraint in identifying the most native-like models from a gallery of models generated by ab initio modeling. The resulting model for the core domain of VP26, based on the 8.5-Å resolution herpes simplex virus type 1 (HSV-1) capsid cryoEM structure and mutational data, exhibited a novel fold. Additionally, the core domain of VP26 appeared to have a complementary interface to the known upper-domain structure of VP5, its cognate binding partner. While this new model provides for a better understanding of the assembly and interactions of VP26 in HSV-1, the approach itself may have broader applications in modeling the components of large macromolecular assemblies.     

Review: electron crystallography: present excitement, a nod to the past, anticipating the future

From a modest beginning with negatively stained samples of the helical T4 bacteriophage tail, electron crystallography has emerged as a powerful tool in structural biology. High-resolution density maps, interpretable in terms of an atomic structure, can be obtained from specimens prepared as well-ordered, two-dimensional crystals, and the resolution achieved with helical specimens and icosahedral viruses is approaching the same goal. A hybrid approach to determining the molecular structure of complex biological assemblies is generating great interest, in which high-resolution structures that have been determined for individual protein components are fitted into a lower resolution envelope of the large complex. With this as background, how much more can be anticipated for the future? Considerable scope still remains to improve the quality of electron microscope images. Automation of data acquisition and data processing, together with the emergence of computational speeds of 10(12) floating point operations per second or higher, will make it possible to extend high-resolution structure determination into the realm of single-particle microscopy. As a result, computational alignment of single particles, i.e., the formation of "virtual crystals," can begin to replace biochemical crystallization. Since single-particle microscopy may remain limited to "large" structures of 200 to 300 kDa or more, however, smaller proteins will continue to be studied as helical assemblies or as two-dimensional crystals. The further development of electron crystallography is thus likely to turn increasingly to the use of single particles and small regions of ordered assemblies, emphasizing more and more the potential for faster, higher throughput.     

Joining forces: integrating proteomics and cross-linking with the mass spectrometry of intact complexes

Protein assemblies are critical for cellular function and understanding their physical organization is the key aim of structural biology. However, applying conventional structural biology approaches is challenging for transient, dynamic, or polydisperse assemblies. There is therefore a growing demand for hybrid technologies that are able to complement classical structural biology methods and thereby broaden our arsenal for the study of these important complexes. Exciting new developments in the field of mass spectrometry and proteomics have added a new dimension to the study of protein-protein interactions and protein complex architecture. In this review, we focus on how complementary mass spectrometry-based techniques can greatly facilitate structural understanding of protein assemblies.     

αB-crystallin polydispersity is a consequence of unbiased quaternary dynamics

The inherent heterogeneity of many protein assemblies complicates characterization of their structure and dynamics, as most biophysical techniques require homogeneous preparations of isolated components. For this reason, quantitative studies of the molecular chaperone αB-crystallin, which populates a range of interconverting oligomeric states, have been difficult, and the physicochemical basis for its polydispersity has remained unknown. Here, we perform mass spectrometry experiments to study αB-crystallin and extract detailed information as to its oligomeric distribution and exchange of subunits under a range of conditions. This allows a determination of the thermodynamic and kinetic parameters that govern the polydisperse ensemble and enables the construction of a simple energy profile for oligomerization. We find that the quaternary structure and dynamics of the protein can be explained using a simple model with just two oligomer-independent interactions (i.e., interactions that are energetically identical in all oligomers from 10mers to 40mers) between constituent monomers. As such, the distribution of oligomers is governed purely by the dynamics of individual monomers. This provides a new means for understanding the polydispersity of αB-crystallin and a framework for interrogating other heterogeneous protein assemblies.     

Protein-protein interaction and quaternary structure

Protein-protein recognition plays an essential role in structure and function. Specific non-covalent interactions stabilize the structure of macromolecular assemblies, exemplified in this review by oligomeric proteins and the capsids of icosahedral viruses. They also allow proteins to form complexes that have a very wide range of stability and lifetimes and are involved in all cellular processes. We present some of the structure-based computational methods that have been developed to characterize the quaternary structure of oligomeric proteins and other molecular assemblies and analyze the properties of the interfaces between the subunits. We compare the size, the chemical and amino acid compositions and the atomic packing of the subunit interfaces of protein-protein complexes, oligomeric proteins, viral capsids and protein-nucleic acid complexes. These biologically significant interfaces are generally close-packed, whereas the non-specific interfaces between molecules in protein crystals are loosely packed, an observation that gives a structural basis to specific recognition. A distinction is made within each interface between a core that contains buried atoms and a solvent accessible rim. The core and the rim differ in their amino acid composition and their conservation in evolution, and the distinction helps correlating the structural data with the results of site-directed mutagenesis and in vitro studies of self-assembly.     

Studying biomolecular complexes with pulsed electron-electron double resonance spectroscopy

The function of biomolecules is intrinsically linked to their structure and the complexes they form during function. Techniques for the determination of structures and dynamics of these nanometre assemblies are therefore important for an understanding on the molecular level. PELDOR (pulsed electron-electron double resonance) is a pulsed EPR method that can be used to reliably and precisely measure distances in the range 1.5-8?nm, to unravel orientations and to determine the number of monomers in complexes. In conjunction with site-directed spin labelling, it can be applied to biomolecules of all sizes in aqueous solutions or membranes. PELDOR is therefore complementary to the methods of X-ray crystallography, NMR and FRET (fluorescence resonance energy transfer) and is becoming a powerful method for structural determination of biomolecules. In the present review, the methods of PELDOR are discussed and examples where PELDOR has been used to obtain structural information on biomolecules are summarized.     

From high symmetry to high resolution in biological electron microscopy: a commentary on Crowther (1971) ‘Procedures for three-dimensional reconstruction of spherical viruses by Fourier synthesis from electron micrographs’

Elucidation of the structure of biological macromolecules and larger assemblies has been essential to understanding the roles they play in living processes. Methods for three-dimensional structure determination of biological assemblies from images recorded in the electron microscope were therefore a key development. In his paper published in Philosophical Transactions B in 1971, Crowther described new computational procedures applied to the first three-dimensional reconstruction of an icosahedral virus from images of virus particles preserved in negative stain. The method for determining the relative orientation of randomly oriented particles and combining their images for reconstruction exploited the high symmetry of the virus particle. Computational methods for image analysis have since been extended to include biological assemblies without symmetry. Further experimental advances, combined with image analysis, have led to the method of cryomicroscopy, which is now used by structural biologists to study the structure and dynamics of biological machines and assemblies in atomic detail. This commentary was written to celebrate the 350th anniversary of the journal Philosophical Transactions of the Royal Society.     

Mapping the structure of macromolecular assemblies by combining chemical modification and separation methods

Several examples will be described in which powerful separation methods are combined with relatively simple chemical modification techniques to provide structural information on complex macromolecular assemblies. Ribosomal RNA structure has been examined by crosslinking, separating individual crosslinked species by gel electrophoresis, and enzymatic methods for determination of crosslink positions in the nucleotide sequence. Chromatin structure has been examined by footprinting the location of individual nucleosomes by a combination of chemical nicking and DNA separations. Virus structure can be examined by using breakable crosslinkers analyzed with diagonal gel electrophoresis. Ultimately such methods may allow structural information to be obtained on systems even as complex as whole chromosomes.     

Visualizing transiently associated catalytic domains in assembly-line biosynthesis using cryo-electron microscopy

While cryo-electron microscopy (cryo-EM) has revolutionized the structure determination of supramolecular protein complexes that are refractory to structure determination by X-ray crystallography, structure determination by cryo-EM can nonetheless be complicated by excessive conformational flexibility or structural heterogeneity resulting from weak or transient protein–protein association. Since such transient complexes are often critical for function, specialized approaches must be employed for the determination of meaningful structure–function relationships. Here, we outline examples in which transient protein–protein interactions have been visualized successfully by cryo-EM in the biosynthesis of fatty acids, polyketides, and terpenes. These studies demonstrate the utility of chemical crosslinking to stabilize transient protein–protein complexes for cryo-EM structural analysis, as well as the use of partial signal subtraction and localized reconstruction to extract useful structural information out of cryo-EM data collected from inherently dynamic systems. While these approaches do not always yield atomic resolution insights on protein–protein interactions, they nonetheless enable direct experimental observation of complexes in assembly-line biosynthesis that would otherwise be too fleeting for structural analysis.     

Structural modeling of heteromeric protein complexes from disassembly pathways and ion mobility-mass spectrometry

Structure determination of macromolecular protein assemblies remains a challenge for well-established methods. Here, we provide an assessment of an emerging structural technique, ion mobility-mass spectrometry (IM-MS), and examine the use of collision cross-sections (CCSs), derived from IM-MS, as restraints for structure characterization of heteromeric protein assemblies. Using 15 complexes selected from the Protein Data Bank, we validate the use of low-resolution models by comparing their CCSs with those calculated for all-atom structures. We then select six heteromeric complexes, disrupting them in solution to form subcomplexes. Experimental and calculated CCSs reveal close similarity for 18 of the 21 (sub)complexes. Exploring the use of CCS as a restraint, we incorporate it into a scoring function and show good correlation between the score and similarity to the native structure for heteromers, especially when an additional symmetry restraint was introduced.     

DynBench3D,a Web-Resource to Dynamically Generate Benchmark Sets of Large Heteromeric Protein Complexes

Multi-protein machines are responsible for most cellular tasks, and many efforts have been invested in the systematic identification and characterization of thousands of these macromolecular assemblies. However, unfortunately, the (quasi) atomic details necessary to understand their function are available only for a tiny fraction of the known complexes. The computational biology community is developing strategies to integrate structural data of different nature, from electron microscopy to X-ray crystallography, to model large molecular machines, as it has been done for individual proteins and interactions with remarkable success. However, unlike for binary interactions, there is no reliable gold-standard set of three-dimensional (3D) complexes to benchmark the performance of these methodologies and detect their limitations. Here, we present a strategy to dynamically generate non-redundant sets of 3D heteromeric complexes with three or more components. By changing the values of sequence identity and component overlap between assemblies required to define complex redundancy, we can create sets of representative complexes with known 3D structure (i.e., target complexes). Using an identity threshold of 20% and imposing a fraction of component overlap of < 0.5, we identify 495 unique target complexes, which represent a real non-redundant set of heteromeric assemblies with known 3D structure. Moreover, for each target complex, we also identify a set of assemblies, of varying degrees of identity and component overlap, that can be readily used as input in a complex modeling exercise (i.e., template subcomplexes). We hope that resources like this will significantly help the development and progress assessment of novel methodologies, as docking benchmarks and blind prediction contests did. The interactive resource is accessible at https://DynBench3D.irbbarcelona.org.     

New developments in protein structure-function analysis by MS and use of hydrogen-deuterium exchange microfluidics

The study of protein structure and function has evolved to become a leading discipline in the biophysical sciences. Although it is not yet possible to determine 3D protein structures from MS data alone, multiple MS-based techniques can be combined to obtain structural and functional data that are complementary to classical protein structure information obtained from NMR or X-ray crystallography. Monitoring gas-phase interactions of noncovalent complexes yields information on binding constants, complex stability, and the nature of interactions. Ion mobility MS and chemical crosslinking strategies can be applied to probe the architecture of macromolecular assemblies and protein-ligand complexes. MS analysis of hydrogen-deuterium exchange can be used to determine the localization of secondary structure elements, binding sites and conformational dynamics of proteins in solution. This minireview focuses first on new strategies that combine these techniques to gain insights into protein structure and function. Using one such strategy, we then demonstrate how a novel hydrogen-deuterium exchange microfluidics tool can be used online with an ESI mass spectrometer to monitor regional accessibility in a peptide, as exemplified with amyloid-β peptide 1-40.     

Direct determination of changes of interdomain orientation on ligation: use of the orientational dependence of 15N NMR relaxation in Abl SH(32).

The relative orientation and motions of domains within many proteins are key to the control of multivalent recognition, or the assembly of protein-based cellular machines. Current methods of structure determination have limited applicability to macromolecular assemblies, characterized by weak interactions between the constituents. Crystal structures of such complexes might be biased by packing forces comparable to the interdomain interactions, while the precision and accuracy of the conventional NMR structural approaches are necessarily limited by the restricted number of NOE contacts and by interdomain flexibility rendering the available NOE information uninterpretable. NMR relaxation studies are capable of providing "long-range" structural information on macromolecules in their native milieu. Here we determine directly the change in domain orientation between unligated and dual ligated subdomains of the SH(32) segment of Abelson kinase in solution, using the orientational dependence of nuclear spin relaxation. These results demonstrate that the change in domain orientation between unligated and ligated forms can be measured directly in solution.     

Transient DNA / RNA-protein interactions

The great pace of biomolecular structure determination has provided a plethora of protein structures, but not as many structures of nucleic acids or of their complexes with proteins. The recognition of DNA and RNA molecules by proteins may produce large and relatively stable assemblies (such as the ribosome) or transient complexes (such as DNA clamps sliding through the DNA). These transient interactions are most difficult to characterize, but even in 'stable' complexes captured in crystal structures, the dynamics of the whole or part of the assembly pose great technical difficulties in understanding their function. The development and refinement of powerful experimental and computational tools have made it possible to learn a great deal about the relevance of these fleeting events for numerous biological processes. We discuss here the most recent findings and the challenges that lie ahead in the quest for a better understanding of protein-nucleic acid interactions.     

Implications of aromatic–aromatic interactions: From protein structures to peptide models

With increasing structural information on proteins, the opportunity to understand physical forces governing protein folding is also expanding. One of the significant non‐covalent forces between the protein side chains is aromatic–aromatic interactions. Aromatic interactions have been widely exploited and thoroughly investigated in the context of folding, stability, molecular recognition, and self‐assembly processes. Through this review, we discuss the contribution of aromatic interactions to the activity and stability of thermophilic, mesophilic, and psychrophilic proteins. Being hydrophobic, aromatic amino acids tend to reside in the protein hydrophobic interior or transmembrane segments of proteins. In such positions, it can play a diverse role in soluble and membrane proteins, and in α‐helix and β‐sheet stabilization. We also highlight here some excellent investigations made using peptide models and several approaches involving aryl–aryl interactions, as an increasingly popular strategy in protein and peptide engineering. A recent survey described the existence of aromatic clusters (trimer, tetramer, pentamer, and higher order assemblies), revealing the self‐associating property of aryl groups, even in folded protein structures. The application of this self‐assembly of aromatics in the generation of modern bionanomaterials is also discussed.