Computational design of a symmetric homodimer using β-strand assembly

Computational design of novel protein–protein interfaces is a test of our understanding of protein interactions and has the potential to allow modification of cellular physiology. Methods for designing high-affinity interactions that adopt a predetermined binding mode have proved elusive, suggesting the need for new strategies that simplify the design process. A solvent-exposed backbone on a β-strand is thought of as “sticky” and β-strand pairing stabilizes many naturally occurring protein complexes. Here, we computationally redesign a monomeric protein to form a symmetric homodimer by pairing exposed β-strands to form an intermolecular β-sheet. A crystal structure of the designed complex closely matches the computational model (rmsd = 1.0 Å). This work demonstrates that β-strand pairing can be used to computationally design new interactions with high accuracy.     

De novo structure prediction and experimental characterization of folded peptoid oligomers

Peptoid molecules are biomimetic oligomers that can fold into unique three-dimensional structures. As part of an effort to advance computational design of folded oligomers, we present blind-structure predictions for three peptoid sequences using a combination of Replica Exchange Molecular Dynamics (REMD) simulation and Quantum Mechanical refinement. We correctly predicted the structure of a N-aryl peptoid trimer to within 0.2 Å rmsd-backbone and a cyclic peptoid nonamer to an accuracy of 1.0 Å rmsd-backbone. X-ray crystallographic structures are presented for a linear N-alkyl peptoid trimer and for the cyclic peptoid nonamer. The peptoid macrocycle structure features a combination of cis and trans backbone amides, significant nonplanarity of the amide bonds, and a unique “basket” arrangement of (S)-N(1-phenylethyl) side chains encompassing a bound ethanol molecule. REMD simulations of the peptoid trimers reveal that well folded peptoids can exhibit funnel-like conformational free energy landscapes similar to those for ordered polypeptides. These results indicate that physical modeling can successfully perform de novo structure prediction for small peptoid molecules.     

Accurate de novo structure prediction of large transmembrane protein domains using fragment-assembly and correlated mutation analysis

A new de novo protein structure prediction method for transmembrane proteins (FILM3) is described that is able to accurately predict the structures of large membrane proteins domains using an ensemble of two secondary structure prediction methods to guide fragment selection in combination with a scoring function based solely on correlated mutations detected in multiple sequence alignments. This approach has been validated by generating models for 28 membrane proteins with a diverse range of complex topologies and an average length of over 300 residues with results showing that TM-scores > 0.5 can be achieved in almost every case following refinement using MODELLER. In one of the most impressive results, a model of mitochondrial cytochrome c oxidase polypeptide I was obtained with a TM-score > 0.75 and an rmsd of only 5.7 Å over all 514 residues. These results suggest that FILM3 could be applicable to a wide range of transmembrane proteins of as-yet-unknown 3D structure given sufficient homologous sequences.     

Lipidic cubic phases: A novel concept for the crystallization of membrane proteins

Understanding the mechanisms of action of membrane proteins requires the elucidation of their structures to high resolution. The critical step in accomplishing this by x-ray crystallography is the routine availability of well-ordered three-dimensional crystals. We have devised a novel, rational approach to meet this goal using quasisolid lipidic cubic phases. This membrane system, consisting of lipid, water, and protein in appropriate proportions, forms a structured, transparent, and complex three-dimensional lipidic array, which is pervaded by an intercommunicating aqueous channel system. Such matrices provide nucleation sites (“seeding”) and support growth by lateral diffusion of protein molecules in the membrane (“feeding”). Bacteriorhodopsin crystals were obtained from bicontinuous cubic phases, but not from micellar systems, implying a critical role of the continuity of the diffusion space (the bilayer) on crystal growth. Hexagonal bacteriorhodopsin crystals diffracted to 3.7 Å resolution, with a space group P6₃, and unit cell dimensions of a = b = 62 Å, c = 108 Å; α = β = 90° and γ = 120°.     

Accurate protein structure modeling using sparse NMR data and homologous structure information

