Symmetric Protein Complexes Research Articles

Mitigating the Blurring Effect of CryoEM Averaging on a Flexible and Highly Symmetric Protein Complex through Sub-Particle Reconstruction.

Many macromolecules are inherently flexible as a feature of their structure and function. During single-particle CryoEM processing, flexible protein regions can be detrimental to high-resolution reconstruction as signals from thousands of particles are averaged together. This "blurring" effect can be difficult to overcome and is possibly more pronounced when averaging highly symmetric complexes. Approaches to mitigating flexibility during CryoEM processing are becoming increasingly critical as the technique advances and is applied to more dynamic proteins and complexes. Here, we detail the use of sub-particle averaging and signal subtraction techniques to precisely target and resolve flexible DARPin protein attachments on a designed tetrahedrally symmetric protein scaffold called DARP14. Particles are first aligned as full complexes, and then the symmetry is reduced by alignment and focused refinement of the constituent subunits. The final reconstructions we obtained were vastly improved over the fully symmetric reconstructions, with observable secondary structure and side-chain placement. Additionally, we were also able to reconstruct the core region of the scaffold to 2.7 Å. The data processing protocol outlined here is applicable to other dynamic and symmetric protein complexes, and our improved maps could allow for new structure-guided variant designs of DARP14.

Improved interface packing and design opportunities revealed by CryoEM analysis of a designed protein nanocage

Symmetric protein assemblies play important roles in nature which makes them an attractive target for engineering. De novo symmetric protein complexes can be created through computational protein design to tailor their properties from first principles, and recently several protein nanocages have been created by bringing together protein components through hydrophobic interactions. Accurate experimental structures of newly-developed proteins are essential to validate their design, improve assembly stability, and tailor downstream applications. We describe the CryoEM structure of the nanocage I3-01, at an overall resolution of 3.5 Å. I3-01, comprising 60 aldolase subunits arranged with icosahedral symmetry, has resisted high-resolution characterization. Some key differences between the refined structure and the original design are identified, such as improved packing of hydrophobic sidechains, providing insight to the resistance of I3-01 to high-resolution averaging. Based on our analysis, we suggest factors important in the design and structural processing of new assemblies.

Computational Design of Single-Peptide Nanocages with Nanoparticle Templating.

Protein complexes perform a diversity of functions in natural biological systems. While computational protein design has enabled the development of symmetric protein complexes with spherical shapes and hollow interiors, the individual subunits often comprise large proteins. Peptides have also been applied to self-assembly, and it is of interest to explore such short sequences as building blocks of large, designed complexes. Coiled-coil peptides are promising subunits as they have a symmetric structure that can undergo further assembly. Here, an α-helical 29-residue peptide that forms a tetrameric coiled coil was computationally designed to assemble into a spherical cage that is approximately 9 nm in diameter and presents an interior cavity. The assembly comprises 48 copies of the designed peptide sequence. The design strategy allowed breaking the side chain conformational symmetry within the peptide dimer that formed the building block (asymmetric unit) of the cage. Dynamic light scattering (DLS) and transmission electron microscopy (TEM) techniques showed that one of the seven designed peptide candidates assembled into individual nanocages of the size and shape. The stability of assembled nanocages was found to be sensitive to the assembly pathway and final solution conditions (pH and ionic strength). The nanocages templated the growth of size-specific Au nanoparticles. The computational design serves to illustrate the possibility of designing target assemblies with pre-determined specific dimensions using short, modular coiled-coil forming peptide sequences.

Geometric description of self-interaction potential in symmetric protein complexes

Proteins can self-associate with copies of themselves to form symmetric complexes called homomers. Homomers are widespread in all kingdoms of life and allow for unique geometric and functional properties, as reflected in viral capsids or allostery. Once a protein forms a homomer, however, its internal symmetry can compound the effect of point mutations and trigger uncontrolled self-assembly into high-order structures. We identified mutation hot spots for supramolecular assembly, which are predictable by geometry. Here, we present a dataset of descriptors that characterize these hot spot positions both geometrically and chemically, as well as computer scripts allowing the calculation and visualization of these properties for homomers of choice. Since the biological relevance of homomers is not readily available from their X-ray crystallographic structure, we also provide reliability estimates obtained by methods we recently developed. These data have implications in the study of disease-causing mutations, protein evolution and can be exploited in the design of biomaterials.

CHDOCK: a hierarchical docking approach for modeling Cn symmetric homo-oligomeric complexes

Protein–protein interactions are crucial in many biological processes. Therefore, determining the complex structure between proteins is valuable for understanding the molecular mechanism and developing drugs. Many proteins like ion channels are formed by symmetric homo-oligomers. In this study, we have proposed a hierarchical docking algorithm to predict the structure of Cn symmetric protein complexes, which is referred to as CHDOCK. The symmetric binding modes were first constructed by an FFT-based docking algorithm and then optimized by our iterative scoring function for protein–protein interactions. When tested on a symmetric protein docking benchmark of 212 homo-oligomeric complexes with Cn symmetry, CHDOCK obtained a significantly better performance in binding mode predictions than three state-of-the-art symmetric docking methods, M-ZDOCK, SAM, and SymmDock. When the top 10 predictions were considered, CHDOCK achieved a success rate of 44.81% and 72.17% for unbound docking and bound docking, respectively in comparison to those of 36.79% and 65.09% for M-ZDOCK, 31.60% and 54.25% for SAM, and 30.66% and 31.60% for SymmDock. CHDOCK is computationally efficient and can normally complete a symmetric docking calculation within 30 min. The CHDOCK can be freely accessed by a web server at http://huanglab.phys.hust.edu.cn/hsymdock/.

