The effect of context on the folding of β-hairpins

Small β-hairpin peptides have been widely used as models for the folding of β-sheets. But how applicable is the folding of such models to β-structure in larger proteins with conventional hydrophobic cores? Here we present multiple unfolding simulations of three such proteins that contain the WW domain double hairpin β-sheet motif: cold shock protein A (CspA), cold shock protein B (CspB) and glucose permease IIA domain. We compare the behavior of the free motif in solution and in the context of proteins of different size and architecture. Both Csp proteins lost contacts between the double-hairpin motif and the protein core as the first step of unfolding and proceeded to unfold with loss of the third β-strand, similar to the isolated WW domain. The glucose permease IIA domain is a larger protein and the contacts between the motif and the core were not lost as quickly. Instead the unfolding pathway of glucose permease IIA followed a different pathway with β1 pulling away from the sheet first. Interestingly, when the double hairpin motif was excised from the glucose permease IIA domain and simulated in isolation in water it unfolded by the same pathway as the WW domain, indicating that it is tertiary interactions with the protein that alter the motif’s unfolding not a sequence dependent effect on its intrinsic unfolding behavior. With respect to the unfolding of the hairpins, there was no consistent order to the loss of hydrogen bonds between the β-strands in the hairpins in any of the systems. Our results show that while the folding behavior of the isolated WW domain is generally consistent with the double hairpin motif’s behavior in the cold shock proteins, it is not the case for the glucose permease IIA domain. So, one must be cautious in extrapolating findings from model systems to larger more complicated proteins where tertiary interactions can overwhelm intrinsic behavior.     

Peptides and peptidomimetics as inhibitors of protein–protein interactions involving β-sheet secondary structures

1. Download : Download high-res image (214KB)
2. Download : Download full-size image

     

Hydroxylysine in the N-terminal regions of the α1- and α2-chains of various collagens

The degree of hydroxylation of the lysine residue located in both alpha(1)- and alpha(2)-chains of collagen in the N-terminal, non-helical telopeptide region of the molecule has been determined in collagen from various sources after isolation of the peptides (alpha(1)- and alpha(2)-CB1) that contain the lysine residue in question and are obtained by cyanogen bromide cleavage of collagen alpha(1)- and alpha(2)-chains respectively. As with collagen from chick tibia, bone collagens from rat tibia and femur and embryonic chick frontal bone, have a high degree of hydroxylation (approx. 50% or more) of the lysine residue in both alpha(1)- and alpha(2)-CB1 peptides. This is in contrast with the lack of hydroxylation of this residue in both alpha(1)- and alpha(2)-chains of all skin collagens so far examined. The presence of hydroxylysine in alpha(1)- and alpha(2)-CB1 peptides from tendon collagen is also indicated. In rat tail tendon collagen the amount of hydroxylation is only slight but in the much less soluble tendon collagen from embryonic chick leg tendons, approximately one-third of the lysine is hydroxylated.     

Solution of Levinthal’s paradox is possible at the level of the formation and assembly of protein secondary structures

The complete volume of a protein’s conformation space is smaller by many orders of magnitude at the level of secondary-structure elements as compared with the conformation of amino-acid residues. According to Levinthal’s estimate, the latter is ~10^2L, with L being the number of residues in the chain, while the former, at the level of secondary structures, increases no faster than ~L^N with N being the number of the secondary-structure elements. N is approximately L/15 according to the statistics of protein structures. This drastic decrease in the exponent (L/15 instead of 2L) considerably reduces the sampling space and explains the reason that sampling of the conformation space at the level of secondary-structure elements does not prevent the protein chain from finding its most stable structure.     

Role of peroxisomes in the biosynthesis and secretion of β-lactams and other secondary metabolites

