Sequence and structural analysis of cellular retinoic acid-binding proteins reveals a network of conserved hydrophobic interactions

Regulation of Absorption and Emission in a Protein/Fluorophore Complex

Human cellular retinol binding protein II (hCRBPII) was used as a protein engineering platform to rationally regulate absorptive and emissive properties of a covalently bound fluorogenic dye. We demonstrate the binding of a thio-dapoxyl analog via formation of a protonated imine between an active site lysine residue and the chromophore's aldehyde. Rational manipulation of the electrostatics of the binding pocket results in a 204 nm shift in absorption and a 131 nm shift in emission. The protein is readily expressed in mammalian systems and binds with exogenously delivered fluorophore as demonstrated by live-cell imaging experiments.

A candidate antibody drug for prevention of malaria

Over 75% of malaria-attributable deaths occur in children under the age of 5 years. However, the first malaria vaccine recommended by the World Health Organization (WHO) for pediatric use, RTS,S/AS01 (Mosquirix), has modest efficacy. Complementary strategies, including monoclonal antibodies, will be important in efforts to eradicate malaria. Here we characterize the circulating B cell repertoires of 45 RTS,S/AS01 vaccinees and discover monoclonal antibodies for development as potential therapeutics. We generated >28,000 antibody sequences and tested 481 antibodies for binding activity and 125 antibodies for antimalaria activity in vivo. Through these analyses we identified correlations suggesting that sequences in Plasmodium falciparum circumsporozoite protein, the target antigen in RTS,S/AS01, may induce immunodominant antibody responses that limit more protective, but subdominant, responses. Using binding studies, mouse malaria models, biomanufacturing assessments and protein stability assays, we selected AB-000224 and AB-007088 for advancement as a clinical lead and backup. We engineered the variable domains (Fv) of both antibodies to enable low-cost manufacturing at scale for distribution to pediatric populations, in alignment with WHO's preferred product guidelines. The engineered clone with the optimal manufacturing and drug property profile, MAM01, was advanced into clinical development.

Function of the Conserved Non-Functional Residues in Apomyoglobin - to Determine and to Preserve Correct Topology of the Protein

In this paper the answer to O. B. Ptitsyn's question "What is the role of conserved non-functional residues in apomyoglobin" is presented, which is based on the research results of three laboratories. The role of conserved non-functional apomyoglobin residues in formation of native topology in the molten globule state of this protein is revealed. This fact allows suggesting that the conserved non-functional residues in this protein are indispensable for fixation and maintaining main elements of the correct topology of its secondary structure in the intermediate state. The correct topology is a native element in the intermediate state of the protein.

Framework Mutations of the 10-1074 bnAb Increase Conformational Stability, Manufacturability, and Stability While Preserving Full Neutralization Activity

The broadly neutralizing anti-HIV antibody, 10-1074, is a highly somatically hypermutated IgG1 being developed for prophylaxis in sub-Saharan Africa. A series of algorithms were applied to identify potentially destabilizing residues in the framework of the Fv region. Of 17 residues defined, a variant was identified encompassing 1 light and 3 heavy chain residues, with significantly increased conformational stability while maintaining full neutralization activity. Central to the stabilization was the replacement of the heavy chain residue T108 with R108 at the base of the CDR3 loop which allowed for the formation of a nascent salt bridge with heavy chain residue D137. Three additional mutations were necessary to confer increased conformational stability as evidenced by differential scanning fluorimetry and isothermal chemical unfolding. In addition, we observed increased stability during low pH incubation in which 40% of the parental monomer aggregated while the combinatorial variant showed no increase in aggregation. Incubation of the variant at 100 mg/mL for 6 weeks at 40°C showed a 9-fold decrease in subvisible particles ≥2 μm relative to the parental molecule. Stability-based designs have also translated to improved pharmacokinetics. Together, these data show that increasing conformational stability of the Fab can have profound effects on the manufacturability and long-term stability of a monoclonal antibody.

Domain-Swapped Dimers of Intracellular Lipid-Binding Proteins: Evidence for Ordered Folding Intermediates

