Skip to main content

Interaction analysis of statistically enriched mutations identified in Cameroon recombinant subtype CRF02_AG that can influence the development of Dolutegravir drug resistance mutations



The Integrase (IN) strand transfer inhibitor (INSTI), Dolutegravir (DTG), has been given the green light to form part of first-line combination antiretroviral therapy (cART) by the World Health Organization (WHO). DTG containing regimens have shown a high genetic barrier against HIV-1 isolates carrying specific resistance mutations when compared with other class of regimens.


We evaluated the HIV-1 CRF02_AG IN gene sequences from Cameroon for the presence of resistance-associated mutations (RAMs) against INSTIs and naturally occurring polymorphisms (NOPs), using study sequences (n = 20) and (n = 287) sequences data derived from HIV Los Alamos National Laboratory database. The possible impact of NOPs on protein structure caused by HIV-1 CRF02_AG variations was addressed within the context of a 3D model of the HIV-1 IN complex and interaction analysis was performed using PyMol to validate DTG binding to the Wild type and seven mutant structures.


We observed 12.8% (37/287) sequences to contain RAMs, with only 1.0% (3/287) of the sequences having major INSTI RAMs: T66A, Q148H, R263K and N155H. Of these,11.8% (34/287) of the sequences contained five different IN accessory mutations; namely Q95K, T97A, G149A, E157Q and D232N. NOPs occurred at a frequency of 66% on the central core domain (CCD) position, 44% on the C-terminal domain (CTD) position and 35% of the N-terminal domain (NTD) position. The interaction analysis revealed that DTG bound to DNA, 2MG ions and DDE motif residues for T66A, T97A, Q148H, N155H and R263K comparable to the WT structure. Except for accessory mutant structure E157Q, only one MG contact was made with DTG, while DTG had no MG ion contacts and no DDE motif residue contacts for structure D232N.


Our analysis indicated that all RAM’s that resulted in a change in the number of interactions with encompassing residues does not affect DTG binding, while accessory mutations E157Q and D232N could affect DTG binding leading to possible DTG resistance. However, further experimental validation is required to validate the in silico findings of our study.

Peer Review reports


Sub-Saharan Africa (SSA) remains one of the regions highly burdened by HIV infection at 70% of the global epidemic. SSA has a particularly high HIV-1 genetic diversity and it is documented that diverse subtypes may affect the clinical treatment outcome in patient management [1]. The HIV-1 CRF02_AG strain continues to be the predominant subtype causing majority of infections in Cameroon, while other strains, including groups N, O and P, account for a minor proportion of infections [2,3,4]. Furthermore, different mutational pathways account for subtype specific differences in drug resistance [5,6,7]. Additionally, other studies have also reported that natural occurring polymorphisms (NOP) which are associated with the occurrence of resistance to Integrase (IN) strand-transfer inhibitors (INSTIs) and IN activity, are subtype-dependent [6,7,8]. These subtype-specific polymorphic mutations in the IN gene have been shown to affect IN DNA binding affinity, in the presence of resistance-associated mutations (RAMs) [6,7,8]. Computational modelling of RAMs against INSTIs, across different HIV-1 subtypes compared to subtype B, showed that the presence of M50I in subtypes A and C, L74I in subtypes A and CRF02_AG, G163R in CRF01_AE, and V165I in subtypes F and CRF01_AE are associated with a lower genetic barrier to resistance in non-B clades [9]. Cameroon has seen a substantial reduction of HIV infection, since the introduction of combination antiretroviral therapy (cART), especially with the rolling-out of programmes like prevention of mother-to-child transmission (PMTCT) and the implementation 90–90-90 strategy to end the AIDS pandemic by 2030 [10]. The ability of the HIV-1 virus to mutate during therapy, can lead to the emergence of HIV-1 drug resistance and this necessitates the need for more effective cART regimens with higher genetic barriers [1]. In Cameroon, the HIV-1 drug resistance rates among cART-initiators stand at approximately, 10% of the Cameroon infected population [11, 12].

