Three-dimensional structure of human protein kinase C interacting protein 1, a member of the HIT family of proteins.

The three-dimensional structure of protein kinase C interacting protein 1 (PKCI-1) has been solved to high resolution by x-ray crystallography using single isomorphous replacement with anomalous scattering. The gene encoding human PKCI-1 was cloned from a cDNA library by using a partial sequence obtained from interactions identified in the yeast two-hybrid system between PKCI-1 and the regulatory domain of protein kinase C-beta. The PKCI-1 protein was expressed in Pichia pastoris as a dimer of two 13.7-kDa polypeptides. PKCI-1 is a member of the HIT family of proteins, shown by sequence identity to be conserved in a broad range of organisms including mycoplasma, plants, and humans. Despite the ubiquity of this protein sequence in nature, no distinct function has been shown for the protein product in vitro or in vivo. The PKCI-1 protomer has an alpha+beta meander fold containing a five-stranded antiparallel sheet and two helices. Two protomers come together to form a 10-stranded antiparallel sheet with extensive contacts between a helix and carboxy terminal amino acids of a protomer with the corresponding amino acids in the other protomer. PKCI-1 has been shown to interact specifically with zinc. The three-dimensional structure has been solved in the presence and absence of zinc and in two crystal forms. The structure of human PKCI-1 provides a model of this family of proteins which suggests a stable fold conserved throughout nature.

Peptides containing a consensus Ras binding sequence from Raf-1 and theGTPase activating protein NF1 inhibit Ras function.

A key event in Ras-mediated signal transduction and transformation involves Ras interaction with its downstream effector targets. Although substantial evidence has established that the Raf-1 serine/threonine kinase is a critical effector of Ras function, there is increasing evidence that Ras function is mediated through interaction with multiple effectors to trigger Raf-independent signaling pathways. In addition to the two Ras GTPase activating proteins (GAPs; p120- and NF1-GAP), other candidate effectors include activators of the Ras-related Ral proteins (RalGDS and RGL) and phosphatidylinositol 3-kinase. Interaction between Ras and its effectors requires an intact Ras effector domain and involves preferential recognition of active Ras-GTP. Surprisingly, these functionally diverse effectors lack significant sequence homology and no consensus Ras binding sequence has been described. We have now identified a consensus Ras binding sequence shared among a subset of Ras effectors. We have also shown that peptides containing this sequence from Raf-1 (RKTFLKLA) and NF1-GAP (RRFFLDIA) block NF1-GAP stimulation of Ras GTPase activity and Ras-mediated activation of mitogen-activated protein kinases. In summary, the identification of a consensus Ras-GTP binding sequence establishes a structural basis for the ability of diverse effector proteins to interact with Ras-GTP. Furthermore, our demonstration that peptides that contain Ras-GTP binding sequences can block Ras function provides a step toward the development of anti-Ras agents.

The human papillomavirus type 16 E7 protein binds human interferon regulatory factor-9 via a novel PEST domain required for transformation

It is critical that viruses are able to avoid the antiviral activities of interferon (IFN). We have shown previously that the human papillomavirus (HPV) is able to avoid IFN-alpha via interaction of the HPV-16 E7 protein with IFN regulatory factor-9 (IRF-9). Here, we investigated the details of the interaction using HPV-16 E7 peptide mapping to show that IRF-9 binds HPV-16 E7 in a domain encompassing amino acids 25-36. A closer examination of this region indicates this is a novel proline, glutamate, serine, and threonine-rich (PEST) domain, with a PEST score of 8.74. We have also mapped the region of interaction within IRF-9 and found that amino acids 354-393 play an important role in binding to HPV-16 E7. This region of IRF-9 encompasses the IRF association domain (IAD), a region important for protein-protein interaction central to IRF function. Finally, we used alanine-scanning mutagenesis to determine if E7-IRF-9 interaction was important for E7-mediated cellular transformation and found that the HPV-16 E7 mutants Y25A, E26A, S31A, S32A, and E35A, but not L28A and N29A, caused loss of transformation ability. Preliminary data suggest loss of IRF-9 interaction with E7 mutants correlated with transformation. Our work suggests E7- IRF- 9 interaction is important for the transforming ability of HPV-16 E7 and that HPV-16 E7 may interact with other IRF proteins that have IAD domains.

The β domain is required for Vps4p oligomerization into a functionally active ATPase

