De novo design of pH-responsive self-assembling helical protein filaments

Biological evolution has led to precise and dynamic nanostructures that reconfigure in response to pH and other environmental conditions. However, designing micrometre-scale protein nanostructures that are environmentally responsive remains a challenge. Here we describe the de novo design of pH-responsive protein filaments built from subunits containing six or nine buried histidine residues that assemble into micrometre-scale, well-ordered fibres at neutral pH. The cryogenic electron microscopy structure of an optimized design is nearly identical to the computational design model for both the subunit internal geometry and the subunit packing into the fibre. Electron, fluorescent and atomic force microscopy characterization reveal a sharp and reversible transition from assembled to disassembled fibres over 0.3 pH units, and rapid fibre disassembly in less than 1 s following a drop in pH. The midpoint of the transition can be tuned by modulating buried histidine-containing hydrogen bond networks. Computational protein design thus provides a route to creating unbound nanomaterials that rapidly respond to small pH changes.

Accurate computational design of three-dimensional protein crystals

Protein crystallization plays a central role in structural biology. Despite this, the process of crystallization remains poorly understood and highly empirical, with crystal contacts, lattice packing arrangements and space group preferences being largely unpredictable. Programming protein crystallization through precisely engineered side-chain-side-chain interactions across protein-protein interfaces is an outstanding challenge. Here we develop a general computational approach for designing three-dimensional protein crystals with prespecified lattice architectures at atomic accuracy that hierarchically constrains the overall number of degrees of freedom of the system. We design three pairs of oligomers that can be individually purified, and upon mixing, spontaneously self-assemble into >100 µm three-dimensional crystals. The structures of these crystals are nearly identical to the computational design models, closely corresponding in both overall architecture and the specific protein-protein interactions. The dimensions of the crystal unit cell can be systematically redesigned while retaining the space group symmetry and overall architecture, and the crystals are extremely porous and highly stable. Our approach enables the computational design of protein crystals with high accuracy, and the designed protein crystals, which have both structural and assembly information encoded in their primary sequences, provide a powerful platform for biological materials engineering.

De novo design of protein structure and function with RFdiffusion

There has been considerable recent progress in designing new proteins using deep-learning methods1-9. Despite this progress, a general deep-learning framework for protein design that enables solution of a wide range of design challenges, including de novo binder design and design of higher-order symmetric architectures, has yet to be described. Diffusion models10,11 have had considerable success in image and language generative modelling but limited success when applied to protein modelling, probably due to the complexity of protein backbone geometry and sequence-structure relationships. Here we show that by fine-tuning the RoseTTAFold structure prediction network on protein structure denoising tasks, we obtain a generative model of protein backbones that achieves outstanding performance on unconditional and topology-constrained protein monomer design, protein binder design, symmetric oligomer design, enzyme active site scaffolding and symmetric motif scaffolding for therapeutic and metal-binding protein design. We demonstrate the power and generality of the method, called RoseTTAFold diffusion (RFdiffusion), by experimentally characterizing the structures and functions of hundreds of designed symmetric assemblies, metal-binding proteins and protein binders. The accuracy of RFdiffusion is confirmed by the cryogenic electron microscopy structure of a designed binder in complex with influenza haemagglutinin that is nearly identical to the design model. In a manner analogous to networks that produce images from user-specified inputs, RFdiffusion enables the design of diverse functional proteins from simple molecular specifications.

Fast and versatile sequence-independent protein docking for nanomaterials design using RPXDock

Computationally designed multi-subunit assemblies have shown considerable promise for a variety of applications, including a new generation of potent vaccines. One of the major routes to such materials is rigid body sequence-independent docking of cyclic oligomers into architectures with point group or lattice symmetries. Current methods for docking and designing such assemblies are tailored to specific classes of symmetry and are difficult to modify for novel applications. Here we describe RPXDock, a fast, flexible, and modular software package for sequence-independent rigid-body protein docking across a wide range of symmetric architectures that is easily customizable for further development. RPXDock uses an efficient hierarchical search and a residue-pair transform (RPX) scoring method to rapidly search through multidimensional docking space. We describe the structure of the software, provide practical guidelines for its use, and describe the available functionalities including a variety of score functions and filtering tools that can be used to guide and refine docking results towards desired configurations.

