Dramatic structural and thermodynamic consequences of repacking a protein's hydrophobic core

Hydrophobic core evolution of major histocompatibility complex class I chain-related protein A for dramatic enhancing binding affinity

Interface residues at sites of protein-protein interaction (PPI) are the focus for affinity optimisation. However, protein hydrophobic cores (HCs) play critical roles and shape the protein surface. We hypothesise that manipulating protein HCs can enhance PPI interaction affinities. A cell stress molecule, major histocompatibility complex class I chain-related protein A (MICA), binds to the natural killer group 2D (NKG2D) homodimer to form three molecule interactions. MICA was used as a study subject to support our hypothesis. We redesigned MICA HCs by directed mutagenesis and isolated high-affinity variants through a newly designed partial-denature panning (PDP) method. A few mutations in MICA HCs increased the NKG2D-MICA interaction affinity by 325-5613-fold. Crystal structures of the NKG2D-MICA variant complexes indicated that mutagenesis of MICA HCs stabilised helical elements for decreasing intermolecular interactive free energy (ΔG) of the NKG2D-MICA heterotrimer. The repacking of MICA HC mutants maintained overall surface residues and the authentic binding specificity of MICA. In conclusion, this study provides a new method for MICA redesign and affinity optimisation through HC manipulation without mutating PPI interface residues. Our study introduces a novel approach to protein manipulation, potentially expanding the toolkit for protein affinity optimisation.

Deploying synthetic coevolution and machine learning to engineer protein-protein interactions

Fine-tuning of protein-protein interactions occurs naturally through coevolution, but this process is difficult to recapitulate in the laboratory. We describe a platform for synthetic protein-protein coevolution that can isolate matched pairs of interacting muteins from complex libraries. This large dataset of coevolved complexes drove a systems-level analysis of molecular recognition between Z domain-affibody pairs spanning a wide range of structures, affinities, cross-reactivities, and orthogonalities, and captured a broad spectrum of coevolutionary networks. Furthermore, we harnessed pretrained protein language models to expand, in silico, the amino acid diversity of our coevolution screen, predicting remodeled interfaces beyond the reach of the experimental library. The integration of these approaches provides a means of simulating protein coevolution and generating protein complexes with diverse molecular recognition properties for biotechnology and synthetic biology.

Design and Selection of Heterodimerizing Helical Hairpins for Synthetic Biology

Synthetic biology applications would benefit from protein modules of reduced complexity that function orthogonally to cellular components. As many subcellular processes depend on peptide-protein or protein-protein interactions, de novo designed polypeptides that can bring together other proteins controllably are particularly useful. Thanks to established sequence-to-structure relationships, helical bundles provide good starting points for such designs. Typically, however, such designs are tested in vitro and function in cells is not guaranteed. Here, we describe the design, characterization, and application of de novo helical hairpins that heterodimerize to form 4-helix bundles in cells. Starting from a rationally designed homodimer, we construct a library of helical hairpins and identify complementary pairs using bimolecular fluorescence complementation in E. coli. We characterize some of the pairs using biophysics and X-ray crystallography to confirm heterodimeric 4-helix bundles. Finally, we demonstrate the function of an exemplar pair in regulating transcription in both E. coli and mammalian cells.

A dimeric proteomimetic prevents SARS-CoV-2 infection by dimerizing the spike protein

Protein tertiary structure mimetics are valuable tools to target large protein-protein interaction interfaces. Here, we demonstrate a strategy for designing dimeric helix-hairpin motifs from a previously reported three-helix-bundle miniprotein that targets the receptor-binding domain (RBD) of severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2). Through truncation of the third helix and optimization of the interhelical loop residues of the miniprotein, we developed a thermostable dimeric helix-hairpin. The dimeric four-helix bundle competes with the human angiotensin-converting enzyme 2 (ACE2) in binding to RBD with 2:2 stoichiometry. Cryogenic-electron microscopy revealed the formation of dimeric spike ectodomain trimer by the four-helix bundle, where all the three RBDs from either spike protein are attached head-to-head in an open conformation, revealing a novel mechanism for virus neutralization. The proteomimetic protects hamsters from high dose viral challenge with replicative SARS-CoV-2 viruses, demonstrating the promise of this class of peptides that inhibit protein-protein interaction through target dimerization.

Rational design of photosynthetic reaction center protein maquettes

New technologies for efficient solar-to-fuel energy conversion will help facilitate a global shift from dependence on fossil fuels to renewable energy. Nature uses photosynthetic reaction centers to convert photon energy into a cascade of electron-transfer reactions that eventually produce chemical fuel. The design of new reaction centers de novo deepens our understanding of photosynthetic charge separation and may one day allow production of biofuels with higher thermodynamic efficiency than natural photosystems. Recently, we described the multi-step electron-transfer activity of a designed reaction center maquette protein (the RC maquette), which can assemble metal ions, tyrosine, a Zn tetrapyrrole, and heme into an electron-transport chain. Here, we detail our modular strategy for rational protein design and show that the intended RC maquette design agrees with crystal structures in various states of assembly. A flexible, dynamic apo-state collapses by design into a more ordered holo-state upon cofactor binding. Crystal structures illustrate the structural transitions upon binding of different cofactors. Spectroscopic assays demonstrate that the RC maquette binds various electron donors, pigments, and electron acceptors with high affinity. We close with a critique of the present RC maquette design and use electron-tunneling theory to envision a path toward a designed RC with a substantially higher thermodynamic efficiency than natural photosystems.

