Variation in Human Cytochrome P-450 Drug-Metabolism Genes: A Gateway to the Understanding of Plasmodium vivax Relapses

Monoamine oxidase-A (MAO-A) low-expression variants and increased risk of Plasmodium vivax malaria relapses

OBJECTIVES

Primaquine is essential for the radical cure of Plasmodium vivax malaria and must be metabolized into its bioactive metabolites. Accordingly, polymorphisms in primaquine-metabolizing enzymes can impact the treatment efficacy. This pioneering study explores the influence of monoamine oxidase-A (MAO-A) on primaquine metabolism and its impact on malaria relapses.

METHODS

Samples from 205 patients with P. vivax malaria were retrospectively analysed by genotyping polymorphisms in MAO-A and cytochrome P450 2D6 (CYP2D6) genes. We measured the primaquine and carboxyprimaquine blood levels in 100 subjects for whom blood samples were available on the third day of treatment. We also examined the relationship between the enzyme variants and P. vivax malaria relapses in a group of subjects with well-documented relapses.

RESULTS

The median carboxyprimaquine level was significantly reduced in individuals carrying low-expression MAO-A alleles plus impaired CYP2D6. In addition, this group experienced significantly more P. vivax relapses. The low-expression MAO-A status was not associated with malaria relapses when CYP2D6 had normal activity. This suggests that the putative carboxyprimaquine contribution is irrelevant when the CYP2D6 pathway is fully active.

CONCLUSIONS

We found evidence that the low-expression MAO-A variants can potentiate the negative impact of impaired CYP2D6 activity, resulting in lower levels of carboxyprimaquine metabolite and multiple relapses. The findings support the hypothesis that carboxyprimaquine may be further metabolized through CYP-mediated pathways generating bioactive metabolites that act against the parasite.

A nanochitosan-D-galactose formulation increases the accumulation of primaquine in the liver

Primaquine is the mainstream antimalarial drug to prevent Plasmodium vivax relapses. However, this drug can induce hemolysis in patients with glucose-6-phosphate dehydrogenase deficiency. Nanostructure formulations of primaquine loaded with D-galactose were used as a strategy to target the drug to the liver and decrease the hemolytic risks. Nanoemulsion (NE-Pq) and nanochitosan (NQ-Pq) formulations of primaquine diphosphate containing D-galactose were prepared and characterized by their physicochemistry properties. Pharmacokinetic and biodistribution studies were conducted using Swiss Webster mice. A single dose of 10 mg/kg of each nanoformulation or free primaquine solution was administered by gavage to the animals, which were killed at 0.5, 1, 2, 4, 8, and 24 hours. Blood samples and tissues were collected, processed, and analyzed by high-performance liquid chromatography. The nanoformulation showed sizes around 200 nm (NE-Pq) and 400 nm (NQ-Pq) and physicochemical stability for over 30 days. Free primaquine solution achieved higher primaquine Cmax in the liver than NE-Pq or NQ-Pq at 0.5 hours. However, the half-life and mean residence time (MRT) of primaquine in the liver were three times higher with the NQ-Pq formulation than with free primaquine, and the volume distribution was four times higher. Conversely, primaquine's half-life, MRT, and volume distribution in the plasma were lower for NQ-Pq than for free primaquine. NE-Pq, on the other hand, accumulated more in the lungs but not in the liver. Galactose-coated primaquine nanochitosan formulation showed increased drug targeting to the liver compared to free primaquine and may represent a promising strategy for a more efficient and safer radical cure for vivax malaria.

Pharmacogenomic variation in the Malagasy population: implications for the antimalarial drug primaquine metabolism

Aim: Antimalarial primaquine (PQ) eliminates liver hypnozoites of Plasmodium vivax. CYP2D6 gene variation contributes to PQ therapeutic failure. Additional gene variation may contribute to PQ efficacy. Information on pharmacogenomic variation in Madagascar, with vivax malaria and a unique population admixture, is scanty. Methods: The authors performed genome-wide genotyping of 55 Malagasy samples and analyzed data with a focus on a set of 28 pharmacogenes most relevant to PQ. Results: Mainly, the study identified 110 coding or splicing variants, including those that, based on previous studies in other populations, may be implicated in PQ response and copy number variation, specifically in chromosomal regions that contain pharmacogenes. Conclusion: With this pilot information, larger genome-wide association analyses with PQ metabolism and response are substantially more feasible.

Individual variation in Plasmodium vivax malaria risk: Are repeatedly infected people just unlucky?

Extensive research has examined why some people have frequent Plasmodium falciparum malaria episodes in sub-Saharan Africa while others remain free of disease most of the time. In contrast, malaria risk heterogeneity remains little studied in regions where P. vivax is the dominant species. Are repeatedly infected people in vivax malaria settings such as the Amazon just unlucky? Here, we briefly review evidence that human genetic polymorphism and acquired immunity after repeated exposure to parasites can modulate the risk of P. vivax infection and disease in predictable ways. One-fifth of the hosts account for 80% or more of the community-wide vivax malaria burden and contribute disproportionally to onward transmission, representing a priority target of more intensive interventions to achieve malaria elimination. Importantly, high-risk individuals eventually develop clinical immunity, even in areas with very low or residual malaria transmission, and may constitute a large but silent parasite reservoir.

Using Graph-Based Signatures to Guide Rational Antibody Engineering

Antibodies are essential experimental and diagnostic tools and as biotherapeutics have significantly advanced our ability to treat a range of diseases. With recent innovations in computational tools to guide protein engineering, we can now rationally design better antibodies with improved efficacy, stability, and pharmacokinetics. Here, we describe the use of the mCSM web-based in silico suite, which uses graph-based signatures to rapidly identify the structural and functional consequences of mutations, to guide rational antibody engineering to improve stability, affinity, and specificity.

Global perspectives on CYP2D6 associations with primaquine metabolism and Plasmodium vivax radical cure

Clinical trial and individual patient treatment outcomes have produced accumulating evidence that effective primaquine (PQ) treatment of Plasmodium vivax and P. ovale liver stage hypnozoites is associated with genetic variation in the human cytochrome P450 gene, CYP2D6. Successful PQ treatment of individual and population-wide infections by the Plasmodium species that generate these dormant liver stage forms is likely to be necessary to reach elimination of malaria caused by these parasites globally. Optimizing safe and effective PQ treatment will require coordination of efforts between the malaria and pharmacogenomics research communities.

Tafenoquine for the treatment of Plasmodium vivax malaria

INTRODUCTION

Plasmodium vivax malaria causes significant disease burden worldwide, especially in Latin America, Southeast Asia, and Oceania. P. vivax is characterized by the production of liver hypnozoites that cause clinical relapses upon periodic activation. Primaquine, an 8-aminoquinoline drug, has been the standard of care for decades to treat liver-stage P. vivax malaria; however, it requires long treatment regimens (one to two weeks) that lead to poor adherence and thus clinical relapses. Tafenoquine (TFQ), a newly available and efficacious single-dose 8-aminoquinoline, aims to address this challenge. Safe administration is possible when paired with the use of glucose-6-phosphate dehydrogenase (G6PD) diagnostics to prevent 8-aminoquinoline-induced hemolysis in patients with underlying G6PD deficiency (G6PDd).

AREAS COVERED

In this review, the authors present the recent literature regarding the pharmacology, efficacy, safety, and tolerability of TFQ and highlight regional differences in these areas. The authors also discuss the potential for TFQ, complemented with primaquine PQ and effective screening for G6PDd, to improve P. vivax clinical management and facilitate targeted mass drug administration in communities to decrease transmission.

EXPERT OPINION

Clinical studies show therapeutic efficacy of TFQ as well as a good performance in terms of safety and tolerability. Additional research regarding the effectiveness and safety TFQ in malaria elimination strategies such as targeted or mass drug administration are needed.

Systematic evaluation of computational tools to predict the effects of mutations on protein stability in the absence of experimental structures

Changes in protein sequence can have dramatic effects on how proteins fold, their stability and dynamics. Over the last 20 years, pioneering methods have been developed to try to estimate the effects of missense mutations on protein stability, leveraging growing availability of protein 3D structures. These, however, have been developed and validated using experimentally derived structures and biophysical measurements. A large proportion of protein structures remain to be experimentally elucidated and, while many studies have based their conclusions on predictions made using homology models, there has been no systematic evaluation of the reliability of these tools in the absence of experimental structural data. We have, therefore, systematically investigated the performance and robustness of ten widely used structural methods when presented with homology models built using templates at a range of sequence identity levels (from 15% to 95%) and contrasted performance with sequence-based tools, as a baseline. We found there is indeed performance deterioration on homology models built using templates with sequence identity below 40%, where sequence-based tools might become preferable. This was most marked for mutations in solvent exposed residues and stabilizing mutations. As structure prediction tools improve, the reliability of these predictors is expected to follow, however we strongly suggest that these factors should be taken into consideration when interpreting results from structure-based predictors of mutation effects on protein stability.

