Background

Understanding the origin and spread of mutations associated with drug resistance, especially in the context of combination therapy, will help guide strategies to halt and prevent the emergence of resistance. Unfortunately, studies have assessed these complex processes when resistance is already highly prevalent. Even further, information on the evolutionary dynamics leading to multidrug-resistant parasites is scattered and limited to areas with low or seasonal malaria transmission. This study describes the dynamics of strong selection for mutations conferring resistance against sulphadoxine-pyrimethamine (SP), a combination therapy, in western Kenya between 1992 and 1999, just before SP became first-line therapy (1999). Importantly, the study is based on longitudinal data, which allows for a comprehensive analysis that contrasts with previous cross-sectional studies carried out in other endemic regions.

Methods

This study used 236 blood samples collected between 1992 and 1999 in the Asembo Bay area of Kenya. Pyrosequencing was used to determine the alleles of dihydrofolate reductase (dhfr) and dihydropterote synthase (dhps) genes. Microsatellite alleles spanning 138 kb around dhfr and dhps, as well as, neutral markers spanning approximately 100 kb on chromosomes 2 and 3 were characterized.

Results

By 1992, the South-Asian dhfr triple mutant was already spreading, albeit in low frequency, in this holoendemic Kenyan population, prior to the use of SP as a first-line therapy. Additionally, dhfr triple mutant alleles that originated independently from the predominant Southeast Asian lineage were present in the sample set. Likewise, dhps double mutants were already present as early as 1992. There is evidence for soft selective sweeps of two dhfr mutant alleles and the possible emergence of a selective sweep of double mutant dhps alleles between 1992 and 1997. The longitudinal structure of the dataset allowed estimation of selection pressures on various dhfr and dhps mutants relative to each other based on a theoretical model tailored to P. falciparum. The data indicate that drug selection acted differently on the resistant alleles of dhfr and dhps, as evidenced by fitness differences. Thus a combination drug therapy such as SP, by itself, does not appear to select for "multidrug"-resistant parasites in areas with high recombination rate.

Conclusions

The complexity of these observations emphasizes the importance of population-based studies to evaluate the effects of strong drug selection on Plasmodium falciparum populations.

The World Health Organization currently does not recommend the use of dried blood spot specimens for drug resistance testing in patients undergoing antiretroviral therapy (ART). Therefore, HIV-1 resistance testing using peripheral blood mononuclear cells (PBMCs) may be of value in resource-limited settings. We compared genotypic resistance profiles in plasma and PBMCs from patients failing ART in two cities of Honduras (Tegucigalpa and San Pedro Sula), a resource-limited country. One hundred patients failing ART were randomly selected from a longitudinal patient monitoring cohort. Plasma and PBMC samples without patient identifier were used for genotypic resistance testing. Sequence data were analyzed, resistance profiles were determined and compared using Stanford HIV Drug Resistance Database algorithm. Specimens with concordant resistance profiles between the two compartments were 88% (95% CI: 80.3% - 94.5 %). Nine specimens (12%, 95% CI: 6.5% - 21.3%) had discordant resistance profiles of clinical significance. Logistic regression analyses indicated that patients on triple therapy were 17.24 times more likely to have concordant drug resistance profile than those on non-triple therapies (OR=17.24, 95% CI: 3.48, 83.33), while patients with increasing number of regimens and years on ART have a decreased rate of concordance (OR = 0.59, 95% CI: 0.32, 1.09 and OR = 0.62, 95% CI: 0.43, 0.88), respectively, than those with less number of regimens and years on ART. Our results show high level of concordance between plasma and PBMC resistance profiles, indicating the possibility of using PBMCs for drug resistance testing in resources-limited settings.

