Mutating the charged residues in the binding pocket of cellular retinoic acid-binding protein simultaneously reduces its binding affinity to retinoic acid and increases its thermostability

A mutational atlas for Parkin proteostasis

Proteostasis can be disturbed by mutations affecting folding and stability of the encoded protein. An example is the ubiquitin ligase Parkin, where gene variants result in autosomal recessive Parkinsonism. To uncover the pathological mechanism and provide comprehensive genotype-phenotype information, variant abundance by massively parallel sequencing (VAMP-seq) is leveraged to quantify the abundance of Parkin variants in cultured human cells. The resulting mutational map, covering 9219 out of the 9300 possible single-site amino acid substitutions and nonsense Parkin variants, shows that most low abundance variants are proteasome targets and are located within the structured domains of the protein. Half of the known disease-linked variants are found at low abundance. Systematic mapping of degradation signals (degrons) reveals an exposed degron region proximal to the so-called "activation element". This work provides examples of how missense variants may cause degradation either via destabilization of the native protein, or by introducing local signals for degradation.

Common mechanism of thermostability in small α- and β-proteins studied by molecular dynamics

Protein thermostability is important to evolution, diseases, and industrial applications. Proteins use diverse molecular strategies to achieve stability at high temperature, yet reducing the entropy of unfolding seems required. We investigated five small α-proteins and five β-proteins with known, distinct structures and thermostability (T_m ) using multi-seed molecular dynamics simulations at 300, 350, and 400 K. The proteins displayed diverse changes in hydrogen bonding, solvent exposure, and secondary structure with no simple relationship to T_m . Our dynamics were in good agreement with experimental B-factors at 300 K and insensitive to force-field choice. Despite the very distinct structures, the native-state (300 + 350 K) free-energy landscapes (FELs) were significantly broader for the two most thermostable proteins and smallest for the three least stable proteins in both the α- and β-group and with both force fields studied independently (tailed t-test, 95% confidence level). Our results suggest that entropic ensembles stabilize proteins at high temperature due to reduced entropy of unfolding, viz., ΔG = ΔH - TΔS. Supporting this mechanism, the most thermostable proteins were also the least kinetically stable, consistent with broader FELs, typified by villin headpiece and confirmed by specific comparison to a mesophilic ortholog of Thermus thermophilus apo-pyrophosphate phosphohydrolase. We propose that molecular strategies of protein thermostabilization, although diverse, tend to converge toward highest possible entropy in the native state consistent with the functional requirements. We speculate that this tendency may explain why many proteins are not optimally structured and why molten-globule states resemble native proteins so much.

Resolving the dark matter of ABCA4 for 1054 Stargardt disease probands through integrated genomics and transcriptomics

PURPOSE

Missing heritability in human diseases represents a major challenge, and this is particularly true for ABCA4-associated Stargardt disease (STGD1). We aimed to elucidate the genomic and transcriptomic variation in 1054 unsolved STGD and STGD-like probands.

METHODS

Sequencing of the complete 128-kb ABCA4 gene was performed using single-molecule molecular inversion probes (smMIPs), based on a semiautomated and cost-effective method. Structural variants (SVs) were identified using relative read coverage analyses and putative splice defects were studied using in vitro assays.

RESULTS

In 448 biallelic probands 14 known and 13 novel deep-intronic variants were found, resulting in pseudoexon (PE) insertions or exon elongations in 105 alleles. Intriguingly, intron 13 variants c.1938-621G>A and c.1938-514G>A resulted in dual PE insertions consisting of the same upstream, but different downstream PEs. The intron 44 variant c.6148-84A>T resulted in two PE insertions and flanking exon deletions. Eleven distinct large deletions were found, two of which contained small inverted segments. Uniparental isodisomy of chromosome 1 was identified in one proband.

CONCLUSION

Deep sequencing of ABCA4 and midigene-based splice assays allowed the identification of SVs and causal deep-intronic variants in 25% of biallelic STGD1 cases, which represents a model study that can be applied to other inherited diseases.

Kinetic versus thermodynamic control of mutational effects on protein homeostasis: A perspective from computational modeling and experiment

The effect of mutations in individual proteins on protein homeostasis, or "proteostasis," can in principle depend on the mutations' effects on the thermodynamics or kinetics of folding, or both. Here, we explore this issue using a computational model of in vivo protein folding that we call FoldEcoSlim. Our model predicts that kinetic versus thermodynamic control of mutational effects on proteostasis hinges on the relationship between how fast a protein's folding reaction reaches equilibrium and a critical time scale that characterizes the lifetime of a protein in its environment: for rapidly dividing bacteria, this time scale is that of cell division; for proteins that are produced in heterologous expression systems, this time scale is the amount of time before the protein is harvested; for proteins that are synthesized in and then exported from the eukaryotic endoplasmic reticulum, this time scale is that of protein secretion, and so forth. This prediction was validated experimentally by examining the expression yields of the wild type and several destabilized mutants of a model protein, the mouse ortholog of cellular retinoic acid-binding protein 1.

Functional Constraint Profiling of a Viral Protein Reveals Discordance of Evolutionary Conservation and Functionality

Viruses often encode proteins with multiple functions due to their compact genomes. Existing approaches to identify functional residues largely rely on sequence conservation analysis. Inferring functional residues from sequence conservation can produce false positives, in which the conserved residues are functionally silent, or false negatives, where functional residues are not identified since they are species-specific and therefore non-conserved. Furthermore, the tedious process of constructing and analyzing individual mutations limits the number of residues that can be examined in a single study. Here, we developed a systematic approach to identify the functional residues of a viral protein by coupling experimental fitness profiling with protein stability prediction using the influenza virus polymerase PA subunit as the target protein. We identified a significant number of functional residues that were influenza type-specific and were evolutionarily non-conserved among different influenza types. Our results indicate that type-specific functional residues are prevalent and may not otherwise be identified by sequence conservation analysis alone. More importantly, this technique can be adapted to any viral (and potentially non-viral) protein where structural information is available.