Human Cellular Retinol Binding Protein II (hCRBPII), a member of the intracellular lipid-binding protein family, is a monomeric protein responsible for the intracellular transport of retinol and retinal. Herein we report that hCRBPII forms an extensive domain-swapped dimer during bacterial expression. The domain-swapped region encompasses almost half of the protein. The dimer represents a novel structural architecture with the mouths of the two binding cavities facing each other, producing a new binding cavity that spans the length of the protein complex. Although wild-type hCRBPII forms the dimer, the propensity for dimerization can be substantially increased via mutation at Tyr60. The monomeric form of the wild-type protein represents the thermodynamically more stable species, making the domain-swapped dimer a kinetically trapped entity. Hypothetically, the wild-type protein has evolved to minimize dimerization of the folding intermediate through a critical hydrogen bond (Tyr60-Glu72) that disfavors the dimeric form.

Energy landscapes of functional proteins are inherently risky

Evolutionary pressure for protein function leads to unavoidable sampling of conformational states that are at risk of misfolding and aggregation. The resulting tension between functional requirements and the risk of misfolding and/or aggregation in the evolution of proteins is becoming more and more apparent. One outcome of this tension is sensitivity to mutation, in which only subtle changes in sequence that may be functionally advantageous can tip the delicate balance toward protein aggregation. Similarly, increasing the concentration of aggregation-prone species by reducing the ability to control protein levels or compromising protein folding capacity engenders increased risk of aggregation and disease. In this Perspective, we describe examples that epitomize the tension between protein functional energy landscapes and aggregation risk. Each case illustrates how the energy landscapes for the at-risk proteins are sculpted to enable them to perform their functions and how the risks of aggregation are minimized under cellular conditions using a variety of compensatory mechanisms.

Delicate balance between functionally required flexibility and aggregation risk in a β-rich protein

Susceptibility to aggregation is general to proteins because of the potential for intermolecular interactions between hydrophobic stretches in their amino acid sequences. Protein aggregation has been implicated in several catastrophic diseases, yet we still lack in-depth understanding about how proteins are channeled to this state. Using a predominantly β-sheet protein whose folding has been explored in detail, cellular retinoic acid-binding protein 1 (CRABP1), as a model, we have tackled the challenge of understanding the links between a protein's natural tendency to fold, 'breathe', and function with its propensity to misfold and aggregate. We identified near-native dynamic species that lead to aggregation and found that inherent structural fluctuations in the native protein, resulting in opening of the ligand-entry portal, expose hydrophobic residues on the most vulnerable aggregation-prone sequences in CRABP1. CRABP1 and related intracellullar lipid-binding proteins have not been reported to aggregate inside cells, and we speculate that the cellular concentration of their open, aggregation-prone conformations is sufficient for ligand binding but below the critical concentration for aggregation. Our finding provides an example of how nature fine-tunes a delicate balance between protein function, conformational variability, and aggregation vulnerability and implies that with the evolutionary requirement for proteins to fold and function, aggregation becomes an unavoidable but controllable risk.

Early folding events protect aggregation-prone regions of a β-rich protein

Protein folding and aggregation inevitably compete with one another. This competition is even keener for proteins with frustrated landscapes, such as those rich in β structure. It is interesting that, despite their rugged energy landscapes and high β sheet content, intracellular lipid-binding proteins (iLBPs) appear to successfully avoid aggregation, as they are not implicated in aggregation diseases. In this study, we used a canonical iLBP, cellular retinoic acid-binding protein 1 (CRABP1), to understand better how folding is favored over aggregation. Analysis of folding kinetics of point mutants reveals that the folding pathway of CRABP1 involves early barrel closure. This folding mechanism protects sequences in CRABP1 that comprise cores of aggregates as identified by nuclear magnetic resonance. The amino acid conservation pattern in other iLBPs suggests that early barrel closure may be a general strategy for successful folding and minimization of aggregation. We suggest that folding mechanisms in general may incorporate steps that disfavor aggregation.

The Role of Aromatic-Aromatic Interactions in Strand-Strand Stabilization of β-Sheets

