NMR assignment of S836: a de novo protein from a designed superfamily

Libraries of de novo proteins provide an opportunity to explore the structural potential of biological macromolecules that have not been biased by billions of years of evolutionary selection. Characterization of individual members of such libraries provides insight into the diversity of structure and dynamics accessible to nascent protein superfamilies in the absence of evolutionary optimization. Here we report the backbone and side chain chemical shifts of protein S836 from a superfamily of designed 4-helix bundles.     

Co-translational domain folding as the structural basis for the rapid de novo folding of firefly luciferase.

The 62 kDa protein firefly luciferase folds very rapidly upon translation on eukaryotic ribosomes. In contrast, the chaperone-mediated refolding of chemically denatured luciferase occurs with significantly slower kinetics. Here we investigate the structural basis for this difference in folding kinetics. We find that an N-terminal domain of luciferase (residues 1-190) folds co-translationally, followed by rapid formation of native protein upon release of the full-length polypeptide from the ribosome. In contrast sequential domain formation is not observed during in vitro refolding. Discrete unfolding steps, corresponding to domain unfolding, are however observed when the native protein is exposed to increasing concentrations of denaturant. Thus, the co-translational folding reaction bears more similarities to the unfolding reaction than to refolding from denaturant. We propose that co-translational domain formation avoids intramolecular misfolding and may be critical in the folding of multidomain proteins.     

Transition-states in protein folding kinetics: the structural interpretation of Phi values

Phi values are experimental measures of the effects of mutations on the folding kinetics of a protein. A central question is what structural information Phi values give about the transition-state of folding. Traditionally, a Phi value is interpreted as representing the "nativeness" of a mutated residue in the transition-state. However, this interpretation is often problematic. We present here a better structural interpretation of Phi values for mutations within a given helix. Our interpretation is based on a simple physical model that distinguishes between secondary and tertiary free energy contributions of helical residues. From a linear fit of the model to experimental data, we obtain two structural parameters: the extent of helix formation in the transition-state, and the nativeness of tertiary interactions in the transition-state. We apply the model to all proteins with well-characterized helices for which more than 10 Phi values are available: protein A, CI2, and protein L. The model is simple to apply to experimental data, captures nonclassical Phi values <0 or >1 in these helices, and explains how different mutations at a given site can lead to different Phi values.     

Folding kinetics of designer proteins. Application of the diffusion-collision model to a de novo designed four-helix bundle.

A folding algorithm is described, based on the diffusion-collision model, combining static and dynamic calculational methods. The algorithm is applied to predict the basic structure and schematic folding pathways of an artificial four-helix bundle.     

Non-linear effects of temperature and urea on the thermodynamics and kinetics of folding and unfolding of hisactophilin

Extensive measurements and analysis of thermodynamic stability and kinetics of urea-induced unfolding and folding of hisactophilin are reported for 5-50 degrees C, at pH 6.7. Under these conditions hisactophilin has moderate thermodynamic stability, and equilibrium and kinetic data are well fit by a two-state transition between the native and the denatured states. Equilibrium and kinetic m values decrease with increasing temperature, and decrease with increasing denaturant concentration. The betaF values at different temperatures and urea concentrations are quite constant, however, at about 0.7. This suggests that the transition state for hisactophilin unfolding is native-like and changes little with changing solution conditions, consistent with a narrow free energy profile for the transition state. The activation enthalpy and entropy of unfolding are unusually low for hisactophilin, as is also the case for the corresponding equilibrium parameters. Conventional Arrhenius and Eyring plots for both folding and unfolding are markedly non-linear, but these plots become linear for constant DeltaG/T contours. The Gibbs free energy changes for structural changes in hisactophilin have a non-linear denaturant dependence that is comparable to non-linearities observed for many other proteins. These non-linearities can be fit for many proteins using a variation of the Tanford model, incorporating empirical quadratic denaturant dependencies for Gibbs free energies of transfer of amino acid constituents from water to urea, and changes in fractional solvent accessible surface area of protein constituents based on the known protein structures. Noteworthy exceptions that are not well fit include amyloidogenic proteins and large proteins, which may form intermediates. The model is easily implemented and should be widely applicable to analysis of urea-induced structural transitions in proteins.     