De novo design of modular peptide-binding proteins by superhelical matching

General approaches for designing sequence-specific peptide-binding proteins would have wide utility in proteomics and synthetic biology. However, designing peptide-binding proteins is challenging, as most peptides do not have defined structures in isolation, and hydrogen bonds must be made to the buried polar groups in the peptide backbone1-3. Here, inspired by natural and re-engineered protein-peptide systems4-11, we set out to design proteins made out of repeating units that bind peptides with repeating sequences, with a one-to-one correspondence between the repeat units of the protein and those of the peptide. We use geometric hashing to identify protein backbones and peptide-docking arrangements that are compatible with bidentate hydrogen bonds between the side chains of the protein and the peptide backbone12. The remainder of the protein sequence is then optimized for folding and peptide binding. We design repeat proteins to bind to six different tripeptide-repeat sequences in polyproline II conformations. The proteins are hyperstable and bind to four to six tandem repeats of their tripeptide targets with nanomolar to picomolar affinities in vitro and in living cells. Crystal structures reveal repeating interactions between protein and peptide interactions as designed, including ladders of hydrogen bonds from protein side chains to peptide backbones. By redesigning the binding interfaces of individual repeat units, specificity can be achieved for non-repeating peptide sequences and for disordered regions of native proteins.

Improving the secretion of designed protein assemblies through negative design of cryptic transmembrane domains

Computationally designed protein nanoparticles have recently emerged as a promising platform for the development of new vaccines and biologics. For many applications, secretion of designed nanoparticles from eukaryotic cells would be advantageous, but in practice, they often secrete poorly. Here we show that designed hydrophobic interfaces that drive nanoparticle assembly are often predicted to form cryptic transmembrane domains, suggesting that interaction with the membrane insertion machinery could limit efficient secretion. We develop a general computational protocol, the Degreaser, to design away cryptic transmembrane domains without sacrificing protein stability. The retroactive application of the Degreaser to previously designed nanoparticle components and nanoparticles considerably improves secretion, and modular integration of the Degreaser into design pipelines results in new nanoparticles that secrete as robustly as naturally occurring protein assemblies. Both the Degreaser protocol and the nanoparticles we describe may be broadly useful in biotechnological applications.

Exploration of Structured Symmetric Cyclic Peptides as Ligands for Metal-Organic Frameworks

Despite remarkable advances in the assembly of highly structured coordination polymers and metal-organic frameworks, the rational design of such materials using more conformationally flexible organic ligands such as peptides remains challenging. In an effort to make the design of such materials fully programmable, we first developed a computational design method for generating metal-mediated 3D frameworks using rigid and symmetric peptide macrocycles with metal-coordinating sidechains. We solved the structures of six crystalline networks involving conformationally constrained 6 to 12 residue cyclic peptides with C2, C3, and S2 internal symmetry and three different types of metals (Zn2+, Co2+, or Cu2+) by single-crystal X-ray diffraction, which reveals how the peptide sequences, backbone symmetries, and metal coordination preferences drive the assembly of the resulting structures. In contrast to smaller ligands, these peptides associate through peptide-peptide interactions without full coordination of the metals, contrary to one of the assumptions underlying our computational design method. The cyclic peptides are the largest peptidic ligands reported to form crystalline coordination polymers with transition metals to date, and while more work is required to develop methods for fully programming their crystal structures, the combination of high chemical diversity with synthetic accessibility makes them attractive building blocks for engineering a broader set of new crystalline materials for use in applications such as sensing, asymmetric catalysis, and chiral separation.

Design of protein-binding proteins from the target structure alone

