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

Opposite Effects of CRABP1 and CRABP2 Homologs on Proliferation of Breast Cancer Cells and Their Sensitivity to Retinoic Acid

Resistance of tumor cells to retinoic acid (RA), a promising therapeutic agent, is the major factor limiting the use of RA in clinical practice. The mechanisms of resistance to RA are still poorly understood. Cellular Retinoic Acid Binding Proteins, CRABP1 and CRABP2, are essential mediators of RA signaling, but role of the two CRABP homologs in regulating cellular sensitivity to RA has not been well studied. In addition, the effects of CRABP1 and CRABP2 on cell proliferation have not been compared. Here, using a broad panel of breast cancer cell lines with different levels of RA sensitivity/resistance, we show for the first time that in the RA-sensitive cells, CRABP1 expression is restricted by methylation, and protein levels are highly variable. In the moderately-RA-resistant cell lines, high level of CRABP1 is observed both at the mRNA and protein levels, unchanged by inhibition of DNA methylation. The cell lines with maximum resistance to RA are characterized by complete repression of CRABP1 expression realized at transcriptional and posttranscriptional levels, and exogenous expression of each of the CRABP homologs has no effect on the studied characteristics. CRABP1 and CRABP2 proteins have opposing effects on proliferation and sensitivity to RA. In particular, CRABP1 stimulates and CRABP2 reduces proliferation and resistance to RA in the initially RA-sensitive cells, while in the more resistant cells the role of each homolog in both of these parameters is reversed. Overall, we have shown for the first time that CRABP proteins exert different effects on the growth and sensitivity to RA of breast cancer cells (stimulation, suppression, or no effect) depending on the baseline level of RA-sensitivity, with the effects of CRABP1 and CRABP2 homologs on the studied properties always being opposite.

Structural and Dynamic Determinants of Molecular Recognition in Bile Acid-Binding Proteins

Disorders in bile acid transport and metabolism have been related to a number of metabolic disease states, atherosclerosis, type-II diabetes, and cancer. Bile acid-binding proteins (BABPs), a subfamily of intracellular lipid-binding proteins (iLBPs), have a key role in the cellular trafficking and metabolic targeting of bile salts. Within the family of iLBPs, BABPs exhibit unique binding properties including positive binding cooperativity and site-selectivity, which in different tissues and organisms appears to be tailored to the local bile salt pool. Structural and biophysical studies of the past two decades have shed light on the mechanism of bile salt binding at the atomic level, providing us with a mechanistic picture of ligand entry and release, and the communication between the binding sites. In this review, we discuss the emerging view of bile salt recognition in intestinal- and liver-BABPs, with examples from both mammalian and non-mammalian species. The structural and dynamic determinants of the BABP-bile–salt interaction reviewed herein set the basis for the design and development of drug candidates targeting the transcellular traffic of bile salts in enterocytes and hepatocytes.

Polyethylene glycol perturbs the unfolding of CRABP I: A correlation between experimental and theoretical approach

The importance of macromolecules paves the way towards a detailed molecular level investigation as all most all cellular processes occurring at the interior of cells in the form of proteins, enzymes, and other biological molecules are significantly affected because of their crowding. Thus, exploring the role of crowding environment on the denaturation and renaturation kinetics of protein molecules is of great importance. Here, CRABP I (cellular retinoic acid binding protein I) is employed as a model protein along with different molecular weights of Polyethylene glycol (PEG) as molecular crowders. The experimental evaluations are done by accessing the protein secondary structure analysis using circular dichroism (CD) spectroscopy and unfolding kinetics using intrinsic fluorescence of CRABP I at 37 °C to mimic the in vivo crowding environment. The unfolding kinetics results indicated that both PEG 2000 and PEG 4000 act as stabilizers by retarding the unfolding kinetic rates. Both kinetic and stability outcomes presented the importance of crowding environment on stability and kinetics of CRABP I. The molecular dynamics (MD) studies revealed that thirteen PEG 2000 molecules assembled during the 500 ns simulation, which increases the stability and percentage of β-sheet. The experimental findings are well supported by the molecular dynamics simulation results.

Long-range Regulation of Partially Folded Amyloidogenic Peptides

Neurodegeneration involves abnormal aggregation of intrinsically disordered amyloidogenic peptides (IDPs), usually mediated by hydrophobic protein-protein interactions. There is mounting evidence that formation of α-helical intermediates is an early event during self-assembly of amyloid-β42 (Aβ42) and α-synuclein (αS) IDPs in Alzheimer’s and Parkinson’s disease pathogenesis, respectively. However, the driving force behind on-pathway molecular assembly of partially folded helical monomers into helical oligomers assembly remains unknown. Here, we employ extensive molecular dynamics simulations to sample the helical conformational sub-spaces of monomeric peptides of both Aβ42 and αS. Our computed free energies, population shifts, and dynamic cross-correlation network analyses reveal a common feature of long-range intra-peptide modulation of partial helical folds of the amyloidogenic central hydrophobic domains via concerted coupling with their charged terminal tails (N-terminus of Aβ42 and C-terminus of αS). The absence of such inter-domain fluctuations in both fully helical and completely unfolded (disordered) states suggests that long-range coupling regulates the dynamicity of partially folded helices, in both Aβ42 and αS peptides. The inter-domain coupling suggests a form of intra-molecular allosteric regulation of the aggregation trigger in partially folded helical monomers. This approach could be applied to study the broad range of amyloidogenic peptides, which could provide a new path to curbing pathogenic aggregation of partially folded conformers into oligomers, by inhibition of sites far from the hydrophobic core.

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.