Aromatic-aromatic interactions have long been believed to play key roles in protein structure, folding, and binding functions. However, we still lack full understanding of the contributions of aromatic-aromatic interactions to protein stability and the timing of their formation during folding. Here, using an aromatic ladder in the β-barrel protein, cellular retinoic acid-binding protein 1 (CRABP1), as a case study, we find that aromatic π stacking plays a greater role in the Phe65-Phe71 cross-strand pair, while in another pair, Phe50-Phe65, hydrophobic interactions are dominant. The Phe65-Phe71 pair spans β-strands 4 and 5 in the β-barrel, which lack interstrand hydrogen bonding, and we speculate that it compensates energetically for the absence of strand-strand backbone interactions. Using perturbation analysis, we find that both aromatic-aromatic pairs form after the transition state for folding of CRABP1, thus playing a role in the final stabilization of the β-sheet rather than in its nucleation as had been earlier proposed. The aromatic interaction between strands 4 and 5 in CRABP1 is highly conserved in the intracellular lipid-binding protein (iLBP) family, and several lines of evidence combine to support a model wherein it acts to maintain barrel structure while allowing the dynamic opening that is necessary for ligand entry. Lastly, we carried out a bioinformatics analysis and found 51 examples of aromatic-aromatic interactions across non-hydrogen-bonded β-strands outside the iLBPs, arguing for the generality of the role played by this structural motif.

Tuning the electronic absorption of protein-embedded all-trans-retinal

Protein-chromophore interactions are a central component of a wide variety of critical biological processes such as color vision and photosynthesis. To understand the fundamental elements that contribute to spectral tuning of a chromophore inside the protein cavity, we redesigned human cellular retinol binding protein II (hCRBPII) to fully encapsulate all-trans-retinal and form a covalent bond as a protonated Schiff base. This system, using rational mutagenesis designed to alter the electrostatic environment within the binding pocket of the host protein, enabled regulation of the absorption maximum of the pigment in the range of 425 to 644 nanometers. With only nine point mutations, the hCRBPII mutants induced a systematic shift in the absorption profile of all-trans-retinal of more than 200 nanometers across the visible spectrum.

Residual interactions in unfolded bile acid-binding protein by 19F NMR

The folding initiation mechanism of human bile acid-binding protein (BABP) has been examined by (19) F NMR. Equilibrium unfolding studies of BABP labeled with fluorine at all eight of its phenylalanine residues showed that at least two sites experience changes in solvent exposure at high denaturant concentrations. Peak assignments were made by site-specific 4FPhe incorporation. The resonances for proteins specifically labeled at Phe17, Phe47, and Phe63 showed changes in chemical shift at denaturant concentrations at which the remaining five phenylalanine residues appear to be fully solvent-exposed. Phe17 is a helical residue that was not expected to participate in a folding initiation site. Phe47 and Phe63 form part of a hydrophobic core region that may be conserved as a site for folding initiation in the intracellular lipid-binding protein family.

Characterizing the regularity of tetrahedral packing motifs in protein tertiary structure

MOTIVATION

While protein secondary structure is well understood, representing the repetitive nature of tertiary packing in proteins remains difficult. We have developed a construct called the relative packing group (RPG) that applies the clique concept from graph theory as a natural basis for defining the packing motifs in proteins. An RPG is defined as a clique of residues, where every member contacts all others as determined by the Delaunay tessellation. Geometrically similar RPGs define a regular element of tertiary structure or tertiary motif (TerMo). This intuitive construct provides a simple approach to characterize general repetitive elements of tertiary structure.

RESULTS

A dataset of over 4 million tetrahedral RPGs was clustered using different criteria to characterize the various aspects of regular tertiary structure in TerMos. Grouping this data within the SCOP classification levels of Family, Superfamily, Fold, Class and PDB showed that similar packing is shared across different folds. Classification of RPGs based on residue sequence locality reveals topological preferences according to protein sizes and secondary structure. We find that larger proteins favor RPGs with three local residues packed against a non-local residue. Classifying by secondary structure, helices prefer mostly local residues, sheets favor at least two local residues, while turns and coil populate with more local residues. To depict these TerMos, we have developed 2 complementary and intuitive representations: (i) Dirichlet process mixture density estimation of the torsion angle distributions and (ii) kernel density estimation of the Cartesian coordinate distribution. The TerMo library and representations software are available upon request.