Endocytic and biosynthetic trafficking pathways to the lysosome/vacuole converge at the prevacuolar endosomal compartment. During transport through this compartment, integral membrane proteins that are destined for delivery to the lysosome/vacuole lumen undergo multivesicular body (MVB) sorting into internal vesicles formed by invagination of the endosomal limiting membrane. Vps4 is an AAA family ATPase which plays a key role in MVB sorting and facilitates transport through endosomes. It possesses an N-terminal microtubule interacting and trafficking domain required for recruitment to endosomes and an AAA domain with an ATPase catalytic site. The recently solved 3D structure revealed a P domain, which protrudes from the AAA domain, and a final C-terminal alpha-helix. However, the in vivo roles of these domains are not known. In this study, we have identified motifs in these domains that are highly conserved between yeast and human Vps4. We have mutated these motifs and studied the effect on yeast Vps4p function in vivo and in vitro. We show that the P domain of the budding yeast Vps4p is not required for recruitment to endosomes, but is essential for all Vps4p endocytic functions in vivo. We also show that the P domain is required for Vps4p homotypic interaction and for full ATPase activity. In addition, it is required for interaction with Vta1p, which works in concert with Vps4p in vivo. Our studies suggest that assembly of a Vps4p oligomeric complex with full ATPase activity that interacts with Vta1p is essential for normal endosome function.

The structural basis of DNA binding by the single-stranded DNA-binding protein from Sulfolobus solfataricus

Canonical single-stranded DNA-binding proteins (SSBs) from the oligosaccharide/oligonucleotide-binding (OB) domain family are present in all known organisms and are critical for DNA replication, recombination and repair. The SSB from the hyperthermophilic crenarchaeote Sulfolobus solfataricus (SsoSSB) has a ‘simple’ domain organization consisting of a single DNA-binding OB fold coupled to a flexible C-terminal tail, in contrast with other SSBs in this family that incorporate up to four OB domains. Despite the large differences in the domain organization within the SSB family, the structure of the OB domain is remarkably similar all cellular life forms. However, there are significant differences in the molecular mechanism of ssDNA binding. We have determined the structure of the SsoSSB OB domain bound to ssDNA by NMR spectroscopy. We reveal that ssDNA recognition is modulated by base-stacking of three key aromatic residues, in contrast with the OB domains of human RPA and the recently discovered human homologue of SsoSSB, hSSB1. We also demonstrate that SsoSSB binds ssDNA with a footprint of five bases and with a defined binding polarity. These data elucidate the structural basis of DNA binding and shed light on the molecular mechanism by which these ‘simple’ SSBs interact with ssDNA.

Elucidating the chromosome 9 association with AS; CARD9 is a candidate gene

Ankylosing spondylitis (AS) is polygenic with contributions from the immunologically relevant genes HLA-B27, ERAP1 and IL23R. A recent genome-wide association screen (GWAS) identified associations (P0.005) with the non-synonymous single-nucleotide polymorphisms (nsSNPs), rs4077515 and rs3812571, in caspase recruitment domain-containing protein 9 (CARD9) and small nuclear RNA-activating complex polypeptide 4 (SNAPC4) on chromosome 9q that had previously been linked to AS. We replicated these associations in a study of 730 AS patients compared with 2879 historic disease controls (rs4077515 P0.0004, odds ratio (OR)1.2, 95% confidence interval (CI)1.1-1.4; rs3812571 P0.0003, OR1.2, 95% CI1.1-1.4). Meta-analysis revealed strong associations of both SNPs with AS, rs4077515 P0.000005, OR1.2, 95% CI1.1-1.3 and rs3812571 P0.000006, OR1.2, 95% CI1.1-1.3. We then typed 1604 AS cases and 1020 controls for 13 tagging SNPs; 6 showed at least nominal association, 5 of which were in CARD9. We imputed genotypes for 13 additional SNPs but none was more strongly associated with AS than the tagging SNPs. Finally, interrogation of an mRNA expression database revealed that the SNPs most strongly associated with AS (or in strong linkage disequilibrium) were those most associated with CARD9 expression. CARD9 is a plausible candidate for AS given its central role in the innate immune response.

Expression of the extracellular domain of the human heat-stable enterotoxin receptor in Escherichia coli and generation of neutralizing antibodies