The design of proteins that bind to a specific site on the surface of a target protein using no information other than the three-dimensional structure of the target remains a challenge1-5. Here we describe a general solution to this problem that starts with a broad exploration of the vast space of possible binding modes to a selected region of a protein surface, and then intensifies the search in the vicinity of the most promising binding modes. We demonstrate the broad applicability of this approach through the de novo design of binding proteins to 12 diverse protein targets with different shapes and surface properties. Biophysical characterization shows that the binders, which are all smaller than 65 amino acids, are hyperstable and, following experimental optimization, bind their targets with nanomolar to picomolar affinities. We succeeded in solving crystal structures of five of the binder-target complexes, and all five closely match the corresponding computational design models. Experimental data on nearly half a million computational designs and hundreds of thousands of point mutants provide detailed feedback on the strengths and limitations of the method and of our current understanding of protein-protein interactions, and should guide improvements of both. Our approach enables the targeted design of binders to sites of interest on a wide variety of proteins for therapeutic and diagnostic applications.

Computational design of mechanically coupled axle-rotor protein assemblies

Natural molecular machines contain protein components that undergo motion relative to each other. Designing such mechanically constrained nanoscale protein architectures with internal degrees of freedom is an outstanding challenge for computational protein design. Here we explore the de novo construction of protein machinery from designed axle and rotor components with internal cyclic or dihedral symmetry. We find that the axle-rotor systems assemble in vitro and in vivo as designed. Using cryo-electron microscopy, we find that these systems populate conformationally variable relative orientations reflecting the symmetry of the coupled components and the computationally designed interface energy landscape. These mechanical systems with internal degrees of freedom are a step toward the design of genetically encodable nanomachines.

Design of multi-scale protein complexes by hierarchical building block fusion

A systematic and robust approach to generating complex protein nanomaterials would have broad utility. We develop a hierarchical approach to designing multi-component protein assemblies from two classes of modular building blocks: designed helical repeat proteins (DHRs) and helical bundle oligomers (HBs). We first rigidly fuse DHRs to HBs to generate a large library of oligomeric building blocks. We then generate assemblies with cyclic, dihedral, and point group symmetries from these building blocks using architecture guided rigid helical fusion with new software named WORMS. X-ray crystallography and cryo-electron microscopy characterization show that the hierarchical design approach can accurately generate a wide range of assemblies, including a 43 nm diameter icosahedral nanocage. The computational methods and building block sets described here provide a very general route to de novo designed protein nanomaterials.

Designed proteins assemble antibodies into modular nanocages

Multivalent display of receptor-engaging antibodies or ligands can enhance their activity. Instead of achieving multivalency by attachment to preexisting scaffolds, here we unite form and function by the computational design of nanocages in which one structural component is an antibody or Fc-ligand fusion and the second is a designed antibody-binding homo-oligomer that drives nanocage assembly. Structures of eight nanocages determined by electron microscopy spanning dihedral, tetrahedral, octahedral, and icosahedral architectures with 2, 6, 12, and 30 antibodies per nanocage, respectively, closely match the corresponding computational models. Antibody nanocages targeting cell surface receptors enhance signaling compared with free antibodies or Fc-fusions in death receptor 5 (DR5)-mediated apoptosis, angiopoietin-1 receptor (Tie2)-mediated angiogenesis, CD40 activation, and T cell proliferation. Nanocage assembly also increases severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pseudovirus neutralization by α-SARS-CoV-2 monoclonal antibodies and Fc-angiotensin-converting enzyme 2 (ACE2) fusion proteins.

Design of biologically active binary protein 2D materials