Using physical features of protein core packing to distinguish real proteins from decoys

The ability to consistently distinguish real protein structures from computationally generated model decoys is not yet a solved problem. One route to distinguish real protein structures from decoys is to delineate the important physical features that specify a real protein. For example, it has long been appreciated that the hydrophobic cores of proteins contribute significantly to their stability. We used two sources to obtain datasets of decoys to compare with real protein structures: submissions to the biennial Critical Assessment of protein Structure Prediction competition, in which researchers attempt to predict the structure of a protein only knowing its amino acid sequence, and also decoys generated by 3DRobot, which have user-specified global root-mean-squared deviations from experimentally determined structures. Our analysis revealed that both sets of decoys possess cores that do not recapitulate the key features that define real protein cores. In particular, the model structures appear more densely packed (because of energetically unfavorable atomic overlaps), contain too few residues in the core, and have improper distributions of hydrophobic residues throughout the structure. Based on these observations, we developed a feed-forward neural network, which incorporates key physical features of protein cores, to predict how well a computational model recapitulates the real protein structure without knowledge of the structure of the target sequence. By identifying the important features of protein structure, our method is able to rank decoy structures with similar accuracy to that obtained by state-of-the-art methods that incorporate many additional features. The small number of physical features makes our model interpretable, emphasizing the importance of protein packing and hydrophobicity in protein structure prediction.

Design and engineering of artificial oxygen-activating metalloenzymes

Many efforts are being made in the design and engineering of metalloenzymes with catalytic properties fulfilling the needs of practical applications. Progress in this field has recently been accelerated by advances in computational, molecular and structural biology. This review article focuses on the recent examples of oxygen-activating metalloenzymes, developed through the strategies of de novo design, miniaturization processes and protein redesign. Considerable progress in these diverse design approaches has produced many metal-containing biocatalysts able to adopt the functions of native enzymes or even novel functions beyond those found in Nature.

Structural plasticity of 4-α-helical bundles exemplified by the puzzle-like molecular assembly of the Rop protein

The dimeric Repressor of Primer (Rop) protein, a widely used model system for the study of coiled-coil 4-α-helical bundles, is characterized by a remarkable structural plasticity. Loop region mutations lead to a wide range of topologies, folding states, and altered physicochemical properties. A protein-folding study of Rop and several loop variants has identified specific residues and sequences that are linked to the observed structural plasticity. Apart from the native state, native-like and molten-globule states have been identified; these states are sensitive to reducing agents due to the formation of nonnative disulfide bridges. Pro residues in the loop are critical for the establishment of new topologies and molten globule states; their effects, however, can be in part compensated by Gly residues. The extreme plasticity in the assembly of 4-α-helical bundles reflects the capacity of the Rop sequence to combine a specific set of hydrophobic residues into strikingly different hydrophobic cores. These cores include highly hydrated ones that are consistent with the formation of interchain, nonnative disulfide bridges and the establishment of molten globules. Potential applications of this structural plasticity are among others in the engineering of bio-inspired materials.

New insights into the interdependence between amino acid stereochemistry and protein structure

To successfully design new proteins and understand the effects of mutations in natural proteins, we must understand the geometric and physicochemical principles underlying protein structure. The side chains of amino acids in peptides and proteins adopt specific dihedral angle combinations; however, we still do not have a fundamental quantitative understanding of why some side-chain dihedral angle combinations are highly populated and others are not. Here we employ a hard-sphere plus stereochemical constraint model of dipeptide mimetics to enumerate the side-chain dihedral angles of leucine (Leu) and isoleucine (Ile), and identify those conformations that are sterically allowed versus those that are not as a function of the backbone dihedral angles ϕ and ψ. We compare our results with the observed distributions of side-chain dihedral angles in proteins of known structure. With the hard-sphere plus stereochemical constraint model, we obtain agreement between the model predictions and the observed side-chain dihedral angle distributions for Leu and Ile. These results quantify the extent to which local, geometrical constraints determine protein side-chain conformations.

Influence of C-H...O interactions on the structural stability of β-lactamases

β-Lactamases produced by pathogenic bacteria cleave β-lactam antibiotics and render them ineffective. Understanding the principles that govern the structural stability of β-lactamases requires elucidation of the nature of the interactions that are involved in stabilization. In the present study, we systematically analyze the influence of CH...O interactions on determining the specificity and stability of β-lactamases in relation to environmental preferences. It is interesting to note that all the residues located in the active site of β-lactamases are involved in CH...O interactions. A significant percentage of CH...O interactions have a higher conservation score and short-range interactions are the predominant type of interactions in β-lactamases. These results will be useful in understanding the stability patterns of β-lactamases.

Energetics of oligomeric protein folding and association

