Origin of the anomalous circular dichroism spectra of many apomyoglobin mutants

Several authors have reported that many sperm whale apomyoglobin mutants show anomalous circular dichroism spectra. These mutants have a low molar ellipticity compared to the wild-type protein but in several cases have the same stability of unfolding. A model in which native apomyoglobin is not folded in the same manner as that in other proteins and in which mutants show progressive reductions in their degree of folding has been suggested to explain this phenomenon. However, nuclear magnetic resonance of the native apomyoglobin conformation has shown that this state is folded and compact, raising the possibility that the anomalous circular dichroism spectra could have another explanation. We studied several mutants with anomalous circular dichroism spectra and found that these proteins were all contaminated with nucleic acid that contributed to the ultraviolet absorption and caused uncertainty in the determination of protein concentration. The resulting overestimation of the concentration of apomyoglobin explains the phenomenon of anomalous circular dichroism spectra. We describe a procedure to remove the contaminant nucleic acid which yields accurate protein concentration measurements and provides the normal circular dichroism spectra. Our findings support a well-structured native conformation for apomyoglobin and may also be of the interest to scientists working with the purification of recombinant proteins.     

Computational protein design with explicit consideration of surface hydrophobic patches

De novo protein design requires the identification of amino-acid sequences that favor the target-folded conformation and are soluble in water. One strategy for promoting solubility is to disallow hydrophobic residues on the protein surface during design. However, naturally occurring proteins often have hydrophobic amino acids on their surface that contribute to protein stability via the partial burial of hydrophobic surface area or play a key role in the formation of protein-protein interactions. A less restrictive approach for surface design that is used by the modeling program Rosetta is to parameterize the energy function so that the number of hydrophobic amino acids designed on the protein surface is similar to what is observed in naturally occurring monomeric proteins. Previous studies with Rosetta have shown that this limits surface hydrophobics to the naturally occurring frequency (～28%), but that it does not prevent the formation of hydrophobic patches that are considerably larger than those observed in naturally occurring proteins. Here, we describe a new score term that explicitly detects and penalizes the formation of hydrophobic patches during computational protein design. With the new term, we are able to design protein surfaces that include hydrophobic amino acids at naturally occurring frequencies, but do not have large hydrophobic patches. By adjusting the strength of the new score term, the emphasis of surface redesigns can be switched between maintaining solubility and maximizing folding free energy.     

Full-sequence computational design and solution structure of a thermostable protein variant

Computational protein design procedures were applied to the redesign of the entire sequence of a 51 amino acid residue protein, Drosophila melanogaster engrailed homeodomain. Various sequence optimization algorithms were compared and two resulting designed sequences were experimentally evaluated. The two sequences differ by 11 mutations and share 22% and 24% sequence identity with the wild-type protein. Both computationally designed proteins were considerably more stable than the naturally occurring protein, with midpoints of thermal denaturation greater than 99 degrees C. The solution structure was determined for one of the two sequences using multidimensional heteronuclear NMR spectroscopy, and the structure was found to closely match the original design template scaffold.     

Circular dichroism studies of an oligo-alpha-thymidylate and of its interactions with complementary sequences.

An octathymidylate synthesized with the alpha anomer of thymidine has been studied using circular dichroism. Its conformation has been compared to that of its analogue containing the naturally occurring beta anomer. In both compounds some degree of intramolecular stacking is present as indicated by the shapes of the circular dichroism spectra and their variations with temperature. As its beta-analogue the alpha-octathymidylate binds to its complementary sequences containing beta-nucleosides. Only complexes with 1A:1T stoechiometry were observed. Binding to ribose- containing oligomers and polymers is much stronger than binding to deoxyribose-containing analogues. Circular dichroism spectra provided evidence for a difference between the geometry of the various complexes formed with an alpha-oligothymidylate and those formed with its beta-anomer-containing analogue.     

The highly cooperative folding of small naturally occurring proteins is likely the result of natural selection