The analysis of sequence conservation is a common approach to identify functional residues within a protein. However, not all functional residues are conserved as natural evolution and species diversification permit continuous innovation of protein functionality through the retention of advantageous mutations. Non-conserved functional residues, which are often species-specific, may not be identified by conventional analysis of sequence conservation despite being biologically important. Here we described a novel approach to identify functional residues within a protein by coupling a high-throughput experimental fitness profiling approach with computational protein modeling. Our methodology is independent of sequence conservation and is applicable to any protein where structural information is available. In this study, we systematically mapped the functional residues on the influenza A PA protein and revealed that non-conserved functional residues are prevalent. Our results not only have significant implication on how functionality evolves during natural evolution, but also highlight the caveats when applying conservation-based approaches to identify functional residues within a protein.

Effect of salt on the formation of salt-bridges in β-hairpin peptides

The formation of salt-bridges in β-hairpin peptides is measured in increasing salt concentrations, indicating a decrease in the salt-bridged population due to charge–charge screening, as well as non-cooperative salt-bridge triads.

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.

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.

Cross-strand split tetra-Cys motifs as structure sensors in a beta-sheet protein

We have designed "split tetra-Cys motifs" that bind the biarsenical fluorescein dye 4',5'-bis(1,3,2-dithioarsolan-2-yl)fluorescein (FlAsH) across strands of a model beta-rich protein. Our strategy was to divide the linear FlAsH binding tetra-Cys sequence such that dye could be fully liganded only when the strands were arranged in space correctly by native protein conformational proximities. We introduced pairs of alternating cysteines on adjacent beta strands of cellular retinoic acid binding protein to create FlAsH binding sites in the native structure. Selective labeling occurred both in vitro and in vivo relative to sites with fewer than four Cys or with inappropriate geometry. Interestingly, two of the split tetra-Cys motif-carrying proteins bound FlAsH whether native or urea unfolded, while one was capable of binding FlAsH only when native. This latter design exemplifies the potential of split motifs as structure sensors.

Common physical basis of macromolecule-binding sites in proteins

Protein–DNA/RNA/protein interactions play critical roles in many biological functions. Previous studies have focused on the different features characterizing the different macromolecule-binding sites and approaches to detect these sites. However, no common unique signature of these sites had been reported. Thus, this work aims to provide a ‘common’ principle dictating the location of the different macromolecule-binding sites founded upon fundamental principles of binding thermodynamics. To achieve this aim, a comprehensive set of structurally nonhomologous DNA-, RNA-, obligate protein- and nonobligate protein-binding proteins, both free and bound to their respective macromolecules, was created and a novel strategy for detecting clusters of residues with electrostatic or steric strain given the protein structure was developed. The results show that regardless of the macromolecule type, the binding strength and conformational changes upon binding, macromolecule-binding sites are energetically less stable than nonmacromolecule-binding sites. They also reveal new energetic features distinguishing DNA- from RNA-binding sites and obligate protein- from nonobligate protein-binding sites in both free/bound protein structures.

Transporter-to-trap conversion: a disulfide bond formation in cellular retinoic acid binding protein I mutant triggered by retinoic acid binding irreversibly locks the ligand inside the protein

Transport proteins must bind their ligands reversibly to enable release at the point of delivery, while irreversible binding is usually associated with the extreme cases of ligand sequestration. Protein conformational dynamics is an important modulator of binding kinetics, as increased flexibility in the regions adjacent to the binding site may facilitate both association and dissociation processes. Ligand entry to, and exit from, the internal binding site of the cellular retinoic acid binding protein I (CRABP I) occurs via a flexible portal region, which functions as a dynamic aperture. We designed and expressed a CRABP I mutant (A35C/T57C), in which a small-scale conformational switch caused by the ligand binding event triggers formation of a disulfide bond in the portal region, thereby arresting structural fluctuations and effectively locking the ligand inside the binding cavity. At the same time, no formation of the disulfide bond is observed in the apo form of the mutant, and most characteristics of the mutant, including protein stability, are very similar to those of the wild-type protein in the absence of retinoic acid. The mutation does not alter the kinetics of retinoic acid binding to the protein, although the disulfide formation makes the binding effectively irreversible, as suggested by the absence of retinoic acid transfer from the holo form of the mutant to lipid vesicles in the absence of a reducing agent. Taken together, these data suggest that the disulfide bond formation in the portal region arrests large-scale structural fluctuations, which are required for retinoic acid release from the protein. The unique properties of the CRABP I mutant described in this work can be used to inspire and guide a design of nanodevices for multiple tasks ranging from sequestering small-molecule toxins in both tissue and circulation to nutrient deprivation of pathogens.

Site-specific fluorescent labeling of poly-histidine sequences using a metal-chelating cysteine

Coupling genetically encoded target sequences with specific and selective labeling strategies has made it possible to utilize fluorescence spectroscopy in complex mixtures to investigate the structure, function, and dynamics of proteins. Thus, there is a growing need for a repertoire of such labeling approaches to deploy based on a given application and to utilize in combination with one another by orthogonal reactivity. We have developed a simple approach to synthesize a fluorescent probe that binds to a poly-histidine sequence. The amino group of cysteine was converted into nitrilotriacetate to create a metal-chelating cysteine molecule, Cys-nitrilotriacetate. Two Cys-nitrilotriacetate molecules were then cross-linked using dibromobimane to generate a fluorophore capable of binding a His-tag on a protein, NTA(2)-BM. NTA(2)-BM is a potential fluorophore for selective tagging of proteins in vivo.

Relating destabilizing regions to known functional sites in proteins

BACKGROUND

Most methods for predicting functional sites in protein 3D structures, rely on information on related proteins and cannot be applied to proteins with no known relatives. Another limitation of these methods is the lack of a well annotated set of functional sites to use as benchmark for validating their predictions. Experimental findings and theoretical considerations suggest that residues involved in function often contribute unfavorably to the native state stability. We examine the possibility of systematically exploiting this intrinsic property to identify functional sites using an original procedure that detects destabilizing regions in protein structures. In addition, to relate destabilizing regions to known functional sites, a novel benchmark consisting of a diverse set of hand-curated protein functional sites is derived.

RESULTS