In nature, proteins most often exist as complexes, with many of these consisting of identical subunits. Understanding of the energetics governing the folding and misfolding of such homooligomeric proteins is central to understanding their function and misfunction, in disease or biotechnology. Much progress has been made in defining the mechanisms and thermodynamics of homooligomeric protein folding. In this review, we outline models as well as calorimetric and spectroscopic methods for characterizing oligomer folding, and describe extensive results obtained for diverse proteins, ranging from dimers to octamers and higher order aggregates. To our knowledge, this area has not been reviewed comprehensively in years, and the collective progress is impressive. The results provide evolutionary insights into the development of subunit interfaces, mechanisms of oligomer folding, and contributions of oligomerization to protein stability, function and regulation. Thermodynamic analyses have also proven valuable for understanding protein misfolding and aggregation mechanisms, suggesting new therapeutic avenues. Successful recent designs of novel, functional proteins demonstrate increased understanding of oligomer folding. Further rigorous analyses using multiple experimental and computational approaches are still required, however, to achieve consistent and accurate prediction of oligomer folding energetics. Modeling the energetics remains challenging but is a promising avenue for future advances.

Increasing sequence diversity with flexible backbone protein design: the complete redesign of a protein hydrophobic core

Protein design tests our understanding of protein stability and structure. Successful design methods should allow the exploration of sequence space not found in nature. However, when redesigning naturally occurring protein structures, most fixed backbone design algorithms return amino acid sequences that share strong sequence identity with wild-type sequences, especially in the protein core. This behavior places a restriction on functional space that can be explored and is not consistent with observations from nature, where sequences of low identity have similar structures. Here, we allow backbone flexibility during design to mutate every position in the core (38 residues) of a four-helix bundle protein. Only small perturbations to the backbone, 1-2 Å, were needed to entirely mutate the core. The redesigned protein, DRNN, is exceptionally stable (melting point >140°C). An NMR and X-ray crystal structure show that the side chains and backbone were accurately modeled (all-atom RMSD = 1.3 Å).

Plant actin-binding protein SCAB1 is dimeric actin cross-linker with atypical pleckstrin homology domain

SCAB1 is a novel plant-specific actin-binding protein that binds, bundles, and stabilizes actin filaments and regulates stomatal movement. Here, we dissected the structure and function of SCAB1 by structural and biochemical approaches. We show that SCAB1 is composed of an actin-binding domain, two coiled-coil (CC) domains, and a fused immunoglobulin and pleckstrin homology (Ig-PH) domain. We determined crystal structures for the CC1 and Ig-PH domains at 1.9 and 1.7 Å resolution, respectively. The CC1 domain adopts an antiparallel helical hairpin that further dimerizes into a four-helix bundle. The CC2 domain also mediates dimerization. At least one of the coiled coils is required for actin binding, indicating that SCAB1 is a bivalent actin cross-linker. The key residues required for actin binding were identified. The PH domain lacks a canonical basic phosphoinositide-binding pocket but can bind weakly to inositol phosphates via a basic surface patch, implying the involvement of inositol signaling in SCAB1 regulation. Our results provide novel insights into the functional organization of SCAB1.

The Nudix hydrolase CDP-chase, a CDP-choline pyrophosphatase, is an asymmetric dimer with two distinct enzymatic activities

A Nudix enzyme from Bacillus cereus (NCBI RefSeq accession no. NP_831800) catalyzes the hydrolysis of CDP-choline to produce CMP and phosphocholine. Here, we show that in addition, the enzyme has a 3'→5' RNA exonuclease activity. The structure of the free enzyme, determined to a 1.8-Å resolution, shows that the enzyme is an asymmetric dimer. Each monomer consists of two domains, an N-terminal helical domain and a C-terminal Nudix domain. The N-terminal domain is placed relative to the C-terminal domain such as to result in an overall asymmetric arrangement with two distinct catalytic sites: one with an "enclosed" Nudix pyrophosphatase site and the other with a more open, less-defined cavity. Residues that may be important for determining the asymmetry are conserved among a group of uncharacterized Nudix enzymes from Gram-positive bacteria. Our data support a model where CDP-choline hydrolysis is catalyzed by the enclosed Nudix site and RNA exonuclease activity is catalyzed by the open site. CDP-Chase is the first identified member of a novel Nudix family in which structural asymmetry has a profound effect on the recognition of substrates.

The crystal structure of the calcium-bound con-G[Q6A] peptide reveals a novel metal-dependent helical trimer

The ability to form and control both secondary structure and oligomerization in short peptides has proven to be challenging owing to the structural instability of such peptides. The conantokin peptides are a family of γ-carboxyglutamic acid containing peptides produced in the venoms of predatory sea snails of the Conus family. They are examples of short peptides that form stable helical structures, especially in the presence of divalent cations. Both monomeric and dimeric conantokin peptides have been identified and represent a new mechanism of helix association, "the metallozipper motif" that is devoid of a hydrophobic interface between monomers. In the present study, a parallel/antiparallel three-helix bundle was identified and its crystal structure determined at high resolution. The three helices are almost perfectly parallel and represent a novel helix-helix association. The trimer interface is dominated by metal chelation between the three helices, and contains no interfacial hydrophobic interactions. It is now possible to produce stable monomeric, dimeric, or trimeric metallozippers depending on the peptide sequence and metal ion. Such structures have important applications in protein design.

An unusual hydrophobic core confers extreme flexibility to HEAT repeat proteins

