Meta-analysis and research on host–parasite interactions: past and future

Host–parasite interactions are characterised by a lack of stable species-specific traits that limits generalisations one can make even about particular host or parasite species. For instance, the virulence, life history traits or transmission mode of a given parasite species can depend on which of its suitable hosts it infects. In the search for general rules or patterns, meta-analysis provides a possible solution to the challenges posed by the highly variable outcomes of host–parasite interactions. It allows an estimate of the overall association between any factor and its biological response that transcends the particulars of given host and parasite taxonomic combinations. In this review, we begin with a historical overview of the use of meta-analysis in research on the ecology and evolution of host–parasite interactions. We then identify several key conceptual advances that were made possible only through meta-analytical synthesis. For example, meta-analysis revealed the predominant association between rates of host and parasite gene flow and local adaptation, as well as an unexpected latitudinal gradient in parasite virulence, or parasite-induced host mortality. Finally, we propose some areas of research on host–parasite interactions that are based on a mature theoretical foundation and for which there now exist sufficient primary results to make them ripe for meta-analysis. The search for the processes causing variability in parasite species richness among host species, and the link between the expression of host resistance and the specificity of parasites, are two such research areas. The main objective of this review is to promote meta-analysis as a synthetic tool overriding the idiosyncrasies of specific host–parasite combinations and capable of uncovering the universal trends, if any, in the evolutionary ecology of parasitism.     

Comparative parasitism of the fish Plagioscion squamosissimus in native and invaded river basins

Biological invasions are considered a major threat to biodiversity around the world, but the role of parasites in this process is still little investigated. Here, we compared parasite infections of a host species in the areas where it originated and where it was introduced, and in native and introduced species in the same environment, using the endoparasites of the fish Plagioscion squamosissimus (Sciaenidae) in 3 Brazilian basins. Samples were taken in 2 rivers where the species is native, i.e., Solim?es River (SO) and Tocantins River (TO), and where the species was introduced, the upper Paraná River (PR). In addition, abundances of diplostomids and larval nematodes were compared between P. squamosissimus and 2 native competitors in the PR, Hoplias malabaricus and Raphiodon vulpinus. In total, 13 species of endoparasites were recorded, but only Austrodiplostomum sp. and cestode cysts were present in all localities. Although infracommunity richness was similar, their species composition was slightly different among localities. General linear models using the relative condition factor of fish as response variables, and abundance of the most prevalent parasites as possible predictors showed that the condition of fish is negatively correlated with parasite abundance only in the native range (TO). Abundance of diplostomid eye flukes was higher in the PR, and in the native species H. malabaricus when compared to the invader, which might present an advantage for P. squamosissimus if they compete for prey. However, although P. squamosissimus may have lost some of its native parasites during its introduction to the PR, it is now possibly acting as a host for native generalist parasites.     

Multiple Conformers in Active Site of Human Dihydrofolate Reductase F31R/Q35E Double Mutant Suggest Structural Basis for Methotrexate Resistance