A procedure for detecting clusters of destabilizing residues in protein structures is presented. Individual residue contributions to protein stability are evaluated using detailed atomic models and a force-field successfully applied in computational protein design. The most destabilizing residues, and some of their closest neighbours, are clustered into destabilizing regions following a rigorous protocol. Our procedure is applied to high quality apo-structures of 63 unrelated proteins. The biologically relevant binding sites of these proteins were annotated using all available information, including structural data and literature curation, resulting in the largest hand-curated dataset of binding sites in proteins available to date. Comparing the destabilizing regions with the annotated binding sites in these proteins, we find that the overlap is on average limited, but significantly better than random. Results depend on the type of bound ligand. Significant overlap is obtained for most polysaccharide- and small ligand-binding sites, whereas no overlap is observed for most nucleic acid binding sites. These differences are rationalised in terms of the geometry and energetics of the binding site.

CONCLUSION

We find that although destabilizing regions as detected here can in general not be used to predict binding sites in protein structures, they can provide useful information, particularly on the location of functional sites that bind polysaccharides and small ligands. This information can be exploited in methods for predicting function in protein structures with no known relatives. Our publicly available benchmark of hand-curated functional sites in proteins should help other workers derive and validate new prediction methods.

Prediction of functional sites based on the fuzzy oil drop model

A description of many biological processes requires knowledge of the 3-D structure of proteins and, in particular, the defined active site responsible for biological function. Many proteins, the genes of which have been identified as the result of human genome sequencing, and which were synthesized experimentally, await identification of their biological activity. Currently used methods do not always yield satisfactory results, and new algorithms need to be developed to recognize the localization of active sites in proteins. This paper describes a computational model that can be used to identify potential areas that are able to interact with other molecules (ligands, substrates, inhibitors, etc.). The model for active site recognition is based on the analysis of hydrophobicity distribution in protein molecules. It is shown, based on the analyses of proteins with known biological activity and of proteins of unknown function, that the region of significantly irregular hydrophobicity distribution in proteins appears to be function related.

Localization of ligand binding site in proteins identified in silico

Knowledge-based models for protein folding assume that the early-stage structural form of a polypeptide is determined by the backbone conformation, followed by hydrophobic collapse. Side chain-side chain interactions, mostly of hydrophobic character, lead to the formation of the hydrophobic core, which seems to stabilize the structure of the protein in its natural environment. The fuzzy-oil-drop model is employed to represent the idealized hydrophobicity distribution in the protein molecule. Comparing it with the one empirically observed in the protein molecule reveals that they are not in agreement. It is shown in this study that the irregularity of hydrophobic distributions is aim-oriented. The character and strength of these irregularities in the organization of the hydrophobic core point to the specificity of a particular protein's structure/function. When the location of these irregularities is determined versus the idealized fuzzy-oil-drop, function-related areas in the protein molecule can be identified. The presented model can also be used to identify ways in which protein-protein complexes can possibly be created. Active sites can be predicted for any protein structure according to the presented model with the free prediction server at http://www.bioinformatics.cm-uj.krakow.pl/activesite. The implication based on the model presented in this work suggests the necessity of active presence of ligand during the protein folding process simulation.

On the derivation of propensity scales for predicting exposed transmembrane residues of helical membrane proteins

Helical membrane proteins (HMPs) play a crucial role in diverse physiological processes. Given the difficulty in determining their structures by experimental techniques, it is desired to develop computational methods for predicting the burial status of transmembrane residues. Deriving a propensity scale for the 20 amino acids to be exposed to the lipid bilayer from known structures is central to developing such methods. A fundamental problem in this regard is what would be the optimal way of deriving propensity scales. Here, we show that this problem can be reformulated such that an optimal scale is straightforwardly obtained in an analytical fashion. The derived scale favorably compares with others in terms of both algorithmic optimality and practical prediction accuracy. It also allows interesting insights into the structural organization of HMPs. Furthermore, the presented approach can be applied to other bioinformatics problems of HMPs, too. All the data sets and programs used in the study and detailed primary results are available upon request.

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.

Point mutations in protein globular domains: contributions from function, stability and misfolding

Several contrasting hypotheses have been formulated about the influence of functional and conformational properties, like stability and avoidance of misfolding, on the evolution of protein globular domains. Selection at functional sites has been suggested to be detrimental to stability or coupled to it. Avoidance of misfolding may be achieved by discarding misfolding-prone sequences or by maintaining a stable native state and thus destabilizing partially or fully unfolded states from which misfolding can take place. We have performed a hierarchical analysis of a large database of point mutations to dissect the relative contributions of function, stability and misfolding in the evolution of natural sequences. We show that at catalytic sites, selection for function overrules selection for stability but find no evidence for an anticorrelation between function and stability. Selection for stability plays a secondary role at binding sites, but is not fully coupled to selection for function. Remarkably, we did not find a selective pressure against misfolding-prone sequences in globular proteins at the level of individual positions. We suggest that such a selection would compromise native-state stability due to a correlation between the stabilities of native and misfolded states. Stabilization of the native state is the most frequent way in which natural proteins avoid misfolding.

Transient structural disorder as a facilitator of protein-ligand binding: native H/D exchange-mass spectrometry study of cellular retinoic acid binding protein I

Binding of all-trans Retinoic Acid (RA) to Cellular Retinoic Acid Binding Protein I (CRABP I) does not result in significant changes of the protein tertiary structure, even though the binding site is inaccessible in a static apo-protein conformation. One of the proposed scenarios for the protein-ligand binding process invokes the notion of a flexible portal region adjacent to the binding site, while another model suggests that the requisite dynamic events are induced by dimerization of the apo-protein in solution. In this work, RA binding to CRABP I is studied in dilute solutions (low micro-molar range), where no dimer and/or oligomer formation occurs. Modulation of backbone dynamics within various segments of the protein by its ligand is assessed using a combination of hydrogen exchange, electrospray ionization mass spectrometry, and collision-induced dissociation of protein ions in the gas phase. Consistent with the portal model of ligand entry, several protein segments (most of them containing residues making hydrophobic contacts to RA in the holo-form of the protein) are flexible in the absence of the ligand. At the same time, the two segments containing arginine residues forming a salt bridge with RA form the least flexible region in the apo-form of the protein. Although the presence of RA in solution reduces flexibility of all protein segments, the largest effect is observed within four strands that form one of the two beta-sheets enveloping a cavity which houses the ligand-binding site. These results are consistent with a model in which ligand binding occurs through a partially unstructured state of the protein with unobstructed access to the ligand-binding site. This intermediate (whose core is formed by the two stable arginine-containing strands) corresponds to a relatively low-energy local minimum on the apo-protein energy surface and is frequently sampled under native conditions.

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