To illuminate the evolutionary pressure acting on the folding free energy landscapes of naturally occurring proteins, we have systematically characterized the folding free energy landscape of Top7, a computationally designed protein lacking an evolutionary history. Stopped-flow kinetics, circular dichroism, and NMR experiments reveal that there are at least three distinct phases in the folding of Top7, that a nonnative conformation is stable at equilibrium, and that multiple fragments of Top7 are stable in isolation. These results indicate that the folding of Top7 is significantly less cooperative than the folding of similarly sized naturally occurring proteins, suggesting that the cooperative folding and smooth free energy landscapes observed for small naturally occurring proteins are not general properties of polypeptide chains that fold to unique stable structures but are instead a product of natural selection.     

Engineering strand-specific DNA nicking enzymes from the type IIS restriction endonucleases BsaI, BsmBI, and BsmAI

To determine the extent to which protein folding rates and free energy landscapes have been shaped by natural selection, we have examined the folding kinetics of five proteins generated using computational design methods and, hence, never exposed to natural selection. Four of these proteins are complete computer-generated redesigns of naturally occurring structures and the fifth protein, called Top7, has a computer-generated fold not yet observed in nature. We find that three of the four redesigned proteins fold much faster than their naturally occurring counterparts. While natural selection thus does not appear to operate on protein folding rates, the majority of the designed proteins unfold considerably faster than their naturally occurring counterparts, suggesting possible selection for a high free energy barrier to unfolding. In contrast to almost all naturally occurring proteins of less than 100 residues but consistent with simple computational models, the folding energy landscape for Top7 appears to be quite complex, suggesting the smooth energy landscapes and highly cooperative folding transitions observed for small naturally occurring proteins may also reflect the workings of natural selection.     

Estimation of the number of alpha-helical and beta-strand segments in proteins using circular dichroism spectroscopy

A simple approach to estimate the number of alpha-helical and beta-strand segments from protein circular dichroism spectra is described. The alpha-helix and beta-sheet conformations in globular protein structures, assigned by DSSP and STRIDE algorithms, were divided into regular and distorted fractions by considering a certain number of terminal residues in a given alpha-helix or beta-strand segment to be distorted. The resulting secondary structure fractions for 29 reference proteins were used in the analyses of circular dichroism spectra by the SELCON method. From the performance indices of the analyses, we determined that, on an average, four residues per alpha-helix and two residues per beta-strand may be considered distorted in proteins. The number of alpha-helical and beta-strand segments and their average length in a given protein were estimated from the fraction of distorted alpha-helix and beta-strand conformations determined from the analysis of circular dichroism spectra. The statistical test for the reference protein set shows the high reliability of such a classification of protein secondary structure. The method was used to analyze the circular dichroism spectra of four additional proteins and the predicted structural characteristics agree with the crystal structure data.     

A tetrapeptide fragment-based design method results in highly stable artificial proteins

Computational protein design has progressed rapidly over the last years. A number of design methods have been proposed and tested. In this paper, we report the successful application of a fragment-based method for protein design. The method uses statistical information on tetrapeptide backbone conformations. The previously published artificial fold of TOP 7 (Kuhlman et al., Science, 2003; 302:1364-1368) was chosen as template. A series of polypeptide sequences were created that were predicted to fold into this target structure. Two of the designed proteins, M5 and M7, were expressed and characterized by fluorescence spectroscopy, circular dichroism and NMR. They showed the hallmarks of well-ordered tertiary structure as well as cooperative folding/unfolding transitions. Furthermore, the two novel proteins were found to be highly stable against temperature and denaturant-induced unfolding.     

Improving computational protein design by using structure‐derived sequence profile

Designing a protein sequence that will fold into a predefined structure is of both practical and fundamental interest. Many successful, computational designs in the last decade resulted from improved understanding of hydrophobic and polar interactions between side chains of amino acid residues in stabilizing protein tertiary structures. However, the coupling between main‐chain backbone structure and local sequence has yet to be fully addressed. Here, we attempt to account for such coupling by using a sequence profile derived from the sequences of five residue fragments in a fragment library that are structurally matched to the five‐residue segments contained in a target structure. We further introduced a term to reduce low complexity regions of designed sequences. These two terms together with optimized reference states for amino‐acid residues were implemented in the RosettaDesign program. The new method, called RosettaDesign‐SR, makes a 12% increase (from 34 to 46%) in fraction of proteins whose designed sequences are more than 35% identical to wild‐type sequences. Meanwhile, it reduces 8% (from 22% to 14%) to the number of designed sequences that are not homologous to any known protein sequences according to psi‐blast. More importantly, the sequences designed by RosettaDesign‐SR have 2–3% more polar residues at the surface and core regions of proteins and these surface and core polar residues have about 4% higher sequence identity to wild‐type sequences than by RosettaDesign. Thus, the proteins designed by RosettaDesign‐SR should be less likely to aggregate and more likely to have unique structures due to more specific polar interactions. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