HSYMDOCK: a docking web server for predicting the structure of protein homo-oligomers with Cn or Dn symmetry.

A major subclass of protein–protein interactions is formed by homo-oligomers with certain symmetry. Therefore, computational modeling of the symmetric protein complexes is important for understanding the molecular mechanism of related biological processes. Although several symmetric docking algorithms have been developed for Cn symmetry, few docking servers have been proposed for Dn symmetry. Here, we present HSYMDOCK, a web server of our hierarchical symmetric docking algorithm that supports both Cn and Dn symmetry. The HSYMDOCK server was extensively evaluated on three benchmarks of symmetric protein complexes, including the 20 CASP11–CAPRI30 homo-oligomer targets, the symmetric docking benchmark of 213 Cn targets and 35 Dn targets, and a nonredundant test set of 55 transmembrane proteins. It was shown that HSYMDOCK obtained a significantly better performance than other similar docking algorithms. The server supports both sequence and structure inputs for the monomer/subunit. Users have an option to provide the symmetry type of the complex, or the server can predict the symmetry type automatically. The docking process is fast and on average consumes 10∼20 min for a docking job. The HSYMDOCK web server is available at http://huanglab.phys.hust.edu.cn/hsymdock/.

The finite number of global motion patterns available to symmetric protein complexes.

In PDB, more than half of the entries are structure complexes and of these complexes, most are symmetric, composed of identical subunits. Complex formation is the way through which larger structures and even molecular machines are assembled and built in nature. In this work, we apply group theory and carry out a comprehensive study of the global motion patterns of protein complexes of various symmetries. The work presents for the first time a comprehensive list of all the symmetric, aesthetically pleasing, global motion patterns available to complexes of cyclic, dihedral, tetrahedral, or octahedral symmetry. Our results clearly demonstrate that complexes with the same symmetry will have the same global motion patterns and thus may function in a similar way, and that there are only a finite number of global motion patterns available to symmetric complexes as the number of protein symmetry groups is effectively finite. The work complements our current understanding of the principle of complex formation that has been mostly structure-based by providing novel dynamics-based insights. Furthermore, as dynamics is closely tied to function, these motion patterns can provide global insights into the general functional mechanisms of protein complexes.

An Integrated Framework Advancing Membrane Protein Modeling and Design.

Membrane proteins are critical functional molecules in the human body, constituting more than 30% of open reading frames in the human genome. Unfortunately, a myriad of difficulties in overexpression and reconstitution into membrane mimetics severely limit our ability to determine their structures. Computational tools are therefore instrumental to membrane protein structure prediction, consequently increasing our understanding of membrane protein function and their role in disease. Here, we describe a general framework facilitating membrane protein modeling and design that combines the scientific principles for membrane protein modeling with the flexible software architecture of Rosetta3. This new framework, called RosettaMP, provides a general membrane representation that interfaces with scoring, conformational sampling, and mutation routines that can be easily combined to create new protocols. To demonstrate the capabilities of this implementation, we developed four proof-of-concept applications for (1) prediction of free energy changes upon mutation; (2) high-resolution structural refinement; (3) protein-protein docking; and (4) assembly of symmetric protein complexes, all in the membrane environment. Preliminary data show that these algorithms can produce meaningful scores and structures. The data also suggest needed improvements to both sampling routines and score functions. Importantly, the applications collectively demonstrate the potential of combining the flexible nature of RosettaMP with the power of Rosetta algorithms to facilitate membrane protein modeling and design.

Conversion of a Dodecahedral Protein Capsid into Pentamers via Minimal Point Mutations

