A large scale test of computational protein design: folding and stability of nine completely redesigned globular proteins

A previously developed computer program for protein design, RosettaDesign, was used to predict low free energy sequences for nine naturally occurring protein backbones. RosettaDesign had no knowledge of the naturally occurring sequences and on average 65% of the residues in the designed sequences differ from wild-type. Synthetic genes for ten completely redesigned proteins were generated, and the proteins were expressed, purified, and then characterized using circular dichroism, chemical and temperature denaturation and NMR experiments. Although high-resolution structures have not yet been determined, eight of these proteins appear to be folded and their circular dichroism spectra are similar to those of their wild-type counterparts. Six of the proteins have stabilities equal to or up to 7kcal/mol greater than their wild-type counterparts, and four of the proteins have NMR spectra consistent with a well-packed, rigid structure. These encouraging results indicate that the computational protein design methods can, with significant reliability, identify amino acid sequences compatible with a target protein backbone.     

Computer-based redesign of a protein folding pathway.

A fundamental test of our current understanding of protein folding is to rationally redesign protein folding pathways. We use a computer-based design strategy to switch the folding pathway of protein G, which normally involves formation of the second, but not the first, beta-turn at the rate limiting step in folding. Backbone conformations and amino acid sequences that maximize the interaction density in the first beta-hairpin were identified, and two variants containing 11 amino acid replacements were found to be approximately 4 kcal mol-1 more stable than wild type protein G. Kinetic studies show that the redesigned proteins fold approximately 100 x faster than wild type protein and that the first beta-turn is formed and the second disrupted at the rate limiting step in folding.     

Microsatellite variation in Drosophila melanogaster and Drosophila simulans: a reciprocal test of the ascertainment bias hypothesis

Interspecific comparisons of microsatellite loci have repeatedly shown that the loci are longer and more variable in the species from which they are derived (the focal species) than are homologous loci in other (nonfocal) species. There is debate as to whether this is due to directional evolution or to an ascertainment bias during the cloning and locus selection processes. This study tests these hypotheses by performing a reciprocal study. Eighteen perfect dinucleotide microsatellite loci identified from a Drosophila simulans library screen and 18 previously identified in an identical Drosophila melanogaster library screen were used to survey natural populations of each species. No difference between focal and nonfocal species was observed for mean PCR fragment length. However, heterozygosity and number of alleles were significantly higher in the focal species than in the nonfocal species. The most common allele in the Zimbabwe population of both species was sequenced for 31 of the 36 loci. The length of the longest stretch of perfect repeat units is, on average, longer in the focal species than in the non-focal species. There is a positive correlation between the length of the longest stretch of perfect repeats and heterozygosity. The difference in heterozygosity can thus be explained by a reduction in the length of the longest stretch of perfect repeats in the nonfocal species. Furthermore, flanking-sequence length difference was noted between the two species at 58% of the loci sequenced. These data do not support the predictions of the directional-evolution hypothesis; however, consistent with the ascertainment bias hypothesis, the lower variability in nonfocal species is an artifact of the microsatellite cloning and isolation process. Our results also suggest that the magnitude of ascertainment bias for repeat unit length is a function of the microsatellite size distribution in the genomes of different species.      

Inhibition of L- and P-selectin by a rationally synthesized novel core 2-like branched structure containing GalNAc-Lewisx and Neu5Acalpha2- 3Galbeta1-3GalNAc sequences

The selectins interact in important normal and pathological situations with certain sialylated, fucosylated glycoconjugate ligands containing sialyl Lewisx(Neu5Acalpha2-3Galbeta1-4(Fucalpha1-3)GlcN Ac). Much effort has gone into the synthesis of sialylated and sulfated Lewisxanalogs as competitive ligands for the selectins. Since the natural selectin ligands GlyCAM-1 and PSGL-1 carry sialyl Lewisxas part of a branched Core 2 O-linked structure, we recently synthesized Galbeta1-4(Fucalpha1-3)GlcNAcbeta1-6(SE-3Galbeta1++ +-3)GalNAc1alphaOMe and found it to be a moderately superior ligand for L and P-selectin (Koenig et al. , Glycobiology 7, 79-93, 1997). Other studies have shown that sulfate esters can replace sialic acid in some selectin ligands (Yeun et al. , Biochemistry, 31, 9126-9131, 1992; Imai et al. , Nature, 361, 555, 1993). Based upon these observations, we hypothesized that Neu5Acalpha2-3Galbeta1-3GalNAc might have the capability of interacting with L- and P-selectin. To examine this hypothesis, we synthesized Galbeta1-4(Fucalpha1-3)GlcNAcbeta1-6(Neu5Acalpha2++ +-3Galbeta1-3)- GalNAc alpha1-OB, which was found to be 2- to 3-fold better than sialyl Lexfor P and L selectin, respectively. We also report the synthesis of an unusual structure GalNAcbeta1-4(Fucalpha1- 3)GlcNAcbeta1-OMe (GalNAc- Lewisx-O-methyl glycoside), which also proved to be a better inhibitor of L- and P-selectin than sialyl Lewisx-OMe. Combining this with our knowledge of Core 2 branched structures, we have synthesized a molecule that is 5- to 6-fold better at inhibiting L- and P-selectin than sialyl Lewisx-OMe, By contrast to unbranched structures, substitution of a sulfate ester group for a sialic acid residue in such a molecule resulted in a considerable loss of inhibition ability. Thus, the combination of a sialic acid residue on the primary (beta1-3) arm, and a modified Lexunit on the branched (beta1-6) arm on an O-linked Core 2 structure generated a monovalent synthetic oliogosaccharide inhibitor superior to SLexfor both L- and P-selectin.      

