Role of local sequence in the folding of cellular retinoic abinding protein I: structural propensities of reverse turns

The experiments described here explore the role of local sequence in the folding of cellular retinoic acid binding protein I (CRABP I). This is a 136-residue, 10-stranded, antiparallel beta-barrel protein with seven beta-hairpins and is a member of the intracellular lipid binding protein (iLBP) family. The relative roles of local and global sequence information in governing the folding of this class of proteins are not well-understood. In question is whether the beta-turns are locally defined by short-range interactions within their sequences, and are thus able to play an active role in reducing the conformational space available to the folding chain, or whether the turns are passive, relying upon global forces to form. Short (six- and seven-residue) peptides corresponding to the seven CRABP I turns were analyzed by circular dichroism and NMR for their tendencies to take up the conformations they adopt in the context of the native protein. The results indicate that two of the peptides, encompassing turns III and IV in CRABP I, have a strong intrinsic bias to form native turns. Intriguingly, these turns are on linked hairpins in CRABP I and represent the best-conserved turns in the iLBP family. These results suggest that local sequence may play an important role in narrowing the conformational ensemble of CRABP I during folding.     

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.     

Protein-DNA interactions: amino acid conservation and the effects of mutations on binding specificity

We investigate the conservation of amino acid residue sequences in 21 DNA-binding protein families and study the effects that mutations have on DNA-sequence recognition. The observations are best understood by assigning each protein family to one of three classes: (i) non-specific, where binding is independent of DNA sequence; (ii) highly specific, where binding is specific and all members of the family target the same DNA sequence; and (iii) multi-specific, where binding is also specific, but individual family members target different DNA sequences. Overall, protein residues in contact with the DNA are better conserved than the rest of the protein surface, but there is a complex underlying trend of conservation for individual residue positions. Amino acid residues that interact with the DNA backbone are well conserved across all protein families and provide a core of stabilising contacts for homologous protein-DNA complexes. In contrast, amino acid residues that interact with DNA bases have variable levels of conservation depending on the family classification. In non-specific families, base-contacting residues are well conserved and interactions are always found in the minor groove where there is little discrimination between base types. In highly specific families, base-contacting residues are highly conserved and allow member proteins to recognise the same target sequence. In multi-specific families, base-contacting residues undergo frequent mutations and enable different proteins to recognise distinct target sequences. Finally, we report that interactions with bases in the target sequence often follow (though not always) a universal code of amino acid-base recognition and the effects of amino acid mutations can be most easily understood for these interactions.     

Folding of intracellular retinol and retinoic acid binding proteins

The folding mechanisms of cellular retinol binding protein II (CRBP II), cellular retinoic acid binding protein I (CRABP I), and cellular retinoic acid binding protein II (CRABP II) were examined. These beta-sheet proteins have very similar structures and higher sequence homologies than most proteins in this diverse family. They have similar stabilities and show completely reversible folding at equilibrium with urea as a denaturant. The unfolding kinetics of these proteins were monitored during folding and unfolding by circular dichroism (CD) and fluorescence. During unfolding, CRABP II showed no intermediates, CRABP I had an intermediate with nativelike secondary structure, and CRBP II had an intermediate that lacked secondary structure. The refolding kinetics of these proteins were more similar. Each protein showed a burst-phase change in intensity by both CD and fluorescence, followed by a single observed phase by both CD and fluorescence and one or two additional refolding phases by fluorescence. The fluorescence spectral properties of the intermediate states were similar and suggested a gradual increase in the amount of native tertiary structure present for each step in a sequential path. However, the rates of folding differed by as much as 3 orders of magnitude and were slower than those expected from the contact order and topology of these proteins. As such, proteins with the same final structure may not follow the same route to the native state.     

Native structural propensity in cellular retinoic acid-binding protein I 64-88: the role of locally encoded structure in the folding of a beta-barrel protein

