Linear B-Cell Epitope Prediction for In Silico Vaccine Design: A Performance Review of Methods Available via Command-Line Interface.

Linear B-cell epitope prediction research has received a steadily growing interest ever since the first method was developed in 1981. B-cell epitope identification with the help of an accurate prediction method can lead to an overall faster and cheaper vaccine design process, a crucial necessity in the COVID-19 era. Consequently, several B-cell epitope prediction methods have been developed over the past few decades, but without significant success. In this study, we review the current performance and methodology of some of the most widely used linear B-cell epitope predictors which are available via a command-line interface, namely, BcePred, BepiPred, ABCpred, COBEpro, SVMTriP, LBtope, and LBEEP. Additionally, we attempted to remedy performance issues of the individual methods by developing a consensus classifier, which combines the separate predictions of these methods into a single output, accelerating the epitope-based vaccine design. While the method comparison was performed with some necessary caveats and individual methods might perform much better for specialized datasets, we hope that this update in performance can aid researchers towards the choice of a predictor, for the development of biomedical applications such as designed vaccines, diagnostic kits, immunotherapeutics, immunodiagnostic tests, antibody production, and disease diagnosis and therapy.

ANTISOMA: A Computational Pipeline for the Reduction of the Aggregation Propensity of Monoclonal Antibodies.

Monoclonal antibodies (mAbs) constitute a promising class of therapeutics, since ca. 25% of all biotech drugs in development are mAbs. Even though their therapeutic value is now well established, human- and murine-derived mAbs do have deficiencies, such as short in vivo lifespan and low stability. However, the most difficult obstacle to overcome, toward the exploitation of mAbs for disease treatment, is the prevention of the formation of protein aggregates. ANTISOMA is a pipeline for the reduction of the aggregation tendency of mAbs through the decrease in their intrinsic aggregation propensity, based on an automated amino acid substitution approach. The method takes into consideration the special features of mAbs and aims at proposing specific point mutations that could lead to the redesign of those promising therapeutics, without affecting their epitope-binding ability. The method is available online at http://bioinformatics.biol.uoa.gr/ANTISOMA .

Revisiting Histone Deacetylases in Human Tumorigenesis: The Paradigm of Urothelial Bladder Cancer.

Urinary bladder cancer is a common malignancy, being characterized by substantial patient mortality and management cost. Its high somatic-mutation frequency and molecular heterogeneity usually renders tumors refractory to the applied regimens. Hitherto, methotrexate-vinblastine-adriamycin-cisplatin and gemcitabine-cisplatin represent the backbone of systemic chemotherapy. However, despite the initial chemosensitivity, the majority of treated patients will eventually develop chemoresistance, which severely reduces their survival expectancy. Since chromatin regulation genes are more frequently mutated in muscle-invasive bladder cancer, as compared to other epithelial tumors, targeted therapies against chromatin aberrations in chemoresistant clones may prove beneficial for the disease. "Acetyl-chromatin" homeostasis is regulated by the opposing functions of histone acetyltransferases (HATs) and histone deacetylases (HDACs). The HDAC/SIRT (super-)family contains 18 members, which are divided in five classes, with each family member being differentially expressed in normal urinary bladder tissues. Since a strong association between irregular HDAC expression/activity and tumorigenesis has been previously demonstrated, we herein attempt to review the accumulated published evidences that implicate HDACs/SIRTs as critical regulators in urothelial bladder cancer. Moreover, the most extensively investigated HDAC inhibitors (HDACis) are also analyzed, and the respective clinical trials are also described. Interestingly, it seems that HDACis should be preferably used in drug-combination therapeutic schemes, including radiation.

Gene-Specific Intron Retention Serves as Molecular Signature that Distinguishes Melanoma from Non-Melanoma Cancer Cells in Greek Patients.