A multinational study of alpha-lactalbumin concentrations in human milk

Alpha-lactalbumin, a 14-kD protein, plays a central biochemical role in the mammary gland as the regulatory subunit of lactose synthase, and also plays a nutritional role for the rapidly growing neonate as the protein in highest concentration in human milk. The current study was undertaken to better characterize alpha-lactalbumin concentrations in human milk from a variety of countries. Mature human milk (lactation duration > or =1 month) was collected from at least 50 women from nine different countries on five continents. Alpha-lactalbumin concentration was determined by HPLC. The mean +/- SD for 452 samples was 2.44 +/- 0.64 g/L. The mean value of the samples from the United States was significantly higher than that from any other country, and the mean in Mexico was significantly lower than that from every country except China and Canada. Alpha-lactalbumin concentration decreased with increasing duration of lactation and was positively correlated with total nitrogen. On average, alpha-lactalbumin contributed 16% of the total nitrogen content of human milk and consequently an important part of the amino acid content.     

Detection of a single nucleotide polymorphism in the human alpha-lactalbumin gene: implications for human milk proteins

Variability in the protein composition of breast milk has been observed in many women and is believed to be due to natural variation of the human population. Single nucleotide polymorphisms (SNPs) are present throughout the entire human genome, but the impact of this variation on human milk composition and biological activity and infant nutrition and health is unclear. The goals of this study were to characterize a variant of human alpha-lactalbumin observed in milk from a Filipino population by determining the location of the polymorphism in the amino acid and genomic sequences of alpha-lactalbumin. Milk and blood samples were collected from 20 Filipino women, and milk samples were collected from an additional 450 women from nine different countries. alpha-Lactalbumin concentration was measured by high-performance liquid chromatography (HPLC), and milk samples containing the variant form of the protein were identified with both HPLC and mass spectrometry (MS). The molecular weight of the variant form was measured by MS, and the location of the polymorphism was narrowed down by protein reduction, alkylation and trypsin digestion. Genomic DNA was isolated from whole blood, and the polymorphism location and subject genotype were determined by amplifying the entire coding sequence of human alpha-lactalbumin by PCR, followed by DNA sequencing. A variant form of alpha-lactalbumin was observed in HPLC chromatograms, and the difference in molecular weight was determined by MS (wild type=14,070 Da, variant=14,056 Da). Protein reduction and digestion narrowed the polymorphism between the 33rd and 77th amino acid of the protein. The genetic polymorphism was identified as adenine to guanine, which translates to a substitution from isoleucine to valine at amino acid 46. The frequency of variation was higher in milk from China, Japan and Philippines, which suggests that this polymorphism is most prevalent in Asia. There are SNPs in the genome for human milk proteins and their implications for protein bioactivity and infant nutrition need to be considered.     

A "solvated rotamer" approach to modeling water-mediated hydrogen bonds at protein-protein interfaces

Water-mediated hydrogen bonds play critical roles at protein-protein and protein-nucleic acid interfaces, and the interactions formed by discrete water molecules cannot be captured using continuum solvent models. We describe a simple model for the energetics of water-mediated hydrogen bonds, and show that, together with knowledge of the positions of buried water molecules observed in X-ray crystal structures, the model improves the prediction of free-energy changes upon mutation at protein-protein interfaces, and the recovery of native amino acid sequences in protein interface design calculations. We then describe a "solvated rotamer" approach to efficiently predict the positions of water molecules, at protein-protein interfaces and in monomeric proteins, that is compatible with widely used rotamer-based side-chain packing and protein design algorithms. Finally, we examine the extent to which the predicted water molecules can be used to improve prediction of amino acid identities and protein-protein interface stability, and discuss avenues for overcoming current limitations of the approach.     

Fab-based bispecific antibody formats with robust biophysical properties and biological activity

A myriad of innovative bispecific antibody (BsAb) platforms have been reported. Most require significant protein engineering to be viable from a development and manufacturing perspective. Single-chain variable fragments (scFvs) and diabodies that consist only of antibody variable domains have been used as building blocks for making BsAbs for decades. The drawback with Fv-only moieties is that they lack the native-like interactions with C_H1/C_L domains that make antibody Fab regions stable and soluble. Here, we utilize a redesigned Fab interface to explore 2 novel Fab-based BsAbs platforms. The redesigned Fab interface designs limit heavy and light chain mixing when 2 Fabs are co-expressed simultaneously, thus allowing the use of 2 different Fabs within a BsAb construct without the requirement of one or more scFvs. We describe the stability and activity of a HER2×HER2 IgG-Fab BsAb, and compare its biophysical and activity properties with those of an IgG-scFv that utilizes the variable domains of the same parental antibodies. We also generated an EGFR × CD3 tandem Fab protein with a similar format to a tandem scFv (otherwise known as a bispecific T cell engager or BiTE). We show that the Fab-based BsAbs have superior biophysical properties compared to the scFv-based BsAbs. Additionally, the Fab-based BsAbs do not simply recapitulate the activity of their scFv counterparts, but are shown to possess unique biological activity.