The entire extracellular domain of the human heat-stable enterotoxin (ST) receptor as well as a truncated N-terminal domain were cloned as glutathione S-transferase fusion proteins and expressed in Escherichia coli. The recombinant fusion proteins were purified from both the cytosol and the inclusion body fractions by selective detergent extraction followed by glutathione-agarose affinity chromatography. The purified protein, corresponding to the entire extracellular domain, bound the stable toxin peptide with an affinity comparable to that of the native receptor characterized from the human colonic T84 cell line. No binding was observed with the N-terminal truncated fragment of the receptor under similar conditions, Polyclonal antibodies were raised to the entire extracellular domain fusion protein as well as the truncated extracellular domain fusion protein, and the antibodies were purified by affinity chromatography. Addition of the purified antibodies to T84 cells inhibited ST binding and abolished ST-mediated cGMP production, indicating that critical epitopes involved in ligand interaction are present in the N-terminal fragment of the receptor, Purified antibodies recognized a single protein of M(r) 160,000 Da on Western blotting with T84 membranes, corresponding to a size of the native glycosylated receptor in T84 cells. These studies are the first report of the expression, purification, and characterization of any member of the guanylyl cyclase family of receptors in E. coli and show that binding of the toxin to the extracellular domain of the receptor is possible in the absence of any posttranslational modifications such as glycosylation. The recombinant fusion proteins as well as the antibodies that we have generated could serve as useful tools in the identification of critical residues of the extracellular domain involved in ligand interaction.

Comparative kinomics of human and chimpanzee reveal unique kinship and functional diversity generated by new domain combinations

Background: Phosphorylation by protein kinases is a common event in many cellular processes. Further, many kinases perform specialized roles and are regulated by non-kinase domains tethered to kinase domain. Perturbation in the regulation of kinases leads to malignancy. We have identified and analysed putative protein kinases encoded in the genome of chimpanzee which is a close evolutionary relative of human. Result: The shared core biology between chimpanzee and human is characterized by many orthologous protein kinases which are involved in conserved pathways. Domain architectures specific to chimp/human kinases have been observed. Chimp kinases with unique domain architectures are characterized by deletion of one or more non-kinase domains in the human kinases. Interestingly, counterparts of some of the multi-domain human kinases in chimp are characterized by identical domain architectures but with kinase-like non-kinase domain. Remarkably, out of 587 chimpanzee kinases no human orthologue with greater than 95% sequence identity could be identified for 160 kinases. Variations in chimpanzee kinases compared to human kinases are brought about also by differences in functions of domains tethered to the catalytic kinase domain. For example, the heterodimer forming PB1 domain related to the fold of ubiquitin/Ras-binding domain is seen uniquely tethered to PKC-like chimpanzee kinase. Conclusion: Though the chimpanzee and human are evolutionary very close, there are chimpanzee kinases with no close counterpart in the human suggesting differences in their functions. This analysis provides a direction for experimental analysis of human and chimpanzee protein kinases in order to enhance our understanding on their specific biological roles.

Structure-Function Studies of GDNF Family Ligand-RET signalling

Glial cell line-derived neurotrophic factor (GDNF) and its family members neurturin (NRTN), artemin (ARTN) and persephin (PSPN) are growth factors, which are involved in the development, differentiation and maintenance of many neuron types. In addition, they function outside of the nervous system, e.g. in the development of kidney, testis and liver. GDNF family ligand (GFL) signalling happens through a tetrameric receptor complex, which includes two glycosylphosphatidylinositol (GPI)-anchored GDNF family receptor (GFRα) molecules and two RET (rearranged during transfection) receptor tyrosine kinases. Each of the ligands binds preferentially one of the four GFRα receptors: GDNF binds to GFRα1, NRTN to GFRα2, ARTN to GFRα3 and PSPN to GFRα4. The signal is then delivered by RET, which cannot bind the GFLs on its own, but can bind the GFL-GFRα complex. Under normal cellular conditions, RET is only phosphorylated on the cell surface after ligand binding. At least the GDNF-GFRα1 complex is believed to recruit RET to lipid rafts, where downstream signalling occurs. In general, GFRαs consist of three cysteine-rich domains, but all GFRα4s except for chicken GFRα4 lack domain 1 (D1). We characterised the biochemical and cell biological properties of mouse PSPN receptor GFRα4 and showed that it has a significantly weaker capacity than GFRα1 to recruit RET to the lipid rafts. In spite of that, it can phosphorylate RET in the presence of PSPN and contribute to neuronal differentiation and survival. Therefore, the recruitment of RET to the lipid rafts does not seem to be crucial for the biological activity of all GFRα receptors. Secondly, we demonstrated that GFRα1 D1 stabilises the GDNF-GFRα1 complex and thus affects the phosphorylation of RET and contributes to the biological activity. This may be important in physiological conditions, where the concentration of the ligand or the soluble GFRα1 receptor is low. Our results also suggest a role for D1 in heparin binding and, consequently, in the biodistribution of released GFRα1 or in the formation of the GFL-GFRα-RET complex. We also presented the crystallographic structure of GDNF in the complex with GFRα1 domains 2 and 3. The structure differs from the previously published ARTN-GFRα3 structure in three significant ways. The biochemical data verify the structure and reveal residues participating in the interactions between GFRα1 and GDNF, and preliminarily also between GFRα1 and RET and heparin. Finally, we showed that, the precursor of the oncogenic MEN 2B (multiple endocrine neoplasia type 2) form of RET gets phosphorylated already during its synthesis in the endoplasmic reticulum (ER). We also demonstrated that it associates with Src homology 2 domain-containing protein (SHC) and growth factor receptor-bound protein (GRB2) in the ER, and has the capacity to activate several downstream signalling molecules.