A central question in protein folding is the relative importance of locally encoded structure and cooperative interactions among residues distant in sequence. We have been exploring this question in a predominantly β-sheet protein, since β-structure formation clearly relies on both local and global sequence information. We present evidence that a 24-residue peptide corresponding to two linked hairpins of cellular retinoic acid-binding protein I (CRABP I) adopts significant native structure in aqueous solution. Prior work from our laboratory showed that the two turns contained in this fragment (turns III and IV) had the highest tendency of any of the eight turns in this anti-parallel β-barrel to fold into native turns. In addition, the primary sequence of these two turns is well conserved throughout the structural family to which CRABP I belongs, and residues in the turns and their associated hairpins participate in a network of conserved long-range interactions. We propose that the strong local-sequence biases within the chain segment comprising turns III and IV favor longer-range interactions that are crucial to the folding and native-state stability of CRABP I, and may play a similar role in related intracellular lipid-binding proteins (iLBPs).     

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.     

Distinguishing structural and functional restraints in evolution in order to identify interaction sites

Structural genomics projects are producing many three-dimensional structures of proteins that have been identified only from their gene sequences. It is therefore important to develop computational methods that will predict sites involved in productive intermolecular interactions that might give clues about functions. Techniques based on evolutionary conservation of amino acids have the advantage over physiochemical methods in that they are more general. However, the majority of techniques neither use all available structural and sequence information, nor are able to distinguish between evolutionary restraints that arise from the need to maintain structure and those that arise from function. Three methods to identify evolutionary restraints on protein sequence and structure are described here. The first identifies those residues that have a higher degree of conservation than expected: this is achieved by comparing for each amino acid position the sequence conservation observed in the homologous family of proteins with the degree of conservation predicted on the basis of amino acid type and local environment. The second uses information theory to identify those positions where environment-specific substitution tables make poor predictions of the overall amino acid substitution pattern. The third method identifies those residues that have highly conserved positions when three-dimensional structures of proteins in a homologous family are superposed. The scores derived from these methods are mapped onto the protein three-dimensional structures and contoured, allowing identification clusters of residues with strong evolutionary restraints that are sites of interaction in proteins involved in a variety of functions. Our method differs from other published techniques by making use of structural information to identify restraints that arise from the structure of the protein and differentiating these restraints from others that derive from intermolecular interactions that mediate functions in the whole organism.     

Localization of protein-binding sites within families of proteins

We address the question of whether or not the positions of protein-binding sites on homologous protein structures are conserved irrespective of the identities of their binding partners. First, for each domain family in the Structural Classification of Proteins (SCOP), protein-binding sites are extracted from our comprehensive database of structurally defined binary domain interactions (PIBASE). Second, the binding sites within each family are superposed using a structural alignment of its members. Finally, the degree of localization of binding sites within each family is quantified by comparing it with localization expected by chance. We found that 72% of the 1847 SCOP domain families in PIBASE have binding sites with localization values greater than expected by chance. Moreover, 554 (30%) of these families have localizations that are statistically significant (i.e., more than four standard deviations away from the mean expected by chance). In contrast, only 144 (8%) families have significantly low localization. The absence of a significant correlation of the binding site localization with the average sequence and structural conservations in a family suggests that localization can be helpful for describing the functional diversity of protein-protein interactions, complementing measures of sequence and structural conservation. Consideration of the binding site localization may also result in spatial restraints for the modeling of protein assembly structures.     

Folding requirements are different between sterol 14alpha-demethylase (CYP51) from Mycobacterium tuberculosis and human or fungal orthologs

Upon sequence alignment of CYP51 sterol 14alpha-demethylase from animals, plants, fungi, and bacteria, arginine corresponding to Arg-448 of CYP51 in Mycobacterium tuberculosis (MT) is conserved near the C terminus of all family members. In MTCYP51 Arg-448 forms a salt bridge with Asp-287, connecting beta-strand 3-2 with helix J. Deletion of the three C-terminal residues of MTCYP51 has little effect on expression of P450 in Escherichia coli. However, truncation of the fourth amino acid (Arg-448) completely abolishes P450 expression. We have investigated whether Arg-448 has other structural or functional roles in addition to folding and whether its conservation reflects conservation of a common folding pathway in the CYP51 family. Characterization of wild type protein and three mutants, R448K, R448I, and R448A, including examination of catalytic activity, secondary and tertiary structure analysis by circular dichroism and tryptophan fluorescence, and studies of both equilibrium and temporal MTCYP51 unfolding behavior, shows that Arg-448 does not play any role in P450 function or maintenance of the native structure. C-terminal truncation of Candida albicans and human CYP51 orthologs reveals that, despite conservation in sequence, the requirement for arginine at the homologous C-terminal position in folding in E. coli is not conserved. Thus, despite similar spatial folds, functionally related but evolutionarily distinct P450s can follow different folding pathways.     