While information from homologous structures plays a central role in X-ray structure determination by molecular replacement, such information is rarely used in NMR structure determination because it can be incorrect, both locally and globally, when evolutionary relationships are inferred incorrectly or there has been considerable evolutionary structural divergence. Here we describe a method that allows robust modeling of protein structures of up to 225 residues by combining , ¹³C, and ¹⁵N backbone and ¹³Cβ chemical shift data, distance restraints derived from homologous structures, and a physically realistic all-atom energy function. Accurate models are distinguished from inaccurate models generated using incorrect sequence alignments by requiring that (i) the all-atom energies of models generated using the restraints are lower than models generated in unrestrained calculations and (ii) the low-energy structures converge to within 2.0 Å backbone rmsd over 75% of the protein. Benchmark calculations on known structures and blind targets show that the method can accurately model protein structures, even with very remote homology information, to a backbone rmsd of 1.2–1.9 Å relative to the conventional determined NMR ensembles and of 0.9–1.6 Å relative to X-ray structures for well-defined regions of the protein structures. This approach facilitates the accurate modeling of protein structures using backbone chemical shift data without need for side-chain resonance assignments and extensive analysis of NOESY cross-peak assignments.     

Ultrahigh resolution protein structures using NMR chemical shift tensors

NMR chemical shift tensors (CSTs) in proteins, as well as their orientations, represent an important new restraint class for protein structure refinement and determination. Here, we present the first determination of both CST magnitudes and orientations for ¹³Cα and ¹⁵N (peptide backbone) groups in a protein, the β1 IgG binding domain of protein G from Streptococcus spp., GB1. Site-specific ¹³Cα and ¹⁵N CSTs were measured using synchronously evolved recoupling experiments in which ¹³C and ¹⁵N tensors were projected onto the ¹H-¹³C and ¹H-¹⁵N vectors, respectively, and onto the ¹⁵N-¹³C vector in the case of ¹³Cα. The orientations of the ¹³Cα CSTs to the ¹H-¹³C and ¹³C-¹⁵N vectors agreed well with the results of ab initio calculations, with an rmsd of approximately 8°. In addition, the measured ¹⁵N tensors exhibited larger reduced anisotropies in α-helical versus β-sheet regions, with very limited variation (18 ± 4°) in the orientation of the z-axis of the ¹⁵N CST with respect to the ¹H-¹⁵N vector. Incorporation of the ¹³Cα CST restraints into structure calculations, in combination with isotropic chemical shifts, transferred echo double resonance ¹³C-¹⁵N distances and vector angle restraints, improved the backbone rmsd to 0.16 Å (PDB ID code 2LGI) and is consistent with existing X-ray structures (0.51 Å agreement with PDB ID code 2QMT). These results demonstrate that chemical shift tensors have considerable utility in protein structure refinement, with the best structures comparable to 1.0-Å crystal structures, based upon empirical metrics such as Ramachandran geometries and χ¹/χ² distributions, providing solid-state NMR with a powerful tool for de novo structure determination.     

Visualizing key hinges and a potential major source of compliance in the lever arm of myosin

We have determined the 2.3-Å-resolution crystal structure of a myosin light chain domain, corresponding to one type found in sea scallop catch (“smooth”) muscle. This structure reveals hinges that may function in the “on” and “off” states of myosin. The molecule adopts two different conformations about the heavy chain “hook” and regulatory light chain (RLC) helix D. This conformational change results in extended and compressed forms of the lever arm whose lengths differ by 10 Å. The heavy chain hook and RLC helix D hinges could thus serve as a potential major and localized source of cross-bridge compliance during the contractile cycle. In addition, in one of the molecules of the crystal, part of the RLC N-terminal extension is seen in atomic detail and forms a one-turn alpha-helix that interacts with RLC helix D. This extension, whose sequence is highly variable in different myosins, may thus modulate the flexibility of the lever arm. Moreover, the relative proximity of the phosphorylation site to the helix D hinge suggests a potential role for conformational changes about this hinge in the transition between the on and off states of regulated myosins.     

Structural dependence of HET-s amyloid fibril infectivity assessed by cryoelectron microscopy