Genetic Variation of G6PD and CYP2D6: Clinical Implications on the Use of Primaquine for Elimination of Plasmodium vivax

Primaquine, an 8-aminoquinoline, is the only medication approved by the World Health Organization to treat the hypnozoite stage of Plasmodium vivax and P. ovale malaria. Relapse, triggered by activation of dormant hypnozoites in the liver, can occur weeks to years after primary infection, and provides the predominant source of transmission in endemic settings. Hence, primaquine is essential for individual treatment and P. vivax elimination efforts. However, primaquine use is limited by the risk of life-threatening acute hemolytic anemia in glucose-6-phosphate dehydrogenase (G6PD) deficient individuals. More recently, studies have demonstrated decreased efficacy of primaquine due to cytochrome P450 2D6 (CYP2D6) polymorphisms conferring an impaired metabolizer phenotype. Failure of standard primaquine therapy has occurred in individuals with decreased or absent CYP2D6 activity. Both G6PD and CYP2D6 are highly polymorphic genes, with considerable geographic and interethnic variability, adding complexity to primaquine use. Innovative strategies are required to overcome the dual challenge of G6PD deficiency and impaired primaquine metabolism. Further understanding of the pharmacogenetics of primaquine is key to utilizing its full potential. Accurate CYP2D6 genotype-phenotype translation may optimize primaquine dosing strategies for impaired metabolizers and expand its use in a safe, efficacious manner. At an individual level the current challenges with G6PD diagnostics and CYP2D6 testing limit clinical implementation of pharmacogenetics. However, further characterisation of the overlap and spectrum of G6PD and CYP2D6 activity may optimize primaquine use at a population level and facilitate region-specific dosing strategies for mass drug administration. This precision public health approach merits further investigation for P. vivax elimination.

Prediction of the CYP2D6 enzymatic activity based on investigating of the CYP2D6 genotypes around the vivax malaria patients in Yunnan Province, China

BACKGROUND

In recent years, the incidence rate of vivax malaria recurrence still had 3.1% in Yunnan Province population after eradication therapy using primaquine (PQ). In order to understand the specific failure reasons for preventing vivax malaria relapses, a preliminary exploration on the CYP2D6 enzyme activity was carried out in the vivax malaria patients in Yunnan Province population by analysing mutational polymorphism in the coding region of CYP2D6 gene.

METHODS

Blood samples were collected from vivax malaria patients with suspected relapse (SR) and non-relapsed (NR) malaria in Yunnan Province. The DNA fragments containing 9 exons regions of human CYP2D6 gene were amplified by performing PCR and sequenced. The sequencing results were aligned by using DNAStar 11.0 to obtain the coding DNA sequence (CDS) of CYP2D6 gene. DnaSP 6.11.01 software was used to identify mutant polymorphisms and haplotypes of the CDS chain. The waterfall function of GenVisR package in R was utilized to visualize the mutational landscape. The alleles of CYP2D6 gene were identified according to the criteria prescribed by Human Cytochrome P450 (CYP) Allele Nomenclature Committee Database and the CYP2D6 enzyme activity was predicted based on diploid genotype.

RESULTS