Ordered two-dimensional arrays such as S-layers^1,2 and designed analogues^3-5 have intrigued bioengineers^6,7, but with the exception of a single lattice formed with flexible linkers⁸, they are constituted from just one protein component. Materials composed of two components have considerable potential advantages for modulating assembly dynamics and incorporating more complex functionality^9-12. Here we describe a computational method to generate co-assembling binary layers by designing rigid interfaces between pairs of dihedral protein building blocks, and use it to design a p6m lattice. The designed array components are soluble at millimolar concentrations, but when combined at nanomolar concentrations, they rapidly assemble into nearly crystalline micrometre-scale arrays nearly identical to the computational design model in vitro and in cells without the need for a two-dimensional support. Because the material is designed from the ground up, the components can be readily functionalized and their symmetry reconfigured, enabling formation of ligand arrays with distinguishable surfaces, which we demonstrate can drive extensive receptor clustering, downstream protein recruitment and signalling. Using atomic force microscopy on supported bilayers and quantitative microscopy on living cells, we show that arrays assembled on membranes have component stoichiometry and structure similar to arrays formed in vitro, and that our material can therefore impose order onto fundamentally disordered substrates such as cell membranes. In contrast to previously characterized cell surface receptor binding assemblies such as antibodies and nanocages, which are rapidly endocytosed, we find that large arrays assembled at the cell surface suppress endocytosis in a tunable manner, with potential therapeutic relevance for extending receptor engagement and immune evasion. Our work provides a foundation for a synthetic cell biology in which multi-protein macroscale materials are designed to modulate cell responses and reshape synthetic and living systems.

Designed proteins assemble antibodies into modular nanocages

Antibodies are widely used in biology and medicine, and there has been considerable interest in multivalent antibody formats to increase binding avidity and enhance signaling pathway agonism. However, there are currently no general approaches for forming precisely oriented antibody assemblies with controlled valency. We describe the computational design of two-component nanocages that overcome this limitation by uniting form and function. One structural component is any antibody or Fc fusion and the second is a designed Fc-binding homo-oligomer that drives nanocage assembly. Structures of 8 antibody nanocages determined by electron microscopy spanning dihedral, tetrahedral, octahedral, and icosahedral architectures with 2, 6, 12, and 30 antibodies per nanocage match the corresponding computational models. Antibody nanocages targeting cell-surface receptors enhance signaling compared to free antibodies or Fc-fusions in DR5-mediated apoptosis, Tie2-mediated angiogenesis, CD40 activation, and T cell proliferation; nanocage assembly also increases SARS-CoV-2 pseudovirus neutralization by α-SARS-CoV-2 monoclonal antibodies and Fc-ACE2 fusion proteins. We anticipate that the ability to assemble arbitrary antibodies without need for covalent modification into highly ordered assemblies with different geometries and valencies will have broad impact in biology and medicine.

Tailored design of protein nanoparticle scaffolds for multivalent presentation of viral glycoprotein antigens

Multivalent presentation of viral glycoproteins can substantially increase the elicitation of antigen-specific antibodies. To enable a new generation of anti-viral vaccines, we designed self-assembling protein nanoparticles with geometries tailored to present the ectodomains of influenza, HIV, and RSV viral glycoprotein trimers. We first de novo designed trimers tailored for antigen fusion, featuring N-terminal helices positioned to match the C termini of the viral glycoproteins. Trimers that experimentally adopted their designed configurations were incorporated as components of tetrahedral, octahedral, and icosahedral nanoparticles, which were characterized by cryo-electron microscopy and assessed for their ability to present viral glycoproteins. Electron microscopy and antibody binding experiments demonstrated that the designed nanoparticles presented antigenically intact prefusion HIV-1 Env, influenza hemagglutinin, and RSV F trimers in the predicted geometries. This work demonstrates that antigen-displaying protein nanoparticles can be designed from scratch, and provides a systematic way to investigate the influence of antigen presentation geometry on the immune response to vaccination.

Efficient flexible backbone protein-protein docking for challenging targets

Motivation

Binding-induced conformational changes challenge current computational docking algorithms by exponentially increasing the conformational space to be explored. To restrict this search to relevant space, some computational docking algorithms exploit the inherent flexibility of the protein monomers to simulate conformational selection from pre-generated ensembles. As the ensemble size expands with increased flexibility, these methods struggle with efficiency and high false positive rates.

Results