The United States of America (USA) Food and Drug Administration (FDA) has approved four HIV-1 INSTIs, including raltegravir (RAL), elvitegravir (EVG), dolutegravir (DTG) and bictegravir (BIC) [13]. However, the high cost of INSTIs, has resulted in restricted access to this class of drugs in resource-limited countries [14]. Despite the cost, the World Health Organization (WHO) has given the green light to include DTG to an alternative 1st-line regimen [14, 15]. The strand transfer reaction catalyzed by HIV-1-expressed IN enzyme is blocked by the activity of INSTIs which bind to the catalytic site in the catalytic core domain (CCD) of the IN protein [16, 17]. Mutations that confer resistance to INSTIs (for example G140S, Q148H and N155H) have been structurally mapped in close proximity to the IN catalytic active site [18, 19]. Primary resistance to INSTIs, along with residues associated with catalytic activity among different subtypes are highly conserved. HIV-1 sequence and structure-based analyses have shown that polymorphic residues can cause subtype-specific effects, which significantly affect the native protein structure, function and activity important for drug-mediated inhibition of enzyme activity [9]. There is limited information available on the IN structure of CRF02_AG [9, 20] and even less on the effect of mutations on the protein structure. There is therefore a need to continue monitoring patients to identify additional RAMs and polymorphic mutations that might affect the genetic barrier to the development of RAMs against INSTI [9]. The goals of this study was to analyse the Cameroonian CRF02_AG IN gene sequences obtained from the Los Alamos National Laboratory ( HIV-1 database to assess the occurrence of mutations and natural occurring polymorphisms (NOPs). NOPs are categorized under secondary mutations which on their own play a limited role in resistance [14]. However, their pre-existence might favour a more rapid evolution towards resistance when additional mutations are selected under therapy [21]. In this study, the possible impact caused by statistically enriched NOPs found in CRF02_AG subtype was modelled within the context of a three-dimensional (3D) protein structure of the HIV-1CRF02-DNA-MG-DTG IN complex. Subsequently, stability predictions was performed using the mutation cut-off scanning matrix server (mCSM) to assess change in Gibbs free energy of mutations on the protein structure followed by interaction analysis to assesses the loss or gain of DTG interactions to Wild type and six mutated HIV-1CRF02-DNA_MG IN structures. Molecular modelling of HIV-1 CRF02 integrase sequences and DTG interaction analysis will help determine which mutations could affect the genetic barrier to the emergence of DTG drug resistance.


Generation of consensus HIV-1CRF02_AG Integrase sequence

To compare our study CRF02_AG treatment naïve sequences (n = 20) available in GenBank under the following accession numbers: MN816445- MN816488 [2], with INSTIs treatment naïve CRF02_AG IN sequences from Cameroon between 1994 and 2010. We performed a search on the HIV Los Alamos National Laboratory database (LANL) database for additional (n = 287) INSTI treatment Naïve patients’ sequences (, completed on 01 February 2020. All Cameroonian HIV-1 subtype CRF02_AG IN sequences for treatment naïve patients, were included in our search criteria [14]. We selected one sequence per patient and every problematic sequence were excluded from further analyses. The consensus sequence was generated using the database-derived HIV-1 CRF02_AG sequences (n = 287) and CRF02_AG cohort sequences from our previous study (n = 20) [20]. An online quality control program that is available on the HIVLANL database ( was employed to screen nucleotide sequences for quality and to verify for stop codons, insertions and/or deletions. The MAFFT version 7 sequence alignment tool, was used to perform a multiple sequence alignment from which the consensus sequence was derived [22]. In an effort to eliminate possible contamination, which was part of our quality control measure, each of the viral sequences were inferred on a phylogenetic tree.

HIV-1 subtyping using online programs

HIV-1 subtyping based on IN sequences was performed using two online available programs; REGA version 3 ( and COMET-HIV (

Drug resistance analysis

Mutations associated with resistance to INSTIs were identified using the Stanford University genotypic resistance interpretation algorithm, HIVdb version 8.3 ( All drug resistance mutations results were classified as either major or minor mutations, last accessed 01 April 2020.

Homology modelling and quality assessment

The crystal structure of the HIV-1B intasome (nucleoprotein complex: containing viral DNA ends and the viral integrase protein) (PDBID: 5U1C) was used as a homologous template to generate a three-dimensional tetrameric structure of HIV-1CRF02_AG IN using the consensus sequence of recombinant subtype CRF02_AG sequence that we generated. The SWISSMODEL webserver was used for model construction by first constructing a pairwise sequence-structure alignment between HIV-1C wild-type (WT) amino acid sequence and template 5U1C [22]. The quality of the resulting model was assessed using SWISSMODEL quality assessment scores; the Qualitative model energy analysis (QMEAN) and Global model quality estimate (GMQE) scores. The QMEAN score is a composite scoring function assessing the major geometrical aspects of protein structures by comparing the predicted protein model to experimental structures of similar size, scores close to zero suggest high compatibility to experimental structures. The GMQE score estimates the quality of the predicted model using the properties from the target-template alignment and the template structure. Values are between 1 and 0, with higher values indicating more reliable models [23]. The Root mean square deviation (RMSD) analysis was done to compare the predicted model to the homologous template (PDBID: 5U1C) and to determine if there were any structural similarity between the two structures, lower RMSD values suggest very little deviation in terms of the backbone between the two structures. We also used publicly available algorithms located at the SAVES webserver ( namely; ERRAT [24] VERIFY3D [25] and PROCHECK [26] to assess the quality of the predicted protein model. ERRAT statistically interrogated the nonbonded atomic interactions of the given target against the interactions of refined structures, with higher scores above 50 indicating higher quality in the protein model. VERIFY3D determines the compatibility of a structure (3D) to its own amino acid sequence (1D), higher values for VERIFY3D indicates high compatibility. PROCHECK on the other hand generates an Ramachandran plot that assesses the stereochemical parameters in a protein structure and if the percentage of phi and psi dihedral angles within the protein structure is more than 90% then the protein model has favourable residue conformations.

Change in free energy predictions and interaction analysis

The predicted 3D structure of HIV-1 CRF02_AG IN was energy minimized using GROMACS software [27] and the resulting structure was used to introduce resistance associated and accessory mutations identified from the Stanford University HIVdb using PyMol mutagenesis wizard. The WT structure and a text file specifying single mutations was used as input to the program, mutation cut-off scanning matrix (mCSM). The change in energy after introduction of a mutation was calculated using mCSM. mCSM uses graph-based distance patterns of neighbouring residues and calculates a Delta-delta G-score for the impact of the mutation on the protein network and provides a phenotypic assessment by annotating a mutation as either being destabilizing (negative value) or stabilizing (positive value). The loss or gain of interactions between the WT and mutant neighbouring residues was calculated using Pymol find polar contacts option. To determine the effect of the mutation on DTG drug binding to the IN structure we energy minimized the WT CRF02-AG-DNA_MG_DTG and mutant complexes using GROMACS. Afterwards, we only selected the mutant structures that showed changes in the number of polar contacts with neighbouring residues and calculated DTG interactions with CRF02-DNA-MG complex IN structures using the PyMol find polar contacts option.


HIV-1 subtyping

HIV-1 subtyping was done using HIV-1 subtyping online-automated tools and all sequences were verified using phylogenetic tree as HIV-1 subtype CRF02_AG (Fig. 1).

Fig. 1

HIV-1 Integrase phylogenetic analysis inferred by ML. The Phylogenetic tree inferred in MEGA contains (n = 20) patient derived sequences in red box and (n = 287) online database sequence indicated in black box. HIV-1 reference sequences were acquired from the HIV-1 LANL database, using the 2010 data set. An ML tree was constructed using Mega version 7.0, with the Kimura 2 parameter. The alignment was based on HXB2 position 4351–5069, of approximately 700 bp in length. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (1000 replicates) are shown next to the branches. The bootstrap values are above 70 indicating significant support for the branches simulated. The reference sequences are the unboxed sequences. All of the sequence’s clusters with HIV-1 subtype CRF02_AG. LANL, Los Alamos National Laboratory; MEGA, Molecular Evolutionary Genetics Analysis; ML, maximum likelihood

Database derived IN sequence resistance analyses

After excluding multiple sequences from a patient to avoid overestimation of the variant calling and problematic sequences, we used 287 sequences collected between 1994 and 2010. These sequences were subsequently screened for the presence of RAMs. We identified 12.8% (37/287) sequences to contain RAMs, with only 1.0% (3/287) having major INSTI RAMs: T66A, Q148H, R263K and N155H. Two mutations, Q148H and R263K, occurred together in one sequence (0.3%), whereas T66A and N155H were present individually in one sequence each. A total of 11.8% (34/287) of the sequences contained five different IN accessory mutations, namely Q95K, T97A, G149A, E157Q and D232N. Mutations G149A and D232N occurred together in one sequence (0.3%). Notably, one sequence dating back from 2010 had two major mutations; Q148H and R263K in combination with two other minor mutations G149A and D232N.

Generation of the consensus sequence for Cameroonian’s HIV-1 CRF02_AG subtype

The consensus sequences generated using the database-derived HIV-1 CRF02_AG sequences (n = 287) and cohort sequences (n = 20), identified 20 naturally occurring polymorphisms (NOPS): E11D,K14R, V31I, M50I, I72V, L74MVI, L101I, T112V, T124A, G134N, I135V, K136K/Q, V201I, T206S, T218I, L234I, A265V, R269K, S283G (Fig. 2). Three of these (E11D, K14R and V31I) belong to the NTD, whereas M50I belongs to the loop region connecting the NTD and CTD. Eleven NOPs (I72V, L74MVI, L101I, T112V, T124A, T124A, G134N, I135V, K136K/Q, V201I and T206S) are part of the CCD, and the remaining five (T218I, L234I, A265V, R269K and S283G) belong to the CTD.

Fig. 2

Prevalence of NOPs in IN genes from CRF02_AG subtypes. The figure shows the distribution of variants among the 287 and 20 CRF02_AG full length integrase sequences. Divided into: N-terminal domain (NTD) (residues 1–50), catalytic core domain (CCD) (residues 50–212) and C-terminal domain (CTD) composes of amino acids 213–288

Molecular modelling and structural quality assessment

The sequence identity between the amino acid sequences of HIV-1 CRF02_AG sequence and the homologous template 5U1C was very high, approximately 93% and the sequence similarity was found to be 60% between the two sequences. The high sequence identity and coverage provides confidence in modelled regions of the protein structure and reduces the occurrence of any problematic or unresolved regions within the final protein model. Figure 3a, shows the 3D tetrameric structure for HIV-1 CRF02_AG IN that consist of 288 amino acids, 10 alpha helices, 9 beta sheets and 19 coil regions. The internal assessment scores calculated for the predicted model of HIV-1 CRF02_AG had an GMQE score of 0.10 and a QMEAN4 score of − 2.23, both scores confirming reliability of the modelled regions within the protein structure. Furthermore, the homology model passed most of the external 3D quality validation checks. The Verify3D score for the model was predicted to be 71.1% (acceptable for crude structures before energy minimization), while ERRAT score for all the chains was 86.0% and higher, the PROCHECK analysis indicated that 98.0% of residues occurred in most favoured and allowed regions of the Ramachandran plot, and the Prosa Z-score was − 6.18 which is in range with proteins of a similar size. Superimposing the template 5U1C onto the energy minimized structure of HIV-1 CRF02_AG indicated an RMSD value of 0.212 Å, suggesting very little backbone deviation in main chain atoms (Fig. 3b). Figure 3c, shows the locations of the 15 mutations relative to the active site.

Fig. 3

3D tetrameric structure for HIV-1 CRF02_AG IN. a Three-dimensional tetrameric structure predicted for HIV-1 CRF_02 AG IN in complex with magnesium ions. Chain A: green, Chain B: cyan, Chain C: Magenta, Chain D: yellow, Magnesium ions shown as spheres coloured in green and Dolutegravir shown as sticks coloured in red. No DNA shown. b Structural superimposition of HIV-1 AG IN onto HIV-1 B (5U1C) in complex with MG ions. CRF02_AG IN: green, 5U1C: blue, MG ions shown as magenta spheres. c Locations of stabilizing and destabilizing mutations on HIV-1 integrase CRF02_AG structure. Mutations that affect the protein structure are labelled and shown as red sticks and mutations with no effects are also labelled and shown as blue sticks. Magnesium ions shown as spheres coloured in green and Dolutegravir shown as sticks coloured by atom type

Gibbs free energy change and interaction analysis

The mCSM predictions indicated that 14 of the 15 mutations (10 RAM’s and five accessory mutations), i.e. the M50I, T66A, L74I, L74M, T97A, G118S, S119R, P145S, Q148H, G149A, N155H, E157Q, D232N and R263K substitutions resulted in destabilizing effects of − 0.582, − 0.703 -1.069, − 0.93, − 1.051, − 0.492, − 0.091, − 0.485, − 0.133, − 0.421, − 0.975, − 1.111, − 0.512 and − 0.455 Kcal/Mol each, respectively. Only substitution Q95K resulted in a slightly stabilizing effect of 0.146Kcal/Mol. Interaction analysis of the single amino acid changes indicated differences in the number and type of interaction between neighbouring residues and the DNA. The T97A showed four polar contacts for T97 compared to the three of A97 (Table 1). This suggests a loss of stable contacts in this region that could destabilize the protein structure. Moreover, T66A, Q148H, N155H, D232N and R263K mutations all indicated a loss of interactions with neighbouring residues after the introduction of substitutions, while only the N155H mutation gained an additional interaction with the 3′ terminal viral DNA Adenine21 (Table 1). Inspection of the E157 residue showed four contacts with neighbouring residues while Q157 revealed five polar contacts of which two were with the 3′ terminal viral DNA (Table 1). In addition, the remaining other six substitutions; M50I, L74I, L74M, Q95K, G118S and P145S showed no changes in the number or type of interactions, implying no strong effect on the protein structure and function (Table 1). Protein drug interaction analysis of energy minimized complexes revealed interesting findings as accessory mutation E157Q made only one MG ion interaction and D232N none, while substitutions T66A, T97A, Q148H, R263K and N155H all had ionic interactions with two MG ions as well as with DDE motif active site residues and the 3′ terminal viral DNA nucleotides (Table 2) and (Figs. 4a-h). Most importantly is to note that MG ions are crucial for DTG coordination to displace viral DNA and thereby preventing HIV viral integration into host DNA.

Table 1 The number of polar contacts observed between WT residue and neighbouring residues before and after the introduction of the RAM’s and Accessory Mutations
Table 2 Summary of all interactions observed between DTG and CRF_02AG IN for the WT and seven mutant structures
Fig. 4

Interactions formed between DTG and the energy minimized WT and seven mutant structures for HIV-1 CRF02_AG Integrase. a WT HIV-1 CRF02_AG IN showing in total six contacts formed between DTG and two IN residues, two DNA nucleotides and two MG ions. b T66A HIV-1 CRF02_AG IN showing in total four contacts formed between DTG and one IN residue, one DNA nucleotide and two MG ions. c T97A HIV-1 CRF02_AG IN showing in total five contacts formed between DTG and two IN residues, one DNA nucleotide and two MG ions. d E157Q HIV-1 CRF02_AG IN showing in total four contacts formed between DTG and one IN residue, two DNA nucleotides and one MG ion. e Q148H HIV-1 CRF02_AG IN showing in total five contacts formed between DTG and one IN residue, two DNA nucleotides and two MG ions. f R263K HIV-1 CRF02_AG IN showing in total four contacts formed between DTG and one IN residue, one DNA nucleotide and two MG ions. g N155H HIV-1 CRF02_AG IN showing in total seven contacts formed between DTG and three IN residues, two DNA nucleotides and two MG ions. h D232N HIV-1 CRF02_AG IN showing in total two contacts formed between DTG and two DNA nucleotides. The drug DTG is shown as sticks and coloured in magenta, the MG ions are shown as spheres coloured in green while the IN protein residues and terminal 3’end viral DNA nucleotides are labelled and shown as sticks. Three letter codes for IN protein residues and numbers are given as well as the DNA nucleotide three letter codes and numbers. The nucleotides represent terminal end viral DNA nucleotides after 3′ processing


Despite INSTIs having an increased genetic barrier against resistance, studies performed from high-income countries shows that the occurrence of RAMS against INSTIs can happen, via acquired drug resistance mutations (DRM) and/or transmitted DRM, leading to reduced susceptibility to INSTIS and possible treatment failure [22, 28]. The IN mutations usually associated with reduced INSTIs susceptibility include both polymorphic mutations and non-polymorphic mutations [29, 30]. Other studies have reported that several NOPs can affect structural stability and flexibility of the IN protein structure [31, 32]. Previous researchers have reported low rates of IN mutations against INSTIs in Cameroon [15, 20]. The WHO has recommended the utilization of DTG as part of first-line regimens [16]. With the approval of INSTIs usage worldwide, it is predicted that approximately 57% of people living with HIV will be receiving DTG based regimens, including new-borns and children [33]. It is therefore imperative to screen for the presence of mutations against INSTIs which can affect treatment outcomes. Currently, there is limited data available for INSTI RAMs from studies that focuses on the SSA region, where over two-thirds of the presently infected individuals reside [34, 35]. In our previous studies we found low level of RAMSs against INSTIs [3].

In a recent study on CRF02_AG IN, we reported that accessory mutations can impact the binding of DTG with or without combination of primary resistance mutations [32, 36]. In this study we analysed the CRF02_AG IN gene sequences for the presence of polymorphic and non-polymorphic mutations. Four major INSTIs mutations were found within the database sequences: T66A, Q148H, N155H and R263K. R263K displayed moderate level resistance against EVG (12-fold) [37] and seems to confer low-level resistance against DTG. Structural analyses have suggested that DTG shares a similar interfacial mechanism of inhibition with EVG and RAL, but is able to make more intimate contacts with the viral DNA [38]. In addition, DTG can adapt its position and conformation in response to structural changes within the active site of EVG- or RAL resistant IN enzymes and in doing so avoid cross-resistance as a result of slower dissociation rates [39, 40]. Two principal mutation pathways identified from our study that reduces susceptibility to RAL are Q148H/K/R and N155H. These mutations are located in close proximity to the Integrase’s active site and each mutation significantly reduces viral fitness by 92-fold for Q148R, 30-fold for N155H [41]. Q148H and N155H mutations are thought to trigger conformational changes within the catalytic pocket that result in lower binding affinity of INSTIs to IN [42]. The variant T66A which is normally selected by EVG treatment, was detected in 0.3% of our sequence cohort. This variant is associated with 5-fold reduced susceptibility to EVG, however, T66A also bears cross-resistance to DTG and is selected by RAL [41]. Abraham et al., 2013, showed that the T66A mutant occurs within the two distal sheet from the DDE triad motif. The close proximity of the T66A/I/K variants to the viral DNA 3′ end and mutation N155H, could sterically hamper viral DNA binding and/or metal ion coordination with DTG [41]. The fact that only 1.0% of sequences analysed contained INSTI primary RAMs suggest that mutations against INSTIs will need to be monitored carefully against Cameroonians living with HIV. This result is in agreement with other studies done in Africa [20, 43,44,45] Asia [46, 47], Europe [48, 49] and South America [50] where studies showed a low frequency of INSTI primary RAMs.

In our study, we observed five IN accessory RAMs; namely Q95K, T97A, G149A, E157Q and D232N. T97A mutation can reduce EVG susceptibility by 3-fold [41] and combination of T97A mutation with other INSTI RAMs lead to reduced susceptibility to RAL [51, 52] and DTG [53, 54]. E157Q acts as a compensatory mutation and individually has a negligible effect on the susceptibility to INSTIs; however, a combination of E157Q with other INSTI RAMs may lead to reduced susceptibility to INSTIs [55, 56]. Individuals containing E157Q mutation in combination with other IN RAMs showed reduced susceptibility to DTG. Moreover, another rare nonpolymorphic accessory resistance mutation Q95K confers little if any effect on drug susceptibility to INSTIs [57]. A study by Axel Fun et al., 2010, showed that this secondary mutation enhances N155H-mediated resistance and partially restores the reduced replication caused by N155H [58]. In our study, we detected L74M mutations at a frequency > 20%, which is not surprising since, nearly 10% of ARV-Naïve patients infected with CRF02_AG viruses harbours L74M mutations [59]. This L74M mutation has minimal if any effect against susceptibly of INSTIs, but in combination with mutations at positions 140 and 148, it reduces susceptibility of DTG [38, 60, 61]. Within the IN CCD, we observed 11 of the reported INSTI NOPs. This IN region is important for recognition of DNA, binding and chromosomal integration of the newly synthesized double-stranded viral DNA into the host genomic DNA [62,63,64]. It contains the endonuclease and polynucleotide transferase site [62,63,64]. While in the CTD, a region that helps stabilize the integrase–viral DNA complex, five other NOP mutations were observed [65]. All of the afore mentioned mutations in either the CCD and/or CTD regions have the potential to affect the IN protein function and interfere with INSTIs binding [65].

We further analysed the effect of NOPs on the stability of the structures and neighbouring residues. Most of the variants noted in our study were shown to destabilise the protein structure, except for one mutation Q95K, that showed to exert a slightly stabilising effect on the protein structure and no changes in the number of polar contacts with neighbouring residues making it unlikely to affect the IN protein structure. It is known that destabilising effects of mutations on the protein structure might reduce drug binding. This was further explored by performing interaction analysis between the drug DTG and energy minimized structures of the WT and mutants T66A, T97A, Q148H, N155H, E157Q, R263K and D232N. The findings revealed accessory mutations E157Q and D232N had the potential to reduce and or prevent DTG binding to HIV-1 CRF02_AG IN structure based on the loss of MG ion interactions, while known RAM’s does not seem to influence DTG drug binding. However, the effect of RAM’s on DTG drug binding needs to be validated using molecular dynamic simulations to calculate the change in free energy of binding between DTG and HIV-1 CRF02_AG IN. Interestingly, the mutation E157Q occurred within the stable alpha-helix secondary structure element and made more contacts with DNA (stabilizing viral DNA complex), while the D232N mutation occurred within the stable Beta-sheet secondary structure element and in close proximity to the flexible G140’s loop region suggesting that these changes can affect the protein conformation and thereby interfere with drug binding leading to resistance.

A limitation of the study is the use of online database sequences that may contain contaminated and truncated sequences leading to spurious phylogenetic tree results and also these databases are not regularly updated. Furthermore, gaps in the aligned regions between the homologous template and target sequence may result in unresolved loop regions within the protein model which is one of the limitations of 3D protein modelling that can result in inaccurate interaction prediction. Furthermore, the interaction analysis was done for only a single static structure of the protein structure and does not consider the dynamic behaviour of the protein structure that might result in the loss and under estimation of crucial interaction partners. An important finding in this study is the fact that sequence diversity amongst different subtypes may affect different folding conformations of the HIV-1 IN subtypes thereby allowing not only RAM’s but accessory mutations to result in less efficacious INSTI binding to HIV-1 IN structures.


Molecular modelling and interaction analysis provided novel insights into the effect of accessory mutations (E157Q and D232N) on HIV-1 CRF02_AG IN drug resistance. This emphasise the need to screen for the presence of INSTIs major RAM’s and accessory mutations in patients on INSTI treatment. This would help identify pathways that contribute to drug resistance and help tailor more effective treatment regimens in INSTI naïve patients.

Availability of data and materials

The datasets generated and/or analysed during the current study are available in the HIV Los Alamos National Laboratory database (LANL) ( repository, and under NCBI database with GenBank accession numbers: MN816445- MN816488. The predicted 3D structural model of HIV-1 CRF02_AG IN is available from the corresponding author on reasonable request because we will use the predicted model in additional simulation studies.





Integrase strand transfer inhibitors


Natural occurring polymorphism


Sub-Saharan Africa


Resistance associated mutations


World Health Organisation


Los Alamos National Laboratory


Central core domain (CCD)


C-terminal domain


N-terminal domain


combination antiretroviral therapy


Prevention of mother-to-child transmission


Food and Drug Administration










Three dimension


Medical Research Council


Mutation cut-off scanning matrix


  1. 1.

    Han Y-S, Mesplède T, Wainberg MA. Differences among HIV-1 subtypes in drug resistance against integrase inhibitors. Infect Genet Evol. 2016;46:286–91.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  2. 2.

    Mikasi SG, Gichana JO, Van der Walt C, Brado D, Obasa AE, Njenda D, et al. HIV-1 Integrase diversity and resistance-associated mutations and polymorphisms among Integrase Strand transfer inhibitor-naive HIV-1 patients from Cameroon. AIDS Res Hum Retrovir. 2020;36(5):450–5.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  3. 3.

    Mikasi SG, Isaacs D, Ikomey GM, Shimba H, Cloete R, Jacobs GB. HIV-1 drug resistance Mutation analyses of Cameroon-derived Integrase sequences. AIDS Res Hum Retrovir. 2020;00(00):1–3.

    Google Scholar 

  4. 4.

    Plantier J, Djemai M, Leme V, Reggiani A, Leoz M, Burc L, et al. Census and analysis of persistent false-negative results in Serological diagnosis of human immunodeficiency virus type 1 group O infections. J Clin Microbiol. 2009;47(9):2906–11.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  5. 5.

    Oliveira MF, Ramalho DB, Abreu CM. Genetic diversity and naturally polymorphisms in HIV type 1 Integrase isolates from Maputo, Mozambique. AIDS Res Hum Retroviruses. 2012;28(12):1788–92.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  6. 6.

    Bar-magen T, Sloan RD, Faltenbacher VH, Donahue DA, Kuhl BD, Oliveira M, et al. Comparative biochemical analysis of HIV-1 subtype B and C integrase enzymes. Retrovirology. 2009;6:2048.

    Article  CAS  Google Scholar 

  7. 7.

    Bar-magen T, Donahue DA, Mcdonough EI, Kuhl D, Faltenbacher VH, Xu H, et al. HIV-1 subtype B and C integrase enzymes exhibit differential patterns of resistance to integrase inhibitors in biochemical assays. Aids. 2010;24:2171–9.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  8. 8.

    Quashie PK, Han Y, Hassounah S, De TM. Structural Studies of the HIV-1 Integrase Protein: Compound Screening and Characterization of a DNA-Binding Inhibitor. Plos one. 2015;10:1–18.

    Article  CAS  Google Scholar 

  9. 9.

    Theys K, Libin PJ, Van Laethem K, Abecasis AB. An evolutionary-based approach to quantify the genetic barrier to drug resistance in fast- evolving viruses: an application to HIV-1 subtypes and integrase inhibitors. Antimicrob Agents Chemother. 2019;63:8.

    Article  Google Scholar 

  10. 10.

    Bain LE, Nkoke C, Noubiap JJN. UNAIDS 90–90 – 90 targets to end the AIDS epidemic by 2020 are not realistic: comment on “ Can the UNAIDS 90–90 – 90 target be achieved ? A systematic analysis of national HIV treatment cascades”. BMJ Global Health. 2020;2:2016–8.

    Google Scholar 

  11. 11.

    Ikomey GM, Claire M, Assoumou O, Gichana JO, Njenda D. Original article observed HIV drug resistance associated mutations amongst naïve immunocompetent children in Yaoundé, Cameroon. Germs. 2017;7(December):178–85.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  12. 12.

    Tchouwa GF, Eymard-Duvernay S, Cournil A, Lamare N, Serrano L, Butel C, et al. Prevalence of pretreatment HIV drug resistance in Cameroon following a nationally representative WHO survey. J Antimicrob Chemother 2018;73(9):2468–2474. Available from:

  13. 13.

    Jacobson JM, Flexner CW. Universal antiretroviral regimens: thinking beyond one-pill-once-a-day. Curr Opin HIV AIDS. 2017;12(4):343–50.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  14. 14.

    Acharya A, Tagny CT, Mbanya D, Fonsah JY, Nchindap E, Kenmogne L, et al. Variability in HIV-1 Integrase gene and 3′-Polypurine tract sequences in Cameroon clinical isolates, and implications for Integrase inhibitors efficacy. Int J Mol Sci. 2020;21(5):1553.

    CAS  Article  PubMed Central  Google Scholar 

  15. 15.

    WHO. Consolidated guidelines on the use of antiretroviral drugs for treating and preventing hiv infection. Recommendations for a public health approach. 2nd ed: World Health Organization; 2019. Available from:

    Google Scholar 

  16. 16.

    Yoshinaga T, Kobayashi M, Seki T, Miki S, Wakasa-morimoto C, Suyama-kagitani A. Antiviral characteristics of GSK1265744, an HIV integrase inhibitor dosed orally or by long-acting injection. Antimicrob Agents Chemother. 2015;59(1):397–406.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  17. 17.

    Hazuda DJ, Felock P, Witmer M, Wolfe A, Stillmock K, Grobler JA, et al. Inhibitors of strand transfer that prevent integration and inhibit HIV-1 replication in cells. Science. 2000;287:646–51.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  18. 18.

    Menéndez-arias L. Molecular basis of human immunodeficiency virus type 1 drug resistance: Overview and recent developments. Antiviral Res. 2013;98:93–120.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  19. 19.

    Hare S, Smith SJ, Me M, Jaxa-chamiec A, Pommier Y, Hughes SH, et al. Structural and functional analyses of the second-generation integrase strand transfer inhibitor dolutegravir (S / GSK1349572). Mol Pharmacol. 2011;80(4):565–72.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  20. 20.

    Mikasi SG, Gichana JO, Walt C, Van Der Brado D, Obasa AE, Njenda D, et al. Transfer Inhibitor-Naive HIV-1 Patients from Cameroon. AIDS Res Hum Retroviruses. 2020;36(5):450–5.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  21. 21.

    Paraschiv S, Otelea D, Dinu M, Maxim D, Tinischi M. Polymorphisms and resistance mutations in the protease and reverse transcriptase genes of HIV-1 F subtype Romanian strains. Int J Infect Dis. 2007;11:123–8.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  22. 22.

    Kantor RS, Wrighton KC, Handley KM, Sharon I, Hug LA, Castelle CJ, et al. Small genomes and sparse metabolisms of sediment-associated bacteria from four candidate phyla. MBio. 2013;4(5):1–11.

    Article  CAS  Google Scholar 

  23. 23.

    Waterhouse A, Bertoni M, Bienert S, Studer G, Tauriello G, Gumienny R, et al. SWISS-MODEL : homology modelling of protein structures and complexes. Nucleic Acids Res. 2018;46(May):296–303.

    Article  CAS  Google Scholar 

  24. 24.

    Colovos C, Yeates TO. Verification of protein structures: patterns of nonbonded atomic interactions. Protein Sci. 1993;2(9):1511–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  25. 25.

    Bowie JU, Lüthy R, Eisenberg D. A method to identify protein sequences that fold into a known three-dimensional structure. Science. 1991;253(5016):164–70.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  26. 26.

    Laskowski RA, Rullmannn JA, MacArthur MW, Kaptein R, Thornton JM. AQUA and PROCHECK-NMR: programs for checking the quality of protein structures solved by NMR. J Biomol NMR. 1996;8(4):477–86.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Oostenbrink C, Villa A, Mark AE, Gunsteren WF. A biomolecular force field based on the free enthalpy of hydration and solvation: The GROMOS Force-Field Parameter sets 53A5 and 53A6. J Comput Chem. 2004;25:1656–76.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  28. 28.

    Lepik KJ, Harrigan PR, Yip B, Wang L, Robbins MA, Zhang WW, et al. Emergent drug resistance with integrase strand transfer inhibitor-based regimens. AIDS. 2017;31(10):1425–34.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  29. 29.

    Rhee S, Liu TF, Kiuchi M, Zioni R, Gifford RJ, Holmes SP, et al. Natural variation of HIV-1 group M integrase : Implications for a new class of antiretroviral inhibitors. Retrovirology. 2008;11:1–11.

    Google Scholar 

  30. 30.

    Anstett K, Brenner B, Mesplede T, Wainberg MA. HIV drug resistance against strand transfer integrase inhibitors. Retrovirology. 2017;14:1–16.

    Article  CAS  Google Scholar 

  31. 31.

    Rogers L, Obasa AE, Jacobs GB, Sarafianos SG, Sönnerborg A, Neogi U, et al. Structural implications of genotypic variations in HIV-1 integrase from diverse subtypes. Front Microbiol. 2018;9:1–9.

    Article  Google Scholar 

  32. 32.

    Chitongo R, Obasa AE, Mikasi SG, Jacobs GB, Cloete R. Molecular dynamic simulations to investigate the structural impact of known drug resistance mutations on HIV-1C Integrase-Dolutegravir binding. PLoS One. 2020;15(5):e0223464.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  33. 33.

    Gupta A, Juneja S, Vitoria M, Habiyambere V, Nguimfack D, Doherty M, et al. Projected uptake of new antiretroviral ( ARV ) medicines in adults in low- and middle- income countries : a forecast analysis 2015–2025. PLoS One. 2016;11:1–18.

    Google Scholar 

  34. 34.

    Mabeya S, Nyamache A, Ngugi C, Nyerere A, Lihana R. Characterization of HIV-1 Integrase gene and Resistance associated mutations prior to roll out of Integrase inhibitors by Kenyan national HIV-treatment program in Kenya. Ethiop J Health Sci. 2020;30(1):37–44.

    Article  PubMed  PubMed Central  Google Scholar 

  35. 35.

    Masoud S, Kamori D, Barabona G, Mahiti M, Sunguya B, Lyamuya E, et al. Circulating HIV-1 integrase genotypes in Tanzania: implication on the introduction of integrase inhibitors-based antiretroviral therapy regimen. AIDS Res Hum Retroviruses. 2020;36(6):539–43.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  36. 36.

    Isaacs D, Mikasi SG, Obasa AE, Ikomey GM. Structural comparison of diverse Hiv-1 subtypes using molecular modelling and docking analyses of integrase inhibitors. Viruses. 2020;12:936.

    CAS  PubMed Central  Article  Google Scholar 

  37. 37.

    Hassounah SA, Mesplède T, Quashie PK, Oliveira M, Sandstrom PA, Wainberg A. Effect of HIV-1 Integrase resistance mutations when introduced into SIVmac239 on susceptibility to Integrase Strand transfer inhibitors. J Virol. 2014;88(17):9683–92.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  38. 38.

    Kobayashi M, Yoshinaga T, Seki T, Wakasa-Morimoto C, Brown KW, Ferris R, et al. Kawauchi- Miki S, Taishi T, Kawasuji T, Johns BA, underwood MR, Garvey EP, Sato a FT. in vitro antiretroviral properties of S/GSK1349572, a next-generation HIV integrase inhibitor. Antimicrob Agents Chemother. 2011;55(2):813–21.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  39. 39.

    Hurt CB, Sebastian J, Hicks CB, Eron JJ. Resistance to HIV integrase strand transfer inhibitors among clinical specimens in the United States, 2009-2012. Clin Infect Dis. 2014;58(3):423–31.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  40. 40.

    Hightower KE, Wang R, Deanda F, Johns BA, Weaver K, Shen Y, et al. Dolutegravir (S/GSK1349572) exhibits significantly slower dissociation than raltegravir and elvitegravir from wild-type and integrase inhibitor-resistant HIV-1 integrase-DNA complexes. Antimicrob Agents Chemother. 2011;55(10):4552–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  41. 41.

    Abram ME, Hluhanich RM, Goodman DD, Andreatta KN, Margot NA, Ye L, et al. Impact of primary elvitegravir resistance-associated mutations in HIV-1 integrase on drug susceptibility and viral replication fitness. Antimicrob Agents Chemother. 2013;57(6):2654–63.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  42. 42.

    Hare S, Gupta SS, Valkov E, Engelman A, Place N. HHS Public Access. 2010;464(7286):232–236.

    Google Scholar 

  43. 43.

    Obasa AE, Mikasi SG, Brado D, Cloete R, Singh K, Neogi U, et al. Drug resistance mutations against protease, reverse transcriptase and integrase inhibitors in people living with HIV-1 receiving boosted protease inhibitors in South Africa. Front Microbiol. 2020;11(March):1–9.

    Google Scholar 

  44. 44.

    Inzaule SC, Hamers RL, Noguera-Julian M, Casadellà M, Parera M, Rinke de Wit TF, et al. Primary resistance to integrase strand transfer inhibitors in patients infected with diverse HIV-1 subtypes in sub-Saharan Africa. J Antimicrob Chemother. 2018;92:1167–72.

    Article  CAS  Google Scholar 

  45. 45.

    Villabona-arenas CJ, Domyeum J, Mouacha F, Delaporte E, Peeters M, Mpoudi-Ngole E, et al. HHS Public Access; 2015. p. 210–6.

    Google Scholar 

  46. 46.

    Kim J-Y, Kim E-J, Choi J-Y, Kwon O-K, Kim GJ, Choi SY, et al. Genetic variation of the HIV-1 integrase region in newly diagnosed anti-retroviral drug-naïve patients with HIV/AIDS in Korea. Clin Microbiol Infect. 2011;17(8):1155–9.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  47. 47.

    Kotaki T, Khairunisa SQ, Sukartiningrum SD, Witaningrum AM, Rusli M, Diansyah MN, et al. Detection of drug resistance-associated mutations in human immunodeficiency virus type 1 integrase derived from drug-naive individuals in Surabaya, Indonesia. AIDS Res Hum Retroviruses. 2014;30(5):489–92.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  48. 48.

    Meixenberger K, Yousef KP, Smith MR, Somogyi S, Fiedler S, Bartmeyer B, et al. Molecular evolution of HIV-1 integrase during the 20 years prior to the first approval of integrase inhibitors. Virol J. 2017;14:1–13.

    Article  CAS  Google Scholar 

  49. 49.

    Tostevin A, White E, Dunn D, Croxford S, Delpech V, Williams I, et al. Recent trends and patterns in HIV-1 transmitted drug resistance in the United Kingdom. HIV Med. 2017; Available from:

  50. 50.

    Passaes CB, Guimaraes ML, Fernandez SL, Lorete Rdos S, Teixeira SL, Fernandez JC, et al. Lack of primary mutations associated with Integrase inhibitors among HIV-1 subtypes B, C, and F Circulating in Brazil. JAIDS. 2009;51(413):7–12.

    CAS  PubMed  PubMed Central  Google Scholar 

  51. 51.

    Eron JJ, Clotet B, Durant J, Katlama C, Kumar P, Lazzarin A, et al. Safety and efficacy of dolutegravir in treatment-experienced subjects with raltegravir-resistant HIV type 1 infection: 24-week results of the VIKING Study. J Infect Dis. 2013;207:740–8.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  52. 52.

    Fransen S, Gupta S, Danovich R, Hazuda D, Miller M, Witmer M, et al. Loss of raltegravir susceptibility by human immunodeficiency virus type 1 is conferred via multiple nonoverlapping genetic pathways. J Virol. 2009;83(22):11440–6.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  53. 53.

    Melikian GL, Rhee S, Varghese V, Porter D, White K, Taylor J, et al. Non-nucleoside reverse transcriptase inhibitor (NNRTI) cross-resistance: implications for preclinical evaluation of novel NNRTIs and clinical genotypic resistance testing. J Antimicrob Chemother. 2014;69:12–20.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  54. 54.

    Naeger LK, Harrington P, Komatsu T, Deming D. Effect of dolutegravir functional monotherapy on HIV-1 virological response in integrase strand transfer inhibitor resistant patients. Antivir Ther. 2016;21(6):481–8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  55. 55.

    Charpentier C, Malet I, Andre-garnier E, Storto A, Bocket L, Amiel C, et al. Phenotypic analysis of HIV-1 E157Q integrase polymorphism and impact on virological outcome in patients initiating an integrase inhibitor-based regimen. J Antimicrob Chemother. 2018;73:1039–44.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  56. 56.

    Anstett K, Cutillas V, Fusco R, Mesplède T, Wainberg MA. Polymorphic substitution E157Q in HIV-1 integrase increases R263K-mediated dolutegravir resistance and decreases DNA binding activity. J Antimicrob Chemother 2016;71(8):2083–2088. Available from:

  57. 57.

    Shimura K, Kodama E, Sakagami Y, Matsuzaki Y, Watanabe W, Yamataka K, et al. Broad antiretroviral activity and resistance profile of the novel human immunodeficiency virus integrase inhibitor elvitegravir (JTK-303/GS-9137). J Virol. 2008;82(2):764–74.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  58. 58.

    Fun A, Van Baelen K, van Lelyveld SFL, Schipper PJ, Stuyver LJ, Wensing AMJ, et al. Mutation Q95K enhances N155H-mediated integrase inhibitor resistance and improves viral replication capacity. J Antimicrob Chemother 2010;65(11):2300–2304. Available from:

  59. 59.

    Blanco J-L, Varghese V, Rhee S-Y, Gatell JM, Shafer RW. HIV-1 integrase inhibitor resistance and its clinical implications. J Infect Dis. 2011;203(9):1204–14.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  60. 60.

    Jones GS, Yu F, Zeynalzadegan A, Hesselgesser J, Chen X, Chen J, et al. Preclinical evaluation of GS-9160, a novel inhibitor of human immunodeficiency virus type 1 integrase. Antimicrob Agents Chemother. 2009;53(3):1194 LP–1203 Available from:

    Article  CAS  Google Scholar 

  61. 61.

    Hachiya A, Kirby KA, Ido Y, Shigemi U, Matsuda M, Okazaki R, et al. Impact of HIV-1 Integrase L74F and V75I mutations in a clinical isolate on resistance to second-generation Integrase Strand transfer inhibitors. Antimicrob Agents Chemother. 2017;61(8):e00315–7 Available from:

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  62. 62.

    Mulu A, Maier M, Liebert UG. Lack of integrase inhibitors associated resistance mutations among HIV - 1C isolates. J Transl Med. 2015;13:1–8.

    CAS  Article  Google Scholar 

  63. 63.

    Li Y, Xuan S, Feng Y, Yan A. Targeting HIV-1 integrase with strand transfer inhibitors. Drug Discov Today. 2015;20(4):435–49.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  64. 64.

    Mouscadet JF, Desmaële D. Chemistry and structure-activity relationship of the styrylquinoline-type HIV integrase inhibitors. Molecules. 2010;15(5):3048–78 Available from:

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  65. 65.

    Cai MI, Zheng RI, Caffrey M, Craigie R, Clore GM, Gronenborn AM. Solution structure of the N-terminal zinc binding domain of HIV-1 integrase. Nat Struct Biol. 1997;4(1):4–5.

    Google Scholar 

Download references


Appreciations goes to the collaborators from South African National Bioinformatics Institute (SANBI), University of the Western Cape, South Africa and University of Yaoundé I, Cameroon. A special acknowledgement goes to the late Dr. Graeme Jacobs who was one of the main driving forces behind this manuscript. His enormous contribution to the HIV drug resistance field will be surely missed. May He rest in peace.


The study was supported by the National Research Foundation (NRF) of South Africa, Poliomyelitis Research Foundation (PRF) of South Africa, Harry Crossley Foundation, South African Research Chairs Initiative of the Department of Science, Technology (DST). South African Medical Research Council (SAMRC) and The Higher Education Department, next Generation of Academic Programme (nGAP), provided support for this study in the form of full time academic positions and salaries to R. Cloete and GBJ. Mr. Darren Isaacs and Miss Rumbidzai Chitongo were funded by the South African Research Chairs Initiative of the Department of Science and Innovation (DSI) and National Research Foundation (NRF) of South Africa, award number UID 64751. The Funders had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information




GBJ and RC2 conceived and designed the study. GMI provided all the sequencing data. SGM conducted sequence analyses and phylogenetic analysis while DI, RC1 and RC2 assisted with the structural analyses. SGM wrote the first draft of the manuscript and reviewed by RC1, GMI, GBJ, DI and RC2. All the authors approved the final version of the manuscript.

Corresponding author

Correspondence to Ruben Cloete.

Ethics declarations

Ethics approval and consent to participate

This study form part of a broader study, for which the ethic was obtained from the Health Research Ethics Committee of Stellenbosch University (N14/10/130 and N15/08/071). The database sequences were conducted according to the ethical guidelines and principles of the international Declaration of Helsinki 2013, South African Guidelines for Good Clinical Practice and the Medical Research Council (MRC) Ethical Guidelines for Research.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Mikasi, S.G., Isaacs, D., Chitongo, R. et al. Interaction analysis of statistically enriched mutations identified in Cameroon recombinant subtype CRF02_AG that can influence the development of Dolutegravir drug resistance mutations. BMC Infect Dis 21, 379 (2021).

Download citation


  • Cameroon
  • Integrase
  • CRF02_AG
  • Molecular modelling
  • Interaction analysis