Proteins in the intracellular lipid-binding protein (iLBP) family show remarkably high structural conservation despite their low-sequence identity. A multiple-sequence alignment using 52 sequences of iLBP family members revealed 15 fully conserved positions, with a disproportionately high number of these (n=7) located in the relatively small helical region. The conserved positions displayed high structural conservation based on comparisons of known iLBP crystal structures. It is striking that the beta-sheet domain had few conserved positions, despite its high structural conservation. This observation prompted us to analyze pair-wise interactions within the beta-sheet region to ask whether structural information was encoded in interacting amino acid pairs. We conducted this analysis on the iLBP family member, cellular retinoic acid-binding protein I (CRABP I), whose folding mechanism is under study in our laboratory. Indeed, an analysis based on a simple classification of hydrophobic and polar amino acids revealed a network of conserved interactions in CRABP I that cluster spatially, suggesting a possible nucleation site for folding. Significantly, a small number of residues participated in multiple conserved interactions, suggesting a key role for these sites in the structure and folding of CRABP I. The results presented here correlate well with available experimental evidence on folding of CRABPs and their family members and suggest future experiments. The analysis also shows the usefulness of considering pair-wise conservation based on a simple classification of amino acids, in analyzing sequences and structures to find common core regions among homologues.

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.

Specificity determinants for lipids bound to beta-barrel proteins

The family of proteins accountable for the intracellular movement of lipids is characterized by a 10-stranded beta-barrel that forms an internalized cavity varying in size and binding preferences. The loop connecting beta-strands E and F (the fifth and sixth strands) is the most striking conformational difference between adipocyte lipid binding protein (ALBP; fatty acids) and cellular retinoic acid binding protein type I (CRABP I). A three-residue mutation was made in wild-type (WT)-ALBP [ALBP with a three-residue mutation (EF-ALBP)] to mimic CRABP I. Crystal structures of ligand-free and EF-ALBP with bound oleic acid were solved to resolutions of 1.5 A and 1.7 A, respectively, and compared with previous studies of WT-ALBP. The changes in three residues of one loop of the protein appear to have altered the positioning of the C18 fatty acid, as observed in the electron density of EF-ALBP. The crystallographic studies made it possible to compare the protein conformation and ligand positioning with those found in the WT protein. Although the cavity binding sites in both the retinoid and fatty acid binding proteins are irregular, the ligand atoms appear to favor a relatively planar region of the cavities. Preliminary chemical characterization of the mutant protein indicated changes in some binding properties and overall protein stability.

Indirect assessment of small hydrophobic ligand binding to a model protein using a combination of ESI MS and HDX/ESI MS

Direct mass spectrometric characterization of interactions between proteins and small hydrophobic ligands often poses a serious problem due to the complex instability in the gas phase. We have developed a method that probes the efficacy of ligand-protein interactions indirectly by monitoring changes in protein flexibility. The latter is assessed quantitatively using a combination of charge state distribution analysis and amide hydrogen exchange under both native and mildly denaturing conditions. The method was used to evaluate binding of a model protein cellular retinoic acid binding protein I to its natural ligand all-trans retinoic acid (RA), isomers 13-cis- and 9-cis-RA, and retinol, yielding the following order of ligand affinities: All-trans RA > 9-cis RA > 13-cis RA, with no detectable binding of retinol. This order is in agreement with the results of earlier fluorimetric titration studies. Furthermore, binding energy of the protein to each of retinoic acid isomers was determined based on the measured hydrogen exchange kinetics data acquired under native conditions.

A simple electrostatic criterion for predicting the thermal stability of proteins

The enhancement of protein thermostability is an important issue for both basic science and biotechnology purposes. We have developed a thermostability criterion for a protein in terms of a quasi-electric dipole moment (contributed by its charged residues) defined for an electric charge distribution whose total charge is not zero. It was found that minimization of the modulus of this dipole moment increased its thermal stability, as demonstrated by surveying these values in pairs of mesostable-thermostable homologous proteins and in mutations described in the literature. The analysis of these observations provides criteria for thermostabilization of a protein, by computing its dipole profile. This profile is obtained by direct substitution of each amino acid of the sequence by either a positive, negative or neutral amino acid, followed by a recalculation of the dipole moment. As an experimental example, these criteria were applied to a beta-glucanase to enhance its thermal stability.

Empirical relationships between protein structure and carboxyl pKa values in proteins

Relationships between protein structure and ionization of carboxyl groups were investigated in 24 proteins of known structure and for which 115 aspartate and 97 glutamate pK(a) values are known. Mean pK(a) values for aspartates and glutamates are < or = 3.4 (+/-1.0) and 4.1 (+/-0.8), respectively. For aspartates, mean pK(a) values are 3.9 (+/-1.0) and 3.1 (+/-0.9) in acidic (pI < 5) and basic (pI > 8) proteins, respectively, while mean pK(a) values for glutamates are approximately 4.2 for acidic and basic proteins. Burial of carboxyl groups leads to dispersion in pK(a) values: pK(a) values for solvent-exposed groups show narrow distributions while values for buried groups range from < 2 to 6.7. Calculated electrostatic potentials at the carboxyl groups show modest correlations with experimental pK(a) values and these correlations are not improved by including simple surface-area-based terms to account for the effects of desolvation. Mean aspartate pK(a) values decrease with increasing numbers of hydrogen bonds but this is not observed at glutamates. Only 10 pK(a) values are > 5.5 and most are found in active sites or ligand-binding sites. These carboxyl groups are buried and usually accept no more than one hydrogen bond. Aspartates and glutamates at the N-termini of helices have mean pK(a) values of 2.8 (+/-0.5) and 3.4 (+/-0.6), respectively, about 0.6 units less than the overall mean values.

Prediction of functionally important residues based solely on the computed energetics of protein structure