Protein Kinase C and Anaplastic Lymphoma Kinase Targeted Compounds

The protein kinases (PKs) belong to the largest single family of enzymes, phosphotransferases, which catalyze the phosphorylation of other enzymes and proteins and function primarily in signal transduction. Consequently, PKs regulate cell mechanisms such as growth, differentiation, and proliferation. Dysfunction of these cellular mechanisms may lead to cancer, a major predicament in health care. Even though there is a range of clinically available cancer-fighting drugs, increasing number of cancer cases and setbacks such as drug resistance, constantly keep cancer research active. At the commencement of this study an isophthalic acid derivative had been suggested to bind to the regulatory domain of protein kinase C (PKC). In order to investigate the biological effects and structure-activity relationships (SARs) of this new chemical entity, a library of compounds was synthesized. The best compounds induced apoptosis in human leukemia HL-60 cells and were not cytotoxic in Swiss 3T3 fibroblasts. In addition, the best apoptosis inducers were neither cytotoxic nor mutagenic. Furthermore, results from binding affinity assays of PKC isoforms revealed the pharmacophores of these isophthalic acid derivatives. The best inhibition constants of the tested compounds were measured to 210 nM for PKCα and to 530 nM for PKCδ. Among natural compounds targeting the regulatory domain of PKC, the target of bistramide A has been a matter of debate. It was initially found to activate PKCδ; however, actin was recently reported as the main target. In order to clarify and to further study the biological effects of bistramide A, the total syntheses of the natural compound and two isomers were performed. Biological assays of the compounds revealed accumulation of 4n polyploid cells as the primary mode of action and the compounds showed similar overall antiproliferative activities. However, each compound showed a distinct distribution of antimitotic effect presumably via actin binding, proapoptotic effect presumably via PKCδ, and pro-differentiation effect as evidenced by CD11b expression. Furthermore, it was shown that the antimitotic and proapoptotic effects of bistramide A were not secondary effects of actin binding but independent effects. The third aim in this study was to synthesize a library of a new class of urea-based type II inhibitors targeted at the kinase domain of anaplastic lymphoma kinase (ALK). The best compounds in this library showed IC50 values as low as 390 nM for ALK while the initial low cellular activities were successfully increased even by more than 70 times for NPM-ALK- positive BaF3 cells. More importantly, selective antiproliferative activity on ALK-positive cell lines was achieved; while the best compound affected the BaF3 and SU-DHL-1 cells with IC50 values of 0.5 and 0.8 μM, respectively, they were less toxic to the NPM-ALK-negative human leukemic cells U937 (IC50 = 3.2 μM) and BaF3 parental cells (IC50 = 5.4 μM). Furthermore, SAR studies of the synthesized compounds revealed functional groups and positions of the scaffold, which enhanced the enzymatic and cellular activities.

Generation of Flat and Curved Membranes by Inverse-BAR Domain Proteins