Dissection of the critical binding determinants of cellular retinoic acid binding protein II by mutagenesis and fluorescence binding assay

The binding of retinoic acid to mutants of Cellular Retinoic Acid Binding Protein II (CRABPII) was evaluated to better understand the importance of the direct protein/ligand interactions. The important role of Arg111 for the correct structure and function of the protein was verified and other residues that directly affect retinoic acid binding have been identified. Furthermore, retinoic acid binding to CRABPII mutants that lack all previously identified interacting amino acids was rescued by providing a carboxylic acid dimer partner in the form of a Glu residue.

Chapter 3: A fluorescent window into protein folding and aggregation in cells

Evolutionary selective pressures have tuned the efficiency of the protein-folding reaction in the crowded complex environment in the cell. Nevertheless, the fidelity of folding is imperfect, leading to off-pathway intermolecular interactions that compete with proper folding and to consequent formation of thermodynamically stable aggregates. Such aggregates constitute the histopathological hallmarks of many neurodegenerative pathologies. Yet, most of the approaches to characterize protein folding and/or misfolding are limited to in vitro conditions. Here, we describe a strategy to directly monitor the behavior of a protein in prokaryotic and eukaryotic cells. The method is based on incorporation of structurally non-perturbing, specific binding motifs for a bis-arsenical fluoroscein dye, FlAsH, in sites that result in distinct dye fluorescence signals for the folded and unfolded states of the protein under study. Our approach has been developed using as a case study the predominantly beta-sheet intracellular lipid-binding protein, cellular retinoic acid-binding protein, alone or as a chimera fused to the exon 1-encoded fragment of huntingtin, which harbors a polyglutamine repeat tract. We have designed protocols to label this protein in vivo and to monitor the resulting fluorescence signal, which reports on any misfolding transition and formation of aggregates, yielding quantitatively interpretable data.

Roles of beta-turns in protein folding: from peptide models to protein engineering

Reverse turns are a major class of protein secondary structure; they represent sites of chain reversal and thus sites where the globular character of a protein is created. It has been speculated for many years that turns may nucleate the formation of structure in protein folding, as their propensity to occur will favor the approximation of their flanking regions and their general tendency to be hydrophilic will favor their disposition at the solvent-accessible surface. Reverse turns are local features, and it is therefore not surprising that their structural properties have been extensively studied using peptide models. In this article, we review research on peptide models of turns to test the hypothesis that the propensities of turns to form in short peptides will relate to the roles of corresponding sequences in protein folding. Turns with significant stability as isolated entities should actively promote the folding of a protein, and by contrast, turn sequences that merely allow the chain to adopt conformations required for chain reversal are predicted to be passive in the folding mechanism. We discuss results of protein engineering studies of the roles of turn residues in folding mechanisms. Factors that correlate with the importance of turns in folding indeed include their intrinsic stability, as well as their topological context and their participation in hydrophobic networks within the protein's structure.

Functional drift of sequence attributes in the FK506-binding proteins (FKBPs)

Diverse members of the FK506-binding proteins (FKBPs) group and their complexes with different macrocyclic ligands of fungal origins such as FK506, rapamycin, ascomycin, and their immunosuppressive and nonimmunosuppressive derivatives display a variety of cellular and biological activities. The functional relatedness of the FKBPs was estimated from the following attributes of their aligned sequences: 1 degrees conservation of the consensus sequence; 2 degrees sequence similarity; 3 degrees pI; 4 degrees hydrophobicity; 5 degrees amino acid hydrophobicity and bulkiness profiles. Analyses of the multiple sequence alignments and intramolecular interaction networks calculated from a series of structures of the FKBPs revealed some variations in the interaction clusters formed by the AA residues that are crucial for sustaining peptidylprolyl cis/trans isomerases (PPIases) activity and binding capacity of the FKBPs. Fine diversification of the sequences of the multiple paralogues and orthologues of the FKBPs encoded in different genomes alter the intramolecular interaction patterns of their structures and allowed them to gain some selectivity in binding to diverse targets (functional drift).