BACKGROUND: Skin cancer represents the most common human malignancy, and it includes BCC, SCC, and melanoma. Since melanoma is one of the most aggressive types of cancer, we have herein attempted to develop a gene-specific intron retention signature that can distinguish BCC and SCC from melanoma biopsy tumors. METHODS: Intron retention events were examined through RT-sqPCR protocols, using total RNA preparations derived from BCC, SCC, and melanoma Greek biopsy specimens. Intron-hosted miRNA species and their target transcripts were predicted via the miRbase and miRDB bioinformatics platforms, respectively. Ιntronic ORFs were recognized through the ORF Finder application. Generation and visualization of protein interactomes were achieved by the IntAct and Cytoscape softwares, while tertiary protein structures were produced by using the I-TASSER online server. RESULTS: c-MYC and Sestrin-1 genes proved to undergo intron retention specifically in melanoma. Interaction maps of proteins encoded by genes being potentially targeted by retained intron-accommodated miRNAs were generated and SRPX2 was additionally delivered to our melanoma-specific signature. Novel ORFs were identified in MCT4 and Sestrin-1 introns, with potentially critical roles in melanoma development. CONCLUSIONS: The property of c-MYC, Sestrin-1, and SRPX2 genes to retain specific introns could be clinically used to molecularly differentiate non-melanoma from melanoma tumors.

MBPpred: Proteome-wide detection of membrane lipid-binding proteins using profile Hidden Markov Models.

A large number of modular domains that exhibit specific lipid binding properties are present in many membrane proteins involved in trafficking and signal transduction. These domains are present in either eukaryotic peripheral membrane or transmembrane proteins and are responsible for the non-covalent interactions of these proteins with membrane lipids. Here we report a profile Hidden Markov Model based method capable of detecting Membrane Binding Proteins (MBPs) from information encoded in their amino acid sequence, called MBPpred. The method identifies MBPs that contain one or more of the Membrane Binding Domains (MBDs) that have been described to date, and further classifies these proteins based on their position in respect to the membrane, either as peripheral or transmembrane. MBPpred is available online at http://bioinformatics.biol.uoa.gr/MBPpred. This method was applied in selected eukaryotic proteomes, in order to examine the characteristics they exhibit in various eukaryotic kingdoms and phyla.

Comparative genomics of the dairy isolate Streptococcus macedonicus ACA-DC 198 against related members of the Streptococcus bovis/Streptococcus equinus complex.

BACKGROUND: Within the genus Streptococcus, only Streptococcus thermophilus is used as a starter culture in food fermentations. Streptococcus macedonicus though, which belongs to the Streptococcus bovis/Streptococcus equinus complex (SBSEC), is also frequently isolated from fermented foods mainly of dairy origin. Members of the SBSEC have been implicated in human endocarditis and colon cancer. Here we compare the genome sequence of the dairy isolate S. macedonicus ACA-DC 198 to the other SBSEC genomes in order to assess in silico its potential adaptation to milk and its pathogenicity status. RESULTS: Despite the fact that the SBSEC species were found tightly related based on whole genome phylogeny of streptococci, two distinct patterns of evolution were identified among them. Streptococcus macedonicus, Streptococcus infantarius CJ18 and Streptococcus pasteurianus ATCC 43144 seem to have undergone reductive evolution resulting in significantly diminished genome sizes and increased percentages of potential pseudogenes when compared to Streptococcus gallolyticus subsp. gallolyticus. In addition, the three species seem to have lost genes for catabolizing complex plant carbohydrates and for detoxifying toxic substances previously linked to the ability of S. gallolyticus to survive in the rumen. Analysis of the S. macedonicus genome revealed features that could support adaptation to milk, including an extra gene cluster for lactose and galactose metabolism, a proteolytic system for casein hydrolysis, auxotrophy for several vitamins, an increased ability to resist bacteriophages and horizontal gene transfer events with the dairy Lactococcus lactis and S. thermophilus as potential donors. In addition, S. macedonicus lacks several pathogenicity-related genes found in S. gallolyticus. For example, S. macedonicus has retained only one (i.e. the pil3) of the three pilus gene clusters which may mediate the binding of S. gallolyticus to the extracellular matrix. Unexpectedly, similar findings were obtained not only for the dairy S. infantarius CJ18, but also for the blood isolate S. pasteurianus ATCC 43144. CONCLUSIONS: Our whole genome analyses suggest traits of adaptation of S. macedonicus to the nutrient-rich dairy environment. During this process the bacterium gained genes presumably important for this new ecological niche. Finally, S. macedonicus carries a reduced number of putative SBSEC virulence factors, which suggests a diminished pathogenic potential.

Beneficial Effects of Sideritis clandestina Extracts and Sideridiol against Amyloid ß Toxicity.