Biological membranes are tightly linked to the evolution of life, because they provide a way to concentrate molecules into partially closed compartments. The dynamic shaping of cellular membranes is essential for many physiological processes, including cell morphogenesis, motility, cytokinesis, endocytosis, and secretion. It is therefore essential to understand the structure of the membrane and recognize the players that directly sculpt the membrane and enable it to adopt different shapes. The actin cytoskeleton provides the force to push eukaryotic plasma membrane in order to form different protrusions or/and invaginations. It has now became evident that actin directly co-operates with many membrane sculptors, including BAR domain proteins, in these important events. However, the molecular mechanisms behind BAR domain function and the differences between the members of this large protein family remain largely unresolved. In this thesis, the structure and functions of the I-BAR domain family members IRSp53 and MIM were thoroughly analyzed. By using several methods such as electron microscopy and systematic mutagenesis, we showed that these I-BAR domain proteins bind to PI(4,5)P2-rich membranes, generate negative membrane curvature and are involved in the formation of plasma membrane protrusions in cells e.g. filopodia. Importantly, we characterized a novel member of the BAR-domain superfamily which we named Pinkbar. We revealed that Pinkbar is specifically expressed in kidney and epithelial cells, and it localizes to Rab13-positive vesicles in intestinal epithelial cells. Remarkably, we learned that the I-BAR domain of Pinkbar does not generate membrane curvature but instead stabilizes planar membranes. Based on structural, mutagenesis and biochemical work we present a model for the mechanism of the novel membrane deforming activity of Pinkbar. Collectively, this work describes the mechanism by which I-BAR domain proteins deform membranes and provides new information about the biological roles of these proteins. Intriguingly, this work also gives evidence that significant functional plasticity exists within the I-BAR domain family. I-BAR proteins can either generate negative membrane curvature or stabilize planar membrane sheets, depending on the specific structural properties of their I-BAR domains. The results presented in this thesis expand our knowledge on membrane sculpting mechanisms and shows for the first time how flat membranes can be generated in cells.

Tethering preferences of domain families co-occurring in multi-domain proteins

Genomic data of several organisms have revealed the presence of a vast repertoire of multi-domain proteins. The role played by individual domains in a multi-domain protein has a profound influence on the overall function of the protein. In the present analysis an attempt has been made to better understand the tethering preferences of domain families that occur in multi-domain proteins. The analysis has been carried out on an exhaustive dataset of 2 961 898 sequences of proteins from 930 organisms, where 741 274 proteins are comprised of at least two domain families. For every domain family, the number of other domain families with which it co-occurs within a protein in this dataset has been enumerated and is referred to as the tethering number of the domain family. It was found that, in the general dataset, the AAA ATPase family and the family of Ser/Thr kinases have the highest tethering numbers of 450 and 444 respectively. Further analysis reveals significant correlation between the number of members in a family and its tethering number. Positive correlation was also observed for the extent of a sequence and functional diversity within a family and the tethering numbers of domain families. Domain families that are present ubiquitously in diverse organisms tend to have large tethering numbers, while organism/kingdom-specific families have low tethering numbers. Thus, the analysis uncovers how domain families recombine and evolve to give rise to multi-domain proteins.

CLAP: A web-server for automatic classification of proteins with special reference to multi-domain proteins

Background: The function of a protein can be deciphered with higher accuracy from its structure than from its amino acid sequence. Due to the huge gap in the available protein sequence and structural space, tools that can generate functionally homogeneous clusters using only the sequence information, hold great importance. For this, traditional alignment-based tools work well in most cases and clustering is performed on the basis of sequence similarity. But, in the case of multi-domain proteins, the alignment quality might be poor due to varied lengths of the proteins, domain shuffling or circular permutations. Multi-domain proteins are ubiquitous in nature, hence alignment-free tools, which overcome the shortcomings of alignment-based protein comparison methods, are required. Further, existing tools classify proteins using only domain-level information and hence miss out on the information encoded in the tethered regions or accessory domains. Our method, on the other hand, takes into account the full-length sequence of a protein, consolidating the complete sequence information to understand a given protein better. Results: Our web-server, CLAP (Classification of Proteins), is one such alignment-free software for automatic classification of protein sequences. It utilizes a pattern-matching algorithm that assigns local matching scores (LMS) to residues that are a part of the matched patterns between two sequences being compared. CLAP works on full-length sequences and does not require prior domain definitions. Pilot studies undertaken previously on protein kinases and immunoglobulins have shown that CLAP yields clusters, which have high functional and domain architectural similarity. Moreover, parsing at a statistically determined cut-off resulted in clusters that corroborated with the sub-family level classification of that particular domain family. Conclusions: CLAP is a useful protein-clustering tool, independent of domain assignment, domain order, sequence length and domain diversity. Our method can be used for any set of protein sequences, yielding functionally relevant clusters with high domain architectural homogeneity. The CLAP web server is freely available for academic use at http://nslab.mbu.iisc.ernet.in/clap/.