"Hot cores" in proteins: comparative analysis of the apolar contact area in structures from hyper/thermophilic and mesophilic organisms

Background

A wide variety of stabilizing factors have been invoked so far to elucidate the structural basis of protein thermostability. These include, amongst the others, a higher number of ion-pairs interactions and hydrogen bonds, together with a better packing of hydrophobic residues. It has been frequently observed that packing of hydrophobic side chains is improved in hyperthermophilic proteins, when compared to their mesophilic counterparts. In this work, protein crystal structures from hyper/thermophilic organisms and their mesophilic homologs have been compared, in order to quantify the difference of apolar contact area and to assess the role played by the hydrophobic contacts in the stabilization of the protein core, at high temperatures.

Results

The construction of two datasets was carried out so as to satisfy several restrictive criteria, such as minimum redundancy, resolution and R-value thresholds and lack of any structural defect in the collected structures. This approach allowed to quantify with relatively high precision the apolar contact area between interacting residues, reducing the uncertainty due to the position of atoms in the crystal structures, the redundancy of data and the size of the dataset. To identify the common core regions of these proteins, the study was focused on segments that conserve a similar main chain conformation in the structures analyzed, excluding the intervening regions whose structure differs markedly. The results indicated that hyperthermophilic proteins underwent a significant increase of the hydrophobic contact area contributed by those residues composing the alpha-helices of the structurally conserved regions.

Conclusion

This study indicates the decreased flexibility of alpha-helices in proteins core as a major factor contributing to the enhanced termostability of a number of hyperthermophilic proteins. This effect, in turn, may be due to an increased number of buried methyl groups in the protein core and/or a better packing of alpha-helices with the rest of the structure, caused by the presence of hydrophobic beta-branched side chains.

Distantly related lipocalins share two conserved clusters of hydrophobic residues: use in homology modeling

Background

Lipocalins are widely distributed in nature and are found in bacteria, plants, arthropoda and vertebra. In hematophagous arthropods, they are implicated in the successful accomplishment of the blood meal, interfering with platelet aggregation, blood coagulation and inflammation and in the transmission of disease parasites such as Trypanosoma cruzi and Borrelia burgdorferi.

The pairwise sequence identity is low among this family, often below 30%, despite a well conserved tertiary structure. Under the 30% identity threshold, alignment methods do not correctly assign and align proteins. The only safe way to assign a sequence to that family is by experimental determination. However, these procedures are long and costly and cannot always be applied. A way to circumvent the experimental approach is sequence and structure analyze. To further help in that task, the residues implicated in the stabilisation of the lipocalin fold were determined. This was done by analyzing the conserved interactions for ten lipocalins having a maximum pairwise identity of 28% and various functions.

Results

It was determined that two hydrophobic clusters of residues are conserved by analysing the ten lipocalin structures and sequences. One cluster is internal to the barrel, involving all strands and the 3₁₀ helix. The other is external, involving four strands and the helix lying parallel to the barrel surface. These clusters are also present in RaHBP2, a unusual "outlier" lipocalin from tick Rhipicephalus appendiculatus. This information was used to assess assignment of LIR2 a protein from Ixodes ricinus and to build a 3D model that helps to predict function. FTIR data support the lipocalin fold for this protein.

Conclusion

By sequence and structural analyzes, two conserved clusters of hydrophobic residues in interactions have been identified in lipocalins. Since the residues implicated are not conserved for function, they should provide the minimal subset necessary to confer the lipocalin fold. This information has been used to assign LIR2 to lipocalins and to investigate its structure/function relationship. This study could be applied to other protein families with low pairwise similarity, such as the structurally related fatty acid binding proteins or avidins.

Folding and stability of the isolated Greek key domains of the long-lived human lens proteins gammaD-crystallin and gammaS-crystallin