Microdrop screening: A rapid method to optimize solvent conditions for NMR spectroscopy of proteins

Determining appropriate solvent conditions is a crucial first step for carrying out NMR spectroscopy of proteins, but rapid and efficient methods for doing so are currently lacking. Microdrop screening examines a large number of different solvent conditions using very small amounts of protein and minimal labor. Starting from one initial buffer condition, small aliquots of protein solution are combined with an array of solutions in which concentration, pH, buffer type, and added stabilizers are systematically varied. The protein concentration of each microliter-sized test drop (microdrop) is gradually changed using vapor diffusion, and the solubility of the protein is determined by visual examination. A variety of analytical techniques may be applied to the contents of the microdrops to monitor enzymatic activity, aggregation, ligand binding, and protein folding.     

De novo folding of membrane proteins: an exploration of the structure and NMR properties of the fd coat protein

De novo folding simulations of the major pVIII coat protein from filamentous fd bacteriophage, using a newly developed implicit membrane generalized Born model and replica-exchange molecular dynamics, are presented and discussed. The quality of the predicted structures, judged by comparison of the root-mean-square deviations of a room temperature ensemble of conformations from the replica-exchange simulations and experimental structures from both solid-state NMR in lipid bilayers and solution-phase NMR on the protein in micelles, was quite good, reinforcing the general quality of the folding simulations. The transmembrane helical segment of the protein was well defined in comparison with experiment and the amphipathic helical fragment remained at the membrane/aqueous phase boundary while undergoing significant conformational flexibility due to the loop connecting the two helical segments of the protein. Additional comparisons of computed solid-state NMR properties, the 15N chemical shift and 15N-1H dipolar coupling constants, showed semi-quantitative agreement with the corresponding measurements. These findings suggest an emerging potential for the de novo investigation of integral membrane peptides and proteins and a mechanism to assist experimental approaches to the characterization and structure determination of these important systems.     

Application of the "codon-shuffling" method. Synthesis and selection of de novo proteins as antibacterials

Library-based methods of non-rational and part-rational designed de novo peptides are worthy beacons in the search for bioactive peptides and proteins of medicinal importance. In this report, we have used a recently developed directed evolution method called "codon shuffling" for the synthesis and selection of bioactive proteins. The selection of such proteins was based on the creation of an inducible library of "codon-shuffled" genes that are constructed from the ligation-based assembly of judiciously designed hexamer DNA duplexes called dicodons. Upon induction with isopropyl 1-thio-beta-D-galactopyranoside, some library members were found to express dicodon-incorporated proteins. Because of this, the host cells, in our case Escherichia coli, were unable to grow any further. The bactereostatic/lytic nature of the dicodon proteins was monitored by growth curves as well as by zone clearance studies. Transmission electron microscopy of the affected cells illustrated the extent of cell damage. The proteins themselves were overexpressed as fusion partners and subsequently purified to homogeneity. One such purified protein was found to strongly bind heparin, an indication that the interaction of the de novo proteins may be with the nucleic acids of the host cell, much like many of the naturally occurring antibacterial peptides, e.g. Buforin. Therefore, our approach may help in generating a multitude of finely tuned antibacterial proteins that can potentially be regarded as lead compounds once the method is extended to pathogenic hosts, such as Mycobacteria, for example.     

Comparison of the folding processes of distantly related proteins. Importance of hydrophobic content in folding