X‐ray vs. NMR structures as templates for computational protein design

Certain protein‐design calculations involve using an experimentally determined high‐resolution structure as a template to identify new sequences that can adopt the same fold. This approach has led to the successful design of many novel, well‐folded, native‐like proteins. Although any atomic‐resolution structure can serve as a template in such calculations, most successful designs have used high‐resolution crystal structures. Because there are many proteins for which crystal structures are not available, it is of interest whether nuclear magnetic resonance (NMR) templates are also appropriate. We have analyzed differences between using X‐ray and NMR templates in side‐chain repacking and design calculations. We assembled a database of 29 proteins for which both a high‐resolution X‐ray structure and an ensemble of NMR structures are available. Using these pairs, we compared the rotamericity, χ₁‐angle recovery, and native‐sequence recovery of X‐ray and NMR templates. We carried out design using RosettaDesign on both types of templates, and compared the energies and packing qualities of the resulting structures. Overall, the X‐ray structures were better templates for use with Rosetta. However, for ～20% of proteins, a member of the reported NMR ensemble gave rise to designs with similar properties. Re‐evaluating RosettaDesign structures with other energy functions indicated much smaller differences between the two types of templates. Ultimately, experiments are required to confirm the utility of particular X‐ray and NMR templates. But our data suggest that the lack of a high‐resolution X‐ray structure should not preclude attempts at computational design if an NMR ensemble is available. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

On the analysis of membrane protein circular dichroism spectra

Analysis of circular dichroism spectra of proteins provides information about protein secondary structure. Analytical methods developed for such an analysis use structures and spectra of a set of reference proteins. The reference protein sets currently in use include soluble proteins with a wide range of secondary structures, and perform quite well in analyzing CD spectra of soluble proteins. The utility of soluble protein reference sets in analyzing membrane protein CD spectra, however, has been questioned in a recent study that found current reference protein sets to be inadequate for analyzing membrane proteins. We have examined the performance of reference protein sets available in the CDPro software package for analyzing CD spectra of 13 membrane proteins with available crystal structures. Our results indicate that the reference protein sets currently available for CD analysis perform reasonably well in analyzing membrane protein CD spectra, with performance indices comparable to those for soluble proteins. Soluble + membrane protein reference sets, which were constructed by combining membrane proteins with soluble protein reference sets, gave improved performance in both soluble and membrane protein CD analysis.     

Isolation, folding and structural investigations of the amino acid transporter OEP16

Membrane proteins compose more than 30% of all proteins in the living cell. However, many membrane proteins have low abundance in the cell and cannot be isolated from natural sources in concentrations suitable for structure analysis. The overexpression, reconstitution, and stabilization of membrane proteins are complex and remain a formidable challenge in membrane protein characterization. Here we describe a novel, in vitro folding procedure for a cation-selective channel protein, the outer envelope membrane protein 16 (OEP16) of pea chloroplast, overexpressed in Escherichia coli in the form of inclusion bodies. The protein is purified and then folded with detergent on a Ni–NTA affinity column. Final concentrations of reconstituted OEP16 of up to 24 mg/ml have been achieved, which provides samples that are sufficient for structural studies by NMR and crystallography. Reconstitution of OEP16 in detergent micelles was monitored by circular dichroism, fluorescence, and NMR spectroscopy. Tryptophan fluorescence spectra of heterologous expressed OEP16 in micelles are similar to spectra of functionally active OEP16 in liposomes, which indicates folding of the membrane protein in detergent micelles. CD spectroscopy studies demonstrate a folded protein consisting primarily of α-helices. ¹⁵N-HSQC NMR spectra also provide evidence for a folded protein. We present here a convenient, effective and quantitative method to screen large numbers of conditions for optimal protein stability by using microdialysis chambers in combination with fluorescence spectroscopy. Recent collection of multidimensional NMR data at 500, 600 and 800 MHz demonstrated that the protein is suitable for structure determination by NMR and stable for weeks during data collection.     