A total of 320 maternal CDS chains, including 63 from SR group and 257 from NR group, were obtained. Twelve mutant loci, including c.31 (rs769259), c.100 (rs1065852), c.271 (rs28371703), c.281 (rs28371704), c.294 (rs28371705), c.297 (rs200269944), c.336 (rs1081003), c.408 (rs1058164), c.505 (rs5030865), c.801 (rs28371718), c.886 (rs16947), and c.1,457 (rs1135840) were observed on the 640 CDS chains (including 320 maternal and 320 paternal chains). The high-frequency mutation at rs1135840 (0.703) and low-frequency mutation, such as rs28371703, were detected only in the SR group. The frequency of mutant rs1058164 and rs1135840 were significantly increased in the SR group ([Formula: see text]= 4.468, 5.889, P < 0.05), as opposed to the NR group. Of the 23 haplotypes (from Hap_1 to Hap_23), the nomenclatures of 11 allelic forms could be found: Hap_3 was non-mutant, Hap_2 accounted for the highest frequency (36.9%, 236/640), and Hap_9 had the most complex sequence structure, containing 7 loci mutations. Allele *10 was the most frequent among these genotypes (0.423). Among the allele *10 standard named genotypes, *1/*10, *1/*1 and *2/*10 were significantly more frequent in the NR group ([Formula: see text]= 3.911, P < 0.05) and all showed uncompromised enzyme activity; the impaired genotype *10/*39 was more frequent in the SR group ([Formula: see text]= 10.050, P < 0.05), and genotype *4/*4was detected only in the SR group.

CONCLUSION

In the patients receiving PQ dosage in Yunnan Province population, both rs1135840 single nucleotide polymorphism and *10 allele form was common in the CYP2D6 gene. Low-frequency mutation sites, such as rs28371703, were only presented in patients with vivax malaria relapse.

Antimalarials for children with Plasmodium vivax infection: Current status, challenges, and research priorities

The aim of this narrative review is to summarise efficacy and pharmacokinetic data for Plasmodium vivax in children. The burden of P. vivax malaria in children continues to remain a significant public health issue, and the need for improved treatment regimens for this vulnerable population is critical. Relapse after re-activation of dormant liver-stage hypnozoites poses additional challenges for treatment, elimination, and control strategies for P. vivax. Whilst it is recognised that paediatric pharmacology may be significantly influenced by anatomical and physiological changes of childhood, dosing regimens often continue to be extrapolated from adult data, highlighting the need for antimalarial dosing in children to be evaluated in early phase clinical trials. This will ensure that globally recommended treatment regimens do not result in suboptimal dosing in children. Furthermore, the development of affordable paediatric formulations to enhance treatment acceptability and widespread G6PD testing to facilitate use of anti-hypnozoite treatment such as primaquine and tafenoquine, should be further prioritised. As the world prepares for malaria elimination, a renewed focus on P. vivax malaria provides an ideal opportunity to harness momentum and ensure that all populations, including children have access to safe, efficacious, and correctly dosed antimalarial therapies.

CYP2D6 Allele Frequency in Five Malaria Vivax Endemic Areas From Brazilian Amazon Region

Genetic variability was linked with individual responses to treatment and susceptibility to malaria by Plasmodium vivax. Polymorphisms in the CYP2D6 gene may modulate enzyme level and activity, thereby affecting individual responses to pharmacological treatment. The aim of the study was to investigate whether or not CYP2D6 single nucleotide polymorphisms rs1065852, rs38920-97, rs16947 and rs28371725 are unequally distributed in malaria by Plasmodium vivax individuals from the Brazilian Amazon region. The blood samples were collected from 220 unrelated Plasmodium vivax patients from five different endemic areas. Genotyping was performed using SNaPshot^® and real-time polymerase chain reaction methods. In all five areas, the rs1065852 (CYP2D6*10, C.100C > T), rs3892097 (CYP2D6*4, 1846C > T) and rs16947 (CYP2D6*2, C.2850G > A), as a homozygous genotype, showed the lowest frequencies. The rs28371725 (CYP2D6*41, 2988G > A) homozygous genotype was not detected, while the allele A was found in a single patient from Macapá region. No deviations from Hardy-Weinberg equilibrium were found, although a borderline p-value was observed (p = 0.048) for the SNP rs3892097 in Goianésia do Pará, Pará state. No significant associations were detected in these frequencies among the five studied areas. For the SNP rs3892097, a higher frequency was observed for the C/T heterozygous genotype in the Plácido de Castro and Macapá, Acre and Amapá states, respectively. The distribution of the CYP2D6 alleles investigated in the different areas of the Brazilian Amazon is not homogeneous. Further investigations are necessary in order to determine which alleles might be informative to assure optimal drug dosing recommendations based on experimental pharmacogenetics.

Simplified and Rapid Determination of Primaquine and 5,6-Orthoquinone Primaquine by UHPLC-MS/MS: Its Application to a Pharmacokinetic Study

The method for the determination of primaquine (PQ) and 5,6-orthoquinone primaquine (5,6-PQ), the representative marker for PQ active metabolites, via CYP2D6 in human plasma and urine has been validated. All samples were extracted using acetonitrile for protein precipitation and analyzed using the ultra-high-performance liquid chromatography–tandem mass spectrometry (UHPLC-MS/MS) system. Chromatography separation was carried out using a Hypersil GOLD^TM aQ C₁₈ column (100 × 2.1 mm, particle size 1.9 μm) with a C₁₈ guard column (4 × 3 mm) flowed with an isocratic mode of methanol, water, and acetonitrile in an optimal ratio at 0.4 mL/min. The retention times of 5,6-PQ and PQ in plasma and urine were 0.8 and 1.6 min, respectively. The method was validated according to the guideline. The linearity of the analytes was in the range of 25–1500 ng/mL. The matrix effect of PQ and 5,6-PQ ranged from 100% to 116% and from 87% to 104% for plasma, and from 87% to 89% and from 86% to 87% for urine, respectively. The recovery of PQ and 5,6-PQ ranged from 78% to 95% and form 80% to 98% for plasma, and from 102% to from 112% to 97% to 109% for urine, respectively. The accuracy and precision of PQ and 5,6-PQ in plasma and urine were within the acceptance criteria. The samples should be kept in the freezer (−80 °C) and analyzed within 7 days due to the metabolite stability. This validated UHPLC-MS/MS method was beneficial for a pharmacokinetic study in subjects receiving PQ.

Impact of CYP2D6 Genetic Variation on Radical Cure of Plasmodium vivax Malaria

Plasmodium vivax (P. vivax) is the most widespread human malaria parasite, with 2.5 billion people at risk of infection worldwide. P. vivax forms liver hypnozoites, which trigger further symptomatic episodes (relapses) weeks or months after the initial episode. Radical cure of vivax malaria requires hypnozoitocide therapy to prevent relapses. The two US Food and Drug Administration (FDA)-approved hypnozoiticides for human use, primaquine, and tafenoquine, are pro-drugs, that require in vivo conversion into metabolites with redox activity. This mini-review focuses on the association between CYP2D6-mediated hydroxylation and hypnozoitocide efficacy of primaquine and tafenoquine. Studies in murine models show that the antimalarial activity of primaquine and tafenoquine is abolished by CYP2D knock-out and partially restored by knock-in of humanized CYP2D6. Human studies explored the impact of CYP2D6 genetic variation and genotype-inferred CYP2D6 phenotype on anti-relapse efficacy. Most, but not all, studies with primaquine report higher rates of relapse in patients with decreased CYP2D6 activity (activity scores (AS) ≤ 1) compared to normal activity (AS ≥ 1.5). Potential factors for discordance among studies include risk of reinfection in endemic areas, adherence to primaquine-treatment, assignment of CYP2D6 phenotypes based on CYP2D6 polymorphism and choice of AS values for dichotomizing the study cohorts. Tafenoquine anti-relapse efficacy did not differ between patients with AS < 1 vs. AS ≥ 1.5 in 2 studies. Absence/small number of poor CYP2D6 metabolizers in AS ≤ 1 groups, combined with lesser dependence of tafenoquine on CYP2D6-mediated conversion into active redox metabolites may account for this result. Additional tafenoquine studies with larger representation of poor CYP2D6 metabolizers are warranted.

CYP2D6 Genetic Variation and Its Implication for Vivax Malaria Treatment in Madagascar

Plasmodium vivax is one of the five human malaria parasite species, which has a wide geographical distribution and can cause severe disease and fatal outcomes. It has the ability to relapse from dormant liver stages (hypnozoites), weeks to months after clearance of the acute blood-stage infection. An 8-aminoquinoline drug primaquine (PQ) can clear the hypnozoites, and thus can be used as an anti-relapse therapeutic agent. Recently, a number of studies have found that its efficacy is compromised by polymorphisms in the cytochrome P450 2D6 (CYP2D6) gene; decreased or absence of CYP2D6 activity contributes to PQ therapeutic failure. The present study sought to characterize CYP2D6 genetic variation in Madagascar, where populations originated from admixture between Asian and African populations, vivax malaria is endemic, and PQ can be deployed soon to achieve national malaria elimination. In a total of 211 samples collected from two health districts, CYP2D6 decreased function alleles CYP2D6*10, *17, *29, *36+*10, and *41 were observed at frequencies of 3.55-17.06%. In addition, nonfunctional alleles were observed, the most common of which were CYP2D6*4 (2.13%), *5 (1.66%), and the *4x2 gene duplication (1.42%). Given these frequencies, 34.6% of the individuals were predicted to be intermediate metabolizers (IM) with an enzyme activity score (AS) ≤ 1.0; both the IM phenotype and AS ≤ 1.0 have been found to be associated with PQ therapeutic failure. Furthermore, the allele and genotype frequency distributions add to the archaeological and genomic evidence of Malagasy populations constituting a unique, Asian-African admixed origin. The results from this exploratory study provide fresh insights about genomic characteristics that could affect the metabolism of PQ into its active state, and may enable optimization of PQ treatment across human genetic diversity, which is critical for achieving P. vivax elimination.

Monitoring Plasmodium vivax resistance to antimalarials: Persisting challenges and future directions

Emerging antimalarial drug resistance may undermine current efforts to control and eliminate Plasmodium vivax, the most geographically widespread yet neglected human malaria parasite. Endemic countries are expected to assess regularly the therapeutic efficacy of antimalarial drugs in use in order to adjust their malaria treatment policies, but proper funding and trained human resources are often lacking to execute relatively complex and expensive clinical studies, ideally complemented by ex vivo assays of drug resistance. Here we review the challenges for assessing in vivo P. vivax responses to commonly used antimalarials, especially chloroquine and primaquine, in the presence of confounding factors such as variable drug absorption, metabolism and interaction, and the risk of new infections following successful radical cure. We introduce a simple modeling approach to quantify the relative contribution of relapses and new infections to recurring parasitemias in clinical studies of hypnozoitocides. Finally, we examine recent methodological advances that may render ex vivo assays more practical and widely used to confirm P. vivax drug resistance phenotypes in endemic settings and review current approaches to the development of robust genetic markers for monitoring chloroquine resistance in P. vivax populations.

The association of CYP2D6 gene polymorphisms in the full-length coding region with higher recurrence rate of vivax malaria in Yunnan Province, China

Background

Accumulating evidence suggest that compromised CYP2D6 enzyme activity caused by gene mutation could contribute to primaquine failure for the radical cure of vivax malaria. The current study aims to preliminarily reveal the association between the recurrence of vivax malaria in Yunnan Province and CYP2D6 gene mutation by analysing polymorphisms in the entire coding region of human CYP2D6 gene.

Methods

Blood samples were collected from patients with vivax malaria, who received "chloroquine and 8-day course of primaquine therapy" in Yunnan Province. The suspected relapsed cases were determined by epidemiological approaches and gene sequence alignment. PCR was conducted to amplify the CYP2D6 gene in the human genome, and the amplified products were then sequenced to compare with the non-mutation “reference” sequence, so as to ensure correct sequencing results and to determine 9 exon regions. Subsequently, the DNA sequences of 9 exons were spliced into the coding DNA sequence (CDS), which, by default, is known as maternal CDS. The paternal CDS was obtained by adjusting the bases according to the sequencing peaks. The mutation loci, haplotypes (star alleles), genotypes and odds ratios (OR) of all the CDSs were analysed.

Results

Of the119 maternal CDS chains in total with 1491 bp in length, 12 mutation sites in the 238 maternal and paternal CDS chains were detected. The c.408G > C mutation was most frequently detected in both suspected relapsed group (SR) and non-relapsed group (NR), reaching 85.2% (75/88) and 76.0% (114/150), respectively. The c.886C > T mutation was most closely related to the recurrence of vivax malaria (OR = 2.167, 95% CI 1.104–4.252, P < 0.05). Among the 23 haplotypes (Hap_1 ~ Hap_23), Hap_3 was non-mutant, and the sequence structure of Hap_9 was the most complicated one. Five star alleles, including *1, *2, *4, *10 and *39, were confirmed by comparison, and CYP2D6*10 allele accounted for the largest percentage (45.4%, 108/238). The frequency of CYP2D6*2 allele in the SR group was significantly higher than that in the NR group (Χ² = 16.177, P < 0.05). Of the defined 24 genotypes, 8 genotypes, including *4/*4, *4/*o, *2/*39, *39/*m, *39/*x, *1/*r, *1/*n, and *v/*10, were detected only in the SR group.

Conclusion

Mutation of CYP2D6*10 allele accounts for the highest proportion of vivax malaria cases in Yunnan Province. The mutations of c. 886C > T and CYP2D6*2 allele, which correspond to impaired PQ metabolizer phenotype, are most closely related to the relapse of vivax malaria. In addition, the genotype *4/*4 with null CYP2D6 enzyme function was only detected in the SR group. These results reveal the risk of defected CYP2D6 enzyme activity that diminishes the therapeutic effect of primaquine on vivax malaria.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-021-03685-3.

Combining structure and genomics to understand antimicrobial resistance

Antimicrobials against bacterial, viral and parasitic pathogens have transformed human and animal health. Nevertheless, their widespread use (and misuse) has led to the emergence of antimicrobial resistance (AMR) which poses a potentially catastrophic threat to public health and animal husbandry. There are several routes, both intrinsic and acquired, by which AMR can develop. One major route is through non-synonymous single nucleotide polymorphisms (nsSNPs) in coding regions. Large scale genomic studies using high-throughput sequencing data have provided powerful new ways to rapidly detect and respond to such genetic mutations linked to AMR. However, these studies are limited in their mechanistic insight. Computational tools can rapidly and inexpensively evaluate the effect of mutations on protein function and evolution. Subsequent insights can then inform experimental studies, and direct existing or new computational methods. Here we review a range of sequence and structure-based computational tools, focussing on tools successfully used to investigate mutational effect on drug targets in clinically important pathogens, particularly Mycobacterium tuberculosis. Combining genomic results with the biophysical effects of mutations can help reveal the molecular basis and consequences of resistance development. Furthermore, we summarise how the application of such a mechanistic understanding of drug resistance can be applied to limit the impact of AMR.

mCSM-membrane: predicting the effects of mutations on transmembrane proteins

Significant efforts have been invested into understanding and predicting the molecular consequences of mutations in protein coding regions, however nearly all approaches have been developed using globular, soluble proteins. These methods have been shown to poorly translate to studying the effects of mutations in membrane proteins. To fill this gap, here we report, mCSM-membrane, a user-friendly web server that can be used to analyse the impacts of mutations on membrane protein stability and the likelihood of them being disease associated. mCSM-membrane derives from our well-established mutation modelling approach that uses graph-based signatures to model protein geometry and physicochemical properties for supervised learning. Our stability predictor achieved correlations of up to 0.72 and 0.67 (on cross validation and blind tests, respectively), while our pathogenicity predictor achieved a Matthew's Correlation Coefficient (MCC) of up to 0.77 and 0.73, outperforming previously described methods in both predicting changes in stability and in identifying pathogenic variants. mCSM-membrane will be an invaluable and dedicated resource for investigating the effects of single-point mutations on membrane proteins through a freely available, user friendly web server at http://biosig.unimelb.edu.au/mcsm_membrane.

mCSM-AB2: guiding rational antibody design using graph-based signatures

MOTIVATION

A lack of accurate computational tools to guide rational mutagenesis has made affinity maturation a recurrent challenge in antibody (Ab) development. We previously showed that graph-based signatures can be used to predict the effects of mutations on Ab binding affinity.

RESULTS

Here we present an updated and refined version of this approach, mCSM-AB2, capable of accurately modelling the effects of mutations on Ab-antigen binding affinity, through the inclusion of evolutionary and energetic terms. Using a new and expanded database of over 1800 mutations with experimental binding measurements and structural information, mCSM-AB2 achieved a Pearson's correlation of 0.73 and 0.77 across training and blind tests, respectively, outperforming available methods currently used for rational Ab engineering.

AVAILABILITY AND IMPLEMENTATION

mCSM-AB2 is available as a user-friendly and freely accessible web server providing rapid analysis of both individual mutations or the entire binding interface to guide rational antibody affinity maturation at http://biosig.unimelb.edu.au/mcsm_ab2.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

DynaMut2: Assessing changes in stability and flexibility upon single and multiple point missense mutations

Predicting the effect of missense variations on protein stability and dynamics is important for understanding their role in diseases, and the link between protein structure and function. Approaches to estimate these changes have been proposed, but most only consider single‐point missense variants and a static state of the protein, with those that incorporate dynamics are computationally expensive. Here we present DynaMut2, a web server that combines Normal Mode Analysis (NMA) methods to capture protein motion and our graph‐based signatures to represent the wildtype environment to investigate the effects of single and multiple point mutations on protein stability and dynamics. DynaMut2 was able to accurately predict the effects of missense mutations on protein stability, achieving Pearson's correlation of up to 0.72 (RMSE: 1.02 kcal/mol) on a single point and 0.64 (RMSE: 1.80 kcal/mol) on multiple‐point missense mutations across 10‐fold cross‐validation and independent blind tests. For single‐point mutations, DynaMut2 achieved comparable performance with other methods when predicting variations in Gibbs Free Energy (ΔΔG) and in melting temperature (ΔT_m). We anticipate our tool to be a valuable suite for the study of protein flexibility analysis and the study of the role of variants in disease. DynaMut2 is freely available as a web server and API at http://biosig.unimelb.edu.au/dynamut2.

Identifying Genotype-Phenotype Correlations via Integrative Mutation Analysis

Mutations in protein-coding regions can lead to large biological changes and are associated with genetic conditions, including cancers and Mendelian diseases, as well as drug resistance. Although whole genome and exome sequencing help to elucidate potential genotype-phenotype correlations, there is a large gap between the identification of new variants and deciphering their molecular consequences. A comprehensive understanding of these mechanistic consequences is crucial to better understand and treat diseases in a more personalized and effective way. This is particularly relevant considering estimates that over 80% of mutations associated with a disease are incorrectly assumed to be causative. A thorough analysis of potential effects of mutations is required to correctly identify the molecular mechanisms of disease and enable the distinction between disease-causing and non-disease-causing variation within a gene. Here we present an overview of our integrative mutation analysis platform, which focuses on refining the current genotype-phenotype correlation methods by using the wealth of protein structural information.

Genetic Polymorphisms Complicate COVID-19 Therapy: Pivotal Role of HO-1 in Cytokine Storm

Coronaviruses are very large RNA viruses that originate in animal reservoirs and include severe acute respiratory distress syndrome (SARS) and Middle East respiratory syndrome (MERS) and other inconsequential coronaviruses from human reservoirs like the common cold. SARS-CoV-2, the virus that causes COVID-19 and is believed to originate from bat, quickly spread into a global pandemic. This RNA virus has a special affinity for porphyrins. It invades the cell at the angiotensin converting enzyme-2 (ACE-2) receptor and binds to hemoproteins, resulting in a severe systemic inflammatory response, particularly in high ACE-2 organs like the lungs, heart, and kidney, resulting in systemic disease. The inflammatory response manifested by increased cytokine levels and reactive oxygen species results in inhibition of heme oxygenase (HO-1), with a subsequent loss of cytoprotection. This has been seen in other viral illness like human immunodeficiency virus (HIV), Ebola, and SARS/MERS. There are a number of medications that have been tried with some showing early clinical promise. This illness disproportionately affects patients with obesity, a chronic inflammatory disease with a baseline excess of cytokines. The majority of the medications used in the treatment of COVID-19 are metabolized by cytochrome P450 (CYP) enzymes, primarily CYP2D6. This is further complicated by genetic polymorphisms of CYP2D6, HO-1, ACE, and ACE-2. There is a potential role for HO-1 upregulation to treat/prevent cytokine storm. Current therapy must focus on antivirals and heme oxygenase upregulation. Vaccine development will be the only magic bullet.

Influence of CYP2C8, CYP3A4, and CYP3A5 Host Genotypes on Early Recurrence of Plasmodium vivax

Cytochrome P450 (CYP) enzymes are involved in the biotransformation of chloroquine (CQ), but the role of the different profiles of metabolism of this drug in relation to Plasmodium vivax recurrences has not been properly investigated. To investigate the influence of the CYP genotypes associated with CQ metabolism on the rates of P. vivax early recurrences, a case-control study was carried out. The cases included patients presenting with an early recurrence (CQ-recurrent individuals), defined as a recurrence during the first 28 days after initial infection and plasma concentrations of CQ plus desethylchloroquine (DCQ; the major CQ metabolite) higher than 100 ng/ml. A control group with no parasite recurrence over the follow-up (the CQ-responsive group) was also included. CQ and DCQ plasma levels were measured on day 28. CQ-metabolizing CYP (CYP2C8, CYP3A4, and CYP3A5) genotypes were determined by real-time PCR. An ex vivo study was conducted to verify the efficacy of CQ and DCQ against P. vivax isolates. The frequency of alleles associated with normal and slow metabolism was similar between the cases and the controls for the CYP2C8 (odds ratio [OR] = 1.45, 95% confidence interval [CI] = 0.51 to 4.14, P = 0.570), CYP3A4 (OR = 2.38, 95% CI = 0.92 to 6.19, P = 0.105), and CYP3A5 (OR = 4.17, 95% CI = 0.79 to 22.04, P = 1.038) genes. DCQ levels were higher than CQ levels, regardless of the genotype. Regarding the DCQ/CQ ratio, there was no difference between groups or between those patients who had a normal genotype and those patients who had a mutant genotype. DCQ and CQ showed similar efficacy ex vivo CYP genotypes had no influence on early recurrence rates. The similar efficacy of CQ and DCQ ex vivo could explain the absence of therapeutic failure, despite the presence of alleles associated with slow metabolism.

Novel Insights into Plasmodium vivax Therapeutic Failure: CYP2D6 Activity and Time of Exposure to Malaria Modulate the Risk of Recurrence

Plasmodium vivax relapse is one of the major causes of sustained global malaria transmission. Primaquine (PQ) is the only commercial drug available to prevent relapses, and its efficacy is dependent on metabolic activation by cytochrome P450 2D6 (CYP2D6). Impaired CYP2D6 function, caused by allelic polymorphisms, leads to the therapeutic failure of PQ as a radical cure for P. vivax malaria. Here, we hypothesized that the host immune response to malaria parasites modulates susceptibility to P. vivax recurrences in association with CYP2D6 activity. We performed a 10-year retrospective study by genotyping CYP2D6 polymorphisms in 261 malaria-exposed individuals from the Brazilian Amazon. The immune responses against a panel of P. vivax blood-stage antigens were evaluated by serological assays. We confirmed our previous findings, which indicated an association between impaired CYP2D6 activity and a higher risk of multiple episodes of P. vivax recurrence (risk ratio, 1.75; 95% confidence interval [CI], 1.2 to 2.6; P = 0.0035). An important finding was a reduction of 3% in the risk of recurrence (risk ratio, 0.97; 95% CI, 0.96 to 0.98; P < 0.0001) per year of malaria exposure, which was observed for individuals with both reduced and normal CYP2D6 activity. Accordingly, subjects with long-term malaria exposure and persistent antibody responses to various antigens showed fewer episodes of malaria recurrence. Our findings have direct implications for malaria control, since it was shown that nonimmune individuals who do not respond adequately to treatment due to reduced CYP2D6 activity may present a significant challenge for sustainable progress toward P. vivax malaria elimination.

Structure guided prediction of Pyrazinamide resistance mutations in pncA

Pyrazinamide plays an important role in tuberculosis treatment; however, its use is complicated by side-effects and challenges with reliable drug susceptibility testing. Resistance to pyrazinamide is largely driven by mutations in pyrazinamidase (pncA), responsible for drug activation, but genetic heterogeneity has hindered development of a molecular diagnostic test. We proposed to use information on how variants were likely to affect the 3D structure of pncA to identify variants likely to lead to pyrazinamide resistance. We curated 610 pncA mutations with high confidence experimental and clinical information on pyrazinamide susceptibility. The molecular consequences of each mutation on protein stability, conformation, and interactions were computationally assessed using our comprehensive suite of graph-based signature methods, mCSM. The molecular consequences of the variants were used to train a classifier with an accuracy of 80%. Our model was tested against internationally curated clinical datasets, achieving up to 85% accuracy. Screening of 600 Victorian clinical isolates identified a set of previously unreported variants, which our model had a 71% agreement with drug susceptibility testing. Here, we have shown the 3D structure of pncA can be used to accurately identify pyrazinamide resistance mutations. SUSPECT-PZA is freely available at: http://biosig.unimelb.edu.au/suspect_pza/.

A Comprehensive Computational Platform to Guide Drug Development Using Graph-Based Signature Methods

High-throughput computational techniques have become invaluable tools to help increase the overall success, process efficiency, and associated costs of drug development. By designing ligands tailored to specific protein structures in a disease of interest, an understanding of molecular interactions and ways to optimize them can be achieved prior to chemical synthesis. This understanding can help direct crucial chemical and biological experiments by maximizing available resources on higher quality leads. Moreover, predicting molecular binding affinity within specific biological contexts, as well as ligand pharmacokinetics and toxicities, can aid in filtering out redundant leads early on within the process. We describe a set of computational tools which can aid in drug discovery at different stages, from hit identification (EasyVS) to lead optimization and candidate selection (CSM-lig, mCSM-lig, Arpeggio, pkCSM). Incorporating these tools along the drug development process can help ensure that candidate leads are chemically and biologically feasible to become successful and tractable drugs.

Pharmacokinetics/pharmacodynamics of chloroquine and artemisinin-based combination therapy with primaquine

Background

Activation of hypnozoites of vivax malaria causes multiple clinical relapses, which contribute to the Plasmodium vivax burden and continuing transmission. Artemisinin-based combination therapy (ACT) is effective against blood-stage P. vivax but requires co-administration with primaquine to achieve radical cure. The therapeutic efficacy of primaquine depends on the generation of a therapeutically active metabolite via cytochrome P450 2D6 (CYP2D6). Impaired CYP2D6 metabolism has been associated with primaquine treatment failure. This study investigated the association between impaired CYP2D6 genotypes, drug-exposure to the long-acting ACT component (schizonticidal drugs) and tolerance and efficacy.

Methods

Adult patients with acute vivax malaria were enrolled in a recently completed trial and treated with artesunate–mefloquine, chloroquine or artemether–lumefantrine. All received concomitant primaquine (0.5 mg/kg/day for 7–9 days). The association between efficacy and safety and drug exposure was explored using area-under-the-curve (AUC) and half-life (t_1/2) estimates obtained by non-compartmental analysis of the long half-life drugs. Parasite recurrences by day 63 were categorized as related relapses or re-infections/unrelated hypnozoite activation by genotyping three microsatellite loci and two polymorphic loci of merozoite surface antigen-1. The CYP2D6 genotype was identified with Taqman assays by real-time PCR to 9 polymorphisms (8 SNPs and one deletion). Impaired CYP2D6 activity was inferred using the Activity Score System.

Results

Most recurrences in the ASMQ (67%), CQ (80%) and AL (85%) groups were considered related relapses. Eight of nine (88.9%) of the patients with impaired CYP2D6 activity relapsed with related parasite compared to 18/25 (72%) with normal activity (RR = 1.23, 0.88; 1.72, p = 0.40). There were no associations between the measured PK parameters and recurrence. Patients with longer chloroquine half-lives had more pruritus (RR = 1.09, 1.03; 1.14, p = 0.001). Higher CQ AUCs were associated with reduced falls in haemoglobin by day 14 (Coef − 0.02, − 0.005; − 0.03, p = 0.01). All regimens were well tolerated.

Conclusion

Genotyping of P. vivax showed that activation of related (homologous) hypnozoites was the most frequent cause of recurrence. The high proportion of the impaired CYP2D6 activity among patients with recurrent infections suggests that slow primaquine metabolism might influence related relapse rates in Brazil among patients receiving primaquine for radical cure, although confirmatory studies are needed. There was no association between drug exposure of the long-acting ACT component (schizonticidal drugs) and risk of related relapse. ACT was well tolerated. These results provide further re-assurance about the safety and efficacy of ACT when combined with short course primaquine to treat uncomplicated malaria vivax in Brazil.

Trial registration RBR-79s56s (http://www.ensaiosclinicos.gov.br/rg/RBR-79s56s/)

mCSM-NA: predicting the effects of mutations on protein-nucleic acids interactions

Over the past two decades, several computational methods have been proposed to predict how missense mutations can affect protein structure and function, either by altering protein stability or interactions with its partners, shedding light into potential molecular mechanisms giving rise to different phenotypes. Effectively and efficiently predicting consequences of mutations on protein–nucleic acid interactions, however, remained until recently a great and unmet challenge. Here we report an updated webserver for mCSM–NA, the only scalable method we are aware of capable of quantitatively predicting the effects of mutations in protein coding regions on nucleic acid binding affinities. We have significantly enhanced the original method by including a pharmacophore modelling and information of nucleic acid properties into our graph-based signatures, considering the reverse mutation and by using a refined, more reliable data set, based on a new release of the ProNIT database, which has significantly improved the reliability and applicability of the methodology. Our new predictive model was capable of achieving a correlation coefficient of up to 0.70 on cross-validation and 0.68 on blind-tests, outperforming its previous version. The server is freely available via a user-friendly web interface at: http://structure.bioc.cam.ac.uk/mcsm_na.

SDM: a server for predicting effects of mutations on protein stability

Here, we report a webserver for the improved SDM, used for predicting the effects of mutations on protein stability. As a pioneering knowledge-based approach, SDM has been highlighted as the most appropriate method to use in combination with many other approaches. We have updated the environment-specific amino-acid substitution tables based on the current expanded PDB (a 5-fold increase in information), and introduced new residue-conformation and interaction parameters, including packing density and residue depth. The updated server has been extensively tested using a benchmark containing 2690 point mutations from 132 different protein structures. The revised method correlates well against the hypothetical reverse mutations, better than comparable methods built using machine-learning approaches, highlighting the strength of our knowledge-based approach for identifying stabilising mutations. Given a PDB file (a Protein Data Bank file format containing the 3D coordinates of the protein atoms), and a point mutation, the server calculates the stability difference score between the wildtype and mutant protein. The server is available at http://structure.bioc.cam.ac.uk/sdm2

Monitoring the Efficacy of Chloroquine-Primaquine Therapy for Uncomplicated Plasmodium vivax Malaria in the Main Transmission Hot Spot of Brazil

Emerging Plasmodium vivax resistance to chloroquine (CQ) may undermine malaria elimination efforts in South America. CQ-resistant P. vivax has been found in the major port city of Manaus but not in the main malaria hot spots across the Amazon Basin of Brazil, where CQ is routinely coadministered with primaquine (PQ) for radical cure of vivax malaria. Here we randomly assigned 204 uncomplicated vivax malaria patients from Juruá Valley, northwestern Brazil, to receive either sequential (arm 1) or concomitant (arm 2) CQ-PQ treatment. Because PQ may synergize the blood schizontocidal effect of CQ and mask low-level CQ resistance, we monitored CQ-only efficacy in arm 1 subjects, who had PQ administered only at the end of the 28-day follow-up. We found adequate clinical and parasitological responses in all subjects assigned to arm 2. However, 2.2% of arm 1 patients had microscopy-detected parasite recrudescences at day 28. When PCR-detected parasitemias at day 28 were considered, response rates decreased to 92.1% and 98.8% in arms 1 and 2, respectively. Therapeutic CQ levels were documented in 6 of 8 recurrences, consistent with true CQ resistance in vivo In contrast, ex vivo assays provided no evidence of CQ resistance in 49 local P. vivax isolates analyzed. CQ-PQ coadministration was not found to potentiate the antirelapse efficacy of PQ over 180 days of surveillance; however, we suggest that larger studies are needed to examine whether and how CQ-PQ interactions, e.g., CQ-mediated inhibition of PQ metabolism, modulate radical cure efficacy in different P. vivax-infected populations. (This study has been registered at ClinicalTrials.gov under identifier NCT02691910.).

Cytochrome P450 2D6 profiles and their relationship with outcomes of primaquine anti-relapse therapy in Australian Defence Force personnel deployed to Papua New Guinea and East Timor

Background

Primaquine, an 8-aminoquinoline with anti-hypnozoite activity against Plasmodium vivax, is metabolized by human cytochrome P450 2D6 (CYP2D6) to its active metabolite. Human CYP2D6 activities may influence the metabolism of primaquine and the risk of experiencing Plasmodium relapses following primaquine anti-relapse therapies (PART). In this study, the CYP2D6 profile and its relationship with outcomes of PART in Australian Defence Force (ADF) personnel is retrospectively investigated.

Methods

Genomic DNA was isolated from stored and de-identified serum or blood samples from ADF personnel deployed on peacekeeping duties to Papua New Guinea (PNG) (1999) and East Timor (1999–2000) who received PART before returning to Australia and after experiencing relapses. CYP2D6 allelic type was determined by PCR and Sanger sequencing. CYP2D6 allele frequency, predicted phenotypes and activity scores were compared among personnel who did not experience P. vivax (ADF-NR, n = 48) and those who experience at least one (ADF-R, n = 109) relapse, as well as between those who experienced 1 (n = 79), 2 (n = 21) and 3–5 (n = 9) relapses within the ADF-R group.

Results

16 CYP2D6 alleles were observed in 157 ADF personnel. Alleles *1, *4, *2 and *41 were major alleles (> 5%). The CYP2D6 allele frequency profile in the ADF-NR group matched that of a European population. There was an increased proportion of non-functional CYP2D6 alleles in the ADF-R group compared to the European population and ADF-NR group. However, frequencies of predicted CYP2D6 phenotype and activity score were not different between the ADF-R and ADF-NR groups, nor among sub-groups experiencing multiple relapses within the ADF-R group.

Conclusions

CYP2D6 phenotype or activity score based on the allele classification was not a major contributor to P. vivax relapse in this ADF cohort. Other factors such as adherence and/or parasite tolerance to primaquine are likely contributing factors to P. vivax relapses in this cohort.

Synthesis and Structure-Activity relationship of 1-(5-isoquinolinesulfonyl)piperazine analogues as inhibitors of Mycobacterium tuberculosis IMPDH

Tuberculosis (TB) is a major infectious disease associated increasingly with drug resistance. Thus, new anti-tubercular agents with novel mechanisms of action are urgently required for the treatment of drug-resistant TB. In prior work, we identified compound 1 (cyclohexyl(4-(isoquinolin-5-ylsulfonyl)piperazin-1-yl)methanone) and showed that its anti-tubercular activity is attributable to inhibition of inosine-5′-monophosphate dehydrogenase (IMPDH) in Mycobacterium tuberculosis. In the present study, we explored the structure–activity relationship around compound 1 by synthesizing and evaluating the inhibitory activity of analogues against M. tuberculosis IMPDH in biochemical and whole-cell assays. X-ray crystallography was performed to elucidate the mode of binding of selected analogues to IMPDH. We establish the importance of the cyclohexyl, piperazine and isoquinoline rings for activity, and report the identification of an analogue with IMPDH-selective activity against a mutant of M. tuberculosis that is highly resistant to compound 1. We also show that the nitrogen in urea analogues is required for anti-tubercular activity and identify benzylurea derivatives as promising inhibitors that warrant further investigation.

•

Forty-eight analogues of 1-(5-isoquinolinesulfonyl)piperazine were synthesized.

•

Biochemical, whole-cell, and X-ray studies were performed to elucidate the IMPDH inhibition.

•

Piperazine and isoquinoline rings were essential for target-selective whole-cell activity.

•

Compound 47 showed improved IC₅₀ against the MtbIMPDH and maintained on-target whole-cell activity.

•

Compound 21 showed activity against IMPDH in both wild type M. tuberculosis and a resistant mutant of compound 1.

Pharmacogenomics research and clinical implementation in Brazil

We searched PubMed entries and the Lattes database of Brazilian Pharmacogenetics Network investigators, for pharmacogenetic/genomic (PGx) studies in the Brazilian population, focusing on the drugs and genes included in the Clinical Pharmacogenetics Implementation Consortium (CPIC) guidelines. Warfarin was the most extensively studied drug in a PGx context: a genomewide association study targeting warfarin stable dose identified significant signals in VKORC1 and CYP2C9, several PGx dosing algorithms were developed based on these and other genes, and the implications of population admixture on extrapolation of dosing recommendations in the CPIC guidelines were examined. A study in renal transplanted patients disclosed association of CYP3A5*6 and CYP3A5*7 with tacrolimus dosing, which led to addition of these variants to CYP3A5*3 in the CPIC tacrolimus guideline. Studies verified predisposition of HIV-positive carriers of UGT1A1*28 to severe atazanavir-induced hyperbilirubinaemia, intolerance to 5-fluorouracyl in gastrointestinal cancer patients with deleterious DPYD variants, failure of HCV-infected carriers of IFNL3 rs12979860 to obtain a sustained viral response to PEG-IFN-α, and hypersensitivity reactions to abacavir in HIV-positive carriers of HLA-B*57:01. No prospective analyses of drug therapy outcomes or cost-effectiveness assessments of PGx-guided therapy were found. In conclusion, the limited adoption of PGx-informed drug prescription in Brazil reflects combination of recognized barriers to PGx implementation worldwide plus factors specific to the Brazilian population. The latter include rarity/absence of genetic variants on which international PGx guidelines are based (eg HLA-B*15.02 for phenytoin and carbamazepine) and the caveat of extrapolating to the admixed Brazilian population, guidelines based on categorical variables, such as continental ancestry (eg warfarin guidelines), "race" or ethnicity.

Single-Dose Tafenoquine to Prevent Relapse of Plasmodium vivax Malaria

BACKGROUND

Treatment of Plasmodium vivax malaria requires the clearing of asexual parasites, but relapse can be prevented only if dormant hypnozoites are cleared from the liver (a treatment termed "radical cure"). Tafenoquine is a single-dose 8-aminoquinoline that has recently been registered for the radical cure of P. vivax.

METHODS

This multicenter, double-blind, double-dummy, parallel group, randomized, placebo-controlled trial was conducted in Ethiopia, Peru, Brazil, Cambodia, Thailand, and the Philippines. We enrolled 522 patients with microscopically confirmed P. vivax infection (>100 to <100,000 parasites per microliter) and normal glucose-6-phosphate dehydrogenase (G6PD) activity (with normal activity defined as ≥70% of the median value determined at each trial site among 36 healthy male volunteers who were otherwise not involved in the trial). All patients received a 3-day course of chloroquine (total dose of 1500 mg). In addition, patients were assigned to receive a single 300-mg dose of tafenoquine on day 1 or 2 (260 patients), placebo (133 patients), or a 15-mg dose of primaquine once daily for 14 days (129 patients). The primary outcome was the Kaplan-Meier estimated percentage of patients who were free from recurrence at 6 months, defined as P. vivax clearance without recurrent parasitemia.

RESULTS

In the intention-to-treat population, the percentage of patients who were free from recurrence at 6 months was 62.4% in the tafenoquine group (95% confidence interval [CI], 54.9 to 69.0), 27.7% in the placebo group (95% CI, 19.6 to 36.6), and 69.6% in the primaquine group (95% CI, 60.2 to 77.1). The hazard ratio for the risk of recurrence was 0.30 (95% CI, 0.22 to 0.40) with tafenoquine as compared with placebo (P<0.001) and 0.26 (95% CI, 0.18 to 0.39) with primaquine as compared with placebo (P<0.001). Tafenoquine was associated with asymptomatic declines in hemoglobin levels, which resolved without intervention.

CONCLUSIONS

Single-dose tafenoquine resulted in a significantly lower risk of P. vivax recurrence than placebo in patients with phenotypically normal G6PD activity. (Funded by GlaxoSmithKline and Medicines for Malaria Venture; DETECTIVE ClinicalTrials.gov number, NCT01376167 .).

Exploring Protein Supersecondary Structure Through Changes in Protein Folding, Stability, and Flexibility

The ability to predict how mutations affect protein structure, folding, and flexibility can elucidate the molecular mechanisms leading to disruption of supersecondary structures, the emergence of phenotypes, as well guiding rational protein engineering. The advent of fast and accurate computational tools has enabled us to comprehensively explore the landscape of mutation effects on protein structures, prioritizing mutations for rational experimental validation.Here we describe the use of two complementary web-based in silico methods, DUET and DynaMut, developed to infer the effects of mutations on folding, stability, and flexibility and how they can be used to explore and interpret these effects on protein supersecondary structures.

Real time PCR detection of common CYP2D6 genetic variants and its application in a Karen population study

BACKGROUND

Plasmodium vivax malaria is characterized by relapses arising from the hypnozoite stages in the liver. The only currently registered drug for radical treatment to prevent relapse is primaquine. Primaquine, a prodrug, requires metabolism through the liver cytochrome CYP2D6 isoenzyme to its active metabolite. Mutations in the CYP2D6 gene may thus affect primaquine efficacy. A SNPs genotyping technique was developed to characterize the CYP2D6 genetic variants and tested this in the patients with Plasmodium vivax infection collected in a Karen population on the Thailand-Myanmar border, where P. vivax malaria is endemic.

METHODS

Direct sequencing of PCR-reamplified products (DSP) was used to uncover exonic CYP2D6 sequence variations. Subsequently, an allele-specific oligonucleotide probe real-time SNPs genotyping (ASO) assay was developed for rapid detection of the four clinically relevant CYP2D6 variants occurring in this population. These two in-house developed assays were used to genotype CYP2D6 mutations in blood samples obtained from 70 Karen adults.

RESULTS

Results showed a high degree of concordance between the DSP and ASO methods. Six CYP2D6 point mutations were identified within the Karen population: C100T, C1039T, G1661C, G1846A, C2850T and G4180C, at frequencies of 0.43, 0.43, 0.76, 0.02, 0.32 and 0.76, respectively. The CYP2D6*2, *4, *5, *10 and *36 allelic frequencies were 0.33, 0.02, 0.03, 0.40 and 0.01, respectively. Alleles conferring an intermediate CYP2D6 metabolizer phenotype comprised 46% of the total number of alleles.

CONCLUSION

The newly developed ASO assay is a reliable and rapid tool for large-scale CYP2D6 genotyping. The high frequency of the CYP2D6*10 allele in the Karen population warrants further assessment of its association with the radical curative efficacy of primaquine.

Frequent transmission of the Mycobacterium tuberculosis Beijing lineage and positive selection for the EsxW Beijing variant in Vietnam

To examine the transmission dynamics of Mycobacterium tuberculosis (Mtb) isolated from tuberculosis patients in Ho Chi Minh City, Vietnam, we sequenced the whole genomes of 1,635 isolates and compared these with 3,144 isolates from elsewhere. The data identify an underlying burden of disease caused by the endemic Mtb lineage 1 associated with the activation of long-term latent infection, and a threefold higher burden associated with the more recently introduced Beijing lineage and lineage 4 Mtb strains. We find that Beijing lineage Mtb is frequently transferred between Vietnam and other countries, and detect higher levels of transmission of Beijing lineage strains within this host population than the endemic lineage 1 Mtb. Screening for parallel evolution of Beijing lineage-associated SNPs in other Mtb lineages as a signal of positive selection, we identify an alteration in the ESX-5 type VII-secreted protein EsxW, which could potentially contribute to the enhanced transmission of Beijing lineage Mtb in Vietnamese and other host populations.

Frequent transmission of the Mycobacterium tuberculosis Beijing lineage and positive selection for the EsxW Beijing variant in Vietnam

To examine the transmission dynamics of Mycobacterium tuberculosis (Mtb) isolated from tuberculosis patients in Ho Chi Minh City, Vietnam, we sequenced the whole genomes of 1,635 isolates and compared these with 3,144 isolates from elsewhere. The data identify an underlying burden of disease caused by the endemic Mtb lineage 1 associated with the activation of long-term latent infection, and a threefold higher burden associated with the more recently introduced Beijing lineage and lineage 4 Mtb strains. We find that Beijing lineage Mtb is frequently transferred between Vietnam and other countries, and detect higher levels of transmission of Beijing lineage strains within this host population than the endemic lineage 1 Mtb. Screening for parallel evolution of Beijing lineage-associated SNPs in other Mtb lineages as a signal of positive selection, we identify an alteration in the ESX-5 type VII-secreted protein EsxW, which could potentially contribute to the enhanced transmission of Beijing lineage Mtb in Vietnamese and other host populations.

Structural Implications of Mutations Conferring Rifampin Resistance in Mycobacterium leprae

The rpoB gene encodes the β subunit of RNA polymerase holoenzyme in Mycobacterium leprae (M. leprae). Missense mutations in the rpoB gene were identified as etiological factors for rifampin resistance in leprosy. In the present study, we identified mutations corresponding to rifampin resistance in relapsed leprosy cases from three hospitals in southern India which treat leprosy patients. DNA was extracted from skin biopsies of 35 relapse/multidrug therapy non-respondent leprosy cases, and PCR was performed to amplify the 276 bp rifampin resistance-determining region of the rpoB gene. PCR products were sequenced, and mutations were identified in four out of the 35 cases at codon positions D441Y, D441V, S437L and H476R. The structural and functional effects of these mutations were assessed in the context of three-dimensional comparative models of wild-type and mutant M. leprae RNA polymerase holoenzyme (RNAP), based on the recently solved crystal structures of RNAP of Mycobacterium tuberculosis, containing a synthetic nucleic acid scaffold and rifampin. The resistance mutations were observed to alter the hydrogen-bonding and hydrophobic interactions of rifampin and the 5' ribonucleotide of the growing RNA transcript. This study demonstrates that rifampin-resistant strains of M. leprae among leprosy patients in southern India are likely to arise from mutations that affect the drug-binding site and stability of RNAP.

Biological concepts in recurrent Plasmodium vivax malaria

A curious aspect of the evolution of the hypnozoite theory of malarial relapse is its transmogrification from theory into 'fact', this being of historical, linguistic, scientific and sociological interest. As far as it goes, the hypnozoite explanation for relapse is almost certainly correct. I contend, however, that many of the genotypically homologous, non-reinfection, relapse-like Plasmodium vivax recurrences that researchers ascribe to hypnozoite activation are probably hypnozoite-independent. Indeed, some malariologists are starting to recognize that homologous P. vivax recurrences have most likely been overattributed to activation of hypnozoites. Hitherto identified, non-hypnozoite, possible plasmodial sources of recurrence that must be considered, besides circulating erythrocytic stages, include parasites in splenic dendritic cells, other cells in the spleen (in addition to infected erythrocytes there), bone marrow (importantly) and the skin. I argue that we need to take into account the possibility of a dual or multiple extra-vascular origin of P. vivax non-reinfection recurrences, not arbitrarily discount it. The existence of a P. vivax reservoir(s) is a topical subject and one of practical importance for malaria eradication. Pertinent drug-associated matters are also discussed, as is the dormancy-related significance of clues provided by blood-stage-induced malarial infection.

CYP2D6 activity and the risk of recurrence of Plasmodium vivax malaria in the Brazilian Amazon: a prospective cohort study

Background

CYP2D6 pathway mediates the activation of primaquine into active metabolite(s) in hepatocytes. CYP2D6 is highly polymorphic, encoding CYP2D6 isoforms with normal, reduced, null or increased activity. It is hypothesized that Plasmodium vivax malaria patients with defective CYP2D6 function would be at increased risk for primaquine failure to prevent recurrence. The aim of this study was to investigate the association of CYP2D6 polymorphisms and inferred CYP2D6 phenotypes with malaria recurrence in patients from the Western Brazilian Amazon, following chloroquine/primaquine combined therapy.

Methods

The prospective cohort consisted of P. vivax malaria patients who were followed for 6 months after completion of the chloroquine/primaquine therapy. Recurrence was defined as one or more malaria episodes, 28–180 days after the initial episode. Genotyping for nine CYP2D6 SNPs and copy number variation was performed using TaqMan assays in a Fast 7500 Real-Time System. CYP2D6 star alleles (haplotypes), diplotypes and CYP2D6 phenotypes were inferred, and the activity score system was used to define the functionality of the CYP2D6 diplotypes. CYP2D6 activity scores (AS) were dichotomized at ≤ 1 (gPM, gIM and gNM-S phenotypes) and ≥ 1.5 (gNM-F and gUM phenotypes).

Results

Genotyping was successfully performed in 190 patients (44 with recurrence and 146 without recurrences). Recurrence incidence was higher in individuals presenting reduced activity CYP2D6 phenotypes (adjusted relative risk = 1.89, 95% CI 1.01–3.70; p = 0.049). Attributable risk and population attributable fraction were 11.5 and 9.9%, respectively. The time elapsed from the first P. vivax malaria episode until the recurrence did not differ between patients with AS of ≤ 1 versus ≥ 1.5 (p = 0.917).

Conclusions

The results suggest that CYP2D6 polymorphisms are associated with increased risk of recurrence of vivax malaria, following chloroquine–primaquine combined therapy. This association is interpreted as the result of reduced conversion of primaquine into its active metabolites in patients with reduced CYP2D6 enzymatic activity.

Correction: Variation in Human Cytochrome P-450 Drug-Metabolism Genes: A Gateway to the Understanding of Plasmodium vivax Relapses

[This corrects the article DOI: 10.1371/journal.pone.0160172.].

Prediction and Optimization of Pharmacokinetic and Toxicity Properties of the Ligand

A crucial factor for the approval and success of any drug is how it behaves in the body. Many drugs, however, do not reach the market due to poor efficacy or unacceptable side effects. It is therefore important to take these into consideration early in the drug development process, both in the prioritization of potential hits, and optimization of lead compounds. In silico approaches offer a cost and time-effective approach to rapidly screen and optimize pharmacokinetic and toxicity properties. Here we demonstrate the use of the comprehensive analysis system pkCSM, to allow early identification of potential problems, prioritization of hits, and optimization of leads.

SLCO1A2, SLCO1B1 and SLCO2B1 polymorphisms influences chloroquine and primaquine treatment in Plasmodium vivax malaria

AIM

The association of transporters gene polymorphisms with chloroquine/primaquine malaria treatment response was investigated in a Brazilian population.

PATIENTS & METHODS

Totally, 164 Plasmodium vivax malaria infected patients were included. Generalized estimating equations were performed to determine gene influences on parasitemia and/or gametocytemia clearance over treatment time.

RESULTS

Significant interaction between SLCO2B1 genotypes and treatment over time for parasitemia clearance rate on day 2 were observed (p FDR = 0.002). SLCO1A2 and SLCO1B1 gene treatment over time interactions were associated with gametocytemia clearance rate (p FDR = 0.018 and p FDR = 0.024). ABCB1, ABCC4 and SLCO1B3 were not associated with treatment response.

CONCLUSION

The present work presents the first pharmacogenetic report of an association between chloroquine/primaquine responses with OATP transporters.

Combating mutations in genetic disease and drug resistance: understanding molecular mechanisms to guide drug design

INTRODUCTION

Mutations introduce diversity into genomes, leading to selective changes and driving evolution. These changes have contributed to the emergence of many of the current major health concerns of the 21st century, from the development of genetic diseases and cancers to the rise and spread of drug resistance. The experimental systematic testing of all mutations in a system of interest is impractical and not cost-effective, which has created interest in the development of computational tools to understand the molecular consequences of mutations to aid and guide rational experimentation. Areas covered: Here, the authors discuss the recent development of computational methods to understand the effects of coding mutations to protein function and interactions, particularly in the context of the 3D structure of the protein. Expert opinion: While significant progress has been made in terms of innovative tools to understand and quantify the different range of effects in which a mutation or a set of mutations can give rise to a phenotype, a great gap still exists when integrating these predictions and drawing causality conclusions linking variants. This often requires a detailed understanding of the system being perturbed. However, as part of the drug development process it can be used preemptively in a similar fashion to pharmacokinetics predictions, to guide development of therapeutics to help guide the design and analysis of clinical trials, patient treatment and public health policy strategies.

Genomes, structural biology and drug discovery: combating the impacts of mutations in genetic disease and antibiotic resistance

For over four decades structural biology has been used to understand the mechanisms of disease, and structure-guided approaches have demonstrated clearly that they can contribute to many aspects of early drug discovery, both computationally and experimentally. Structure can also inform our understanding of impacts of mutations in human genetic diseases and drug resistance in cancers and infectious diseases. We discuss the ways that structural insights might be useful in both repurposing off-licence drugs and guide the design of new molecules that might be less susceptible to drug resistance in the future.

Assessment of copy number variation in genes related to drug resistance in Plasmodium vivax and Plasmodium falciparum isolates from the Brazilian Amazon and a systematic review of the literature

Background

Parasite resistance to anti-malarials represents a great obstacle for malaria elimination. The majority of studies have investigated the association between single-nucleotide polymorphisms (SNPs) and drug resistance; however, it is becoming clear that the copy number variation (CNV) is also associated with this parasite phenotype. To provide a baseline for molecular surveillance of anti-malarial drug resistance in the Brazilian Amazon, the present study characterized the genetic profile of both markers in the most common genes associated with drug resistance in Plasmodium falciparum and Plasmodium vivax isolates. Additionally, these data were compared to data published elsewhere applying a systematic review of the literature published over a 20-year time period.

Methods

The genomic DNA of 67 patients infected by P. falciparum and P. vivax from three Brazilian States was obtained between 2002 and 2012. CNV in P. falciparum multidrug resistance gene-1 (pfmdr1), GTP cyclohydrolase 1 (pfgch1) and P. vivax multidrug resistance gene-1 (pvmdr1) were assessed by real-time PCR assays. SNPs in the pfmdr1 and pfcrt genes were assessed by PCR–RFLP. A literature search for studies that analysed CNP in the same genes of P. falciparum and P. vivax was conducted between May 2014 and March 2017 across four databases.

Results

All analysed samples of P. falciparum carried only one copy of pfmdr1 or pfgch1. Although the pfcrt K76T polymorphism, a determinant of CQ resistance, was present in all samples genotyped, the pfmdr1 N86Y was absent. For P. vivax isolates, an amplification rate of 20% was found for the pvmdr1 gene. The results of the study are in agreement with the low amplification rates for pfmdr1 gene evidenced in the Americas and Africa, while higher rates have been described in Southeast Asia. For P. vivax, very low rates of amplification for pvmdr1 have been described worldwide, with exceptions in French Guiana, Cambodia, Thailand and Brazil.

Conclusions

The present study was the first to evaluate gch1 CNV in P. falciparum isolates from Brazil, showing an absence of amplification of this gene more than 20 years after the withdrawal of the Brazilian antifolates therapeutic scheme. Furthermore, the rate of pvmdr1 amplification was significantly higher than that previously reported for isolates circulating in Northern Brazil.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-017-1806-z) contains supplementary material, which is available to authorized users.