Peroxisomes are eukaryotic organelles surrounded by a single bilayer membrane, containing a variety of proteins depending on the organism; they mainly perform degradation reactions of toxic metabolites (detoxification), catabolism of linear and branched-chain fatty acids, and removal of H₂O₂ (formed in some oxidative processes) by catalase. Proteins named peroxins are involved in recruiting, transporting, and introducing the peroxisomal matrix proteins into the peroxisomes. The matrix proteins contain the peroxisomal targeting signals PTS1 and/or PTS2 that are recognized by the peroxins Pex5 and Pex7, respectively. Initial evidence indicated that the penicillin biosynthetic enzyme isopenicillin N acyltransferase (IAT) of Penicillium chrysogenum is located inside peroxisomes. There is now solid evidence (based on electron microscopy and/or biochemical data) confirming that IAT and the phenylacetic acid- and fatty acid-activating enzymes are also located in peroxisomes. Similarly, the Acremonium chrysogenum CefD1 and CefD2 proteins that perform the central reactions (activation and epimerization of isopenicillin N) of the cephalosporin pathway are targeted to peroxisomes. Growing evidence supports the conclusion that some enzymes involved in the biosynthesis of mycotoxins (e.g., AK-toxin), and the biosynthesis of signaling molecules in plants (e.g., jasmonic acid or auxins) occur in peroxisomes. The high concentration of substrates (in many cases toxic to the cytoplasm) and enzymes inside the peroxisomes allows efficient synthesis of metabolites with interesting biological or pharmacological activities. This compartmentalization poses additional challenges to the cell due to the need to import the substrates into the peroxisomes and to export the final products; the transporters involved in these processes are still very poorly known. This article focuses on new aspects of the metabolic processes occurring in peroxisomes, namely the degradation and detoxification processes that lead to the biosynthesis and secretion of secondary metabolites.     

Analysis of secondary structural and physicochemical changes in protein–protein complexes

Conformation switching in protein–protein complexes is considered important for the molecular recognition process. Overall analysis of 123 protein–protein complexes in a benchmark data-set showed that 6.8% of residues switched over their secondary structure conformation upon complex formation. Amino acid residue-wise preference for conformation change has been analyzed in binding and non-binding site residues separately. In this analysis, residues such as Ser, Leu, Glu, and Lys had higher frequency of secondary structural conformation change. The change of helix to coil and sheet to coil conformation and vice versa has been observed frequently, whereas the conformation change of helix to extended sheet occurred rarely in the studied complexes. Influence of conformation change toward the N and C terminal on either side of the binding site residues has been analyzed. Further, analysis on φ and ψ angle variation, conservation, stability, and solvent accessibility have been performed on binding site residues. Knowledge obtained from the present study could be effectively employed in the protein–protein modeling and docking studies.     

How well are protein structures annotated in secondary databases?

We investigated to what extent Protein Data Bank (PDB) entries are annotated with second-party information based on existing cross-references between PDB and 15 other databases. We report 2 interesting findings. First, there is a clear "annotation gap" for structures less than 7 years old for secondary databases that are manually curated. Second, the examined databases overlap with each other quite well, dividing the PDB into 2 well-annotated thirds and one poorly annotated third. Both observations should be taken into account in any study depending on the selection of protein structures by their annotation.     

Steered molecular dynamics simulation of the binding of the β2 and β3 regions in domain-swapped human cystatin C dimer

The crystal structure of the human cystatin C (hCC) dimer revealed that a stable twofold-symmetric dimer was formed via 3D domain swapping. Domain swapping with the need for near-complete unfolding has been proposed as a possible route for amyloid fibril initiation. Thus, the interesting interactions that occur between the two molecules may be important for the further aggregation of the protein. In this work, we performed steered molecular dynamics (SMD) simulations to investigate the dissociation of the β2 and β3 strands in the hCC dimer. The energy changes observed during the SMD simulations showed that electrostatic interactions were the dominant interactions involved in stabilizing the two parts of the dimer during the early stages of SMD simulation, whereas van der Waals (VDW) interactions and electrostatic interactions were equally matched during the latter stages. Furthermore, our data indicated that the two parts of the dimer are stabilized by intermolecular hydrogen bonds among the residues Arg51 (β2), Gln48 (β2), Asp65 (β3), and Glu67 (β3), salt bridges among the residues Arg53 (β2), Arg51 (β2), and Asp65 (β3), and VDW interactions among the residues Gln48 (β2), Arg51 (β2), Glu67 (β3), Asp65 (β3), Phe63 (β3), and Asn61 (β3). The residues Gln48 (β2), Arg51 (β2), Asp65 (β3) and Glu67 (β3) appear to be crucial, as they play important roles in both electrostatic and VDW interactions. Thus, the present study determined the key residues involved in the stabilization of the domain-swapped dimer structure, and also provided molecular-level insights into the dissociation process of the hCC dimer.     