The N-terminal domain of HypF from Escherichia coli (HypF-N) is a 91 residue protein module sharing the same folding topology and a significant sequence identity with two extensively studied human proteins, muscle and common-type acylphosphatases (mAcP and ctAcP). With the aim of learning fundamental aspects of protein folding from the close comparison of so similar proteins, the folding process of HypF-N has been studied using stopped-flow fluorescence. While mAcP and ctAcP fold in a two-state fashion, HypF-N was found to collapse into a partially folded intermediate before reaching the fully folded conformation. Formation of a burst-phase intermediate is indicated by the roll over in the Chevron plot at low urea concentrations and by the large jump of intrinsic and 8-anilino-1-naphtalenesulphonic acid-derived fluorescence immediately after removal of denaturant. Furthermore, HypF-N was found to fold rapidly with a rate constant that is approximately two and three orders of magnitudes faster than ctAcP and mAcP, respectively. Differences between the bacterial protein and the two human counterparts were also found as to the involvement of proline isomerism in their respective folding processes. The results clearly indicate that features that are often thought to be relevant in protein folding are not highly conserved in the evolution of the acylphosphatase superfamily. The large difference in folding rate between mAcP and HypF-N cannot be entirely accounted for by the difference in relative contact order or related topological metrics. The analysis shows that the higher folding rate of HypF-N is in part due to the relatively high hydrophobic content of this protein. This conclusion, which is also supported by the highly significant correlation found between folding rate and hydrophobic content within a group of proteins displaying the topology of HypF-N and AcPs, suggests that the average hydrophobicity of a protein sequence is an important determinant of its folding rate.     

Computer simulations of the properties of the alpha2, alpha2C, and alpha2D de novo designed helical proteins

Reduced lattice models of the three de novo designed helical proteins alpha2, alpha2C, and alpha2D were studied. Low temperature stable folds were obtained for all three proteins. In all cases, the lowest energy folds were four-helix bundles. The folding pathway is qualitatively the same for all proteins studied. The energies of various topologies are similar, especially for the alpha2 polypeptide. The simulated crossover from molten globule to native-like behavior is very similar to that seen in experimental studies. Simulations on a reduced protein model reproduce most of the experimental properties of the alpha2, alpha2C, and alpha2D proteins. Stable four-helix bundle structures were obtained, with increasing native-like behavior on-going from alpha2 to alpha2D that mimics experiment.     

Sequence optimization for native state stability determines the evolution and folding kinetics of a small protein

Investigating the relative importance of protein stability, function, and folding kinetics in driving protein evolution has long been hindered by the fact that we can only compare modern natural proteins, the products of the very process we seek to understand, to each other, with no external references or baselines. Through a large-scale all-atom simulation of protein evolution, we have created a large diverse alignment of SH3 domain sequences which have been selected only for native state stability, with no other influencing factors. Although the average pairwise identity between computationally evolved and natural sequences is only 17%, the residue frequency distributions of the computationally evolved sequences are similar to natural SH3 sequences at 86% of the positions in the domain, suggesting that optimization for the native state structure has dominated the evolution of natural SH3 domains. Additionally, the positions which play a consistent role in the transition state of three well-characterized SH3 domains (by phi-value analysis) are structurally optimized for the native state, and vice versa. Indeed, we see a specific and significant correlation between sequence optimization for native state stability and conservation of transition state structure.     

NMR spectroscopic filtration of polypeptides and proteins in complex mixtures

Due to the inherent complexity of the natural biological environment, most studies on polypeptides, proteins and nucleic acids have so far been performed in vitro, away from physiologically relevant conditions. Nuclear magnetic resonance is an ideal technique to extend the in vitro analysis of simple model systems to the more complex biological context. This work shows how diffusion-based spectroscopic selection can be combined with isotopic labeling to tackle and optimize the NMR analysis of specific macromolecules in multicomponent mixtures. Typical media include cell-free systems containing overexpressed proteins, lysates and proteolytic mixtures. We present a few variants of diffusion-edited HSQC pulse sequences for the selective spectroscopic detection of protein and polypeptide resonances within complex mixtures containing undesired species of smaller molecular weight. Due to diffusion-based filtering, peak intensities of fast diffusing small molecules are attenuated more than peaks due to large molecules. The basic sequence, denoted as PFGSTE-HSQC, combines translational diffusion-ordering with two dimensional heteronuclear single quantum correlation spectroscopy. The GCSTE-HSQC and BPPSTE-HSQC sequences include bipolar gradients and are therefore suitable for both diffusion-based filtering and determination of diffusion coefficients of individual mixture components. Practical applications range from protein stability/folding investigations in physiologically relevant contexts to prescreening of tertiary fold and resonance assignments in structural genomics studies. A few applications of diffusion-edited HSQC to an E. coli cell lysate containing the (15)N-labeled B domain of streptococcal protein G (GB1), and to a (15)N-labeled N-acetylglycine/apomyoglobin mixture are presented. In addition, we provide specific guidelines for experimental setup and parameter optimization.     