Methotrexate is a slow, tight-binding, competitive inhibitor of human dihydrofolate reductase (hDHFR), an enzyme that provides key metabolites for nucleotide biosynthesis. In an effort to better characterize ligand binding in drug resistance, we have previously engineered hDHFR variant F31R/Q35E. This variant displays a >650-fold decrease in methotrexate affinity, while maintaining catalytic activity comparable to the native enzyme. To elucidate the molecular basis of decreased methotrexate affinity in the doubly substituted variant, we determined kinetic and inhibitory parameters for the simple variants F31R and Q35E. This demonstrated that the important decrease of methotrexate affinity in variant F31R/Q35E is a result of synergistic effects of the combined substitutions. To better understand the structural cause of this synergy, we obtained the crystal structure of hDHFR variant F31R/Q35E complexed with methotrexate at 1.7-Å resolution. The mutated residue Arg-31 was observed in multiple conformers. In addition, seven native active-site residues were observed in more than one conformation, which is not characteristic of the wild-type enzyme. This suggests that increased residue disorder underlies the observed methotrexate resistance. We observe a considerable loss of van der Waals and polar contacts with the p-aminobenzoic acid and glutamate moieties. The multiple conformers of Arg-31 further suggest that the amino acid substitutions may decrease the isomerization step required for tight binding of methotrexate. Molecular docking with folate corroborates this hypothesis.Human dihydrofolate reductase (hDHFR)⁶ catalyzes the reduction of 7,8-dihydrofolate (DHF) to 5,6,7,8-tetrahydrofolate in a NADPH-dependent manner. 5,6,7,8-Tetrahydrofolate is a cofactor in purine and thymidylate biosynthesis, which are essential metabolites in cell division and proliferation. As a consequence of its essential role in nucleoside biosynthesis, hDHFR has been extensively exploited as a drug target. Inhibition with folate antagonists, or antifolates, arrests cell proliferation. The most effective clinical antifolate to date is methotrexate (MTX (Fig. 1)), a slow, tight-binding competitive inhibitor that displays high affinity for hDHFR (K_i^MTX = 3.4 pm). MTX is currently used to treat a variety of diseases, including cancer (1–3), and autoimmune diseases such as juvenile idiopathic arthritis (4). A number of resistance mechanisms to MTX have been observed in cancer patients, including impaired transport of MTX to the cytoplasm (5) and decreased retention of MTX in the cell (6). Numerous ex vivo studies have reported mutations in the hDHFR gene resulting in an enzyme variant with decreased affinity for MTX (7, 8). These have contributed to increase our understanding of the molecular basis for active-site discrimination between the substrate, DHF, and its competitive inhibitor, MTX. Understanding the molecular interactions that affect tight binding of MTX to the active site of DHFR will contribute to our understanding of antifolate binding to DHFR, which can in turn contribute to the design of more efficient inhibitors.Open in a separate windowFIGURE 1.Chemical structures of hDHFR ligands. Atom numbering is shown on DHF.A considerable number of DHFR active-site variants have been identified in MTX-resistant cancer cell lines (although never in patients) (9) or engineered in vitro to elucidate the role of active site residues in the binding of MTX. Amino acid substitutions at residues Ile-7 (10), Leu-22 (11, 12), Phe-31 (13), Phe-34 (14), Arg-70 (15), and Val-115 (16) have yielded MTX-resistant variants. These residues are all present in the folate-binding pocket (17). Because MTX and DHF bind to the active site of hDHFR in a similar manner, all known substitutions causing a decrease in MTX affinity also decrease DHF affinity and overall catalytic efficiency (7, 16, 18). However, the loss of DHF affinity and catalytic efficiency is generally smaller than the loss of MTX affinity. This is often attributed to formation of different contacts with either ligand due to the 180° inversion of the pterin ring of bound DHF relative to MTX (17, 19).Crystal structures of MTX-resistant point mutants have offered insight into the causes of decreased binding of MTX or other antifolates (17, 20–24). To this day, crystal structures of MTX-resistant hDHFR variants L22F, L22R, and L22Y (12), as well as F31G and F31S (25), complexed to various antifolates, have been reported. Only the L22Y variant has been co-crystallized with MTX. Despite its decreased affinity for MTX (L22Y K_i^MTX = 11 nm versus WT K_i^MTX < 31 pm (18)), the inhibitor in the variant structure was bound in the same way as in the native enzyme, making interpretation of decreased affinity difficult to assess. Nonetheless, the low probability conformation of residue Tyr-22 suggested that the presence of a bulky aromatic residue in this area of the folate-binding pocket generated unfavorable hydrophobic interactions with the 2,4-diaminopterin moiety of the inhibitor (12). This is also expected to reduce DHF substrate binding. Structures of MTX-resistant variants F31G and F31S were obtained complexed to N-[4-[(2,4-diaminofuro[2,3-d]pyrimidin-5-yl)methyl]methylamino]benzoyl]-l-glutamate (MTXO) (25), a MTX analog in which the 2,4–2,4-diaminopterin moiety is replaced by a 2,4-diaminofuropyrimidine moiety. Superposition of MTXO-bound variants with MTX-bound WT hDHFR revealed that the ligands bind to the active site in an analogous manner. It was suggested that decreased MTX binding in the substituted variants resulted from the loss of van der Waals and hydrophobic contacts established between the native Phe-31 and the p-ABA and 2,4-diaminopterin moieties of MTX. F31G and F31S display a 10-fold decrease in affinity for MTX relative to WT hDHFR (K_i^MTX < 31 pm (18)). Further Phe-31 variants (i.e. F31R; K_i^MTX = 7 nm, 200-fold decrease in MTX affinity) (10) display larger decreases in affinity relative to F31G and F31S. This cannot be rationalized by reduction of side-chain contacts with the inhibitor due to the presence of a smaller side chain.These results illustrate the difficulty of gaining insight into the molecular causes for altered MTX binding. This may be partly attributed to the very tight binding of MTX to the native enzyme, such that binding to resistant variants often remains in the sub-nanomolar or low nanomolar range, where the general mode of ligand binding has not changed appreciably relative to the native enzyme. Combining active-site mutations in hDHFR by protein engineering has been shown to generate variants with greatly decreased affinity to MTX (18, 26). Studying the molecular interactions in highly MTX-resistant hDHFR variants offers the possibility of capturing more important changes in enzyme-ligand interactions.Here, we report detailed observations for the mode of MTX resistance in the combinatorial variant F31R/Q35E. Variant F31R/Q35E is a relevant candidate for better understanding the specific interactions that govern ligand recognition in the folate binding site, because it displays a >650-fold decrease in MTX affinity (K_i^MTX = 21 nm) accompanied by a modest, 9-fold decrease of affinity for the substrate DHF relative to WT hDHFR (18). In addition, we have recently shown that this variant is an efficient selectable marker for various mammalian cell types, including murine hematopoietic stem cells (18).⁷ Because mutations giving rise to MTX resistance are not observed in mammals, and because MTX is approved for human treatment, engineered resistant DHFRs offer great potential as human selective markers ex vivo or in vivo (10, 27, 28). To better understand the effect of either amino acid substitution on each ligand, a kinetic double mutant cycle was constructed with the simple variants F31R and Q35E. The crystal structure of the F31R/Q35E variant was obtained with bound MTX at 1.7-Å resolution, to elucidate the structural basis of MTX resistance in this variant. In addition, molecular docking was performed with the F31R/Q35E structure to evaluate the role of the two substitutions toward folate binding. Overall, the results reveal synergistic effects of the combined substitutions toward loss of MTX binding, characterized by increased disorder of specific residues throughout the active site of the highly MTX-resistant F31R/Q35E variant.     