HET-s is a prion protein of the fungus Podospora anserina which, in the prion state, is active in a self/nonself recognition process called heterokaryon incompatibility. Its prionogenic properties reside in the C-terminal “prion domain.” The HET-s prion domain polymerizes in vitro into amyloid fibrils whose properties depend on the pH of assembly; above pH 3, infectious singlet fibrils are produced, and below pH 3, noninfectious triplet fibrils. To investigate the correlation between structure and infectivity, we performed cryo-EM analyses. Singlet fibrils have a helical pitch of approximately 410 Å and a left-handed twist. Triplet fibrils have three protofibrils whose lateral dimensions (36 × 25 Å) and axial packing (one subunit per 9.4 Å) match those of singlets but differ in their supercoiling. At 8.5-Å resolution, the cross-section of the singlet fibril reconstruction is largely consistent with that of a β-solenoid model previously determined by solid-state NMR. Reconstructions of the triplet fibrils show three protofibrils coiling around a common axis and packed less tightly at pH 3 than at pH 2, eventually peeling off. Taken together with the earlier observation that fragmentation of triplet fibrils by sonication does not increase infectivity, these observations suggest a novel mechanism for self-propagation, whereby daughter fibrils nucleate on the lateral surface of singlet fibrils. In triplets, this surface is occluded, blocking nucleation and thereby explaining their lack of infectivity.     

Chemical synthesis and X-ray structure of a heterochiral {D-protein antagonist plus vascular endothelial growth factor} protein complex by racemic crystallography

Total chemical synthesis was used to prepare the mirror image (D-protein) form of the angiogenic protein vascular endothelial growth factor (VEGF-A). Phage display against D-VEGF-A was used to screen designed libraries based on a unique small protein scaffold in order to identify a high affinity ligand. Chemically synthesized D- and L- forms of the protein ligand showed reciprocal chiral specificity in surface plasmon resonance binding experiments: The L-protein ligand bound only to D-VEGF-A, whereas the D-protein ligand bound only to L-VEGF-A. The D-protein ligand, but not the L-protein ligand, inhibited the binding of natural VEGF₁₆₅ to the VEGFR1 receptor. Racemic protein crystallography was used to determine the high resolution X-ray structure of the heterochiral complex consisting of {D-protein antagonist + L-protein form ofVEGF-A}. Crystallization of a racemic mixture of these synthetic proteins in appropriate stoichiometry gave a racemic protein complex of more than 73 kDa containing six synthetic protein molecules. The structure of the complex was determined to a resolution of 1.6 Å. Detailed analysis of the interaction between the D-protein antagonist and the VEGF-A protein molecule showed that the binding interface comprised a contact surface area of approximately 800 Å² in accord with our design objectives, and that the D-protein antagonist binds to the same region of VEGF-A that interacts with VEGFR1-domain 2.     

Structural basis for substrate activation and regulation by cystathionine beta-synthase (CBS) domains in cystathionine {beta}-synthase

The catalytic potential for H₂S biogenesis and homocysteine clearance converge at the active site of cystathionine β-synthase (CBS), a pyridoxal phosphate-dependent enzyme. CBS catalyzes β-replacement reactions of either serine or cysteine by homocysteine to give cystathionine and water or H₂S, respectively. In this study, high-resolution structures of the full-length enzyme from Drosophila in which a carbanion (1.70 Å) and an aminoacrylate intermediate (1.55 Å) have been captured are reported. Electrostatic stabilization of the zwitterionic carbanion intermediate is afforded by the close positioning of an active site lysine residue that is initially used for Schiff base formation in the internal aldimine and later as a general base. Additional stabilizing interactions between active site residues and the catalytic intermediates are observed. Furthermore, the structure of the regulatory “energy-sensing” CBS domains, named after this protein, suggests a mechanism for allosteric activation by S-adenosylmethionine.     

Structure of a type III secretion needle at 7-Å resolution provides insights into its assembly and signaling mechanisms

Type III secretion systems of Gram-negative bacteria form injection devices that deliver effector proteins into eukaryotic cells during infection. They span both bacterial membranes and the extracellular space to connect with the host cell plasma membrane. Their extracellular portion is a needle-like, hollow tube that serves as a secretion conduit for effector proteins. The needle of Shigella flexneri is approximately 50-nm long and 7-nm thick and is made by the helical assembly of one protein, MxiH. We provide a 7-Å resolution 3D image reconstruction of the Shigella needle by electron cryomicroscopy, which resolves α-helices and a β-hairpin that has never been observed in the crystal and solution structures of needle proteins, including MxiH. An atomic model of the needle based on the 3D-density map, in comparison with that of the bacterial-flagellar filament, provides insights into how such a thin tubular structure is stably assembled by intricate intermolecular interactions. The map also illuminates how the needle-length control protein functions as a ruler within the central channel during export of MxiH for assembly at the distal end of the needle, and how the secretion-activation signal may be transduced through a conformational change of the needle upon host-cell contact.     