Alzheimer's disease (AD) is the most common form of dementia. Given the link between oxidative stress and AD, many studies focus on the identification of natural antioxidants against AD. Although their antioxidant capacity is important, increasing data suggest that additional activities are related to their beneficial effects, including properties against amyloid beta (Aß) aggregation. Sideritis spp. (mountain tea) extracts possess not only antioxidant activity but also other bioactivities that confer neuroprotection. Although various Sideritis spp. extracts have been extensively studied, there are scarce data on S. clandestina subsp. peloponnesiaca (SCP) phytochemical composition and neuroprotective potential, while nothing is known of the responsible compounds. Given that SCP is a weaker antioxidant compared to other Sideritis spp., here, we investigated its potential beneficial properties against Aß aggregation. We characterized different SCP extracts and revealed their anti-aggregation activity by taking advantage of established C. elegans AD models. Importantly, we identified two pure compounds, namely, sideridiol and verbascoside, being responsible for the beneficial effects. Furthermore, we have revealed a potential anti-Aß aggregation mechanism for sideridiol. Our results support the use of mountain tea in the elderly against dementia and demonstrate the activity of sideridiol against Aß aggregation that could be exploited for drug development.

Bacterial Lectin FimH and Its Aggregation Hot-Spots: An Alternative Strategy against Uropathogenic Escherichia coli.

Type I fimbriae are the main adhesive organelles of uropathogenic Escherichia coli (UPEC), consisting of four different subunits. Their component with the most important role in establishing bacterial infections is the FimH adhesin located at the fimbrial tip. This two-domain protein mediates adhesion to host epithelial cells through interaction with terminal mannoses on epithelial glycoproteins. Here, we propose that the amyloidogenic potential of FimH can be exploited for the development of therapeutic agents against Urinary Tract Infections (UTIs). Aggregation-prone regions (APRs) were identified via computational methods, and peptide-analogues corresponding to FimH lectin domain APRs were chemically synthesized and studied with the aid of both biophysical experimental techniques and molecular dynamic simulations. Our findings indicate that these peptide-analogues offer a promising set of antimicrobial candidate molecules since they can either interfere with the folding process of FimH or compete for the mannose-binding pocket.

A Multilevel Study of Eupatorin and Scutellarein as Anti-Amyloid Agents in Alzheimer's Disease.

Today, Alzheimer's disease (AD)-the most common neurodegenerative disorder, which affects 50 million people-remains incurable. Several studies suggest that one of the main pathological hallmarks of AD is the accumulation of abnormal amyloid beta (Aß) aggregates; therefore, many therapeutic approaches focus on anti-Aß aggregation inhibitors. Taking into consideration that plant-derived secondary metabolites seem to have neuroprotective effects, we attempted to assess the effects of two flavones-eupatorin and scutellarein-on the amyloidogenesis of Aß peptides. Biophysical experimental methods were employed to inspect the aggregation process of Aß after its incubation with each natural product, while we monitored their interactions with the oligomerized Aß through molecular dynamics simulations. More importantly, we validated our in vitro and in silico results in a multicellular organismal model-namely, Caenorhabditis elegans-and we concluded that eupatorin is indeed able to delay the amyloidogenesis of Aß peptides in a concentration-dependent manner. Finally, we propose that further investigation could lead to the exploitation of eupatorin or its analogues as potential drug candidates.

SARS-CoV-2 Amino Acid Mutations Detection in Greek Patients Infected in the First Wave of the Pandemic.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a novel virus that belongs to the Coronoviridae family, emerged in December 2019, causing the COVID-19 pandemic in March 2020. Unlike previous SARS and Middle East respiratory syndrome (MERS) outbreaks, this virus has a higher transmissibility rate, albeit a lower case fatality rate, which results in accumulation of a significant number of mutations and a faster evolution rate. Genomic studies on the mutation rate of the virus, as well as the identification of mutations that prevail and their impact on disease severity, are of great importance for pandemic surveillance and vaccine and drug development. Here, we aim to identify mutations on the SARS-CoV-2 viral genome and their effect on the proteins they are located in, in Greek patients infected in the first wave of the pandemic. To this end, we perform SARS-CoV-2 amplicon-based NGS sequencing on nasopharyngeal swab samples from Greek patients and bioinformatic analysis of the results. Although SARS-CoV-2 is considered genetically stable, we discover a variety of mutations on the viral genome. In detail, 18 mutations are detected in total on 10 SARS-CoV-2 isolates. The mutations are located on ORF1ab, S protein, M protein, ORF3a and ORF7a. Sixteen are also detected in patients from other regions around the world, and two are identified for the first time in the present study. Most of them result in amino acid substitutions. These substitutions are analyzed using computational tools, and the results indicate minor or major impact on the proteins' structural stability, which could probably affect viral transmissibility and pathogenesis. The correlation of these variations with the viral load levels is examined, and their implication for disease severity and the biology of the virus are discussed.