Here, we develop and benchmark RosettaDock 4.0, which efficiently samples large conformational ensembles of flexible proteins and docks them using a novel, six-dimensional, coarse-grained score function. A strong discriminative ability allows an eight-fold higher enrichment of near-native candidate structures in the coarse-grained phase compared to RosettaDock 3.2. It adaptively samples 100 conformations each of the ligand and the receptor backbone while increasing computational time by only 20-80%. In local docking of a benchmark set of 88 proteins of varying degrees of flexibility, the expected success rate (defined as cases with ≥50% chance of achieving 3 near-native structures in the 5 top-ranked ones) for blind predictions after resampling is 77% for rigid complexes, 49% for moderately flexible complexes and 31% for highly flexible complexes. These success rates on flexible complexes are a substantial step forward from all existing methods. Additionally, for highly flexible proteins, we demonstrate that when a suitable conformer generation method exists, the method successfully docks the complex.

Availability and implementation

As a part of the Rosetta software suite, RosettaDock 4.0 is available at https://www.rosettacommons.org to all non-commercial users for free and to commercial users for a fee.

Supplementary information

Supplementary data are available at Bioinformatics online.

Multi-input chemical control of protein dimerization for programming graded cellular responses

Chemical and optogenetic methods for post-translationally controlling protein function have enabled modulation and engineering of cellular functions. However, most of these methods only confer single-input, single-output control. To increase the diversity of post-translational behaviors that can be programmed, we built a system based on a single protein receiver that can integrate multiple drug inputs, including approved therapeutics. Our system translates drug inputs into diverse outputs using a suite of engineered reader proteins to provide variable dimerization states of the receiver protein. We show that our single receiver protein architecture can be used to program a variety of cellular responses, including graded and proportional dual-output control of transcription and mammalian cell signaling. We apply our tools to titrate the competing activities of the Rac and Rho GTPases to control cell morphology. Our versatile tool set will enable researchers to post-translationally program mammalian cellular processes and to engineer cell therapies.

Enhancing and shaping the immunogenicity of native-like HIV-1 envelope trimers with a two-component protein nanoparticle

The development of native-like HIV-1 envelope (Env) trimer antigens has enabled the induction of neutralizing antibody (NAb) responses against neutralization-resistant HIV-1 strains in animal models. However, NAb responses are relatively weak and narrow in specificity. Displaying antigens in a multivalent fashion on nanoparticles (NPs) is an established strategy to increase their immunogenicity. Here we present the design and characterization of two-component protein NPs displaying 20 stabilized SOSIP trimers from various HIV-1 strains. The two-component nature permits the incorporation of exclusively well-folded, native-like Env trimers into NPs that self-assemble in vitro with high efficiency. Immunization studies show that the NPs are particularly efficacious as priming immunogens, improve the quality of the Ab response over a conventional one-component nanoparticle system, and are most effective when SOSIP trimers with an apex-proximate neutralizing epitope are displayed. Their ability to enhance and shape the immunogenicity of SOSIP trimers make these NPs a promising immunogen platform.

De novo design of self-assembling helical protein filaments

We describe a general computational approach to designing self-assembling helical filaments from monomeric proteins and use this approach to design proteins that assemble into micrometer-scale filaments with a wide range of geometries in vivo and in vitro. Cryo-electron microscopy structures of six designs are close to the computational design models. The filament building blocks are idealized repeat proteins, and thus the diameter of the filaments can be systematically tuned by varying the number of repeat units. The assembly and disassembly of the filaments can be controlled by engineered anchor and capping units built from monomers lacking one of the interaction surfaces. The ability to generate dynamic, highly ordered structures that span micrometers from protein monomers opens up possibilities for the fabrication of new multiscale metamaterials.

De novo design of a fluorescence-activating β-barrel