Alpha-solenoid proteins are suggested to constitute highly flexible macromolecules, whose structural variability and large surface area is instrumental in many important protein-protein binding processes. By equilibrium and nonequilibrium molecular dynamics simulations, we show that importin-beta, an archetypical alpha-solenoid, displays unprecedentedly large and fully reversible elasticity. Our stretching molecular dynamics simulations reveal full elasticity over up to twofold end-to-end extensions compared to its bound state. Despite the absence of any long-range intramolecular contacts, the protein can return to its equilibrium structure to within 3 A backbone RMSD after the release of mechanical stress. We find that this extreme degree of flexibility is based on an unusually flexible hydrophobic core that differs substantially from that of structurally similar but more rigid globular proteins. In that respect, the core of importin-beta resembles molten globules. The elastic behavior is dominated by nonpolar interactions between HEAT repeats, combined with conformational entropic effects. Our results suggest that alpha-solenoid structures such as importin-beta may bridge the molecular gap between completely structured and intrinsically disordered proteins.

Cysteine-free Rop: a four-helix bundle core mutant has wild-type stability and structure but dramatically different unfolding kinetics

Cysteine residues can complicate the folding and storage of proteins due to improper formation of disulfide bonds or oxidation of residues that are natively reduced. Wild-type Rop is a homodimeric four-helix bundle protein and an important model for protein design in the understanding of protein stability, structure and folding kinetics. In the native state, Rop has two buried, reduced cysteine residues in its core, but these are prone to oxidation in destabilized variants, particularly upon extended storage. To circumvent this problem, we designed and characterized a Cys-free variant of Rop, including solving the 2.3 A X-ray crystal structure. We show that the C38A C52V variant has similar structure, stability and in vivo activity to wild-type Rop, but that it has dramatically faster unfolding kinetics like virtually every other mutant of Rop that has been characterized. This cysteine-free Rop has already proven useful for studies on solution topology and on the relationship of core mutations to stability. It also suggests a general strategy for removal of pairs of Cys residues in proteins, both to make them more experimentally tractable and to improve their storage properties for therapeutic or industrial purposes.

Structure of the GLD-1 homodimerization domain: insights into STAR protein-mediated translational regulation

Posttranscriptional regulation of gene expression is an important mechanism for modulating protein levels in eukaryotes, especially in developmental pathways. The highly conserved homodimeric STAR/GSG proteins play a key role in regulating translation by binding bipartite consensus sequences in the untranslated regions of target mRNAs, but the exact mechanism remains unknown. Structures of STAR protein RNA binding subdomains have been determined, but structural information is lacking for the homodimerization subdomain. Here, we present the structure of the C. elegans GLD-1 homodimerization domain dimer, determined by a combination of X-ray crystallography and NMR spectroscopy, revealing a helix-turn-helix monomeric fold with the two protomers stacked perpendicularly. Structure-based mutagenesis demonstrates that the dimer interface is not easily disrupted, but the structural integrity of the monomer is crucial for GLD-1 dimerization. Finally, an improved model for STAR-mediated translational regulation of mRNA, based on the GLD-1 homodimerization domain structure, is presented.

Approaches to the assignment of (19)F resonances from 3-fluorophenylalanine labeled calmodulin using solution state NMR

Traditional single site replacement mutations (in this case, phenylalanine to tyrosine) were compared with methods which exclusively employ (15)N and (19)F-edited two- and three-dimensional NMR experiments for purposes of assigning (19)F NMR resonances from calmodulin (CaM), biosynthetically labeled with 3-fluorophenylalanine (3-FPhe). The global substitution of 3-FPhe for native phenylalanine was tolerated in CaM as evidenced by a comparison of (1)H-(15)N HSQC spectra and calcium binding assays in the presence and absence of 3-FPhe. The (19)F NMR spectrum reveals six resolved resonances, one of which integrates to three 3-FPhe species, making for a total of eight fluorophenylalanines. Single phenylalanine to tyrosine mutants of five phenylalanine positions resulted in (19)F NMR spectra with significant chemical shift perturbations of the remaining resonances, and provided only a single definitive assignment. Although (1)H-(19)F heteronucleclear NOEs proved weak, (19)F-edited (1)H-(1)H NOESY connectivities were relatively easy to establish by making use of the (3)J(FH) coupling between the fluorine nucleus and the adjacent fluorophenylalanine delta proton. (19)F-edited NOESY connectivities between the delta protons and alpha and beta nuclei in addition to (15)N-edited (1)H, (1)H NOESY crosspeaks proved sufficient to assign 4 of 8 (19)F resonances. Controlled cleavage of the protein into two fragments using trypsin, and a repetition of the above 2D and 3D techniques resulted in unambiguous assignments of all 8 (19)F NMR resonances. Our studies suggest that (19)F-edited NOESY NMR spectra are generally adequate for complete assignment without the need to resort to mutational analysis.

Rapid model building of alpha-helices in electron-density maps