Femtoliter compartment in liposomes for in vitro selection of proteins

The aqueous compartment in liposomes provides a reaction resembling the cell and therefore is used as a microcompartment in which to study enzymatic reactions. However, regardless of their method of preparation, the heterogeneity in size of cell-size liposomes limits their potential uses. We established a strategy to estimate the internal aqueous volume of cell-size liposomes using a fluorescence-activated cell sorter (FACS). Reactions inside individual liposomes can be measured in a high-throughput format provided that the encapsulated proteins give rise to a fluorescent signal such as by exhibiting fluorescence themselves or by catalyzing production of a fluorescent compound. The strategy of volume estimation was applied to in vitro selection experiments. The green fluorescent protein (GFP) gene was encapsulated into liposomes together with an in vitro translation system. Here liposomes carrying a single copy of the gene were identified using the internal aqueous volume information of individual liposomes, and those exhibiting higher green fluorescence intensity were sorted by the FACS machine. This system was able to enrich those encoding GFP with higher fluorescence intensity over those with lower intensity. These results suggest the possibility of performing evolutionary experiments in an environment that mimics the cell.     

Effects of domain dissection on the folding and stability of the 43 kDa protein PGK probed by NMR

The characterization of early folding intermediates is key to understanding the protein folding process. Previous studies of the N-domain of phosphoglycerate kinase (PGK) from Bacillus stearothermophilus combined equilibrium amide exchange data with a kinetic model derived from stopped-flow kinetics. Together, these implied the rapid formation of an intermediate with extensive native-like hydrogen bonding. However, there was an absence of protection in the region proximal to the C-domain in the intact protein. We now report data for the intact PGK molecule, which at 394 residues constitutes a major extension to the protein size for which such data can be acquired. The methods utilised to achieve the backbone assignment are described in detail, including a semi-automated protocol based on a simulated annealing Monte Carlo technique. A substantial increase in the stability of the contact region is observed, allowing protection to be inferred on both faces of the beta-sheet in the intermediate. Thus, the entire N-domain acts concertedly in the formation of the kinetic refolding intermediate rather than there existing a distinct local folding nucleus.     

Male cognitive performance declines in the absence of sexual selection

Sexual selection is responsible for the evolution of male ornaments and armaments, but its role in the evolution of cognition—the ability to process, retain and use information—is largely unexplored. Because successful courtship is likely to involve processing information in complex, competitive sexual environments, we hypothesized that sexual selection contributes to the evolution and maintenance of cognitive abilities in males. To test this, we removed mate choice and mate competition from experimental populations of Drosophila melanogaster by enforcing monogamy for over 100 generations. Males evolved under monogamy became less proficient than polygamous control males at relatively complex cognitive tasks. When faced with one receptive and several unreceptive females, polygamous males quickly focused on receptive females, whereas monogamous males continued to direct substantial courtship effort towards unreceptive females. As a result, monogamous males were less successful in this complex setting, despite being as quick to mate as their polygamous counterparts with only one receptive female. This diminished ability to use past information was not limited to the courtship context: monogamous males (but not females) also showed reduced aversive olfactory learning ability. Our results provide direct experimental evidence that the intensity of sexual selection is an important factor in the evolution of male cognitive ability.     