The relationship between conservation, thermodynamic stability, and function in the SH3 domain hydrophobic core

To investigate the relationships between sequence conservation, protein stability, and protein function, we have measured the thermodynamic stability, folding kinetics, and in vitro peptide-binding activity of a large number of single-site substitutions in the hydrophobic core of the Fyn SH3 domain. Comparison of these data to that derived from an analysis of a large alignment of SH3 domain sequences revealed a very good correlation between the distinct pattern of conservation observed at each core position and the thermodynamic stability of mutants. Conservation was also found to correlate well with the unfolding rates of mutants, but not to the folding rates, suggesting that evolution selects more strongly for optimal native state packing interactions than for maximal folding rates. Structural analysis suggests that residue-residue core packing interactions are very similar in all SH3 domains, which provides an explanation for the correlation between conservation and mutant stability effects studied in a single SH3 domain. We also demonstrate a correlation between stability and the in vivo activity of mutants, and between conservation and activity. However, the relationship between conservation and activity was very strong only for the three most conserved hydrophobic core positions. The weaker correlation between activity and conservation seen at the other seven core positions indicates that maintenance of protein stability is the dominant selective pressure at these positions. In general, the pattern of conservation at hydrophobic core positions appears to arise from conserved packing constraints, and can be effectively utilized to predict the destabilizing effects of amino acid substitutions.     

Fluorescent dyes as probes to study lipid-binding proteins

We studied the equilibrium binding of two hydrophobic fluorescent dyes, ANS and bisANS, to four members of a family of intracellular lipid-binding proteins: IFABP, CRABP I, CRABP II, and ILBP. The spectral and binding parameters for the probes bound to the proteins were determined. Typically, there was a single binding site on each protein for the ligands. However, IFABP cooperatively bound a second bisANS molecule in the binding pocket. Comparative analysis of affinities and spectral characteristics for the two probes allowed us to examine the contributions of electrostatic and hydrophobic interactions to the binding process, and to address some aspects of the internal structure of the studied proteins.     

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.     

Expression and functional analysis of the cellular retinoic acid binding protein from silkworm pupae (Bombyx mori)

Cellular retinoic acid binding protein (CRABP) is a member of intracellular lipid-binding protein (iLBP), and closely associated with retinoic acid (RA) activity. We have cloned the CRABP gene from silkworm pupae and studied the interaction between Bombyx mori CRABP (BmCRABP) and all-trans retinoic acid (atRA). The MTT assay data indicated that when BmCRABP is overexpressed in Bm5 cells, the cells dramatically resisted to atRA-induced growth inhibition. Conversely, the cells were sensitive to atRA treatment upon knocking down the BmCRABP expression. Subcellular localization revealed that BmCRABP is a cytoplasm protein, even when treated with atRA, the CRABP still remained in the cytoplasm. These data demonstrated that the function of BmCRABP have an effect on the physiological function of atRA.     

Fold change in evolution of protein structures

Typically, protein spatial structures are more conserved in evolution than amino acid sequences. However, the recent explosion of sequence and structure information accompanied by the development of powerful computational methods led to the accumulation of examples of homologous proteins with globally distinct structures. Significant sequence conservation, local structural resemblance, and functional similarity strongly indicate evolutionary relationships between these proteins despite pronounced structural differences at the fold level. Several mechanisms such as insertions/deletions/substitutions, circular permutations, and rearrangements in beta-sheet topologies account for the majority of detected structural irregularities. The existence of evolutionarily related proteins that possess different folds brings new challenges to the homology modeling techniques and the structure classification strategies and offers new opportunities for protein design in experimental studies.     