Analyzing the pathways enriched in genes associated with nicotine dependence in the context of human protein–protein interaction network

Nicotine dependence is the primary addictive stage of cigarette smoking. Although a lot of studies have been performed to explore the molecular mechanism underlying nicotine dependence, our understanding on this disorder is still far from complete. Over the past decades, an increasing number of candidate genes involved in nicotine dependence have been identified by different technical approaches, including the genetic association analysis. In this study, we performed a comprehensive collection of candidate genes reported to be genetically associated with nicotine dependence. Then, the biochemical pathways enriched in these genes were identified by considering the gene’s propensity to be related to nicotine dependence. One of the most widely used pathway enrichment analysis approach, over-representation analysis, ignores the function non-equivalence of genes in candidate gene set and may have low discriminative power in identifying some dysfunctional pathways. To overcome such drawbacks, we constructed a comprehensive human protein–protein interaction network, and then assigned a function weighting score to each candidate gene based on their network topological features. Evaluation indicated the function weighting score scheme was consistent with available evidence. Finally, the function weighting scores of the candidate genes were incorporated into pathway analysis to identify the dysfunctional pathways involved in nicotine dependence, and the interactions between pathways was detected by pathway crosstalk analysis. Compared to conventional over-representation-based pathway analysis tool, the modified method exhibited improved discriminative power and detected some novel pathways potentially underlying nicotine dependence. In summary, we conducted a comprehensive collection of genes associated with nicotine dependence and then detected the biochemical pathways enriched in these genes using a modified pathway enrichment analysis approach with function weighting score of candidate genes integrated. Our results may provide insight into the molecular mechanism underlying nicotine dependence.     

The genetic defect in the various types of human β-galactosidase deficiency

Summary Gene localization studies revealed the presence of two structural -galactosidase (GAL) loci on the human chromosomes 3 and 22 (de Wit et al., 1979). To determine the function of these genes, proliferating hybrid cell lines were isolated following fusion of fibroblasts from two different patients with a GAL deficiency and Chinese hamster cells. The hybrids were analyzed electrophoretically and immunologically.Fibroblasts from a patient with an adult type of GAL deficiency associated with a neuraminidase deficiency were used for the first fusion. No evidence for a structural GAL mutation was found in these hybrids. The absence of a structural GAL mutation is consistent with a primary defect in neuraminidase in this adult patient.Fibroblasts from a patient with the infantile type 1 G_M1-gangliosidosis were used for the second fusion. It is concluded that the human determinants present in the isolated hybrid lines occur in heteropolymeric man-Chinese hamster molecules. The heteropolymeric isoenzyme in (+3–22) hybrids is very labile and is sensitive to neuraminidase treatment. Therefore it is concluded that the infantile type 1 patient is mutated in the structural GAL gene on chromosome 3. Because this patient has a primary defect in G_M1-GAL, the GAL gene on chromosome 3 is apparently a G _M1-GAL gene. Interaction of the two GAL loci results in an additional band of GAL activity on electrophoresis. This suggests that the gene on chromosome 22 is also a structural G _M1-GAL gene.     

Intracellular trafficking of the β-secretase and processing of amyloid precursor protein

The main component of the amyloid plaques found in the brains of those with Alzheimer's disease (AD) is a polymerized form of the β-amyloid peptide (Aβ) and is considered to play a central role in the pathogenesis of this neurodegenerative disorder. Aβ is derived from the proteolytic processing of the amyloid precursor protein (APP). Beta site APP-cleaving enzyme, BACE1 (also known as β-secretase) is a membrane-bound aspartyl protease responsible for the initial step in the generation of Aβ peptide and is thus a prime target for therapeutic intervention. Substantive evidence now indicates that the processing of APP by BACE1 is regulated by the intracellular sorting of the enzyme and, moreover, perturbations in these intracellular trafficking pathways have been linked to late-onset AD. In this review, we highlight the recent advances in the understanding of the regulation of the intracellular sorting of BACE1 and APP and illustrate why the trafficking of these cargos represent a key issue for understanding the membrane-mediated events associated with the generation of the neurotoxic Aβ products in AD.     