Clusterin in Alzheimer's disease: An amyloidogenic inhibitor of amyloid formation?

Clusterin is a heterodimeric glycoprotein (α- and ß-chain), which has been described as an extracellular molecular chaperone. In humans, clusterin is an amyloid-associated protein, co-localizing with fibrillar deposits in several amyloidoses, including Alzheimer's disease. To clarify its potential implication in amyloid formation, we located aggregation-prone regions within the sequence of clusterin α-chain, via computational methods. We had peptide-analogues, which correspond to each of these regions, chemically synthesized and experimentally demonstrated that all of them can form amyloid-like fibrils. We also provide evidence that the same peptide-analogues can inhibit amyloid-ß fibril formation, potentially making them appropriate drug candidates for Alzheimer's disease. At the same time, our findings hint that the respective aggregation-prone clusterin regions may be implicated in the molecular mechanism in which clusterin inhibits amyloid formation. Furthermore, we suggest that molecular chaperones with amyloidogenic properties might have a role in the regulation of amyloid formation, essentially acting as functional amyloids.

From Proteomic Mapping to Invasion-Metastasis-Cascade Systemic Biomarkering and Targeted Drugging of Mutant BRAF-Dependent Human Cutaneous Melanomagenesis.

Melanoma is classified among the most notoriously aggressive human cancers. Despite the recent progress, due to its propensity for metastasis and resistance to therapy, novel biomarkers and oncogenic molecular drivers need to be promptly identified for metastatic melanoma. Hence, by employing nano liquid chromatography-tandem mass spectrometry deep proteomics technology, advanced bioinformatics algorithms, immunofluorescence, western blotting, wound healing protocols, molecular modeling programs, and MTT assays, we comparatively examined the respective proteomic contents of WM115 primary (n = 3955 proteins) and WM266-4 metastatic (n = 6681 proteins) melanoma cells. It proved that WM115 and WM266-4 cells have engaged hybrid epithelial-to-mesenchymal transition/mesenchymal-to-epithelial transition states, with TGF-ß controlling their motility in vitro. They are characterized by different signatures of SOX-dependent neural crest-like stemness and distinct architectures of the cytoskeleton network. Multiple signaling pathways have already been activated from the primary melanoma stage, whereas HIF1α, the major hypoxia-inducible factor, can be exclusively observed in metastatic melanoma cells. Invasion-metastasis cascade-specific sub-routines of activated Caspase-3-triggered apoptosis and LC3B-II-dependent constitutive autophagy were also unveiled. Importantly, WM115 and WM266-4 cells exhibited diverse drug response profiles, with epirubicin holding considerable promise as a beneficial drug for metastatic melanoma clinical management. It is the proteome navigation that enables systemic biomarkering and targeted drugging to open new therapeutic windows for advanced disease.

Targeting of copper-trafficking chaperones causes gene-specific systemic pathology in Drosophila melanogaster: prospective expansion of mutational landscapes that regulate tumor resistance to cisplatin.

Copper, a transition metal, is an essential component for normal growth and development. It acts as a critical co-factor of many enzymes that play key roles in diverse cellular processes. The present study attempts to investigate the regulatory functions decisively controlling copper trafficking during development and aging of the Drosophila model system. Hence, through engagement of the GAL4/UAS genetic platform and RNAi technology, we herein examined the in vivo significance of Atox1 and CCS genes, products of which pivotally govern cellular copper trafficking in fly tissue pathophysiology. Specifically, we analyzed the systemic effects of their targeted downregulation on the eye, wing, neuronal cell populations and whole-body tissues of the fly. Our results reveal that, in contrast to the eye, suppression of their expression in the wing leads to a notable increase in the percentage of malformed organs observed. Furthermore, we show that Atox1 or CCS gene silencing in either neuronal or whole-body tissues can critically affect the viability and climbing capacity of transgenic flies, while their double-genetic targeting suggests a rather synergistic mode of action of the cognate protein products. Interestingly, pharmacological intervention with the anti-cancer drug cisplatin indicates the major contribution of CCS copper chaperone to cisplatin's cellular trafficking, and presumably to tumor resistance often acquired during chemotherapy. Altogether, it seems that Atox1 and CCS proteins serve as tissue/organ-specific principal regulators of physiological Drosophila development and aging, while their tissue-dependent downregulation can provide important insights for Atox1 and CCS potential exploitation as predictive gene biomarkers of cancer-cell chemotherapy responses.