The transparency of the eye lens depends on the high solubility and stability of the lens crystallin proteins. The monomeric gamma-crystallins and oligomeric beta-crystallins have paired homologous double Greek key domains, presumably evolved through gene duplication and fusion. Prior investigation of the refolding of human gammaD-crystallin revealed that the C-terminal domain folds first and nucleates the folding of the N-terminal domain. This result suggested that the human N-terminal domain might not be able to fold on its own. We constructed and expressed polypeptide chains corresponding to the isolated N- and C-terminal domains of human gammaD-crystallin, as well as the isolated domains of human gammaS-crystallin. Both circular dichroism and fluorescence spectroscopy indicated that the isolated domains purified from Escherichia coli were folded into native-like monomers. After denaturation, the isolated domains refolded efficiently at pH 7 and 37 degrees C into native-like structures. The in vitro refolding of all four domains revealed two kinetic phases, identifying partially folded intermediates for the Greek key motifs. When subjected to thermal denaturation, the isolated N-terminal domains were less stable than the full-length proteins and less stable than the C-terminal domains, and this was confirmed in equilibrium unfolding/refolding experiments. The decrease in stability of the N-terminal domain of human gammaD-crystallin with respect to the complete protein indicated that the interdomain interface contributes of 4.2 kcal/mol to the overall stability of this very long-lived protein.

Evolutionary coupling of structural and functional sequence information in the intracellular lipid-binding protein family