Ligand-binding specificity of an invertebrate (Manduca sexta) putative cellular retinoic acid binding protein

Intracellular lipid-binding proteins (iLBPs) are small cytoplasmic proteins that specifically interact with hydrophobic ligands. Fatty acid-binding proteins (FABPs), cellular retinoic acid-binding proteins (CRABPs) and cellular retinol-binding proteins (CRBPs) belong to the iLBP family. A recently identified insect (Manduca sexta) iLBP has been reported to possibly represent an invertebrate CRABP mimicking the role of CRABPs in vertebrate organisms. The presence in this protein of the characteristic binding triad residues involved in the interaction with ligand carboxylate head groups, a feature pertaining to several FABPs and to CRABPs, and the close phylogenetic relationships with both groups of vertebrate heart-type FABPs and CRBPs/CRABPs, makes it difficult to assign it to either FABPs or CRABPs. However, its negligible interaction with retinoic acid and high affinity (K(d) values in the 10(-8) M range) for fatty acids have been established by means of direct and competitive binding assays. As shown by phylogenetic analysis, the M. sexta iLBP belongs to a wide group of invertebrate iLBPs, which, besides being closely related phylogenetically, share distinctive features, such as the conservation of chemically distinct residues in their amino acid sequences and the ability to bind fatty acids. Our results are in keeping with the lack of cellular retinoid-binding proteins in invertebrates and with their later appearance during the course of chordate evolution.     

Evolution of the family of intracellular lipid binding proteins in vertebrates

Members of the family of intracellular lipid binding proteins (iLBPs) have been implicated in cytoplasmic transport of lipophilic ligands, such as long-chain fatty acids and retinoids. iLBPs are low molecular mass proteins (14–16 kDa) sharing a common structural fold. The iLBP family likely arose through duplication and diversification of an ancestral iLBP gene. Phylogenetic analysis undertaken in the present study indicates that the ancestral iLBP gene arose after divergence of animals from fungi and plants. The first gene duplication was dated around 930 millions of years ago, and subsequent duplications in the succeeding 550 millions of years gave rise to the 16 iLBP types currently recognized in vertebrates. Four clusters of proteins, each binding a characteristic range of ligands, are evident from the phylogenetic tree. Evolution of different binding properties probably allowed cytoplasmic trafficking of distinct ligands. It is speculated that recruitment of an iLBP during evolution of animals enabled the mitochondrial oxidation of long-chain fatty acids.     

Context-Dependent Remodeling of Structure in Two Large Protein Fragments

Protein folding involves the formation of secondary structural elements from the primary sequence and their association with tertiary assemblies. The relation of this primary sequence to a specific folded protein structure remains a central question in structural biology. An increasing body of evidence suggests that variations in homologous sequence ranging from point mutations to substantial insertions or deletions can yield stable proteins with markedly different folds. Here we report the structural characterization of domain IV (D4) and ΔD4 (polypeptides with 222 and 160 amino acids, respectively) that differ by virtue of an N-terminal deletion of 62 amino acids (28% of the overall D4 sequence). The high-resolution crystal structures of the monomeric D4 and the dimeric ΔD4 reveal substantially different folds despite an overall conservation of secondary structure. These structures show that the formation of tertiary structures, even in extended polypeptide sequences, can be highly context dependent, and they serve as a model for structural plasticity in protein isoforms.     

Molecular evolution of class A β-lactamases: phylogeny and patterns of sequence conservation

We present a multiple alignment of the amino acid sequences of eight class A beta-lactamases and utilized it to propose a phylogeny, based on the nucleotide sequences of their corresponding genes. We have also used the alignment, together with the alpha-carbon co-ordinates of the Staphylococcus aureus protein, to search systematically for neighbouring residues that share the same pattern of conservation among the different members of the protein family. The distribution of invariant residues and of groups of residues with co-ordinate changes map, predominantly, at the region of the active site and at interfaces between structural elements, respectively. We have also contrasted the distribution of conserved residues with the positions which are known to differ in mutants and variants of class A beta-lactamases.