Catalytic and other functionally important residues in proteins can often be mutated to yield more stable proteins. Many of these residues are charged residues that are located in electrostatically unfavorable environments. Here it is demonstrated that because continuum electrostatics methods can identify these destabilizing residues, the same methods can also be used to identify functionally important residues in otherwise uncharacterized proteins. To establish this point, detailed calculations are performed on six proteins for which good structural and mutational data are available from experiments. In all cases it is shown that functionally important residues known to be destabilizing experimentally are among the most destabilizing residues found in the calculations. A larger scale analysis performed on 216 different proteins demonstrates the existence of a general relationship between the calculated electrostatic energy of a charged residue and its degree of evolutionary conservation. This relationship becomes obscured when electrostatic energies are calculated using Coulomb's law instead of the more complete continuum electrostatics method. Finally, in a first predictive application of the method, calculations are performed on three proteins whose structures have recently been reported by a structural genomics consortium.

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.

Direct photoaffinity labeling of cellular retinoic acid-binding protein I (CRABP-I) with all-trans-retinoic acid: identification of amino acids in the ligand binding site

Cellular retinoic acid-binding proteins I and II (CRABP-I and -II, respectively) are transport proteins for all-trans-retinoic acid (RA), an active metabolite of vitamin A (retinol), and have been reported to be directly involved in the metabolism of RA. In this study, direct photoaffinity labeling with [11,12-(3)H]RA was used to identify amino acids comprising the ligand binding site of CRABP-I. Photoaffinity labeling of CRABP-I with [(3)H]RA was light- and concentration-dependent and was protected by unlabeled RA and various retinoids, indicating that the labeling was directed to the RA-binding site. Photolabeled CRABP-I was hydrolyzed with endoproteinase Lys-C to yield radioactive peptides, which were separated by reversed-phase HPLC for analysis by Edman degradation peptide sequencing. This method identified five modified amino acids from five separate HPLC fractions: Trp7, Lys20, Arg29, Lys38, and Trp109. All five amino acids are located within one side of the "barrel" structure in the area indicated by the reported crystal structure as the ligand binding site. This is the first direct identification of specific amino acids in the RA-binding site of CRABPs by photoaffinity labeling. These results provide significant information about the ligand binding site of the CRABP-I molecule in solution.

Multiple roles of prolyl residues in structure and folding

To explore the ways that proline residues may influence the conformational options of a polypeptide backbone, we have characterized Pro-->Ala mutants of cellular retinoic acid-binding protein I (CRABP I). While all three Xaa-Pro bonds are in the trans conformation in the native protein and the equilibrium stability of each mutant is similar to that of the parent protein, each has distinct effects on folding and unfolding kinetics. The mutation of Pro105 does not alter the kinetics of folding of CRABP I, which indicates that the flexible loop containing this residue is passive in the folding process. By contrast, replacement of Pro85 by Ala abolishes the observable slow phase of folding, revealing that correct configuration of the 84-85 peptide bond is prerequisite to productive folding. Substitution of Pro39 by Ala yields a protein that folds and unfolds more slowly. Removal of the conformational constraint imposed by the proline ring likely raises the transition state barrier by increasing the entropic cost of narrowing the conformational ensemble. Additionally, the Pro-->Ala mutation removes a helix-termination signal that is important for efficient folding to the native state.

Contribution of the hydrogen-bond network involving a tyrosine triad in the active site to the structure and function of a highly proficient ketosteroid isomerase from Pseudomonas putida biotype B

Delta(5)-3-Ketosteroid isomerase from Pseudomonas putida biotype B is one of the most proficient enzymes catalyzing an allylic isomerization reaction at rates comparable to the diffusion limit. The hydrogen-bond network (Asp99... Wat504...Tyr14...Tyr55...Tyr30) which links the two catalytic residues, Tyr14 and Asp99, to Tyr30, Tyr55, and a water molecule in the highly apolar active site has been characterized in an effort to identify its roles in function and stability. The DeltaG(U)(H2O) determined from equilibrium unfolding experiments reveals that the elimination of the hydroxyl group of Tyr14 or Tyr55 or the replacement of Asp99 with leucine results in a loss of conformational stability of 3.5-4.4 kcal/mol, suggesting that the hydrogen bonds of Tyr14, Tyr55, and Asp99 contribute significantly to stability. While decreasing the stability by about 6.5-7.9 kcal/mol, the Y55F/D99L or Y30F/D99L double mutation also reduced activity significantly, exhibiting a synergistic effect on k(cat) relative to the respective single mutations. These results indicate that the hydrogen-bond network is important for both stability and function. Additionally, they suggest that Tyr14 cannot function efficiently alone without additional support from the hydrogen bonds of Tyr55 and Asp99. The crystal structure of Y55F as determined at 1.9 A resolution shows that Tyr14 OH undergoes an alteration in orientation to form a new hydrogen bond with Tyr30. This observation supports the role of Tyr55 OH in positioning Tyr14 properly to optimize the hydrogen bond between Tyr14 and C3-O of the steroid substrate. No significant structural changes were observed in the crystal structures of Y30F and Y30F/Y55F, which allowed us to estimate approximately the interaction energies mediated by the hydrogen bonds Tyr30...Tyr55 and Tyr14...Tyr55. Taken together, our results demonstrate that the hydrogen-bond network provides the structural support that is needed for the enzyme to maintain the active-site geometry optimized for both function and stability.

Unfolding dynamics of a beta-sheet protein studied by mass spectrometry

The unfolding dynamics of cellular retinoic acid-binding protein I (CRABP I), an 18 kDa predominantly beta-sheet protein, were studied by monitoring the hydrogen-deuterium (H-D) exchange reaction under various solution conditions. A bimodal charge state distribution was observed when a denaturing agent was added to the protein aqueous solution. These two populations exhibit different kinetics of H-D exchange, with the high charge state ions undergoing very rapid isotope exchange, while the low charge state protein ions exchange cooperatively but at much slower rates. Transiently populated intermediate states were detected indirectly using hydrogen exchange measurement in aqueous solution at various pHs. At pH 2.5 and room temperature, three distinct populations of CRABP I ions exist over an extended period of time, each corresponding to a specific degree of backbone amide hydrogen atom protection. Mass spectral data are complementary to hydrogen exchange measurements by NMR, since the former samples a much faster time-scale of dynamic events in solution.