Two structurally atypical HEAT domains in the C-terminal portion of human eIF4G support binding to eIF4A and Mnk1

The X-ray structure of the C-terminal region of human eukaryotic translation initiation factor 4G (eIF4G) has been determined at 2.2 A resolution, revealing two atypical HEAT-repeat domains. eIF4G recruits various translation factors and the 40S ribosomal subunit to the mRNA 5' end. In higher eukaryotes, the C terminus of eIF4G (4G/C) supports translational regulation by recruiting eIF4A, an RNA helicase, and Mnk1, the kinase responsible for phosphorylating eIF4E. Structure-guided surface mutagenesis and protein-protein interaction assays were used to identify binding sites for eIF4A and Mnk1 within the HEAT-repeats of 4G/C. p97/DAP5, a translational modulator homologous to eIF4G, lacks an eIF4A binding site in the corresponding region. The second atypical HEAT domain of the 4G/C binds Mnk1 using two conserved aromatic/acidic-box (AA-box) motifs. Within the first AA-box, the aromatic residues contribute to the hydrophobic core of the domain, while the acidic residues form a negatively charged surface feature suitable for electrostatic interactions with basic residues in Mnk1.     

Parasitological Consequences of Overcrowding in Protected Areas

For the past several years, there has been growing interest in understanding the dynamics of parasites in ecosystems, as well as the diversity of ways in which they interfere with conservation and health preoccupations. Although it is widely recognized that many conservation practices (e.g., wildlife translocations, species removal, food supplementation) may be associated with parasite-related problems, less attention has been devoted to exploring the parasitological consequences of the overcrowding of animals in protected wildlife areas. Here, we discuss this important ecological/epidemiological problem, presenting at the same time an overview of the main questions and challenges in this field. Using empirical and theoretical examples chosen from the literature, we focus particularly on the interactions between the overcrowding of free living species and parasite population dynamics, the evolution of parasite virulence, the indirect effects on the structure of invertebrate communities, as well as the nutritional value of prey species. We argue that conservation policies should be aware more than ever of this problem, especially given the serious health risks currently posed by the spread of virulent viruses (e.g., avian influenza).