A method for the identification of alpha-helices in electron-density maps at low resolution followed by interpretation at moderate to high resolution is presented. Rapid identification is achieved at low resolution, where alpha-helices appear as tubes of density. The positioning and direction of the alpha-helices is obtained at moderate to high resolution, where the positions of side chains can be seen. The method was tested on a set of 42 experimental electron-density maps at resolutions ranging from 1.5 to 3.8 A. An average of 63% of the alpha-helical residues in these proteins were built and an average of 76% of the residues built matched helical residues in the refined models of the proteins. The overall average r.m.s.d. between main-chain atoms in the modeled alpha-helices and the nearest atom with the same name in the refined models of the proteins was 1.3 A.

Rapid model building of beta-sheets in electron-density maps

A method for rapidly building beta-sheets into electron-density maps is presented. beta-Strands are identified as tubes of high density adjacent to and nearly parallel to other tubes of density. The alignment and direction of each strand are identified from the pattern of high density corresponding to carbonyl and C(beta) atoms along the strand averaged over all repeats present in the strand. The beta-strands obtained are then assembled into a single atomic model of the beta-sheet regions. The method was tested on a set of 42 experimental electron-density maps at resolutions ranging from 1.5 to 3.8 A. The beta-sheet regions were nearly completely built in all but two cases, the exceptions being one structure at 2.5 A resolution in which a third of the residues in beta-sheets were built and a structure at 3.8 A in which under 10% were built. The overall average r.m.s.d. of main-chain atoms in the residues built using this method compared with refined models of the structures was 1.5 A.

Rapid chain tracing of polypeptide backbones in electron-density maps

A method for rapid chain tracing of polypeptide backbones at moderate resolution is presented.

A method for the rapid tracing of polypeptide backbones has been developed. The method creates an approximate chain tracing that is useful for visual evaluation of whether a structure has been solved and for use in scoring the quality of electron-density maps. The essence of the method is to (i) sample candidate C^α positions at spacings of approximately 0.6 Å along ridgelines of high electron density, (ii) list all possible nonapeptides that satisfy simple geometric and density criteria using these candidate C^α positions, (iii) score the nonapeptides and choose the highest scoring ones, and (iv) find the longest chains that can be made by connecting nonamers. An indexing and storage scheme that allows a single calculation of most distances and density values is used to speed up the process. The method was applied to 42 density-modified electron-density maps at resolutions from 1.5 to 3.8 Å. A total of 21 428 residues in these maps were traced in 24 CPU min with an overall r.m.s.d. of 1.61 Å for C^α atoms compared with the known refined structures. The method appears to be suitable for rapid evaluation of electron-density map quality.

Robustness and generalization of structure-based models for protein folding and function

Functional dynamics of native proteins share the energy landscape that guides folding into the native state. Folding simulations of structure-based protein models, using an minimally frustrated energy landscape dominated by native interactions, can describe the geometrical aspects of the folding mechanism. Technical limitations imposed by the fixed shape of conventional contact potentials are a key obstacle toward advanced applications of structure-based models like allostery or ligand binding, which require multiple stable conformations. Generalizations of existing models, commonly using Lennard-Jones-like potentials, lead to inevitable clashes between their repulsive branches. To resolve these challenges, a new contact potential is developed that combines an attractive part based on Gaussians with a separate repulsive term allowing flexibility for adjustments of the potential shape. With this new model multiple minima for studies of functional transitions can be introduced easily and consistently. A sensitivity analysis for five small proteins confirms the robust behavior of structure-based models with our adaptable potential and explores their capacity for quantitative adjustment of the folding thermodynamics. We demonstrate its ability to incorporate alternative contact distances in simulations of structural transitions for the well-studied ROP dimer. Individual contact pairs can switch between distinct states to match the competing syn and anti structures. The flexibility of the new potential facilitates advanced uses of structure-based models. Depending on the application, features can be chosen from physical considerations or to match experiments. Generalized models can be built from multiple structures to study structural transitions or effects of disorder.

Direct single-molecule observation of a protein living in two opposed native structures

Biological activity in proteins requires them to share the energy landscape for folding and global conformational motions, 2 key determinants of function. Although most structural studies to date have focused on fluctuations around a single structural basin, we directly observe the coexistence of 2 symmetrically opposed conformations for a mutant of the Rop-homodimer (Repressor of Primer) in single-molecule fluorescence resonance energy transfer (smFRET) measurements. We find that mild denaturing conditions can affect the sensitive balance between the conformations, generating an equilibrium ensemble consisting of 2 equally occupied structural basins. Despite the need for large-scale conformational rearrangement, both native structures are dynamically and reversibly adopted for the same paired molecules without separation of the constituent monomers. Such an ability of some proteins or protein complexes to switch between conformations by thermal fluctuations and/or minor environmental changes could be central to their ability to control biological function.

Decision-making in structure solution using Bayesian estimates of map quality: the PHENIX AutoSol wizard

Ten measures of experimental electron-density-map quality are examined and the skewness of electron density is found to be the best indicator of actual map quality. A Bayesian approach to estimating map quality is developed and used in the PHENIX AutoSol wizard to make decisions during automated structure solution.