Protein self-assembly relies upon the formation of stabilizing noncovalent interactions across subunit interfaces. Identifying the determinants of self-assembly is crucial for understanding structure-function relationships in symmetric protein complexes and for engineering responsive nanoscale architectures for applications in medicine and biotechnology. Lumazine synthases (LS's) comprise a protein family that forms diverse quaternary structures, including pentamers and 60-subunit dodecahedral capsids. To improve our understanding of the basis for this difference in assembly, we attempted to convert the capsid-forming LS from Aquifex aeolicus (AaLS) into pentamers through a small number of rationally designed amino acid substitutions. Our mutations targeted side chains at ionic (R40), hydrogen bonding (H41), and hydrophobic (L121 and I125) interaction sites along the interfaces between pentamers. We found that substitutions at two or three of these positions could reliably generate pentameric variants of AaLS. Biophysical characterization indicates that this quaternary structure change is not accompanied by substantial changes in secondary or tertiary structure. Interestingly, previous homology-based studies of the assembly determinants in LS's had identified only one of these four positions. The ability to control assembly state in protein capsids such as AaLS could aid efforts in the development of new systems for drug delivery, biocatalysis, or materials synthesis.

A Geometric Arrangement Algorithm for Structure Determination of Symmetric Protein Homo-Oligomers from NOEs and RDCs

Nuclear magnetic resonance (NMR) spectroscopy is a primary tool to perform structural studies of proteins in physiologically-relevant solution conditions. Restraints on distances between pairs of nuclei in the protein, derived from the nuclear Overhauser effect (NOE), provide information about the structure of the protein in its folded state. NMR studies of symmetric protein homo-oligomers present a unique challenge. Using X-filtered NOESY experiments, it is possible to determine whether an NOE restrains a pair of protons across different subunits or within a single subunit, but current experimental techniques are unable to determine in which subunits the restrained protons lie. Consequently, it is difficult to assign NOEs to particular pairs of subunits with certainty, thus hindering the structural analysis of the oligomeric state. Computational approaches are needed to address this subunit ambiguity, but traditional solutions often rely on stochastic search coupled with simulated annealing and simulations of simplified molecular dynamics, which have many tunable parameters that must be chosen carefully and can also fail to report structures consistent with the experimental restraints. In addition, these traditional approaches rarely provide guarantees on running time or solution quality. We reduce the structure determination of homo-oligomers with cyclic symmetry to computing geometric arrangements of unions of annuli in a plane. Our algorithm, disco, runs in expected O(n²) time, where n is the number of distance restraints, potentially assigned ambiguously. disco is guaranteed to report the exact set of oligomer structures consistent with the distance restraints and also with orientational restraints from residual dipolar couplings (RDCs). We demonstrate our method using two symmetric protein complexes: the trimeric E. coli diacylglycerol kinase (DAGK) and a dimeric mutant of the immunoglobulin-binding domain B1 of streptococcal protein G (GB1). In both cases, disco computes oligomer structures with high precision and also finds distance restraints that are either mutually inconsistent or inconsistent with the RDCs. The entire protocol DISCO has been completely automated in a software package that is freely available and open-source at www.cs.duke.edu/donaldlab/software.php.

The 2010 Rosetta Developers Meeting: Macromolecular Prediction and Design Meets Reproducible Publishing

The 2010 Rosetta Developers Meeting: Macromolecular Prediction and Design Meets Reproducible Publishing

SymmRef: a flexible refinement method for symmetric multimers.

Symmetric protein complexes are abundant in the living cell. Predicting their atomic structure can shed light on the mechanism of many important biological processes. Symmetric docking methods aim to predict the structure of these complexes given the unbound structure of a single monomer, or its model. Symmetry constraints reduce the search-space of these methods and make the prediction easier compared to asymmetric protein-protein docking. However, the challenge of modeling the conformational changes that the monomer might undergo is a major obstacle. In this article, we present SymmRef, a novel method for refinement and reranking of symmetric docking solutions. The method models backbone and side-chain movements and optimizes the rigid-body orientations of the monomers. The backbone movements are modeled by normal modes minimization and the conformations of the side-chains are modeled by selecting optimal rotamers. Since solved structures of symmetric multimers show asymmetric side-chain conformations, we do not use symmetry constraints in the side-chain optimization procedure. The refined models are re-ranked according to an energy score. We tested the method on a benchmark of unbound docking challenges. The results show that the method significantly improves the accuracy and the ranking of symmetric rigid docking solutions. SymmRef is available for download at http:// bioinfo3d.cs.tau.ac.il/SymmRef/download.html.

An Automatic Particle Pickup Method Using a Neural Network Applicable to Low-Contrast Electron Micrographs

Three-dimensionalreconstruction from electron micrographs requires the selection of many single-particle projection images; more than 10 000 are generally required to obtain 5- to 10-Å structural resolution. Consequently, various automatic detection algorithms have been developed and successfully applied to large symmetric protein complexes. This paper presents a new automated particle recognition and pickup procedure based on the three-layer neural network that has a large application range than other automated procedures. Its use for both faint and noisy electron micrographs is demonstrated. The method requires only 200 selected particles as learning data and is able to detect images of proteins as small as 200 kDa.

Statistical significance of DNA sequence symmetries.

SOME regions of DNA involved in regulatory or catalytic interactions with proteins have sites with twofold rotation symmetry in the primary nucleotide sequence1–3. It would be worthwhile to determine whether certain of these symmetries occur so rarely by chance that they probably have a function, for example, as recognition sites for symmetric protein complexes. Likewise, it would be useful to be able to dismiss other apparent symmetries as probable statistical accidents of no functional importance. The following quantitative method can be used to evaluate the statistical significance of hyphenated (that is, incomplete) symmetries, where a more intuitive approach may be misleading. Since the probability of random occurrence of a given subsequence increases in almost direct proportion to the length of the sequence considered, a rigorous approach will become increasingly important in the evaluation of the significance of subsequence anomalies occurring in the much longer DNA sequences that will be determined in the future.