Commercially available HIV-1 drug resistance (HIVDR) genotyping assays are expensive and have limitations in detecting non-B subtypes and circulating recombinant forms that are co-circulating in resource-limited settings (RLS). This study aimed to optimize a low cost and broadly sensitive in-house assay in detecting HIVDR mutations in the protease (PR) and reverse transcriptase (RT) regions of pol gene. The overall plasma genotyping sensitivity was 95.8% (N = 96). Compared to the original in-house assay and two commercially available genotyping systems, TRUGENE® and ViroSeq®, the optimized in-house assay showed a nucleotide sequence concordance of 99.3%, 99.6% and 99.1%, respectively. The optimized in-house assay was more sensitive in detecting mixture bases than the original in-house (N = 87, P<0.001) and TRUGENE® and ViroSeq® assays. When the optimized in-house assay was applied to genotype samples collected for HIVDR surveys (N = 230), all 72 (100%) plasma and 69 (95.8%) of the matched dried blood spots (DBS) in the Vietnam transmitted HIVDR survey were genotyped and nucleotide sequence concordance was 98.8%; Testing of treatment-experienced patient plasmas with viral load (VL) ≥ and <3 log10 copies/ml from the Nigeria and Malawi surveys yielded 100% (N = 46) and 78.6% (N = 14) genotyping rates, respectively. Furthermore, all 18 matched DBS stored at room temperature from the Nigeria survey were genotyped. Phylogenetic analysis of the 236 sequences revealed that 43.6% were CRF01_AE, 25.9% subtype C, 13.1% CRF02_AG, 5.1% subtype G, 4.2% subtype B, 2.5% subtype A, 2.1% each subtype F and unclassifiable, 0.4% each CRF06_CPX, CRF07_BC and CRF09_CPX.

Conclusions

The optimized in-house assay is broadly sensitive in genotyping HIV-1 group M viral strains and more sensitive than the original in-house, TRUGENE® and ViroSeq® in detecting mixed viral populations. The broad sensitivity and substantial reagent cost saving make this assay more accessible for RLS where HIVDR surveillance is recommended to minimize the development and transmission of HIVDR.

As antiretroviral therapy (ART) is scaled up in resource-limited countries, surveillance for HIV drug resistance (DR) is vital to ensure sustained effectiveness of first-line ART. We have developed and applied a broadly sensitive dried-blood-spot (DBS)-based genotyping assay for surveillance of HIV-1 DR in international settings. In 2005 and 2006, 171 DBS samples were collected under field conditions from newly diagnosed HIV-1-infected individuals from Malawi (n = 58), Tanzania (n = 60), and China (n =53). In addition, 30 DBS and 40 plasma specimens collected from ART patients in China and Cameroon, respectively, were also tested. Of the 171 DBS analyzed at the protease and RT regions, 149 (87.1%) could be genotyped, including 49 (81.7%) from Tanzania, 47 (88.7%) from China, and 53 (91.4%) from Malawi. Among the 70 ART patient samples analyzed, 100% (30/30) of the Chinese DBS and 90% (36/40) of the Cameroonian plasma specimens were genotyped, including 8 samples with a viral load of <400 copies/ml. The results of phylogenetic analyses indicated that the subtype, circulating recombinant form (CRF), and unique recombinant form (URF) distribution was as follows: 73 strains were subtype C (34%), 37 were subtype B (17.2%), 24 each were CRF01_AE or CRF02_AG (11.2% each), 22 were subtype A1 (10.2%), and 9 were unclassifiable (UC) (4.2%). The remaining samples were minor strains comprised of 6 that were CRF07_BC (2.8%), 5 that were CRF10_CD (2.3%), 3 each that were URF_A1C and CRF08_BC (1.4%), 2 each that were G, URF_BC, and URF_D/UC (0.9%), and 1 each that were subtype F1, subtype F2, and URF_A1D (0.5%). Our results indicate that this broadly sensitive genotyping assay can be used to genotype DBS collected from areas with diverse HIV-1 group M subtypes and CRFs. Thus, the assay is likely to become a useful screening tool in the global resistance surveillance and monitoring of HIV-1 where multiple subtypes and CRFs are found.

Virtual screening by molecular docking has become a widely used approach to lead discovery in the pharmaceutical industry when a high resolution structure of the biological target of interest is available. The performance of three widely-used docking programs (Glide, GOLD, and DOCK) for virtual database screening is studied when they are applied to the same protein target and ligand set. Comparisons of the docking programs and scoring functions using a large and diverse data set of pharmaceutically interesting targets and active compounds are carried out. We focus on the problem of docking and scoring flexible compounds which are sterically capable of docking into a rigid conformation of the receptor. The Glide XP methodology is shown to consistently yield enrichments superior to the two alternative methods, while GOLD outperforms DOCK on average. The study also shows that docking into multiple receptor structures can decrease the docking error in screening a diverse set of active compounds.