Simultaneous NMR assignment of backbone and side chain amides in large proteins with IS-TROSY

A new strategy for the simultaneous NMR assignment of both backbone and side chain amides in large proteins with isotopomer-selective transverse-relaxation-optimized spectroscopy (IS-TROSY) is reported. The method considers aspects of both the NMR sample preparation and the experimental design. First, the protein is dissolved in a buffer with 50%H₂O/50%D₂O in order to promote the population of semideuterated NHD isotopomers in side chain amides of Asn/Gln residues. Second, a ¹³C′-coupled 2D ¹⁵N–¹H IS-TROSY spectrum provides a stereospecific distinction between the geminal protons in the E and Z configurations of the carboxyamide group. Third, a suite of IS-TROSY-based triple-resonance NMR experiments, e.g. 3D IS-TROSY-HNCA and 3D IS-TROSY-HNCACB, are designed to correlate aliphatic carbon atoms with backbone amides and, for Asn/Gln residues, at the same time with side chain amides. The NMR assignment procedure is similar to that for small proteins using conventional 3D HNCA/3D HNCACB spectra, in which, however, signals from NH₂ groups are often very weak or even missing due to the use of broad-band proton decoupling schemes and NOE data have to be used as a remedy. For large proteins, the use of conventional TROSY experiments makes resonances of side chain amides not observable at all. The application of IS-TROSY experiments to the 35-kDa yeast cytosine deaminase has established a complete resonance assignment for the backbone and stereospecific assignment for side chain amides, which otherwise could not be achieved with existing NMR experiments. Thus, the development of IS-TROSY-based method provides new opportunities for the NMR study of important structural and biological roles of carboxyamides and side chain moieties of arginine and lysine residues in large proteins as well as amino moieties in nucleic acids.Electronic Supplementary Material Supplementary material is available to authorised users in the online version of this article at .     

Prospects for resonance assignments in multidimensional solid-state NMR spectra of uniformly labeled proteins

Summary The feasibility of assigning the backbone ¹⁵N and ¹³C NMR chemical shifts in multidimensional magic angle spinning NMR spectra of uniformly isotopically labeled proteins and peptides in unoriented solid samples is assessed by means of numerical simulations. The goal of these simulations is to examine how the upper limit on the size of a peptide for which unique assignments can be made depends on the spectral resolution, i.e., the NMR line widths. Sets of simulated three-dimensional chemical shift correlation spectra for artificial peptides of varying length are constructed from published liquid-state NMR chemical shift data for ubiquitin, a well-characterized soluble protein. Resonance assignments consistent with these spectra to within the assumed spectral resolution are found by a numerical search algorithm. The dependence of the number of consistent assignments on the assumed spectral resolution and on the length of the peptide is reported. If only three-dimensional chemical shift correlation data for backbone ¹⁵N and ¹³C nuclei are used, and no residue-specific chemical shift information, information from amino acid side-chain signals, and proton chemical shift information are available, a spectral resolution of 1 ppm or less is generally required for a unique assignment of backbone chemical shifts for a peptide of 30 amino acid residues.     

A simple and effective NMR cell for studies of encapsulated proteins dissolved in low viscosity solvents

Application of triple-resonance and isotope-edited-NOE methods to the study of increasingly larger macromolecules and their complexes remains a central goal of solution NMR spectroscopy. The slow reorientational motion of larger molecules leads to rapid transverse relaxation and results in losses in both resolution and sensitivity of multidimensional-multinuclear solution NMR experiments. A recently described technique employs a physical approach to increase the tumbling rate of macromolecules in an attempt to preserve access to the full range of structural restraints available to studies of smaller systems. This technique involves encapsulation of a hydrated protein in a surfactant shell which is subsequently solubilized in a low viscosity solvent. A simple, efficient and cost effective NMR cell that accommodates the moderate liquefaction pressures required in the encapsulation method is described. Application of the method to the 56 kD triose phosphate isomerase homodimer is demonstrated.