Equilibrium collapse and the kinetic 'foldability' of proteins.

An important element of protein folding theory has been the identification of equilibrium parameters that might uniquely distinguish rapidly folding polypeptide sequences from those that fold slowly. One such parameter, termed sigma, is a dimensionless, equilibrium measure of the coincidence of chain compaction and folding that is predicted to be an important determinant of relative folding kinetics. To test this prediction and improve our understanding of the putative relationship between nonspecific compaction of the unfolded state and protein folding kinetics, we have used small-angle X-ray scattering and circular dichroism spectroscopy to measure the sigma of five well-characterized proteins. Consistent with theoretical predictions, we find that near-perfect coincidence of the unfolded state contraction and folding (sigma approximately 0) is associated with the rapid kinetics of these naturally occurring proteins. We do not, however, observe any significant correlation between sigma and either the relative folding rates of these proteins or the presence or absence of well-populated kinetic intermediates. Thus, while sigma approximately 0 may be a necessary condition to ensure rapid folding, differences in sigma do not account for the wide range of rates and mechanisms with which naturally occurring proteins fold.     

Circular dichroism analyses of membrane proteins: an examination of differential light scattering and absorption flattening effects in large membrane vesicles and membrane sheets

The circular dichroism spectra of membrane suspensions are distorted by differential light scattering and absorption flattening effects, which arise as a consequence of the large size of the membrane particles relative to the wavelength of light and the high concentration of proteins in the membranes. In this paper, the consequences of these phenomena on the protein spectra of large membrane particles are discussed, and methods for eliminating them are examined. The distortions due to differential light scattering are relatively small in membrane systems, and can be compensated for by use of a large detector acceptance angle geometry. Several methods for correcting for differential flattening, which introduces a substantial distortion, have been evaluated, and a new method, the flattening quotient approach, which produces by far the best results, is described. Since the secondary structures calculated from circular dichroism spectra are highly dependent on accurate spectral shape and magnitude, this method for correcting the spectra may find general application in circular dichroism studies of membrane proteins.     

Chloranthones A – D: Minor and Unprecedented Dinor‐Eudesmenes from Chloranthus elatior

Four novel naturally occurring diastereoisomers of dinor‐eudesmenes, named chloranthones A–D ( 1 – 4 , resp.), were isolated as minor components from the EtOH extract of the aerial parts of Chloranthus elatior. The unprecedented framework was established using extensive 2D‐NMR techniques. Their absolute configurations were deduced from the observed Cotton effects in their circular dichroism (CD) spectra. A plausible biosynthetic pathway to the dinor‐eudesmenes is proposed.     

Induction of protein-like molecular architecture by self-assembly processes

One of the most intriguing self-assembly processes is the folding of peptide chains into native protein structures. We have developed a method for building protein-like structural motifs that incorporate sequences of biological interest. A lipophilic moiety is attached onto an N(alpha)-amino group of a peptide chain, resulting in a 'peptide-amphiphile'. The alignment of amphiphilic compounds at the lipid solvent interface is used to facilitate peptide alignment and structure initiation and propagation. Peptide-amphiphiles containing potentially triple-helical structural motifs have been synthesized. The resultant head group structures have been characterized by circular dichroism and NMR spectroscopies. Evidence for a self-assembly process of peptide-amphiphiles has been obtained from: (a) circular dichroism spectra and melting curves characteristic of triple-helices, (b) one- and two-dimensional NMR spectra indicative of stable triple-helical structure at low temperatures and melted triple-helices at high temperatures, and (c) pulsed-field gradient NMR experiments demonstrating different self-diffusion coefficients between proposed triple-helical and non-triple-helical species. The peptide-amphiphiles described here provide a simple approach for building stable protein structural motifs using peptide head groups.     