The regular arrangements of β-strands around a central axis in β-barrels and of α-helices in coiled coils contrast with the irregular tertiary structures of most globular proteins, and have fascinated structural biologists since they were first discovered. Simple parametric models have been used to design a wide range of α-helical coiled-coil structures, but to date there has been no success with β-barrels. Here we show that accurate de novo design of β-barrels requires considerable symmetry-breaking to achieve continuous hydrogen-bond connectivity and eliminate backbone strain. We then build ensembles of β-barrel backbone models with cavity shapes that match the fluorogenic compound DFHBI, and use a hierarchical grid-based search method to simultaneously optimize the rigid-body placement of DFHBI in these cavities and the identities of the surrounding amino acids to achieve high shape and chemical complementarity. The designs have high structural accuracy and bind and fluorescently activate DFHBI in vitro and in Escherichia coli, yeast and mammalian cells. This de novo design of small-molecule binding activity, using backbones custom-built to bind the ligand, should enable the design of increasingly sophisticated ligand-binding proteins, sensors and catalysts that are not limited by the backbone geometries available in known protein structures.

Accurate computational design of multipass transmembrane proteins

The computational design of transmembrane proteins with more than one membrane-spanning region remains a major challenge. We report the design of transmembrane monomers, homodimers, trimers, and tetramers with 76 to 215 residue subunits containing two to four membrane-spanning regions and up to 860 total residues that adopt the target oligomerization state in detergent solution. The designed proteins localize to the plasma membrane in bacteria and in mammalian cells, and magnetic tweezer unfolding experiments in the membrane indicate that they are very stable. Crystal structures of the designed dimer and tetramer-a rocket-shaped structure with a wide cytoplasmic base that funnels into eight transmembrane helices-are very close to the design models. Our results pave the way for the design of multispan membrane proteins with new functions.

Computational design of self-assembling cyclic protein homo-oligomers

Self-assembling cyclic protein homo-oligomers play important roles in biology and the ability to generate custom homo-oligomeric structures could enable new approaches to probe biological function. Here we report a general approach to design cyclic homo-oligomers that employs a new residue pair transform method for assessing the design ability of a protein-protein interface. This method is sufficiently rapid to enable systematic enumeration of cyclically docked arrangements of a monomer followed by sequence design of the newly formed interfaces. We use this method to design interfaces onto idealized repeat proteins that direct their assembly into complexes that possess cyclic symmetry. Of 96 designs that were experimentally characterized, 21 were found to form stable monodisperse homo-oligomers in solution, and 15 (4 homodimers, 6 homotrimers, 6 homotetramers and 1 homopentamer) had solution small angle X-ray scattering data consistent with the design models. X-ray crystal structures were obtained for five of the designs and each of these were shown to be very close to their design model.

Accurate design of megadalton-scale two-component icosahedral protein complexes

Nature provides many examples of self- and co-assembling protein-based molecular machines, including icosahedral protein cages that serve as scaffolds, enzymes, and compartments for essential biochemical reactions and icosahedral virus capsids, which encapsidate and protect viral genomes and mediate entry into host cells. Inspired by these natural materials, we report the computational design and experimental characterization of co-assembling, two-component, 120-subunit icosahedral protein nanostructures with molecular weights (1.8 to 2.8 megadaltons) and dimensions (24 to 40 nanometers in diameter) comparable to those of small viral capsids. Electron microscopy, small-angle x-ray scattering, and x-ray crystallography show that 10 designs spanning three distinct icosahedral architectures form materials closely matching the design models. In vitro assembly of icosahedral complexes from independently purified components occurs rapidly, at rates comparable to those of viral capsids, and enables controlled packaging of molecular cargo through charge complementarity. The ability to design megadalton-scale materials with atomic-level accuracy and controllable assembly opens the door to a new generation of genetically programmable protein-based molecular machines.

Accelerated electron transport from photosystem I to redox partners by covalently linked ferredoxin