De-DUFing the DUFs: Deciphering distant evolutionary relationships of Domains of Unknown Function using sensitive homology detection methods

Background: In the post-genomic era where sequences are being determined at a rapid rate, we are highly reliant on computational methods for their tentative biochemical characterization. The Pfam database currently contains 3,786 families corresponding to ``Domains of Unknown Function'' (DUF) or ``Uncharacterized Protein Family'' (UPF), of which 3,087 families have no reported three-dimensional structure, constituting almost one-fourth of the known protein families in search for both structure and function. Results: We applied a `computational structural genomics' approach using five state-of-the-art remote similarity detection methods to detect the relationship between uncharacterized DUFs and domain families of known structures. The association with a structural domain family could serve as a start point in elucidating the function of a DUF. Amongst these five methods, searches in SCOP-NrichD database have been applied for the first time. Predictions were classified into high, medium and low-confidence based on the consensus of results from various approaches and also annotated with enzyme and Gene ontology terms. 614 uncharacterized DUFs could be associated with a known structural domain, of which high confidence predictions, involving at least four methods, were made for 54 families. These structure-function relationships for the 614 DUF families can be accessed on-line at http://proline.biochem.iisc.ernet.in/RHD_DUFS/. For potential enzymes in this set, we assessed their compatibility with the associated fold and performed detailed structural and functional annotation by examining alignments and extent of conservation of functional residues. Detailed discussion is provided for interesting assignments for DUF3050, DUF1636, DUF1572, DUF2092 and DUF659. Conclusions: This study provides insights into the structure and potential function for nearly 20 % of the DUFs. Use of different computational approaches enables us to reliably recognize distant relationships, especially when they converge to a common assignment because the methods are often complementary. We observe that while pointers to the structural domain can offer the right clues to the function of a protein, recognition of its precise functional role is still `non-trivial' with many DUF domains conserving only some of the critical residues. It is not clear whether these are functional vestiges or instances involving alternate substrates and interacting partners. Reviewers: This article was reviewed by Drs Eugene Koonin, Frank Eisenhaber and Srikrishna Subramanian.

Dietary animal and plant protein intakes and their associations with obesity and cardio-metabolic indicators in European adolescents: the HELENA cross-sectional study

Background: Previous studies suggest that dietary protein might play a beneficial role in combating obesity and its related chronic diseases. Total, animal and plant protein intakes and their associations with anthropometry and serum biomarkers in European adolescents using one standardised methodology across European countries are not well documented. Objectives: To evaluate total, animal and plant protein intakes in European adolescents stratified by gender and age, and to investigate their associations with cardio-metabolic indicators (anthropometry and biomarkers). Methods: The current analysis included 1804 randomly selected adolescents participating in the HELENA study (conducted in 2006-2007) aged 12.5-17.5 y (47% males) who completed two non-consecutive computerised 24-h dietary recalls. Associations between animal and plant protein intakes, and anthropometry and serum biomarkers were examined with General linear Model multivariate analysis. Results: Average total protein intake exceeded the recommendations of World Health Organization and European Food Safety Authority. Mean total protein intake was 96 g/d (59% derived from animal protein). Total, animal and plant protein intakes (g/d) were significantly lower in females than in males and total and plant protein intakes were lower in younger participants (12.5-14.9 y). Protein intake was significantly lower in underweight subjects and higher in obese ones; the direction of the relationship was reversed after adjustments for body weight (g/(kg.d)). The inverse association of plant protein intakes was stronger with BMI z-score and body fat percentage (BF%) compared to animal protein intakes. Additionally, BMI and BF% were positively associated with energy percentage of animal protein. Conclusions: This sample of European adolescents appeared to have adequate total protein intake. Our findings suggest that plant protein intakes may play a role in preventing obesity among European adolescents. Further longitudinal studies are needed to investigate the potential beneficial effects observed in this study in the prevention of obesity and related chronic diseases.