Estimates of the quality of experimental maps are important in many stages of structure determination of macromolecules. Map quality is defined here as the correlation between a map and the corresponding map obtained using phases from the final refined model. Here, ten different measures of experimental map quality were examined using a set of 1359 maps calculated by re-analysis of 246 solved MAD, SAD and MIR data sets. A simple Bayesian approach to estimation of map quality from one or more measures is presented. It was found that a Bayesian estimator based on the skewness of the density values in an electron-density map is the most accurate of the ten individual Bayesian estimators of map quality examined, with a correlation between estimated and actual map quality of 0.90. A combination of the skewness of electron density with the local correlation of r.m.s. density gives a further improvement in estimating map quality, with an overall correlation coefficient of 0.92. The PHENIX AutoSol wizard carries out automated structure solution based on any combination of SAD, MAD, SIR or MIR data sets. The wizard is based on tools from the PHENIX package and uses the Bayesian estimates of map quality described here to choose the highest quality solutions after experimental phasing.

Organic ligand binding by a hydrophobic cavity in a designed tetrameric coiled-coil protein

The design and characterization of a hydrophobic cavity in de novo designed proteins provides a wide range of information about the functions of de novo proteins. We designed a de novo tetrameric coiled-coil protein with a hydrophobic pocketlike cavity. Tetrameric coiled coils with hydrophobic cavities have previously been reported. By replacing one Leu residue at the a position with Ala, hydrophobic cavities that did not flatten out due to loose peptide chains were reliably created. To perform a detailed examination of the ligand-binding characteristics of the cavities, we originally designed two other coiled-coil proteins: AM2, with eight Ala substitutions at the adjacent a and d positions at the center of a bundled structure, and AM2W, with one Trp and seven Ala substitutions at the same positions. To increase the association of the helical peptides, each helical peptide was connected with flexible linkers, which resulted in a single peptide chain. These proteins exhibited CD spectra corresponding to superhelical structures, despite weakened hydrophobic packing. AM2W exhibited binding affinity for size-complementary organic compounds. The dissociation constants, K(d), of AM2W were 220 nM for adamantane, 81 microM for 1-adamantanol, and 294 microM for 1-adamantaneacetic acid, as measured by fluorescence titration analyses. Although it was contrary to expectations, AM2 did not exhibit any binding affinity, probably due to structural defects around the designed hydrophobic cavity. Interestingly, AM2W exhibited incremental structure stability through ligand binding. Plugging of structural defects with organic ligands would be expected to facilitate protein folding.

Structure and dynamics of de novo proteins from a designed superfamily of 4-helix bundles

Libraries of de novo proteins provide an opportunity to explore the structural and functional potential of biological molecules that have not been biased by billions of years of evolutionary selection. Given the enormity of sequence space, a rational approach to library design is likely to yield a higher fraction of folded and functional proteins than a stochastic sampling of random sequences. We previously investigated the potential of library design by binary patterning of hydrophobic and hydrophilic amino acids. The structure of the most stable protein from a binary patterned library of de novo 4-helix bundles was solved previously and shown to be consistent with the design. One structure, however, cannot fully assess the potential of the design strategy, nor can it account for differences in the stabilities of individual proteins. To more fully probe the quality of the library, we now report the NMR structure of a second protein, S-836. Protein S-836 proved to be a 4-helix bundle, consistent with design. The similarity between the two solved structures reinforces previous evidence that binary patterning can encode stable, 4-helix bundles. Despite their global similarities, the two proteins have cores that are packed at different degrees of tightness. The relationship between packing and dynamics was probed using the Modelfree approach, which showed that regions containing a high frequency of chemical exchange coincide with less well-packed side chains. These studies show (1) that binary patterning can drive folding into a particular topology without the explicit design of residue-by-residue packing, and (2) that within a superfamily of binary patterned proteins, the structures and dynamics of individual proteins are modulated by the identity and packing of residues in the hydrophobic core.

Mutations as trapdoors to two competing native conformations of the Rop-dimer

Conformational transitions play a central role in regulating protein function. Structure-based models with multiple basins have been used to understand the mechanisms governing these transitions. A model able to accommodate multiple folding basins is proposed to explore the mutational effects in the folding of the Rop-dimer (Rop). In experiments, Rop mutants show unusually strong increases in folding rates with marginal effects on stability. We investigate the possibility of two competing conformations representing a parallel (P) and the wild-type antiparallel (AP) arrangement of the monomers as possible native conformations. We observe occupation of both distinct states and characterize the transition pathways. An interesting observation from the simulations is that, for equivalent energetic bias, the transition to the P basin (non-wild-type basin) shows a lower free-energy barrier. Thus, the rapid kinetics observed in experiments appear to be the result of two competing states with different kinetic behavior, triggered upon mutation by the opening of a trapdoor arising from Rop's symmetric structure. The general concept of having competing conformations for the native state goes beyond explaining Rop's mutational behaviors and can be applied to other systems. A switch between competing native structures might be triggered by external factors to allow, for example, allosteric control or signaling.

Nanoscale dewetting transition in protein complex folding