Nonplanar peptide bonds in proteins are common and conserved but not biased toward active sites

The planarity of peptide bonds is an assumption that underlies decades of theoretical modeling of proteins. Peptide bonds strongly deviating from planarity are considered very rare features of protein structure that occur for functional reasons. Here, empirical analyses of atomic-resolution protein structures reveal that trans peptide groups can vary by more than 25° from planarity and that the true extent of nonplanarity is underestimated even in 1.2 Å resolution structures. Analyses as a function of the φ,ψ-backbone dihedral angles show that the expected value deviates by ± 8° from planar as a systematic function of conformation, but that the large majority of variation in planarity depends on tertiary effects. Furthermore, we show that those peptide bonds in proteins that are most nonplanar, deviating by over 20° from planarity, are not strongly associated with active sites. Instead, highly nonplanar peptides are simply integral components of protein structure related to local and tertiary structural features that tend to be conserved among homologs. To account for the systematic φ,ψ-dependent component of nonplanarity, we present a conformation-dependent library that can be used in crystallographic refinement and predictive protein modeling.     

Structure of a bacterial cell surface decaheme electron conduit

Some bacterial species are able to utilize extracellular mineral forms of iron and manganese as respiratory electron acceptors. In Shewanella oneidensis this involves decaheme cytochromes that are located on the bacterial cell surface at the termini of trans-outer-membrane electron transfer conduits. The cell surface cytochromes can potentially play multiple roles in mediating electron transfer directly to insoluble electron sinks, catalyzing electron exchange with flavin electron shuttles or participating in extracellular intercytochrome electron exchange along “nanowire” appendages. We present a 3.2-Å crystal structure of one of these decaheme cytochromes, MtrF, that allows the spatial organization of the 10 hemes to be visualized for the first time. The hemes are organized across four domains in a unique crossed conformation, in which a staggered 65-Å octaheme chain transects the length of the protein and is bisected by a planar 45-Å tetraheme chain that connects two extended Greek key split β-barrel domains. The structure provides molecular insight into how reduction of insoluble substrate (e.g., minerals), soluble substrates (e.g., flavins), and cytochrome redox partners might be possible in tandem at different termini of a trifurcated electron transport chain on the cell surface.     

Functional modulation of a protein folding landscape via side-chain distortion

Ultrahigh-resolution (< 1.0 Å) structures have revealed unprecedented and unexpected details of molecular geometry, such as the deformation of aromatic rings from planarity. However, the functional utility of such energetically costly strain is unknown. The 0.83 Å structure of α-lytic protease (αLP) indicated that residues surrounding a conserved Phe side-chain dictate a rotamer which results in a ∼6° distortion along the side-chain, estimated to cost 4 kcal/mol. By contrast, in the closely related protease Streptomyces griseus Protease B (SGPB), the equivalent Phe adopts a different rotamer and is undistorted. Here, we report that the αLP Phe side-chain distortion is both functional and conserved in proteases with large pro regions. Sequence analysis of the αLP serine protease family reveals a bifurcation separating those sequences expected to induce distortion and those that would not, which correlates with the extent of kinetic stability. Structural and folding kinetics analyses of family members suggest that distortion of this side-chain plays a role in increasing kinetic stability within the αLP family members that use a large Pro region. Additionally, structural and kinetic folding studies of mutants demonstrate that strain alters the folding free energy landscape by destabilizing the transition state (TS) relative to the native state (N). Although side-chain distortion comes at a cost of foldability, it suppresses the rate of unfolding, thereby enhancing kinetic stability and increasing protein longevity under harsh extracellular conditions. This ability of a structural distortion to enhance function is unlikely to be unique to αLP family members and may be relevant in other proteins exhibiting side-chain distortions.     

A designed protein as experimental model of primordial folding

