Search

Accurate computational identification of B-cell epitopes is crucial for the development of vaccines, therapies, and diagnostic tools. However, current structure-based prediction methods face limitations due to the dependency on experimentally solved structures. Here, we introduce DiscoTope-3.0, a markedly improved B-cell epitope prediction tool that innovatively employs inverse folding structure representations and a positive-unlabelled learning strategy, and is adapted for both solved and predicted structures. Our tool demonstrates a considerable improvement in performance over existing methods, accurately predicting linear and conformational epitopes across multiple independent datasets. Most notably, DiscoTope-3.0 maintains high predictive performance across solved, relaxed and predicted structures, alleviating the need for experimental structures and extending the general applicability of accurate B-cell epitope prediction by 3 orders of magnitude. DiscoTope-3.0 is made widely accessible on two web servers, processing over 100 structures per submission, and as a downloadable package. In addition, the servers interface with RCSB and AlphaFoldDB, facilitating large-scale prediction across over 200 million cataloged proteins. DiscoTope-3.0 is available at: https://services.healthtech.dtu.dk/service.php?DiscoTope-3.0.

The aim of this study was to investigate the potential of potato proteins and peptides as emulsifiers in the microencapsulation of fish oil by spray-drying. Microcapsules were produced using a potato protein extract, and fractions enriched in patatin and protease inhibitors. Furthermore, bioinformatically predicted emulsifier peptides from abundant potato proteins and a hydrolysate, obtained through targeted proteolysis of the extract, were investigated. During 28 days of storage at 25 °C, peptides and the hydrolysate exhibited better emulsifying properties and higher encapsulation efficiencies compared to native proteins and sodium caseinate. Significant differences (p < 0.05) were observed in the peroxide value (PV) and secondary volatile oxidation products between the microcapsules produced with peptides and native proteins. Microcapsules produced with peptides and hydrolysate showed the highest oxidative stability, not exceeding a PV of 10 meq/kg oil, and with concentrations of volatiles below the odor threshold in oil for five of the six studied compounds. These results show the emulsifying potential of potato peptides and hydrolysate for use in microencapsulation of hydrophobic bioactive ingredients such as fish oil.

The TCRDivER algorithm incorporates both clone size and sequence similarity of T cell receptors (TCR) to reveal greater complexity in TCR repertoires and novel TCR repertoire traits due to a nonlinear correlation of the two.

Abstract Background In a range of human disorders such as multiple myeloma (MM), immunoglobulin light chains (IgLCs) can be produced at very high concentrations. This can lead to pathological aggregation and deposition of IgLCs in different tissues, which in turn leads to severe and potentially fatal organ damage. However, IgLCs can also be highly soluble and non-toxic. It is generally thought that the cause for this differential solubility behaviour is solely found within the IgLC amino acid sequences, and a variety of individual sequence-related biophysical properties (e.g. thermal stability, dimerisation) have been proposed in different studies as major determinants of the aggregation in vivo. Here, we investigate biophysical properties underlying IgLC amyloidogenicity. Results We introduce a novel and systematic workflow, Thermodynamic and Aggregation Fingerprinting (ThAgg-Fip), for detailed biophysical characterisation, and apply it to nine different MM patient-derived IgLCs. Our set of pathogenic IgLCs spans the entire range of values in those parameters previously proposed to define in vivo amyloidogenicity; however, none actually forms amyloid in patients. Even more surprisingly, we were able to show that all our IgLCs are able to form amyloid fibrils readily in vitro under the influence of proteolytic cleavage by co-purified cathepsins. Conclusions We show that (I) in vivo aggregation behaviour is unlikely to be mechanistically linked to any single biophysical or biochemical parameter and (II) amyloidogenic potential is widespread in IgLC sequences and is not confined to those sequences that form amyloid fibrils in patients. Our findings suggest that protein sequence, environmental conditions and presence and action of proteases all determine the ability of light chains to form amyloid fibrils in patients.