Haloadaptation: insights from comparative modeling studies of halophilic archaeal DHFRs.

Proteins of halophilic archaea function in high-salt concentrations that inactivate or precipitate homologous proteins from non-halophilic species. Haloadaptation and the mechanism behind the phenomenon are not yet fully understood. In order to obtain useful information, homology modeling studies of dihydrofolate reductases (DHFRs) from halophilic archaea were performed that led to the construction of structural models. These models were subjected to energy minimization, structural evaluation and analysis. Complementary approaches concerning calculations of the amino acid composition and visual inspection of the surfaces and cores of the models, as well as calculations of electrostatic surface potentials, in comparison to non-halophilic DHFRs were also performed. The results provide evidence that sheds some light on the phenomenon of haloadaptation: DHFRs from halophilic archaea may maintain their fold, in high-salt concentrations, by sharing highly negatively charged surfaces and weak hydrophobic cores.

Acquisition through horizontal gene transfer of plasmid pSMA198 by Streptococcus macedonicus ACA-DC 198 points towards the dairy origin of the species.

BACKGROUND: Streptococcus macedonicus is an intriguing streptococcal species whose most frequent source of isolation is fermented foods similarly to Streptococcus thermophilus. However, S. macedonicus is closely related to commensal opportunistic pathogens of the Streptococcus bovis/Streptococcus equinus complex. METHODOLOGY/PRINCIPAL FINDINGS: We analyzed the pSMA198 plasmid isolated from the dairy strain Streptococcus macedonicus ACA-DC 198 in order to provide novel clues about the main ecological niche of this bacterium. pSMA198 belongs to the narrow host range pCI305/pWV02 family found primarily in lactococci and to the best of our knowledge it is the first such plasmid to be reported in streptococci. Comparative analysis of the pSMA198 sequence revealed a high degree of similarity with plasmids isolated from Lactococcus lactis strains deriving from milk or its products. Phylogenetic analysis of the pSMA198 Rep showed that the vast majority of closely related proteins derive from lactococcal dairy isolates. Additionally, cloning of the pSMA198 ori in L. lactis revealed a 100% stability of replication over 100 generations. Both pSMA198 and the chromosome of S. macedonicus exhibit a high percentage of potential pseudogenes, indicating that they have co-evolved under the same gene decay processes. We identified chromosomal regions in S. macedonicus that may have originated from pSMA198, also supporting a long co-existence of the two replicons. pSMA198 was also found in divergent biotypes of S. macedonicus and in strains isolated from dispersed geographic locations (e.g. Greece and Switzerland) showing that pSMA198's acquisition is not a recent event. CONCLUSIONS/SIGNIFICANCE: Here we propose that S. macedonicus acquired plasmid pSMA198 from L. lactis via an ancestral genetic exchange event that took place most probably in milk or dairy products. We provide important evidence that point towards the dairy origin of this species.

CutProtFam-Pred: detection and classification of putative structural cuticular proteins from sequence alone, based on profile hidden Markov models.