In a previous study, a surprising drying transition was observed to take place inside the nanoscale hydrophobic channel in the tetramer of the protein melittin. The goal of this paper is to determine if there are other protein complexes capable of displaying a dewetting transition during their final stage of folding. We searched the entire protein data bank (PDB) for all possible candidates, including protein tetramers, dimers, and two-domain proteins, and then performed the molecular dynamics (MD) simulations on the top candidates identified by a simple hydrophobic scoring function based on aligned hydrophobic surface areas. Our large scale MD simulations found several more proteins, including three tetramers, six dimers, and two two-domain proteins, which display a nanoscale dewetting transition in their final stage of folding. Even though the scoring function alone is not sufficient (i.e., a high score is necessary but not sufficient) in identifying the dewetting candidates, it does provide useful insights into the features of complex interfaces needed for dewetting. All top candidates have two features in common: (1) large aligned (matched) hydrophobic areas between two corresponding surfaces, and (2) large connected hydrophobic areas on the same surface. We have also studied the effect on dewetting of different water models and different treatments of the long-range electrostatic interactions (cutoff vs PME), and found the dewetting phenomena is fairly robust. This work presents a few proteins other than melittin tetramer for further experimental studies of the role of dewetting in the end stages of protein folding.

Structural plasticity of an aminoacyl-tRNA synthetase active site

Recently, tRNA aminoacyl-tRNA synthetase pairs have been evolved that allow one to genetically encode a large array of unnatural amino acids in both prokaryotic and eukaryotic organisms. We have determined the crystal structures of two substrate-bound Methanococcus jannaschii tyrosyl aminoacyl-tRNA synthetases that charge the unnatural amino acids p-bromophenylalanine and 3-(2-naphthyl)alanine (NpAla). A comparison of these structures with the substrate-bound WT synthetase, as well as a mutant synthetase that charges p-acetylphenylalanine, shows that altered specificity is due to both side-chain and backbone rearrangements within the active site that modify hydrogen bonds and packing interactions with substrate, as well as disrupt the alpha8-helix, which spans the WT active site. The high degree of structural plasticity that is observed in these aminoacyl-tRNA synthetases is rarely found in other mutant enzymes with altered specificities and provides an explanation for the surprising adaptability of the genetic code to novel amino acids.

Exploring amino-acid radical chemistry: protein engineering and de novo design

Amino-acid radical enzymes are often highly complex structures containing multiple protein subunits and cofactors. These properties have in many cases hampered the detailed characterization of their amino-acid redox cofactors. To address this problem, a range of approaches has recently been developed in which a common strategy is to reduce the complexity of the radical-containing system. This work will be reviewed and it includes the light-induced generation of aromatic radicals in small-molecule and peptide systems. Natural redox proteins, including the blue copper protein azurin and a bacterial photosynthetic reaction center, have been engineered to introduce amino-acid radical chemistry. The redesign strategies to achieve this remarkable change in the properties of these proteins will be described. An additional approach to gain insights into the properties of amino-acid radicals is to synthesize de novo designed model proteins in which the redox chemistry of these species can be studied. Here we describe the design, synthesis and characteristics of monomeric three-helix bundle and four-helix bundle proteins designed to study the redox chemistry of tryptophan and tyrosine. This work demonstrates that de novo protein design combined with structural, electrochemical and quantum chemical analyses can provide detailed information on how the protein matrix tunes the thermodynamic properties of tryptophan.

Symmetry and frustration in protein energy landscapes: a near degeneracy resolves the Rop dimer-folding mystery

Protein folding has become one of the best understood biochemical reactions from a kinetic viewpoint. The funneled energy landscape, a consequence of the minimal frustration achieved by evolution in sequences, explains how most proteins fold efficiently and robustly to their functional structure and allows robust prediction of folding kinetics. The folding of Rop (repressor of primer) dimer is exceptional because some of its mutants with a redesigned hydrophobic core both fold and unfold much faster than the WT protein, which seems to conflict with a simple funneled energy landscape for which topology mainly determines the kinetics. We propose that the mystery of Rop folding can be unraveled by assuming a double-funneled energy landscape on which there are two basins that correspond to distinct but related topological structures. Because of the near symmetry of the molecule, mutations can cause a conformational switch to a nearly degenerate yet distinct topology or lead to a mixture of both topologies. The topology predicted to have the lower free-energy barrier height for folding was further found by all-atom modeling to give a better structural fit for those mutants with the extreme folding and unfolding rates. Thus, the non-Hammond effects can be understood within energy-landscape theory if there are in fact two different but nearly degenerate structures for Rop. Mutations in symmetric and regular structures may give rise to frustration and thus result in degeneracy.

Hepatitis B Virus Capsid Assembly Is Enhanced by Naturally Occurring Mutation F97L

In chronic hepatitis B virus (HBV) infections, one of the most common mutations to the virus occurs at amino acid 97 of the core protein, where leucine replaces either phenylalanine or isoleucine, depending on strain. This mutation correlates with changes in viral nucleic acid metabolism and/or secretion. We hypothesize that this phenotype is due in part to altered core assembly, a process required for DNA synthesis. We examined in vitro assembly of empty HBV capsids from wild-type and F97L core protein assembly domains. The mutation enhanced both the rate and extent of assembly relative to those for the wild-type protein. The difference between the two proteins was most obvious in the temperature dependence of assembly, which was dramatically stronger for the mutant protein, indicating a much more positive enthalpy. Since the structures of the mutant and wild-type capsids are essentially the same and the mutation is not involved in the contact between dimers, we suggest that the F97L mutation affects the dynamic behavior of dimer and capsid.

Combinatorial approaches to protein stability and structure