NMR study of the binding of all-trans-retinoic acid to type II human cellular retinoic acid binding protein

Cellular RA binding proteins are thought to play important roles in the (RA), a hormonally active metabolite of vitamin A that has profound effects on cell growth, + differentiation and morphogenesis. Binding of RA to type II human cellular RA binding proteins (CRABPII) has been investigated by NMR spectroscopy. The sequential resonance assignments of +CRABPII in the presence of RA were established by heteronuclear three-dimensional NMR at pH 7.3. The resonance assignments of the bound RA were achieved by homonucl NMR. The secondary structures of holo-CRABPII determined by NMR were ess as revealed by the crystal structure of holo-CRABPII. Most of the nuclear Overhauser effects (NOEs) between CRABPII and the bound RA were consistent with those predicted crystal structure of holo-CRABPII. The results suggested that the conformations in solution and in the crystalline state are highly similar. Compared to the ligand binding pocket, especially the ligand entrance, was stabilize Ser12-Leu18, one of the structure elements that constitute the ligand binding pocket, became more mobile upon binding of RA. Intramolecular NOEs between protons of the bo the carboxylate end of the bound RA is well fixed but the &beta;-ionone

Probing the folding pathway of a beta-clam protein with single-tryptophan constructs

BACKGROUND

Cellular retinoic acid binding protein I (CRABPI) is a small, predominantly beta-sheet protein with a simple architecture and no disulfides or cofactors. Folding of mutants containing only one of the three native tryptophans has been examined using stopped-flow fluorescence and circular dichroism at multiple wavelengths.

RESULTS

Within 10 ms, the tryptophan fluorescence of all three mutants shows a blue shift, and stopped-flow circular dichroism shows significant secondary structure content. The local environment of Trp7, a completely buried residue located near the intersection of the N and C termini, develops on a 100 ms time scale. Spectral signatures of the other two tryptophan residues (87 and 109) become native-like in a 1 s kinetic phase.

CONCLUSIONS

Formation of the native beta structure of CRABPI is initiated by rapid hydrophobic collapse, during which local segments of chain adopt significant secondary structure. Subsequently, transient yet specific interactions of amino acid residues restrict the arrangement of the chain topology and initiate long-range associations such as the docking of the N and C termini. The development of native tertiary environments, including the specific packing of the beta-sheet sidechains, occurs in a final, highly cooperative step simultaneous with stable interstrand hydrogen bonding.

Folding mechanism of three structurally similar beta-sheet proteins

The folding mechanism of cellular retinoic acid binding protein I (CRABP I), cellular retinol binding protein II (CRBP II), and intestinal fatty acid binding protein (IFABP) were investigated to determine if proteins with similar native structures have similar folding mechanisms. These mostly beta-sheet proteins have very similar structures, despite having as little as 33% sequence similarity. The reversible urea denaturation of these proteins was characterized at equilibrium by circular dichroism and fluorescence. The data were best fit by a two-state model for each of these proteins, suggesting that no significant population of folding intermediates were present at equilibrium. The native states were of similar stability with free energies (linearly extrapolated to 0 M urea, deltaGH2O) of 6.5, 8.3, and 5.5 kcal/mole for CRABP I, CRBP II, and IFABP, respectively. The kinetics of the folding and unfolding processes for these proteins was monitored by stopped-flow CD and fluorescence. Intermediates were observed during both the folding and unfolding of all of these proteins. However, the overall rates of folding and unfolding differed by nearly three orders of magnitude. Further, the spectroscopic properties of the intermediate states were different for each protein, suggesting that different amounts of secondary and/or tertiary structure were associated with each intermediate state for each protein. These data show that the folding path for proteins in the same structural family can be quite different, and provide evidence for different folding landscapes for these sequences.

NMR study suggests a major role for Arg111 in maintaining the structure and dynamical properties of type II human cellular retinoic acid binding protein

The solution structure of a site-directed mutant of type-II human cellular retinoic acid binding protein (CRABPII) with Arg111 replaced by methionine (R111M) has been determined by NMR spectroscopy. The sequential assignments of the 1H and 15N resonances of apo-R111M were established by multinuclear multidimensional NMR. The solution structure was calculated from 2302 distance restraints and 77 phi dihedral restraints derived from the NMR data. The root-mean-square deviation of the ensemble of 28 refined conformers that represent the structure from the mean coordinate set derived from them was 0.54 +/- 0.26 and 0.98 +/- 0.23 A for the backbone atoms and all heavy atoms, respectively. The solution structure of apo-R111M is similar to that of wild-type apo-CRABPII. However, there are significant conformational differences between the two proteins, localized mainly to three segments (Leu19-Ala36, Glu73-Cys81, and Leu99-Pro105) clustered around the ligand entrance more than 17 A away from the point mutation. In apo-R111M, all the three segments move toward the center of the ligand entrance so that the opening of the ligand-binding pocket in apo-R111M is much smaller than that in wild-type apo-CRABPII. Furthermore, the ligand-binding pocket of apo-R111M, especially the ligand entrance, is much less flexible than that of apo-CRABPII. Surprisingly, apo-R111M is more similar to holo-CRABPII than to apo-CRABPII in both structure and dynamical properties. The conformational and dynamical changes caused by the mutation are similar to those induced by binding of RA, although the magnitudes of the changes caused by the mutation are smaller than those induced by binding of RA. The results suggest that Arg111 plays a critical role in determining the structure and dynamical properties of CRABPII.

Cellular retinaldehyde-binding protein ligand interactions. Gln-210 and Lys-221 are in the retinoid binding pocket