Structural and functional analysis of the flexible regions of the catalytic α-subunit of protein kinase CK2

CK2 is a Ser/Thr protein kinase essential for cell viability. Its activity is anomalously high in several solid (prostate, mammary gland, lung, kidney and head and neck) and haematological tumours (AML, CML and PML), creating conditions favouring the onset of cancer. Cancer cells become addicted to high levels of CK2 activity and therefore this kinase is a remarkable example of "non-oncogene addiction". CK2 is a validated target for cancer therapy with one inhibitor in phase I clinical trials. Several crystal structures of CK2 are available, many in complex with ATP-competitive inhibitors, showing the presence of regions with remarkable flexibility. We present the structural characterisation of these regions by means of seven new crystal structures, in the apo form and in complex with inhibitors. We confirm previous findings about the unique flexibility of the CK2α catalytic subunit in the hinge/αD region, the p-loop and the β4β5 loop, and show here that there is no clear-cut correlation between the conformations of these flexible zones. Our findings challenge some of the current interpretations on the functional role of these regions and dispute the hypothesis that small ligands stabilize an inactive state. The mobility of the hinge/αD region in the human enzyme is unique among protein kinases, and this can be exploited for the development of more selective ATP-competitive inhibitors. The identification of different ligand binding modes to a secondary site can provide hints for the design of non-ATP-competitive inhibitors targeting the interaction between the α catalytic and the β regulatory subunits.     

Synthesis of the tri- and tetra-saccharides related to the fine structures of lichenan and cereal β-d-glucans

Syntheses, based on silver trifluoromethanesulfonate-promoted Koenigs-Knorr type condensations, are described of the d-glucotrioses, β-d-Glcp-(1→3)-β-d-Glcp-(1→4)-d-Glcp and β-d-Glcp-(1→4)-β-d-Glcp-(1→3)-d-Glcp, and the d-Glucotetraoses, β-d-Glcp-(1→3)-β-d-Glcp-(1→4)-β-d-Glcp-(1→4)-d-Glcp, β-d-Glcp-(1→4)-β-d-Glcp-(1→3)-β-d-Glcp-(1→4)-d-Glcp, and β-d-Glcp-(1→4)-β-d-Glcp-(1→4)-β-d-Glcp-(1→3)-d-Glcp, corresponding to the tri- and tetra-saccharide units in the linear chains of (1→4)- and (1→3)-linked β-d-glucopyranosyl residues of lichenan, and of oat and barley β-d-glucans.     

Computational design and experimental testing of the fastest-folding β-sheet protein

One of the most important and elusive goals of molecular biology is the formulation of a detailed, atomic-level understanding of the process of protein folding. Fast-folding proteins with low free-energy barriers have proved to be particularly productive objects of investigation in this context, but the design of fast-folding proteins was previously driven largely by experiment. Dramatic advances in the attainable length of molecular dynamics simulations have allowed us to characterize in atomic-level detail the folding mechanism of the fast-folding all-β WW domain FiP35. In the work reported here, we applied the biophysical insights gained from these studies to computationally design an even faster-folding variant of FiP35 containing only naturally occurring amino acids. The increased stability and high folding rate predicted by our simulations were subsequently validated by temperature-jump experiments. The experimentally measured folding time was 4.3 μs at 80 °C—about three times faster than the fastest previously known protein with β-sheet content and in good agreement with our prediction. These results provide a compelling demonstration of the potential utility of very long molecular dynamics simulations in redesigning proteins well beyond their evolved stability and folding speed.     

Chronic effects of β2 agonists on body composition and protein synthesis in the rat

Chronic treatment of rats with the ₂-adrenergic agonists clenbuterol and fenoterol over 16–19 d raised energy intake, expenditure, and body weight gain but did not affect fat or energy deposition, and body protein gain was increased by 50 and 18%, respectively. Both drugs increased the protein content and mitochondrial GDP-binding capacity of brown adipose tissue. Clenbuterol did not affect plasma insulin, growth hormone, or triiodothyronine levels, although insulin levels were reduced by fenoterol. Both drugs caused hypertrophy of skeletal muscle (gastrocnemius), and muscle protein synthesis in vivo (fractional rate) was elevated by 34 and 26% in clenbuterol and fenoteroltreated rats, respectively.