The frequency of alleles with triple mutations conferring sulfadoxine-pyrimethamine (SP) resistance in the Peruvian Amazon Basin has declined (16.9% for dhfr and 0% for dhps compared to 47% for both alleles in 1997) 5 years after SP was replaced as the first-line treatment for Plasmodium falciparum malaria. Microsatellite analysis showed that the dhfr and dhps alleles are of common origin.

Torsion-angle sampling, as implemented in the Protein Local Optimization Program (PLOP), is used to generate multiple structurally variable single-conformer models which are in good agreement with X-ray data. An ensemble-refinement approach to differentiate between positional uncertainty and conformational heterogeneity is proposed.

Modeling structural variability is critical for understanding protein function and for modeling reliable targets for in silico docking experiments. Because of the time-intensive nature of manual X-ray crystallographic refinement, automated refinement methods that thoroughly explore conformational space are essential for the systematic construction of structurally variable models. Using five proteins spanning resolutions of 1.0–2.8 Å, it is demonstrated how torsion-angle sampling of backbone and side-chain libraries with filtering against both the chemical energy, using a modern effective potential, and the electron density, coupled with minimization of a reciprocal-space X-ray target function, can generate multiple structurally variable models which fit the X-ray data well. Torsion-angle sampling as implemented in the Protein Local Optimization Program (PLOP) has been used in this work. Models with the lowest R free values are obtained when electrostatic and implicit solvation terms are included in the effective potential. HIV-1 protease, calmodulin and SUMO-conjugating enzyme illustrate how variability in the ensemble of structures captures structural variability that is observed across multiple crystal structures and is linked to functional flexibility at hinge regions and binding interfaces. An ensemble-refinement procedure is proposed to differentiate between variability that is a consequence of physical con­formational heterogeneity and that which reflects uncertainty in the atomic coordinates.

A pyrosequencing protocol was developed as a rapid and reliable method to identify the mutations of the dhfr and dhps genes of Plasmodium falciparum that are associated with antifolate resistance. The accuracy and specificity of this method were tested using six laboratory-cultured P. falciparum isolates harboring known single nucleotide polymorphisms (SNPs) in the genes dhfr (codons 50, 51, 59, 108, and 164) and dhps (codons 436, 437, 540, 581, and 613). The lowest threshold for detection of all the SNPs tested by pyrosequencing was the equivalent of two to four parasite genomes. Also, this method was highly specific for P. falciparum, as it did not amplify any DNA products from the other species of human malaria parasites. We also mixed wild-type and mutant-type parasite DNAs in various proportions to determine how pyrosequencing, restriction fragment length polymorphism (RFLP), and direct conventional sequencing (for dhfr) compared with each other in detecting different SNPs in the mixture. In general, pyrosequencing and RFLP showed comparable sensitivities in detecting most of the SNPs in dhfr except for the 164L mutation, which required at least twice the amount of DNA for pyroseqencing as for RFLP. For detecting SNPs in dhps, pyrosequencing was slightly more sensitive than RFLP and direct sequencing. Overall, pyrosequencing was faster and less expensive than either RFLP or direct sequencing. Thus, pyrosequencing is a practical alternative method that can be used in a high-throughput format for molecular surveillance of antimalarial-drug resistance.

To determine the effect of placental malaria (PM) infection on the development of antibody responses to malaria in infants, we measured immunoglobulin G levels to seven different Plasmodium falciparum epitopes by using plasma samples collected at monthly intervals from infants born to mothers with and without PM. Overall, PM was associated with diminished antibody levels to all of the epitopes tested, especially with infants aged ≥4 to 12 months, and the difference was statistically significant for four of the seven epitopes (P < 0.0035). These findings suggest that PM can negatively influence the development of immune responses to malaria in infants.