Local and non-local topological information in the denatured state ensemble of a β-barrel protein

The folding of predominantly β-sheet proteins is complicated by the presence of a large number of non-local interactions in their native states, which increase the ruggedness of their folding energy landscapes. However, forming non-local contacts early in folding or even in the unfolded state can smooth the energy landscape and facilitate productive folding. We report that several sequence regions of a β-barrel protein, cellular retinoic acid-binding protein 1 (CRABP1), populate native-like secondary structure to a significant extent in the denatured state in 8 M urea. In addition, we provide evidence for both local and non-local interactions in the denatured state of CRABP1. NMR chemical shift perturbations (CSPs) under denaturing conditions upon substitution of single residues by mutation support the presence of several non-local interactions in topologically key sites, arguing that the denatured state is conformationally restricted and contains topological information for the native fold. Among the most striking non-local interactions are those between the N- and C-terminal regions, which are involved in closure of the native β-barrel. In addition, CSPs support the presence of two features in the denatured state: a major hydrophobic cluster involving residues from various parts of the sequence and a native-like interaction similar to one identified in previous studies as forming early in folding (Budyak et al., Structure 21, 476 [2013]). Taken together, our data support a model in which transient structures involving nonlocal interactions prime early folding interactions in CRABP1, determine its barrel topology, and may protect this predominantly β-sheet protein against aggregation.

Protein Domain-Swapping Can Be a Consequence of Functional Residues

Monomer topology has been implicated in domain-swapping, a potential first step on the route to disease-causing protein aggregation. Despite having the same topology (β1-α1-β2-β3-β4-β5), the cysteine protease inhibitor stefin-B domain swaps more readily than a single-chain variant of the heterodimeric sweet protein monellin (scMn). Here, we computationally study the folding of stefin-B and scMn in order to understand the molecular basis for the difference in their domain-swapping propensities. In agreement with experiments, our structure-based simulations show that scMn folds cooperatively without the population of an intermediate while stefin-B populates an equilibrium intermediate state. Since the simulation intermediate has only one domain structured (β3-β4-β5), it can directly lead to domain-swapping. Using computational variants of stefin-B, we show that the population of this intermediate is caused by regions of stefin-B that have been implicated in protease inhibition. We also find that the protease-binding regions are located on two structural elements and localized in space. In contrast, the residues that contribute to the sweetness of monellin are not localized to a few structural elements but are distributed over the protein fold. We conclude that the distributed functional residues of monellin do not induce large local perturbations in the protein structure, eliminating the formation of folding intermediates and in turn domain-swapping. On the other hand, the localized protease-binding regions of stefin-B promote the formation of a folding intermediate which can lead to domain-swapping. Thus, domain-swapping can be a direct consequence of the constraints that function imposes on the protein structure.

Topological and sequence information predict that foldons organize a partially overlapped and hierarchical structure

It has been suggested that proteins have substructures, called foldons, which can cooperatively fold into the native structure. However, several prior investigations define foldons in various ways, citing different foldon characteristics, thereby making the concept of a foldon ambiguous. In this study, we perform a Gō model simulation and analyze the characteristics of substructures that cooperatively fold into the native-like structure. Although some results do not agree well with the experimental evidence due to the simplicity of our coarse-grained model, our results strongly suggest that cooperatively folding units sometimes organize a partially overlapped and hierarchical structure. This view makes us easy to interpret some different proposal about the foldon as a difference of the hierarchical structure. On the basis of this finding, we present a new method to assign foldons and their hierarchy, using structural and sequence information. The results show that the foldons assigned by our method correspond to the intermediate structures identified by some experimental techniques. The new method makes it easy to predict whether a protein folds sequentially into the native structure or whether some foldons fold into the native structure in parallel.

Individual and collective contributions of chaperoning and degradation to protein homeostasis in E. coli

The folding fate of a protein in vivo is determined by the interplay between a protein's folding energy landscape and the actions of the proteostasis network, including molecular chaperones and degradation enzymes. The mechanisms of individual components of the E. coli proteostasis network have been studied extensively, but much less is known about how they function as a system. We used an integrated experimental and computational approach to quantitatively analyze the folding outcomes (native folding versus aggregation versus degradation) of three test proteins biosynthesized in E. coli under a variety of conditions. Overexpression of the entire proteostasis network benefited all three test proteins, but the effect of upregulating individual chaperones or the major degradation enzyme, Lon, varied for proteins with different biophysical properties. In sum, the impact of the E. coli proteostasis network is a consequence of concerted action by the Hsp70 system (DnaK/DnaJ/GrpE), the Hsp60 system (GroEL/GroES), and Lon.

Energy landscapes of functional proteins are inherently risky

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

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

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

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.