The arthropod cuticle is a composite, bipartite system, made of chitin filaments embedded in a proteinaceous matrix. The physical properties of cuticle are determined by the structure and the interactions of its two major components, cuticular proteins (CPs) and chitin. The proteinaceous matrix consists mainly of structural cuticular proteins. The majority of the structural proteins that have been described to date belong to the CPR family, and they are identified by the conserved R&R region (Rebers and Riddiford Consensus). Two major subfamilies of the CPR family RR-1 and RR-2, have also been identified from conservation at sequence level and some correlation with the cuticle type. Recently, several novel families, also containing characteristic conserved regions, have been described. The package HMMER v3.0 (http://hmmer.janelia.org/) was used to build characteristic profile Hidden Markov Models based on the characteristic regions for 8 of these families, (CPF, CPAP3, CPAP1, CPCFC, CPLCA, CPLCG, CPLCW, Tweedle). In brief, these families can be described as having: CPF (a conserved region with 44 amino acids); CPAP1 and CPAP-3 (analogous to peritrophins, with 1 and 3 chitin-binding domains, respectively); CPCFC (2 or 3 C-x(5)-C repeats); and four of five low complexity (LC) families, each with characteristic domains. Using these models, as well as the models previously created for the two major subfamilies of the CPR family, RR-1 and RR-2 (Karouzou et al., 2007), we developed CutProtFam-Pred, an on-line tool (http://bioinformatics.biol.uoa.gr/CutProtFam-Pred) that allows one to query sequences from proteomes or translated transcriptomes, for the accurate detection and classification of putative structural cuticular proteins. The tool has been applied successfully to diverse arthropod proteomes including a crustacean (Daphnia pulex) and a chelicerate (Tetranychus urticae), but at this taxonomic distance only CPRs and CPAPs were recovered.

A consensus method for the prediction of 'aggregation-prone' peptides in globular proteins.

The purpose of this work was to construct a consensus prediction algorithm of 'aggregation-prone' peptides in globular proteins, combining existing tools. This allows comparison of the different algorithms and the production of more objective and accurate results. Eleven (11) individual methods are combined and produce AMYLPRED2, a publicly, freely available web tool to academic users (http://biophysics.biol.uoa.gr/AMYLPRED2), for the consensus prediction of amyloidogenic determinants/'aggregation-prone' peptides in proteins, from sequence alone. The performance of AMYLPRED2 indicates that it functions better than individual aggregation-prediction algorithms, as perhaps expected. AMYLPRED2 is a useful tool for identifying amyloid-forming regions in proteins that are associated with several conformational diseases, called amyloidoses, such as Altzheimer's, Parkinson's, prion diseases and type II diabetes. It may also be useful for understanding the properties of protein folding and misfolding and for helping to the control of protein aggregation/solubility in biotechnology (recombinant proteins forming bacterial inclusion bodies) and biotherapeutics (monoclonal antibodies and biopharmaceutical proteins).

LepChorionDB, a database of Lepidopteran chorion proteins and a set of tools useful for the identification of chorion proteins in Lepidopteran proteomes.

Chorion proteins of Lepidoptera have a tripartite structure, which consists of a central domain and two, more variable, flanking arms. The central domain is highly conserved and it is used for the classification of chorion proteins into two major classes, A and B. Annotated and unreviewed Lepidopteran chorion protein sequences are available in various databases. A database, named LepChorionDB, was constructed by searching 5 different protein databases using class A and B central domain-specific profile Hidden Markov Models (pHMMs), developed in this work. A total of 413 Lepidopteran chorion proteins from 9 moths and 1 butterfly species were retrieved. These data were enriched and organised in order to populate LepChorionDB, the first relational database, available on the web, containing Lepidopteran chorion proteins grouped in A and B classes. LepChorionDB may provide insights in future functional and evolutionary studies of Lepidopteran chorion proteins and thus, it will be a useful tool for the Lepidopteran scientific community and Lepidopteran genome annotators, since it also provides access to the two pHMMs developed in this work, which may be used to discriminate A and B class chorion proteins. LepChorionDB is freely available at http://bioinformatics.biol.uoa.gr/LepChorionDB.

A possible structural model of members of the CPF family of cuticular proteins implicating binding to components other than chitin.

The physical properties of cuticle are determined by the structure of its two major components, cuticular proteins (CPs) and chitin, and, also, by their interactions. A common consensus region (extended R&R Consensus) found in the majority of cuticular proteins, the CPRs, binds to chitin. Previous work established that beta-pleated sheet predominates in the Consensus region and we proposed that it is responsible for the formation of helicoidal cuticle. Remote sequence similarity between CPRs and a lipocalin, bovine plasma retinol binding protein (RBP), led us to suggest an antiparallel beta-sheet half-barrel structure as the basic folding motif of the R&R Consensus. There are several other families of cuticular proteins. One of the best defined is CPF. Its four members in Anopheles gambiae are expressed during the early stages of either pharate pupal or pharate adult development, suggesting that the proteins contribute to the outer regions of the cuticle, the epi- and/or exo-cuticle. These proteins did not bind to chitin in the same assay used successfully for CPRs. Although CPFs are distinct in sequence from CPRs, the same lipocalin could also be used to derive homology models for one A. gambiae and one Drosophila melanogaster CPF. For the CPFs, the basic folding motif predicted is an eight-stranded, antiparallel beta-sheet, full-barrel structure. Possible implications of this structure are discussed and docking experiments were carried out with one possible Drosophila ligand, 7(Z),11(Z)-heptacosadiene.