We have mined the evolutionary record for the large family of intracellular lipid-binding proteins (iLBPs) by calculating the statistical coupling of residue variations in a multiple sequence alignment using methods developed by Ranganathan and coworkers (Lockless and Ranganathan, Science 1999:286;295-299). The 213 sequences analyzed have a wide range of ligand-binding functions as well as highly divergent phylogenetic origins, assuring broad sampling of sequence space. Emerging from this analysis were two major clusters of coupled residues, which when mapped onto the structure of a representative iLBP under study in our laboratory, cellular retinoic-acid binding protein I, are largely contiguous and provide useful points of comparison to available data for the folding of this protein. One cluster comprises a predominantly hydrophobic core away from the ligand-binding site and likely represents key structural information for the iLBP fold. The other cluster includes the portal region where ligand enters its binding site, regions of the ligand-binding cavity, and the region where the 10-stranded beta-barrel characteristic of this family closes (between strands 1' and 10). Linkages between these two clusters suggest that evolutionary pressures on this family constrain structural and functional sequence information in an interdependent fashion. The necessity of the structure to wrap around a hydrophobic ligand confounds the typical sequestration of hydrophobic side chains. Additionally, ligand entry and exit require these structures to have a capacity for specific conformational change during binding and release. We conclude that an essential and structurally apparent separation of local and global sequence information is conserved throughout the iLBP family.

The Structure of Apo-wild-type Cellular Retinoic Acid Binding Protein II at 1.4 Å and its Relationship to Ligand Binding and Nuclear Translocation

CRABPII is a small, cytosolic protein that solubilizes and transfers retinoic acid (RA) to the nucleus while also enhancing its transcriptional activity. We have determined the first high-resolution structure of apo-wild type (WT) CRABPII at 1.35 A. Using three different data sets collected on apo-WT CRABPII we have shown that apo- and holo-CRABPII share very similar structures. Binding of RA appears to increase the overall rigidity of the structure, although the induced structural changes are not as pronounced as previously thought. The enhanced structural rigidity may be an important determinant for the enhanced nuclear localization of the RA-bound protein. Comparison of our apo-WT with a mutant apo-CRABPII structure shows that mutation of Arg111, a conserved residue of CRABPII and a key residue in RA binding, causes structural changes in the molecule. We further investigated the structural importance of conserved residues by determining the structure of the F15W mutant CRABPII (F15W-CRABPII). Our structures also demonstrate structural changes induced by crystal packing and show that a crystal can harbor demonstrative structural differences in the asymmetric unit.

Determinants of protein stability and folding: comparative analysis of beta-lactoglobulins and liver basic fatty acid binding protein

A new energy decomposition approach, aimed at identifying residues playing a folding key role, has been applied here to three homologous proteins, belonging to the calycin superfamily, namely bovine and porcine beta-lactoglobulins and Liver basic fatty acid binding protein, sharing the same beta-barrel fold and different degree of sequence identities. All-atom, explicit solvent molecular dynamics simulations around the native conformation were used to generate, for each of the three proteins, energy maps which were further simplified through eigenvalue decomposition. Analysis of the components of the eigenvector associated with the lowest eigenvalue singled out those residues (hot sites) behaving as strongly interacting and possible nucleation centers. The results fit well with experimental folding data and, especially, with the analysis of side chain-side chain interaction conservation.

Extended polyglutamine tracts cause aggregation and structural perturbation of an adjacent beta barrel protein

Formation of fibrillar intranuclear inclusions and related neuropathologies of the CAG-repeat disorders are linked to the expansion of a polyglutamine tract. Despite considerable effort, the etiology of these devastating diseases remains unclear. Although polypeptides with glutamine tracts recapitulate many of the observed characteristics of the gene products with CAG repeats, such as in vitro and in vivo aggregation and toxicity in model organisms, extended polyglutamine segments have also been reported to structurally perturb proteins into which they are inserted. Additionally, the sequence context of a polyglutamine tract has recently been shown to modulate its propensity to aggregate. These findings raise the possibility that indirect influences of the repeat tract on adjacent protein domains are contributory to pathologies. Destabilization of an adjacent domain may lead to loss of function, as well as favoring non-native structures in the neighboring domain causing them to be prone to intermolecular association and consequent aggregation. To explore these phenomena, we have used chimeras of a well studied globular protein and exon 1 of huntingtin. We find that expansion of the polyglutamine segment beyond the pathological threshold (>35 glutamines) results in structural perturbation of the neighboring protein whether the huntingtin exon is N- or C-terminal. Elongation of the polyglutamine region also substantially increases the propensity of the chimera to aggregate, both in vitro and in vivo, and in vitro aggregation kinetics of a chimera with a 53-glutamine repeat follow a nucleation polymerization mechanism with a monomeric nucleus.

Characterization of the molten globule state of retinol-binding protein using a molecular dynamics simulation approach

Retinol-binding protein transports retinol, and circulates in the plasma as a macromolecular complex with the protein transthyretin. Under acidic conditions retinol-binding protein undergoes a transition to the molten globule state, and releases the bound retinol ligand. A biased molecular dynamics simulation method has been used to generate models for the ensemble of conformers populated within this molten globule state. Simulation conformers, with a radius of gyration at least 1.1 A greater than that of the native state, contain on average 37%beta-sheet secondary structure. In these conformers the central regions of the two orthogonal beta-sheets that make up the beta-barrel in the native protein are highly persistent. However, there are sizable fluctuations for residues in the outer regions of the beta-sheets, and large variations in side chain packing even in the protein core. Significant conformational changes are seen in the simulation conformers for residues 85-104 (beta-strands E and F and the E-F loop). These changes give an opening of the retinol-binding site. Comparisons with experimental data suggest that the unfolding in this region may provide a mechanism by which the complex of retinol-binding protein and transthyretin dissociates, and retinol is released at the cell surface.

Porcine beta-lactoglobulin chemical unfolding: Identification of a non-native α-helical intermediate

The chemical unfolding behavior of porcine beta-lactoglobulin (PLG) has been followed at pH 2 and 6 in the presence of guanidinium hydrochloride. The PLG unfolding transition, monitored by tryptophan fluorescence, far and near UV circular dichroism and 1D-NMR, can be described by a three-state transition suggesting the presence of at least one intermediate state that appears to display an excess of non-native alpha-helical structures. The thermodynamic parameters, as determined through a global analysis fitting procedure, give estimates of the free energy differences of the transitions connecting the native, the intermediate and the unfolded state: DeltaG(NI) (0) = 2.8 +/- 0.7 kcal mol(-1) (pH 2) and 4.2 +/- 0.5 kcal mol(-1) (pH 6) and DeltaG(NU) (0) = 7.2 +/- 0.6 kcal mol(-1) (pH 2) and 6.9 +/- 0.6 kcal mol(-1) (pH 6). CD unfolding data of the bovine species (BLG) have been collected here under the same experimental conditions of PLG to allow a careful comparison of the two beta-lactoglobulins. Intermediates with different characteristics have been identified for BLG and PLG, and their nature has been discussed on a structural analysis basis. The thermodynamic data reported here for PLG and BLG and the comparative analysis with data reported for equine beta lactoglobulin, show that homologous beta-barrel proteins, belonging to the same family and displaying high sequence identity (52-64%) populate unfolding intermediates to different extents, even though a common tendency to the formation of non-native alpha-helical intermediates, can be envisaged. The present results provide a prerequisite foundation of knowledge for the design and interpretation of future folding kinetic studies.

Evolutionarily conserved regions and hydrophobic contacts at the superfamily level: The case of the fold-type I, pyridoxal-5'-phosphate-dependent enzymes

The wealth of biological information provided by structural and genomic projects opens new prospects of understanding life and evolution at the molecular level. In this work, it is shown how computational approaches can be exploited to pinpoint protein structural features that remain invariant upon long evolutionary periods in the fold-type I, PLP-dependent enzymes. A nonredundant set of 23 superposed crystallographic structures belonging to this superfamily was built. Members of this family typically display high-structural conservation despite low-sequence identity. For each structure, a multiple-sequence alignment of orthologous sequences was obtained, and the 23 alignments were merged using the structural information to obtain a comprehensive multiple alignment of 921 sequences of fold-type I enzymes. The structurally conserved regions (SCRs), the evolutionarily conserved residues, and the conserved hydrophobic contacts (CHCs) were extracted from this data set, using both sequence and structural information. The results of this study identified a structural pattern of hydrophobic contacts shared by all of the superfamily members of fold-type I enzymes and involved in native interactions. This profile highlights the presence of a nucleus for this fold, in which residues participating in the most conserved native interactions exhibit preferential evolutionary conservation, that correlates significantly (r = 0.70) with the extent of mean hydrophobic contact value of their apolar fraction.

Recognition of functional sites in protein structures

Recognition of regions on the surface of one protein, that are similar to a binding site of another is crucial for the prediction of molecular interactions and for functional classifications. We first describe a novel method, SiteEngine, that assumes no sequence or fold similarities and is able to recognize proteins that have similar binding sites and may perform similar functions. We achieve high efficiency and speed by introducing a low-resolution surface representation via chemically important surface points, by hashing triangles of physico-chemical properties and by application of hierarchical scoring schemes for a thorough exploration of global and local similarities. We proceed to rigorously apply this method to functional site recognition in three possible ways: first, we search a given functional site on a large set of complete protein structures. Second, a potential functional site on a protein of interest is compared with known binding sites, to recognize similar features. Third, a complete protein structure is searched for the presence of an a priori unknown functional site, similar to known sites. Our method is robust and efficient enough to allow computationally demanding applications such as the first and the third. From the biological standpoint, the first application may identify secondary binding sites of drugs that may lead to side-effects. The third application finds new potential sites on the protein that may provide targets for drug design. Each of the three applications may aid in assigning a function and in classification of binding patterns. We highlight the advantages and disadvantages of each type of search, provide examples of large-scale searches of the entire Protein Data Base and make functional predictions.

Monitoring protein stability and aggregation in vivo by real-time fluorescent labeling

In vivo fluorescent labeling of an expressed protein has enabled the observation of its stability and aggregation directly in bacterial cells. Mammalian cellular retinoic acid-binding protein I (CRABP I) was mutated to incorporate in a surface-exposed omega loop the sequence Cys-Cys-Gly-Pro-Cys-Cys, which binds specifically to a biarsenical fluorescein dye (FlAsH). Unfolding of labeled tetra-Cys CRABP I is accompanied by enhancement of FlAsH fluorescence, which made it possible to determine the free energy of unfolding of this protein by urea titration in cells and to follow in real time the formation of inclusion bodies by a slow-folding, aggregationprone mutant (FlAsH-labeled P39A tetra-Cys CRABP I). Aggregation in vivo displayed a concentration-dependent apparent lag time similar to observations of protein aggregation in purified in vitro model systems.