Cellular retinaldehyde-binding protein (CRALBP) carries 11-cis-retinal and/or 11-cis-retinol as endogenous ligands in the retinal pigment epithelium (RPE) and Müller cells of the retina and has been linked with autosomal recessive retinitis pigmentosa. Ligand interactions determine the physiological role of CRALBP in the RPE where the protein is thought to function as a substrate carrier for 11-cis-retinol dehydrogenase in the synthesis of 11-cis-retinal for visual pigment regeneration. However, CRALBP is also present in optic nerve and brain where its natural ligand and function are not yet known. We have characterized the interactions of retinoids with native bovine CRALBP, human recombinant CRALBP (rCRALBP) and five mutant rCRALBPs. Efforts to trap and/or identify a Schiff base in the dark, under a variety of reducing, denaturing, and pH conditions were unsuccessful, suggesting the lack of covalent interactions between CRALBP and retinoid. Buried and solvent-exposed lysine residues were identified in bovine CRALBP by reductive methylation of the holoprotein followed by denaturation and reaction with [3H]acetic anhydride. Radioactive lysine residues were identified by Edman degradation and electrospray mass spectrometry following proteolysis and purification of modified peptides. Human rCRALBP mutants K152A, K221A, and K294A were prepared to investigate possible retinoid interactions with buried or partially buried lysines. Two other rCRALBP mutants, I162V and Q210R, were also prepared to identify substitutions altering the retinoid binding properties of a random mutant. The structures of all the mutants were verified by amino acid and mass spectral analyses and retinoid binding properties evaluated by UV-visible and fluorescence spectroscopy. All of the mutants bound 11-cis-retinal essentially like the wild type protein, indicating that the proteins were not grossly misfolded. Three of the mutants bound 9-cis-retinal like the wild type protein; however, Q210R and K221A bound less than stoichiometric amounts of the 9-cis-isomer and exhibited lower affinity for this retinoid relative to wild type rCRALBP. Residues Gln-210 and Lys-221 are located within a region of CRALBP exhibiting sequence homology with the ligand binding cavity of yeast phosphatidylinositol-transfer protein. The data implicate Gln-210 and Lys-221 as components of the CRALBP retinoid binding cavity and are discussed in the context of ligand interactions in structurally or functionally related proteins with known crystallographic structures.

Molecular cloning and characterization of an invertebrate cellular retinoic acid binding protein

We have cloned a cDNA and gene from the tobacco hornworm, Manduca sexta, which is related to the vertebrate cellular retinoic acid binding proteins (CRABPs). CRABPs are members of the superfamily of lipid binding proteins (LBPs) and are thought to mediate the effects of retinoic acid (RA) on morphogenesis, differentiation, and homeostasis. This discovery of a Manduca sexta CRABP (msCRABP) demonstrates the presence of a CRABP in invertebrates. Compared with bovine/murine CRABP I, the deduced amino acid sequence of msCRABP is 71% homologous overall and 88% homologous for the ligand binding pocket. The genomic organization of msCRABP is conserved with other CRABP family members and the larger LBP superfamily. Importantly, the promoter region contains a motif that resembles an RA response element characteristic of the promoter region of most CRABPs analyzed. Three-dimensional molecular modeling based on postulated structural homology with bovine/murine CRABP I shows msCRABP has a ligand binding pocket that can accommodate RA. The existence of an invertebrate CRABP has significant evolutionary implications, suggesting CRABPs appeared during the evolution of the LBP superfamily well before vertebrate/invertebrate divergence, instead of much later in evolution in selected vertebrates.

Crystal structure of apo-cellular retinoic acid-binding protein type II (R111M) suggests a mechanism of ligand entry

The crystal structure of unliganded mutant R111M of human cellular retinoic acid-binding protein type II (apo-CRABPII (R111M)) has been determined at 2.3 A and refined to a crystallographic R-factor of 0. 18. Although the mutant protein has lower affinity for all-trans-retinoic acid (RA) than the wild-type, it is properly folded, and its conformation is very similar to the wild-type. apo-CRABPII (R111M) crystallizes in space group P1 with two molecules in the unit cell. The two molecules have high structural similarity except that their alpha2 helices differ strikingly. Analyses of the molecular conformation and crystal packing environment suggest that one of the two molecules assumes a conformation compatible with RA entry. Three structural elements encompassing the opening of the binding pocket exhibit large conformational changes, when compared with holo-CRABPII, which include the alpha2 helix and the betaC-betaD and betaE-betaF hairpin loops. The alpha2 helix is unwound at its N terminus, which appears to be essential for the opening of the RA-binding pocket. Three arginine side-chains (29, 59, and 132) are found with their guanidino groups pointing into the RA-binding pocket. A three-step mechanism of RA entry has been proposed, addressing the opening of the RA entrance, the electrostatic potential that directs entry of RA into the binding pocket, and the intramolecular interactions that stabilize the RA.CRABPII complex via locking the three flexible structural elements when RA is bound.

Coupling protein stability and protein function in Escherichia coli CspA

BACKGROUND

CspA is a small protein that binds single-stranded RNA and DNA. The binding site of CspA consists of a cluster of aromatic amino acids, which form an unusually large nonpolar patch on the surface of the protein. Because nonpolar residues are generally found in the interiors of proteins, this cluster may have evolved to bind nucleic acids at the expense of protein stability.

RESULTS

Three neighboring phenylalanines have been mutated singly and in combination to leucine and to serine. All mutations adversely affect DNA binding. Surprisingly, all mutations, and especially those to serine, are destabilizing.

CONCLUSIONS

The aromatic cluster in CspA is required not only for protein function but also for protein stability. This result is pertinent to the design of beta-sheet proteins and single-stranded nucleic acid binding proteins, whose binding mode is proposed to be of aromatic-aromatic intercalation.

Soluble expression in E. coli of a functional interphotoreceptor retinoid-binding protein module fused to thioredoxin: correlation of vitamin A binding regions with conserved domains of C-terminal processing proteases