Why do proteins adopt the conformations that they do, and what determines their stabilities? While we have come to some understanding of the forces that underlie protein architecture, a precise, predictive, physicochemical explanation is still elusive. Two obstacles to addressing these questions are the unfathomable vastness of protein sequence space, and the difficulty in making direct physical measurements on large numbers of protein variants. Here, we review combinatorial methods that have been applied to problems in protein biophysics over the last 15 years. The effects of hydrophobic core composition, the most important determinant of structure and stability, are still poorly understood. Particular attention is given to core composition as addressed by library methods. Increasingly useful screens and selections, in combination with modern high-throughput approaches borrowed from genomics and proteomics efforts, are making the empirical, statistical correlation between sequence and structure a tractable problem for the coming years.

Local complexity of amino acid interactions in a protein core

Atomic resolution structures of proteins indicate that the core is typically well packed, suggesting a densely connected network of interactions between amino acid residues. The combinatorial complexity of energetic interactions in such a network could be enormous, a problem that limits our ability to relate structure and function. Here, we report a case study of the complexity of amino acid interactions in a localized region within the core of the GFP, a particularly stable and tightly packed molecule. Mutations at three sites within the chromophore-binding pocket display an overlapping pattern of conformational change and are thermodynamically coupled, seemingly consistent with the dense network model. However, crystallographic and energetic analyses of coupling between mutations paint a different picture; pairs of mutations couple through independent "hotspots" in the region of structural overlap. The data indicate that, even in highly stable proteins, the core contains sufficient plasticity in packing to uncouple high-order energetic interactions of residues, a property that is likely general in proteins.

X-ray structure of a maquette scaffold

Maquettes are de novo designed mimicries of nature used to test the construction and engineering criteria of oxidoreductases. One type of scaffold used in maquette construction is a four-alpha-helical bundle. The sequence of the four-alpha-helix bundle maquettes follows a heptad repeat pattern typical of left-handed coiled-coils. Initial designs were molten globular due partly to the minimalist approach taken by the designers. Subsequent iterative redesign generated several structured scaffolds with similar heme binding properties. Variant [I(6)F(13)](2), a structured scaffold, was partially resolved with NMR spectroscopy and found to have a set of mobile inter-helical packing interfaces. Here, the X-ray structure of a similar peptide ([I(6)F(13)M(31)](2) i.e. ([CGGG EIWKL HEEFLKK FEELLKL HEERLKKM](2))(2) which we call L31M), has been solved using MAD phasing and refined to 2.8A resolution. The structure shows that the maquette scaffold is an anti-parallel four-helix bundle with "up-up-down-down" topology. No pre-formed heme-binding pocket exists in the protein scaffold. We report unexpected inter-helical crossing angles, residue positions and translations between the helices. The crossing angles between the parallel helices are -5 degrees rather than the expected +20 degrees for typical left-handed coiled-coils. Deviation of the scaffold from the design is likely due to the distribution and size of hydrophobic residues. The structure of L31M points out that four identical helices may interact differently in a bundle and heptad repeats with an alternating [HPPHHPP]/[HPPHHPH] (H: hydrophobic, P: polar) pattern are not a sufficient design criterion to generate left-hand coiled-coils.

Rotamer libraries in the 21st century

Rotamer libraries are widely used in protein structure prediction, protein design, and structure refinement. As the size of the structure data base has increased rapidly in recent years, it has become possible to derive well-refined rotamer libraries using strict criteria for data inclusion and for studying dependence of rotamer populations and dihedral angles on local structural features.

Heat capacity changes upon burial of polar and nonpolar groups in proteins

In this paper we address the question of whether the burial of polar and nonpolar groups in the protein locale is indeed accompanied by the heat capacity changes, DeltaC(p), that have an opposite sign, negative for nonpolar groups and positive for polar groups. To accomplish this, we introduced amino acid substitutions at four fully buried positions of the ubiquitin molecule (Val5, Val17, Leu67, and Gln41). We substituted Val at positions 5 and 17 and Leu at position 67 with a polar residue, Asn. As a control, Ala was introduced at the same three positions. We also replaced the buried polar Gln41 with Val and Leu, nonpolar residues that have similar size and shape as Gln. As a control, Asn was introduced at Gln41 as well. The effects of these amino acid substitutions on the stability, and in particular, on the heat capacity change upon unfolding were measured using differential scanning calorimetry. The effect of the amino acid substitutions on the structure was also evaluated by comparing the (1)H-(15)N HSQC spectra of the ubiquitin variants. It was found that the Ala substitutions did not have a considerable effect on the heat capacity change upon unfolding. However, the substitutions of aliphatic side chains (Val or Leu) with a polar residue (Asn) lead to a significant (> 30%) decrease in the heat capacity change upon unfolding. The decrease in heat capacity changes does not appear to be the result of significant structural perturbations as seen from the HSQC spectra of the variants. The substitution of a buried polar residue (Gln41) to a nonpolar residue (Leu or Val) leads to a significant (> 25%) increase in heat capacity change upon unfolding. These results indicate that indeed the heat capacity change of burial of polar and nonpolar groups has an opposite sign. However, the observed changes in DeltaC(p) are several times larger than those predicted, based on the changes in water accessible surface area upon substitution.