How do proteins accomplish folding during early evolution? Theoretically the mechanism involves the selective stabilization of the native structure against all other competing compact conformations in a process that involves cumulative changes in the amino acid sequence along geological timescales. Thus, an evolved protein folds into a single structure at physiological temperature, but the conformational competition remains latent. For natural proteins such competition should emerge only near cryogenic temperatures, which places it beyond experimental testing. Here, we introduce a designed monomeric miniprotein (FSD-1ss) that within biological temperatures (330–280 K) switches between simple fast folding and highly complex conformational dynamics in a structurally degenerate compact ensemble. Our findings demonstrate the physical basis for protein folding evolution in a designed protein, which exhibits poorly evolved or primordial folding. Furthermore, these results open the door to the experimental exploration of primitive folding and the switching between alternative protein structures that takes place in evolutionary branching points and prion diseases, as well as the benchmarking of de novo design methods.     

Control over overall shape and size in de novo designed proteins

We recently described general principles for designing ideal protein structures stabilized by completely consistent local and nonlocal interactions. The principles relate secondary structure patterns to tertiary packing motifs and enable design of different protein topologies. To achieve fine control over protein shape and size within a particular topology, we have extended the design rules by systematically analyzing the codependencies between the lengths and packing geometry of successive secondary structure elements and the backbone torsion angles of the loop linking them. We demonstrate the control afforded by the resulting extended rule set by designing a series of proteins with the same fold but considerable variation in secondary structure length, loop geometry, β-strand registry, and overall shape. Solution NMR structures of four designed proteins for two different folds show that protein shape and size can be precisely controlled within a given protein fold. These extended design principles provide the foundation for custom design of protein structures performing desired functions.Protein design holds promise for applications ranging from therapeutics to biomaterials, with recent progress in designing small molecule binding proteins (1, 2), inhibitors of protein–protein interactions (3, 4), and self-assembling nanomaterials (5–7). Most of these efforts have repurposed naturally occurring scaffolds, which are likely not optimal starting points for creating new functions because they generally contain sequence and structural idiosyncrasies that arose during evolutionary optimization for their natural functions (8). Robust design of new functional proteins would be considerably enabled by the capability of precisely designing from scratch arbitrary protein structures.We previously described general principles that allowed the de novo design of ideal protein structures with five different folds (9). In this paper, we focus on the “variations on a theme” problem of precisely controlling structural variation within the same fold. To achieve such control, we begin by characterizing the coupling between loop backbone geometry and the packing of the flanking secondary elements. We then use the resulting extended set of design principles to systematically vary structure for two different folds and describe the experimental characterization of five of these de novo designed proteins.     

Structure of the 26S proteasome from Schizosaccharomyces pombe at subnanometer resolution

The structure of the 26S proteasome from Schizosaccharomyces pombe has been determined to a resolution of 9.1 Å by cryoelectron microscopy and single particle analysis. In addition, chemical cross-linking in conjunction with mass spectrometry has been used to identify numerous residue pairs in close proximity to each other, providing an array of spatial restraints. Taken together these data clarify the topology of the AAA-ATPase module in the 19S regulatory particle and its spatial relationship to the α-ring of the 20S core particle. Image classification and variance analysis reveal a belt of high “activity” surrounding the AAA-ATPase module which is tentatively assigned to the reversible association of proteasome interacting proteins and the conformational heterogeneity among the particles. An integrated model is presented which sheds light on the early steps of protein degradation by the 26S complex.     

Accessing protein conformational ensembles using room-temperature X-ray crystallography

Modern protein crystal structures are based nearly exclusively on X-ray data collected at cryogenic temperatures (generally 100 K). The cooling process is thought to introduce little bias in the functional interpretation of structural results, because cryogenic temperatures minimally perturb the overall protein backbone fold. In contrast, here we show that flash cooling biases previously hidden structural ensembles in protein crystals. By analyzing available data for 30 different proteins using new computational tools for electron-density sampling, model refinement, and molecular packing analysis, we found that crystal cryocooling remodels the conformational distributions of more than 35% of side chains and eliminates packing defects necessary for functional motions. In the signaling switch protein, H-Ras, an allosteric network consistent with fluctuations detected in solution by NMR was uncovered in the room-temperature, but not the cryogenic, electron-density maps. These results expose a bias in structural databases toward smaller, overpacked, and unrealistically unique models. Monitoring room-temperature conformational ensembles by X-ray crystallography can reveal motions crucial for catalysis, ligand binding, and allosteric regulation.