The exchange of all-trans retinol and 11-cis retinal between the photoreceptors and retinal pigmented epithelium is mediated by interphotoreceptor retinoid-binding protein (IRBP). IRBP contains binding sites for retinoids, docosahexaenoic acid and probably cell surface and matrix receptors. IRBP arose through the quadruplication of an ancient protein, represented by its carboxy-terminal module (module 4 in amphibians and mammals). Module 4 has retinol binding activity and is composed of regions coded for by each of IRBP's four exons. Determining the function of the exons has been hampered by insoluble expression of module 4 in Escherichia coli. Here, we found that module 4 of Xenopus IRBP (X4IRBP), as well as its exon segments, can be expressed in a soluble form as thioredoxin fusion proteins. The recombinant proteins were purified by ion exchange and arsenical-based affinity chromatography. Liquid chromatography/mass spectrometry confirmed that the sequence of X4IRBP is correct. All-trans retinol binding was characterized by monitoring enhancement of retinol fluorescence, quenching of intrinsic protein fluorescence, and transfer of energy to the bound retinol. Retinol bound to X4IRBP at 2.20+/-0.29 sites with a KD=1.25+/-0.39. One of the two sites was localized to Exons(2+3) and had a KD=0.26+/-0.13 micron. This site, which supported protein quenching and energy transfer, probably contains at least one of the two conserved tryptophans present in this segment. The second site was localized to Exon 4. This site supported the enhancement of retinol fluorescence but not protein quenching or energy transfer and had a KD=1.94+/-0.20 micron. Exon 1 had no retinol binding activity. The location of the retinol binding regions correlated with the distribution of domains conserved between IRBPs and the newly recognized family of C-terminal processing proteases (CtpAs), proteins which bind and cleave non-polar carboxy termini.

Cavity formation before stable hydrogen bonding in the folding of a beta-clam protein

The time course of folding of a small beta-sheet protein reveals formation of a central ligand binding cavity before the consolidation of the native hydrogen bonding network. These results suggest that side chain interactions and not stable hydrogen bonding determine the beta-sheet architecture and play crucial roles in the overall chain topology.

Lipocalin-type prostaglandin D synthase (beta-trace) is a newly recognized type of retinoid transporter

Lipocalin-type prostaglandin D synthase is responsible for the biosynthesis of prostaglandin D2 in the central nervous system and the genital organs and is secreted into the cerebrospinal fluid and the seminal plasma as beta-trace. Here we analyzed retinoids binding of the enzyme by monitoring the fluorescence quenching of an intrinsic tryptophan residue, and appearance of circular dichroism around 330 nm, and a red shift of the UV absorption spectra of retinoids. We found that the enzyme binds all-trans- or 9-cis-retinoic acid and all-trans- or 13-cis-retinal, but not all-trans-retinol, with affinities (Kd of 70-80 nM) sufficient for function as a retinoid transporter. All-trans-retinoic acid inhibited the enzyme activity in a noncompetitive manner, suggesting that it binds to the same hydrophobic pocket as prostaglandin H2, the substrate for prostaglandin D synthase, but at a different site in this pocket. It is likely that this enzyme is a bifunctional protein that acts as both retinoid transporter and prostaglandin D2-producing enzyme.

Biochemical and crystallographic analyses of a portal mutant of the adipocyte lipid-binding protein

A number of crystallographic studies of the adipocyte lipid-binding protein have established that the fatty acid-binding site is within an internalized water-filled cavity. The same studies have also suggested the existence of a region physically distinct from the fatty acid-binding site which connects the cavity of the protein with the external solvent, hereafter referred to as the portal. In an effort to examine the portal region, we have used site-directed mutagenesis to introduce the mutations V32D/F57H into the murine ALBP cDNA. Mutant protein has been isolated, crystallized, and its stability and binding properties studied by biochemical methods. As assessed by guanidine-HCl denaturation, the mutant form exhibited a slight overall destabilization relative to the wild-type protein under both acid and alkaline conditions. Accessibility to the cavity in both the mutant and wild-type proteins was observed by stopped-flow analysis of the modification of a cavity residue, Cys117, by the sulfhydryl reactive agent 5, 5'-dithiobis(2-nitrobenzoic acid) at pH 8.5. Cys117 of V32D/F57H ALBP was modified 7-fold faster than the wild-type protein. The ligand binding properties of both the V32D/F57H mutant and wild-type proteins were analyzed using a fluorescent probe at pH 6.0 and 8.0. The apparent dissociation constants for 1-anilinonaphthalene-8-sulfonic acid were approximately 9-10-fold greater than the wild-type protein, independent of pH. In addition, there is a 6-fold increase in the Kd for oleic acid for the portal mutant relative to the wild-type at pH 8.0. To study the effect of pH on the double mutant, it was crystallized and analyzed in two distinct space groups at pH 4.5 and 6.4. While in general the differences in the overall main chain conformations are negligible, changes were observed in the crystallographic structures near the site of the mutations. At both pH values, the mutant side chains are positioned somewhat differently than in wild-type protein. To ensure that the mutations had not altered ionic conditions near the binding site, the crystallographic coordinates were used to monitor the electrostatic potentials from the head group site to the positions near the portal region. The differences in the electrostatic potentials were small in all regions, and did not explain the differences in ligand affinity. We present these results within the context of fatty acid binding and suggest lipid association is more complex than that described within a single equilibrium event.

Structure-function relationships of cellular retinoic acid-binding proteins. Quantitative analysis of the ligand binding properties of the wild-type proteins and site-directed mutants

It has been suggested that electrostatic interactions are critical for binding of retinoic acid by cellular retinoic acid-binding proteins (CRABP-I and CRABP-II). However, the roles of two conserved arginine residues (Arg-111 and Arg-131 in CRABP-I; Arg-111 and Arg-132 in CRABP-II) that interact with the carboxyl group of retinoic acid have not been evaluated. A novel competitive binding assay has been developed for measuring the relative dissociation constants of the site-directed mutants of CRABPs. Arg-111 and Arg-132 of CRABP-II were replaced with methionine by site-directed mutagenesis. The relative dissociation constants of R111M and R132M (Kd (R111M)/Kd (CRABP-II) and Kd (R132M)/Kd(CRABP-II)) were determined to be 40-45 and 6-8, respectively. The ring protons of the aromatic residues of the wild-type CRABP-II and the two mutants were sequentially assigned by two-dimensional homonuclear NMR in conjunction with three-dimensional heteronuclear NMR. Detailed analysis of the nuclear Overhauser effect spectroscopy spectra of the proteins indicated that the conformations of the two mutants are highly similar to that of the wild-type CRABP-II. These results taken together showed that Arg-111 and Arg-132 are important for binding retinoic acid but contribute to the binding energy only by approximately 2.2 and 1.2 kcal/mol, respectively. In addition, the relative dissociation constant of CRABP-II and CRABP-I (Kd (CRABP-II)/Kd (CRABP-I)) was determined to be 2-3, in close agreement with that calculated using the apparent Kd values determined under the same conditions by fluorometric titrations.