Lactic acid bacteria (LAB) have an extracellular proteolytic system that includes a multi-domain, cell envelope protease (CEP) with a subtilisin homologous protease domain. These CEPs have different proteolytic activities despite having similar protein sequences. Structural characterization has previously been limited to CEP homologs of dairy- and human-derived LAB strains, excluding CEPs of plant-derived LAB strains. CEP structures are a challenge to determine experimentally due to their large size and attachment to the cell envelope. This study aims to clarify the prevalence and structural diversity of CEPs by using the structure prediction software AlphaFold 2. Domain boundaries are clarified based on a comparative analysis of 21 three-dimensional structures, revealing novel domain architectures of CEP homologs that are not necessarily restricted to specific LAB species or ecological niches. The C-terminal flanking region of the protease domain is divided into fibronectin type-III-like domains with various structural traits. The analysis also emphasizes the existence of two distinct domains for cell envelope attachment that are preceded by an intrinsically disordered cell wall spanning domain. The domain variants and their combinations provide CEPs with different stability, proteolytic activity, and potentially adhesive properties, making CEPs targets for steering proteolytic activity with relevance for both food development and human health.

Abstract Dietary antioxidants are an important preservative in food and have been suggested to help in disease prevention. With consumer demands for less synthetic and safer additives in food products, the food industry is searching for antioxidants that can be marketed as natural. Peptides derived from natural proteins show promise, as they are generally regarded as safe and potentially contain other beneficial bioactivities. Antioxidative peptides are usually obtained by testing various peptides derived from hydrolysis of proteins by a selection of proteases. This slow and cumbersome trial-and-error approach to identify antioxidative peptides has increased interest in developing computational approaches for prediction of antioxidant activity and thereby reduce laboratory work. A few antioxidant predictors exist, however, no tool predicting the antioxidative properties of peptides is, to the best of our knowledge, currently available as a web-server. We here present the AnOxPePred tool and web-server ( http://services.bioinformatics.dtu.dk/service.php?AnOxPePred-1.0 ) that uses deep learning to predict the antioxidant properties of peptides. Our model was trained on a curated dataset consisting of experimentally-tested antioxidant and non-antioxidant peptides. For a variety of metrics our method displays a prediction performance better than a k-NN sequence identity-based approach. Furthermore, the developed tool will be a good benchmark for future predictors of antioxidant peptides.

There has been increasing interest in the role of T cells and their involvement in cancer, autoimmune and infectious diseases. However, the nature of T cell receptor (TCR) epitope recognition at a repertoire level is not yet fully understood. Due to technological advances a plethora of TCR sequences from a variety of disease and treatment settings has become readily available. Current efforts in TCR specificity analysis focus on identifying characteristics in immune repertoires which can explain or predict disease outcome or progression, or can be used to monitor the efficacy of disease therapy. In this context, clustering of TCRs by sequence to reflect biological similarity, and especially to reflect antigen specificity have become of paramount importance. We review the main TCR sequence clustering methods and the different similarity measures they use, and discuss their performance and possible improvement. We aim to provide guidance for non-specialists who wish to use TCR repertoire sequencing for disease tracking, patient stratification or therapy prediction, and to provide a starting point for those aiming to develop novel techniques for TCR annotation through clustering.

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disorder associated to deteriorating motor and cognitive functions, and short survival. The disease is caused by neuronal death which results in progressive muscle wasting and weakness, ultimately leading to lethal respiratory failure. The misbehaviour of a specific protein, TDP-43, which aggregates and becomes toxic in ALS patient’s neurons, is supposed to be one of the causes. TDP-43 is a DNA/RNA-binding protein involved in several functions related to nucleic acid metabolism. Sequestration of TDP-43 aggregates is a possible therapeutic strategy that could alleviate or block pathology. Here, we describe the selection and characterization of a new intracellular antibody (intrabody) against TDP-43 from a llama nanobody library. The structure of the selected intrabody was predicted in silico and the model was used to suggest mutations that enabled to improve its expression yield, facilitating its experimental validation. We showed how coupling experimental methodologies with in silico design may allow us to obtain an antibody able to recognize the RNA binding regions of TDP-43. Our findings illustrate a strategy for the mitigation of TDP-43 proteinopathy in ALS and provide a potential new tool for diagnostics.

T cells play a crucial role in controlling and driving the immune response with their ability to discriminate peptides derived from healthy as well as pathogenic proteins. In this review, we focus on the currently available computational tools for epitope prediction, with a particular focus on tools aimed at identifying neoepitopes, i.e. cancer-specific peptides and their potential for use in immunotherapy for cancer treatment. This review will cover how these tools work, what kind of data they use, as well as pros and cons in their respective applications.

Background Methods We evaluated the impact of a high percentage of leukemia blast contamination in patient-derived starting material (SM) on CAR T-cell drug product (DP) manufacturing. In vitro as well as in vivo models were employed to identify characteristics of the construct associated with better profile of safety in case of inadvertent B-cell leukemia transduction during CAR T-cell manufacturing.Results The presence of large amounts of CD19+ cells in SM did not affect the transduction level of DPs, as well as the CAR T-cell rate of expansion at the end of standard production of 14 days. DPs were deeply characterized by flow cytometry and molecular biology for Ig-rearrangements, showing that the level of B-cell contamination in DPs did not correlate with the percentage of CD19+ cells in SM, in the studied patient cohort. Moreover, we investigated whether CAR design may affect the control of CAR+ leukemia cells. We provided evidences that CAR.CD19 short linker (SL) prevents complete epitope masking in CD19+CAR+ leukemia cells and we demonstrated in vitro and in vivo that CD19 +CAR(SL)+leukemic cells are killed by CAR.CD19 T-cells.Conclusions Taken together, these data suggest that a VL-VH SL may result in a safe CAR-T product, even when manufacturing starts from biological materials characterized by heavy contamination of leukemia blasts.

Abstract Background The development of accurate epitope prediction tools is important in facilitating disease diagnostics, treatment and vaccine development. The advent of new approaches making use of antibody and TCR sequence information to predict receptor-specific epitopes have the potential to transform the epitope prediction field. Development and validation of these new generation of epitope prediction methods would benefit from regularly updated high-quality receptor-antigen complex datasets. Results To address the need for high-quality datasets to benchmark performance of these new generation of receptor-specific epitope prediction tools, a webserver called SCEptRe (Structural Complexes of Epitope-Receptor) was created. SCEptRe extracts weekly updated 3D complexes of antibody-antigen, TCR-pMHC and MHC-ligand from the Immune Epitope Database and clusters them based on antigen, receptor and epitope features to generate benchmark datasets. SCEptRe also provides annotated information such as CDR sequences and VDJ genes on the receptors. Users can generate custom datasets based by selecting thresholds for structural quality and clustering parameters (e.g. resolution, R-free factor, antigen or epitope sequence identity) based on their need. Conclusions SCEptRe provides weekly updated, user-customized comprehensive benchmark datasets of immune receptor-epitope structural complexes. These datasets can be used to develop and benchmark performance of receptor-specific epitope prediction tools in the future. SCEptRe is freely accessible at http://tools.iedb.org/sceptre.

The patent literature should reflect the past 30 years of engineering efforts directed toward developing monoclonal antibody therapeutics. Such information is potentially valuable for rational antibody design. Patents, however, are designed not to convey scientific knowledge, but to provide legal protection. It is not obvious whether antibody information from patent documents, such as antibody sequences, is useful in conveying engineering know-how, rather than as a legal reference only. To assess the utility of patent data for therapeutic antibody engineering, we quantified the amount of antibody sequences in patents destined for medicinal purposes and how well they reflect the primary sequences of therapeutic antibodies in clinical use. We identified 16,526 patent families covering major jurisdictions (e.g., US Patent and Trademark Office (USPTO) and World Intellectual Property Organization) that contained antibody sequences. These families held 245,109 unique antibody chains (135,397 heavy chains and 109,712 light chains) that we compiled in our Patented Antibody Database (PAD, http://naturalantibody.com/pad). We find that antibodies make up a non-trivial proportion of all patent amino acid sequence depositions (e.g., 11% of USPTO Full Text database). Our analysis of the 16,526 families demonstrates that the volume of patent documents with antibody sequences is growing, with the majority of documents classified as containing antibodies for medicinal purposes. We further studied the 245,109 antibody chains from patent literature to reveal that they very well reflect the primary sequences of antibody therapeutics in clinical use. This suggests that the patent literature could serve as a reference for previous engineering efforts to improve rational antibody design.

The cell envelope proteinase (CEP) of Lactococcus lactis is a large extracellular protease covalently linked to the peptidoglycan of the cell wall. Strains of L. lactis are typically auxotrophic for several amino acids and in order to grow to high cell densities in milk they need an extracellular protease. The structure of the entire CEP enzyme is difficult to determine experimentally due to the large size and due to the attachment to the cell surface. We here describe the use of a combination of structure prediction tools to create a structural model for the entire CEP enzyme of Lactococcus lactis. The model has implications for how the bacterium interacts with casein micelles during growth in milk, and it has implications regarding the energetics of the proteolytic system. Our model for the CEP indicates that the catalytic triad is activated through a structural change caused by interaction with the substrate. The CEP of L. lactis might become a useful model for the mode of action for enzymes belonging to the large class of S8 proteinases with a PA (protease associated) domain and a downstream fibronectin like domain.

Abstract Background Although the etiology of chronic lymphocytic leukemia (CLL), the most common type of adult leukemia, is still unclear, strong evidence implicates antigen involvement in disease ontogeny and evolution. Primary and 3D structure analysis has been utilised in order to discover indications of antigenic pressure. The latter has been mostly based on the 3D models of the clonotypic B cell receptor immunoglobulin (BcR IG) amino acid sequences. Therefore, their accuracy is directly dependent on the quality of the model construction algorithms and the specific methods used to compare the ensuing models. Thus far, reliable and robust methods that can group the IG 3D models based on their structural characteristics are missing. Results Here we propose a novel method for clustering a set of proteins based on their 3D structure focusing on 3D structures of BcR IG from a large series of patients with CLL. The method combines techniques from the areas of bioinformatics, 3D object recognition and machine learning. The clustering procedure is based on the extraction of 3D descriptors, encoding various properties of the local and global geometrical structure of the proteins. The descriptors are extracted from aligned pairs of proteins. A combination of individual 3D descriptors is also used as an additional method. The comparison of the automatically generated clusters to manual annotation by experts shows an increased accuracy when using the 3D descriptors compared to plain bioinformatics-based comparison. The accuracy is increased even more when using the combination of 3D descriptors. Conclusions The experimental results verify that the use of 3D descriptors commonly used for 3D object recognition can be effectively applied to distinguishing structural differences of proteins. The proposed approach can be applied to provide hints for the existence of structural groups in a large set of unannotated BcR IG protein files in both CLL and, by logical extension, other contexts where it is relevant to characterize BcR IG structural similarity. The method does not present any limitations in application and can be extended to other types of proteins.

T-cell receptors (TCR) mediate immune responses recognizing peptides in complex with major histocompatibility complexes (pMHC) displayed on the surface of cells. Resolving the challenge of predicting the cognate pMHC target of a TCR would benefit many applications in the field of immunology, including vaccine design/discovery and the development of immunotherapies. Here, we developed a model for prediction of TCR targets based on similarity to a database of TCRs with known targets. Benchmarking the model on a large set of TCRs with known target, we demonstrated how the predictive performance is increased (i) by focusing on CDRs rather than the full length TCR protein sequences, (ii) by incorporating information from paired α and β chains, and (iii) integrating information for all 6 CDR loops rather than just CDR3. Finally, we show how integration of the structure of CDR loops, as obtained through homology modeling, boosts the predictive power of the model, in particular in situations where no high-similarity TCRs are available for the query. These findings demonstrate that TCRs that bind to the same target also share, to a very high degree, sequence, and structural features. This observation has profound impact for future development of prediction models for TCR-pMHC interactions and for the use of such models for the rational design of T cell based therapies.

B-cells can neutralize pathogenic molecules by targeting them with extreme specificity using receptors secreted or expressed on their surface (antibodies). This is achieved via molecular interactions between the paratope (i.e., the antibody residues involved in the binding) and the interacting region (epitope) of its target molecule (antigen). Discerning the rules that define this specificity would have profound implications for our understanding of humoral immunogenicity and its applications. The aim of this work is to produce improved, antibody-specific epitope predictions by exploiting features derived from the antigens and their cognate antibodies structures, and combining them using statistical and machine learning algorithms. We have identified several geometric and physicochemical features that are correlated in interacting paratopes and epitopes, used them to develop a Monte Carlo algorithm to generate putative epitopes-paratope pairs, and train a machine-learning model to score them. We show that, by including the structural and physicochemical properties of the paratope, we improve the prediction of the target of a given B-cell receptor. Moreover, we demonstrate a gain in predictive power both in terms of identifying the cognate antigen target for a given antibody and the antibody target for a given antigen, exceeding the results of other available tools.

Antibodies are important for immunity and exist in several classes (IgM, IgD, IgA, IgG, IgE). They are composed of symmetric dimeric molecules with two antigen binding regions (Fab) and a constant part (Fc), usually depicted as Y-shaped molecules. Rheumatoid factors found in patients with rheumatoid arthritis are autoantibodies binding to IgG and paradoxically appear to circulate in blood alongside with their antigen (IgG) without reacting with it. Here, it is shown that rheumatoid factors do not react with native IgG in solution, and that their epitopes only become accessible upon certain physico-chemical treatments (e.g. heat treatment at 57 °C), by physical adsorption on a hydrophobic surface or by antigen binding. Moreover, chemical cross-linking in combination with mass spectrometry showed that the native state of IgG is a compact (closed) form and that the Fab parts of IgG shield the Fc region and thereby control access of rheumatoid factors and presumably also some effector functions. It can be inferred that antibody binding to pathogen surfaces induces a conformational change, which exposes the Fc part with its effector sites and rheumatoid factor epitopes. This has strong implications for understanding antibody structure and physiology and necessitates a conceptual reformulation of IgG models.

Protein hydrolysates show great promise as bioactive food and feed ingredients and for valorization of side-streams from e.g., the fish processing industry. We present a novel approach for hydrolysate characterization that utilizes proteomics data for calculation of weighted mean peptide properties (length, molecular weight, and charge) and peptide-level abundance estimation. Using a novel bioinformatic approach for subsequent prediction of biofunctional properties of identified peptides, we are able to provide an unprecedented, in-depth characterization. The study further characterizes bulk emulsifying, foaming, and in vitro antioxidative properties of enzymatic hydrolysates derived from cod frame by application of Alcalase and Neutrase, individually and sequentially, as well as the influence of heat pre-treatment. All hydrolysates displayed comparable or higher emulsifying activity and stability than sodium caseinate. Heat-treatment significantly increased stability but showed a negative effect on the activity and degree of hydrolysis. Lower degrees of hydrolysis resulted in significantly higher chelating activity, while the opposite was observed for radical scavenging activity. Combining peptide abundance with bioinformatic prediction, we identified several peptides that are likely linked to the observed differences in bulk emulsifying properties. The study highlights the prospects of applying proteomics and bioinformatics for hydrolysate characterization and in food protein science.

The Immune Epitope Database (IEDB) is a free public resource which catalogs experiments characterizing immune epitopes. To accommodate data from next generation repertoire sequencing experiments, we recently updated how we capture and query epitope specific antibodies and T cell receptors. Specifically, we are now storing partial receptor sequences sufficient to determine CDRs and VDJ gene usage which are commonly identified by repertoire sequencing. For previously captured full length receptor sequencing data, we have calculated the corresponding CDR sequences and gene usage information using IMGT numbering and VDJ gene nomenclature format. To integrate information from receptors defined at different levels of resolution, we grouped receptors based on their host species, receptor type and CDR3 sequence. As of August 2018, we have cataloged sequence information for more than 22,510 receptors in 18,292 receptor groups, shown to bind to more than 2,241 distinct epitopes. These data are accessible as full exports and through a new dedicated query interface. The later combines the new ability to search by receptor characteristics with previously existing capability to search by epitope characteristics such as the infectious agent the epitope is derived from, or the kind of immune response involved in its recognition. We expect that this comprehensive capture of epitope specific immune receptor information will provide new insights into receptor-epitope interactions, and facilitate the development of novel tools that help in the analysis of receptor repertoire data.

Identification of epitopes targeted by antibodies (B cell epitopes) is of critical importance for the development of many diagnostic and therapeutic tools. For clinical usage, such epitopes must be extensively characterized in order to validate specificity and to document potential cross-reactivity. B cell epitopes are typically classified as either linear epitopes, i.e. short consecutive segments from the protein sequence or conformational epitopes adapted through native protein folding. Recent advances in high-density peptide microarrays enable high-throughput, high-resolution identification and characterization of linear B cell epitopes. Using exhaustive amino acid substitution analysis of peptides originating from target antigens, these microarrays can be used to address the specificity of polyclonal antibodies raised against such antigens containing hundreds of epitopes. However, the interpretation of the data provided in such large-scale screenings is far from trivial and in most cases it requires advanced computational and statistical skills. Here, we present an online application for automated identification of linear B cell epitopes, allowing the non-expert user to analyse peptide microarray data. The application takes as input quantitative peptide data of fully or partially substituted overlapping peptides from a given antigen sequence and identifies epitope residues (residues that are significantly affected by substitutions) and visualize the selectivity towards each residue by sequence logo plots. Demonstrating utility, the application was used to identify and address the antibody specificity of 18 linear epitope regions in Human Serum Albumin (HSA), using peptide microarray data consisting of fully substituted peptides spanning the entire sequence of HSA and incubated with polyclonal rabbit anti-HSA (and mouse anti-rabbit-Cy3). The application is made available at: www.cbs.dtu.dk/services/ArrayPitope.

The immune system has developed a number of distinct complex mechanisms to shape and control the antibody repertoire. One of these mechanisms, the affinity maturation process, works in an evolutionary-like fashion: after binding to a foreign molecule, the antibody-producing B-cells exhibit a high-frequency mutation rate in the genome region that codes for the antibody active site. Eventually, cells that produce antibodies with higher affinity for their cognate antigen are selected and clonally expanded. Here, we propose a new statistical approach based on maximum entropy modeling in which a scoring function related to the binding affinity of antibodies against a specific antigen is inferred from a sample of sequences of the immune repertoire of an individual. We use our inference strategy to infer a statistical model on a data set obtained by sequencing a fairly large portion of the immune repertoire of an HIV-1 infected patient. The Pearson correlation coefficient between our scoring function and the IC50 neutralization titer measured on 30 different antibodies of known sequence is as high as 0.77 (p-value 10-6), outperforming other sequence- and structure-based models.

Rheumatoid arthritis (RA) is a chronic systemic autoimmune disorder of unknown etiology, which is characterized by inflammation in the synovium and joint damage. Although the pathogenesis of RA remains to be determined, a combination of environmental (e.g., viral infections) and genetic factors influence disease onset. Especially genetic factors play a vital role in the onset of disease, as the heritability of RA is 50–60%, with the human leukocyte antigen (HLA) alleles accounting for at least 30% of the overall genetic risk. Some HLA-DR alleles encode a conserved sequence of amino acids, referred to as the shared epitope (SE) structure. By analyzing the structure of a HLA-DR molecule in complex with Epstein-Barr virus (EBV), the SE motif is suggested to play a vital role in the interaction of MHC II with the viral glycoprotein (gp) 42, an essential entry factor for EBV. EBV has been repeatedly linked to RA by several lines of evidence and, based on several findings, we suggest that EBV is able to induce the onset of RA in predisposed SE-positive individuals, by promoting entry of B-cells through direct contact between SE and gp42 in the entry complex.

Antigen stimulation may be important for splenic marginal zone lymphoma pathogenesis. To address this hypothesis, the occurrence of stereotyped B-cell receptors was investigated in 133 SMZL (26 HCV+) compared with 4,414 HCDR3 sequences from public databases. Sixteen SMZL (12%) showed stereotyped BCR; 7 of 86 (8%) SMZL sequences retrieved from public databases also belonged to stereotyped HCDR3 subsets. Three categories of subsets were identified: i) "SMZL-specific subsets" (n=5), composed only of 12 SMZL (9 HCV-from our series); ii) "Non-Hodgkin’s lymphoma-like subsets" (n=5), comprising 5 SMZL (4 from our series) clustering with other indolent lymphomas; iii) “CLL-like subsets” (n=6), comprising 6 SMZL (3 from our series) that belonged to known CLL subsets (n=4) or clustered with public CLL sequences. Immunoglobulin 3D modeling of 3 subsets revealed similarities in antigen binding regions not limited to HCDR3. Overall, data suggest that the pathogenesis of splenic marginal zone lymphoma may involve also HCV-unrelated epitopes or an antigenic trigger common to other indolent lymphomas.

Resistance to chloroquine of malaria strains is known to be associated with a parasite protein named PfCRT, the mutated form of which is able to reduce chloroquine accumulation in the digestive vacuole of the pathogen. Whether the protein mediates extrusion of the drug acting as a channel or as a carrier and which is the protonation state of its chloroquine substrate is the subject of a scientific debate. We present here an analytical approach that explores which combination of hypotheses on the mechanism of transport and the protonation state of chloroquine are consistent with available equilibrium experimental data. We show that the available experimental data are not, by themselves, sufficient to conclude whether the protein acts as a channel or as a transporter, which explains the origin of their different interpretation by different authors. Interestingly, though, each of the two models is only consistent with a subset of hypotheses on the protonation state of the transported molecule. The combination of these results with a sequence and structure analysis of PfCRT, which strongly suggests that the molecule is a carrier, indicates that the transported species is either or both the mono and di-protonated forms of chloroquine. We believe that our results, besides shedding light on the mechanism of chloroquine resistance in P. falciparum, have implications for the development of novel therapies against resistant malaria strains and demonstrate the usefulness of an approach combining systems biology strategies with structural bioinformatics and experimental data.