Relation between structure and function in some partially synthetic ribonucleases S'. Enzymic and spectroscopic investigation on [Orn10, Asn14]-RNase S' and 1epsilon, 7epsilon, 10delta-triguanidino-[Orn10, Asn14]-RNase S'.

Some analogues have been prepared of S-peptide, the peptide obtained together with S-protein from subtilisn-modified beef pancreatic R Nase A. The syntheses are described of [Orn10, Asn14]-S-peptide and 1epsilon, 7epsilon, 10delta-triguanidino-[Orn10, Asn14]-S-peptide. The S-peptide analogues are able to activate S-protein at the level of the parent [Orn10]-S-peptide and 1epsilon, 7epsilon-diguanidino-S-peptide respectively, although at high peptide-to-protein molar ratios. After their recombination with S-protein the buried character of Tyr-25 was restored, as judged from difference absorption and circular dichroism spectra in the near-ultraviolet region. These findings indicate that the asparaginyl residue is a possible naturally occurring substituent in the R Nase A sequences whose state of amidation in position 14 has not yet been defined.     

A de novo redesign of the WW domain

We have used a sequence prediction algorithm and a novel sampling method to design protein sequences for the WW domain, a small beta-sheet motif. The procedure, referred to as SPANS, designs sequences to be compatible with an ensemble of closely related polypeptide backbones, mimicking the inherent flexibility of proteins. Two designed sequences (termed SPANS-WW1 and SPANS-WW2), using only naturally occurring L-amino acids, were selected for study and the corresponding polypeptides were prepared in Escherichia coli. Circular dichroism data suggested that both purified polypeptides adopted secondary structure features related to those of the target without the aid of disulfide bridges or bound cofactors. The structure exhibited by SPANS-WW2 melted cooperatively by raising the temperature of the solution. Further analysis of this polypeptide by proton nuclear magnetic resonance spectroscopy demonstrated that at 5 degrees C, it folds into a structure closely resembling a natural WW domain. This achievement constitutes one of a small number of successful de novo protein designs through fully automated computational methods and highlights the feasibility of including backbone flexibility in the design strategy.     

Solution NMR structures of IgG binding domains with artificially evolved high levels of sequence identity but different folds

We describe here the solution NMR structures of two IgG binding domains with highly homologous sequences but different three-dimensional structures. The proteins, G311 and A219, are derived from the IgG binding domains of their wild-type counterparts, protein G and protein A, respectively. Through a series of site-directed mutations and phage display selections, the sequences of G311 and A219 were designed to converge to a point of high-level sequence identity while keeping their respective wild-type tertiary folds. Structures of both artificially evolved sequences were determined by NMR spectroscopy. The main chain fold of G311 can be superimposed on the wild-type alpha/beta protein G structure with a backbone rmsd of 1.4 A, and the A219 structure can be overlaid on the wild-type three-alpha-helix protein A fold also with a backbone rmsd of 1.4 A. The structure of G311, in particular, accommodates a large number of mutational changes without undergoing a change in the overall fold of the main chain. The structural differences are maintained despite a high level (59%) of sequence identity. These proteins serve as starting points for further experiments that will probe basic concepts of protein folding and conformational switching.     

NMR structure of the bovine prion protein isolated from healthy calf brains

NMR structures of recombinant prion proteins from various species expressed in Escherichia coli have been solved during the past years, but the fundamental question of the relevancy of these data relative to the naturally occurring forms of the prion protein has not been directly addressed. Here, we present a comparison of the cellular form of the bovine prion protein isolated and purified from healthy calf brains without use of detergents, so that it contains the two carbohydrate moieties and the part of the GPI anchor that is maintained after enzymatic cleavage of the glycerolipid moiety, with the recombinant bovine prion protein expressed in E. coli. We show by circular dichroism and (1)H-NMR spectroscopy that the three-dimensional structure and the thermal stability of the natural glycoprotein and the recombinant polypeptide are essentially identical. This result indicates possible functional roles of the glycosylation of prion proteins in healthy organisms, and provides a platform and validation for future work on the structural biology of prion proteins, which will have to rely primarily on the use of recombinant polypeptides.