Photosystem I is a highly efficient and potent light-induced reductase that is considered to be an appealing target for integration into hybrid solar fuel production systems. However, rapid transport of multiple electrons from the reducing end of photosystem I to downstream processes in vivo is limited by the diffusion of its native redox partner ferredoxin that is a single electron carrier. Here, we describe the design and construction of a faster electron transfer interface based on anchoring ferredoxin to the reducing end of photosystem I thereby confining the diffusion space of ferredoxin to the near vicinity of its photosystem I binding and reduction site. This was achieved by fusing ferredoxin to the PsaE subunit of photosystem I by a flexible peptide linker and reconstituting PSI in vitro with the new fusion protein. A computational algorithm was developed in order to determine the optimal linker length that will confine ferredoxin to the vicinity of photosystem I's reducing end without restricting the formation of electron transfer complexes. According to the calculation, we reconstituted photosystem I with three fusion proteins comprising PsaE and ferredoxin separated by linkers of different lengths, namely 14, 19, and 25 amino acids, and tested their effect on electron transfer rates from photosystem I to downstream processes. Indeed, we found a significant enhancement of light dependent NADPH synthesis using photosystems containing the PsaE-ferredoxin fusion proteins, equivalent to a ten-fold increase in soluble ferredoxin concentration. We propose that such a system could be used for other ferredoxin dependent redox reactions, such as the enzymatic production of hydrogen, a promising alternative fuel. As the system is comprised entirely of natural amino acids and biological cofactors, it could be integrated into the energy conversion apparatus of photosynthetic organisms by genetic engineering.

Computational design of self-assembling protein nanomaterials with atomic level accuracy

We describe a general computational method for designing proteins that self-assemble to a desired symmetric architecture. Protein building blocks are docked together symmetrically to identify complementary packing arrangements, and low-energy protein-protein interfaces are then designed between the building blocks in order to drive self-assembly. We used trimeric protein building blocks to design a 24-subunit, 13-nm diameter complex with octahedral symmetry and a 12-subunit, 11-nm diameter complex with tetrahedral symmetry. The designed proteins assembled to the desired oligomeric states in solution, and the crystal structures of the complexes revealed that the resulting materials closely match the design models. The method can be used to design a wide variety of self-assembling protein nanomaterials.

Modeling disordered regions in proteins using Rosetta

Protein structure prediction methods such as Rosetta search for the lowest energy conformation of the polypeptide chain. However, the experimentally observed native state is at a minimum of the free energy, rather than the energy. The neglect of the missing configurational entropy contribution to the free energy can be partially justified by the assumption that the entropies of alternative folded states, while very much less than unfolded states, are not too different from one another, and hence can be to a first approximation neglected when searching for the lowest free energy state. The shortcomings of current structure prediction methods may be due in part to the breakdown of this assumption. Particularly problematic are proteins with significant disordered regions which do not populate single low energy conformations even in the native state. We describe two approaches within the Rosetta structure modeling methodology for treating such regions. The first does not require advance knowledge of the regions likely to be disordered; instead these are identified by minimizing a simple free energy function used previously to model protein folding landscapes and transition states. In this model, residues can be either completely ordered or completely disordered; they are considered disordered if the gain in entropy outweighs the loss of favorable energetic interactions with the rest of the protein chain. The second approach requires identification in advance of the disordered regions either from sequence alone using for example the DISOPRED server or from experimental data such as NMR chemical shifts. During Rosetta structure prediction calculations the disordered regions make only unfavorable repulsive contributions to the total energy. We find that the second approach has greater practical utility and illustrate this with examples from de novo structure prediction, NMR structure calculation, and comparative modeling.

RosettaHoles2: a volumetric packing measure for protein structure refinement and validation

We present an improved version of RosettaHoles, a methodology for quantitative and visual characterization of protein core packing. RosettaHoles2 features a packing measure more rapidly computable, accurate and physically transparent, as well as a new validation score intended for structures submitted to the Protein Data Bank. The differential packing measure is parameterized to maximize the gap between computationally generated and experimentally determined X-ray structures, and can be used in refinement of protein structure models. The parameters of the model provide insight into components missing in current force fields, and the validation score gives an upper bound on the X-ray resolution of Protein Data Bank structures; a crystal structure should have a validation score as good as or better than its resolution.