Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Advertisement
Scientific Reports volume 12, Article number: 16936 (2022)
314
Metrics details
Single nucleotide polymorphisms detected in the solute carrier member family-22 has been shown to result in a variable response in the treatment of type 2 diabetes mellitus with Metformin. This study predicted a three-dimensional protein structure for the SLC22A2 protein sequence using AlphaFold 2 and modelled five haplotypes within SLC22A2 protein structure observed in the Xhosa population of South Africa. The protein models were used to determine the effect(s) of haplotype variations on the transport function of Metformin and 10 other drugs by the SLC22A2 protein. Molecular dynamic simulation studies, molecular docking and interaction analysis of the five SLC22A2 haplotypes were performed in complex with the ligand 5RE in a POPC lipid bilayer to understand the mechanism of drug binding. Weakest binding free energy was found between 5RE and haplotype 1. Molecular docking studies indicated the top binding ligands as well as Metformin to bind inside the transport channel in all haplotypes increasing the probability of Metformin inhibition during co-administration of drugs. Metformin showed reduced binding affinity and number of interactions compared to the top four binding molecules. Molecular dynamic simulation analysis indicated that haplotypes 1, 3 and 4 were less stable than 2 and 5. The findings suggest haplotypes 4 and 5 having stronger preference for large inhibitor molecule binding in the active site and this could result in haplotypes 4 and 5 demonstrating reduced Metformin clearance via the SLC22A2 transporter during co-administration of drugs. The current study is the first to investigate the potential effect(s) of haplotype variation on the protein structure of SLC22A2 to assess its ability to transport Metformin in an indigenous South African population.
The incidence of diabetes mellitus (DM) is constantly on the rise across the world. In the African continent, South Africa contributes 7% of DM burden1. Type 2 diabetes mellitus (T2DM) accounts for ~ 90% of all DM cases and the preferred first line treatment for the management of T2DM is Metformin. Nonetheless, variability in drug response has resulted in several other classes of drugs being prescribed to treat T2DM, including: sulfonylureas, meglitinides, thiazolidinediones, α-glucosidase inhibitors, dipeptidyl peptidase-4 inhibitors, glucagon-like peptide-1 agonist, sodium glucose cotransporter-2 inhibitors, insulin and its analogues2,3,4.
The variability in response to the treatment of T2DM has been linked to genetic and environmental factors5,6. Genetic polymorphisms in candidate genes involved in drug absorption, transportation, distribution, and excretion (ADME) contribute to the observed heterogeneity in the treatment and management of T2DM2,7. Genetic variants in the solute carrier family (SLC) are amongst the drug related genes that have altered response in the treatment of T2DM and have therefore been selected for further examination8.
Members of the SLC22 family share a common structure of 12 transmembrane domains (TMDs), a large extracellular loop with glycosylation sites between domains 1 and 2, and a intracellular loop with phosphorylation sites between domains 6 and 79. The SLC22A2 gene encodes for the organic cation transporter 2 (OCT2) isoform which acts as a renal uptake transporter10. OCT2 is expressed on the basolateral side of the proximal tubule cells and plays a crucial role in the disposition and renal clearance of cationic drugs and endogenous compounds11. OCT2 is a poly-specific transporter that facilitates the bi-directional transport of structurally unrelated small organic cations down their electrochemical gradients12,13. It transports positively charged drugs and since most prescribed drugs belong to the organic cations which are positively charged, the function of it as a transporter has pharmacological consequences.
Many drugs (ligands) have been identified as substrates and/or inhibitors of SLC22A2 including the oral anti-diabetic drug Metformin and the chemotherapeutic agent Cisplatin14,15,16. Metformin (C4H11N5) has a molecular weight of 129.16 g/mol, whilst Cisplatin (Cl2H6N2Pt) which is larger, has a molecular weight of 300 g/mol17. Larger ligands are more likely to interact with more residues at the same binding pocket of a protein than smaller ligands inferring that the protein would preferentially bind larger ligands. However, smaller ligands can have more options in binding pockets on a protein structure. Furthermore, the binding of ligands to proteins are further affected by amino acid (aa) changes18. Amino acid (aa) residue(s) substitutions/mutations can affect protein conformation and the function of the protein regardless of their location18. Since single nucleotide polymorphisms (SNPs) are most likely to occur in combinations, the incorporation of haplotypes in pharmacogenetic studies is becoming more relevant in the practice of genetic medicine at both an individual and population level19.
A computational approach, using molecular docking, was previously used to determine the effect(s) of the common non-synonymous SNP (nsSNP) rs316019 (A270S) on the interaction of Metformin and other drugs within the SLC22A2 protein20. The study determined that the WT protein was more efficient at ligand binding than the variant since it had a more open and wider space in its binding pocket20. In addition, the study also explored the two haplotypes, i.e. WT and the A270S variant. The current study aims to explore these haplotypes including additional haplotypes which have been identified in a sub-Saharan African population, i.e. the Xhosa population of South Africa. This study will determine if the haplotypes identified in the Xhosa population of a South African cohort could alter the transport function of SLC22A2. This will be achieved by building structural models for each of the five haplotype sequences that have been observed in the Xhosa population. Molecular dynamic simulation studies were used to understand the difference in the binding of inhibitor 5RE to the five haplotype structures in a POPC lipid bilayer. Furthermore, molecular docking and interaction analysis of SLC22A2 substrates &/or inhibitors (i.e. Cimetidine, Creatinine, Dolutegravir, Isavuconazole, Metformin, Ranitidine, Ranolazine, Trimethoprim and Vandetanib) will provide novel insights into the mechanism of binding of these substrates and inhibitors to SLC22A2.
Previously, Jacobs et al. (2015) characterised the frequency of various haplotypes in the SLC22A2 gene that exists within the South African Xhosa population21. Three nsSNP, i.e. rs316019 results in aa substitution (Ala270Ser), rs8177516 in aa substitution (Arg400Cys) and rs8177517 in aa substitution (Lys432Gln), were detected. These three nsSNPs characterises the five identified haplotypes which were found in the study cohort.
The phenotypic effects of three aa substitutions, (i.e. Ala270Ser, Arg400Cys and Lys432Gln), on the structure and function of SLC22A2 protein sequence was investigated using a variety of algorithms that included: the Sorting Intolerant from Tolerant (SIFT)22,23,24, the Polymorphism Phenotyping 2 (PolyPhen-2)25,26,27 and the I-mutant algorithms28. The SIFT program uses sequence homology to predict the potential effect(s) an aa substitution could have on protein function22. It calculates a tolerance index score ranging from 0 (deleterious) to 1 (neutral) based on multiple sequence alignments. Variants with scores ranging from 0.00 to 0.049 are considered deleterious, with scores closer to 0.00 being more confidently predicted as deleterious. Variant scores ranging between 0.05 and 1.00 are predicted as neutral or tolerated with scores closer to 1.00 more confidently predicted as neutral. PolyPhen-2 predicts the possible impact of an aa substitution on protein structure and function by physical and evolutionary comparisons27. It categorizes aa substitutions from 0 (neutral) to 1 (deleterious). In addition, it also provides the functional significance of each aa substitution as: benign (0.00–0.14), possibly damaging (0.15–0.84) and probably damaging (0.85–1). The I-mutant algorithm, predicts the effect(s) of single point mutations on protein folding using a support vector machine-based approach that includes either sequence or structural information28. The predicted I-mutant score is a reliability index ranging from 0 (least likely to result in the altered protein stability as predicted) to 8 (most likely to result in the altered protein stability as predicted)28,29.
The SLC22A2 primary aa sequence was retrieved and downloaded from the UniProtKB database with the accession number O15244. The protein structure of SLC22A2 has not been experimentally resolved. Therefore homology modelling methods were employed to predict the three-dimensional (3D) structure of SLC22A2 using the online 3D structure prediction webserver AlphaFold 2 Protein Structure Database30,31. AlphaFold 2 is an artificial intelligence system developed to predict the 3D structure of a protein from its aa sequence30. It constructs a multiple sequence alignment by comparing the input amino acid sequence to several protein sequence databases. AlphaFold 2 identifies similar sequences to the input aa sequence and extracts the information to produce an accurate 3D structure model30.
The generated protein model was further assessed using internal scores to validate the quality of the predicted model. Since the SLC22A2 protein is a transmembrane protein, the local model quality estimation method for membrane proteins (QMEANBrane) score was determined32. The QMEANBrane score is a version of QMEAN developed to assess the local quality of alpha-helical transmembrane protein models32. It takes into the consideration the three different segments in a transmembrane protein (membrane, interface and soluble) model. Values are normally between 1 and 0, with higher values indicating more reliable models32. To determine if the correct fold/conformation was assigned to the protein sequence a structural alignment was performed to calculate a root mean square deviation (RMSD). The lower the RMSD value, the closer the structures are in terms of backbone deviation. External quality analysis was done using the Structure Assessment and Verification Server (SAVES) webserver and included tools such as ERRAT and Procheck. ERRAT statistically interrogates the non-bonded atomic interactions of the given protein model against the interactions of well-refined structures33. Higher overall quality factors of 80% and more suggest very few errors within the protein model, suggesting it is highly reliable. Procheck generates a Ramachandran plot and determines if the phi and psi dihedral angle distribution of protein residues satisfy conformations found in crystal structures of high resolution34,35. Usually, more than 90% of residues in a protein structure needs to be located in the most favourable regions of the plot to be considered an accurate protein structure.
Five simulation systems (i.e. haplotypes 1 to 5) were prepared in total for SLC22A2 in complex with inhibitor (2~{S}-3-(4-fluorophenyl-2-[2-(3-hydroxyphenyl)ethanoylamino]-~{N}-[(1~{S}-1-phenylethyl] propenamide (PDBID: 5RE) using Pymol by aligning the template structure and SLC22A2 predicted structure and extracting the compound 5RE36. The variant haplotype structures were prepared using the PyMol mutagenesis wizard where: haplotype 1 (WT) had no aa change, haplotype 2 (A270S), haplotype 3 (A270S and R400C), haplotype 4 (A270S and K432Q) and haplotype 5 (R400C). The five haplotype complex systems were embedded into a heterogeneous lipid bilayer consisting of 120 1-palmitoyl-2-oleoylphosphatidylcholine (POPC) lipids. The lipid bilayer consisted of 60 lipids in the upper leaflet and 60 in the lower leaflet with a water thickness of 22.5 Å from the top leaflet and the bottom leaflet. The POPC lipid bilayer is one of the most abundant lipids found in human transmembrane proteins. The CHARMM-GUI membrane builder was used to orient the protein relative to the POPC lipid bilayer membrane37. The proteins of each system were aligned along the first principal axis the Z-axis of the membrane with the protein in the center (Z = 0). All five simulations were performed using the GROMACS-2019 package38 along with the CHARMM36M all-atom force field39. The accurate topologies for the ligand was generated using ParamChem/CHARMM general force field40. All the systems were solvated with TIP3 water molecules in a cubic box of at least 10 Å of water between the protein and edges of the box at a concentration of 0.15M. To neutralize the positive and negative charges of the systems for haplotype one and two, 63 potassium (K) ions and 70 chloride (Cl) ions were added to neutralize the charge of the system. While, haplotype three, four and five had 62 K ions and 68 Cl ions, respectively.
Each system underwent 50,000 steps of steepest descents energy minimization to remove close Van der Waals force contacts. Subsequently, all systems were subjected to a two-step equilibration phase namely; NVT (constant number of particles, Volume and Temperature) for 500 ps to stabilize the temperature of the system and a short position restraint NPT (constant number of particles, Pressure and Temperature) for 500 ps to stabilize the pressure of the system by relaxing the system and keeping the protein restrained. For the NVT simulation the system was gradually heated by switching on the water bath and the V-rescale temperature-coupling method was used, with constant coupling of 0.1 ps at 300 K under a random sampling seed. While for NPT the Parrinello-Rahman pressure coupling41 was turned on with constant semi-isotropic pressure coupling for uniform scaling of x, y box vectors of 0.1 ps at 300 K under conditions of position restraints (all-bonds). For both NVT and NPT, electrostatic forces were calculated using the Particle Mesh Ewald method42 and all systems were subjected to a full 200 ns simulation.
The analysis of the trajectory files was done using GROMACS utilities. The root-mean square deviation (RMSD) was calculated using gmx rmsd for the protein back bone atoms and the root mean square fluctuation (RMSF) analysis for the protein residues was calculated using gmx rms. The average number of hydrogen bonds was calculated using the gmx hbond tool, while the distance between the drug and the protein was calculated using gmx pairdist. The non-bonded interaction energy was calculated between the protein and the drug using (gmx energy) and the free energy of binding was calculated between the protein and the drug 5RE using the Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) protocol implemented in g_mmpbsa package over the last 300 frames of the simulation trajectory43. We characterized the simulated POPC lipid bilayer in complex with the protein by calculating several parameters. This lipid bilayer analysis comprised of different calculations with the corresponding GROMACS utility indicated in parentheses: the area per lipid (gmx energy) to analyze the X and Y box dimensions, the bilayer thickness (gmx density) for the phosphate atom headgroups, the lateral diffusion coefficients (gmx mean square displacement (msd)) for the lipid head groups, the deuterium order parameters for all carbons in the acyl chains was calculated using g_lomepro44. The average number of hydrogen bonds were also calculated between the protein and POPC lipids (gmx hbond), and the non-bonded interaction energy (gmx energy) between the protein and the POPC lipid. The final snapshot of the simulation systems was extracted at 200 ns and were visually illustrated using PyMol.
The equilibrated structural conformations were extracted at 40 ns simulation time period for each of the simulation trajectories using gmx trjconv tool. Eleven compounds, including 5RE were docked to the conformational snapshot structures. The eleven drug structures (ligands) were retrieved from Pubchem17. They were selected because they were identified as known cationic compounds to be transported by SLC22A211. Metformin (a SLC22A2 substrate), Glibenclamide and Gliclazide (substrates of CYP2C9) are commonly prescribed for the treatment and management of T2DM45,46. The remaining ligands, (i.e. Dolutegravir, Isovuconazole, Vandetanib, Ranolazine, Trimethoprime, Cimetidine and Creatinine) are drugs prescribed in the treatment of a range of diseases and illnesses, and are known to be substrates &/or inhibitors of SLC22A2. Co-administration of drugs can influence the efficacy and action of the anti-diabetic drug Metformin.
Molecular docking of the drug substrate 5RE and the selected 11 ligands to the five haplotype variant structures of the SLC22A2 protein was performed using Autodock Vina47. All haplotype receptors were prepared using AutoDock tools, which add polar hydrogens, Gasteiger charges and saves the output file in pdbqt format. The twelve ligands were prepared with automated bash and python scripts namely, split_multi_mol2_file.py and prepare_ligand4.py48. This was done to correct for errors such as missing atoms, added H2O, more than one molecule chain break and alternate locations. After receptor and ligand preparations, the centre of mass was calculated for the ligand 5RE in complex with each of the receptor structures to determine the X, Y and Z coordinates of the docking grid space using a python script center_of_mass.py. Various parameters for the docking process were stored in a configuration file. The configuration file contained all the input parameters for all docking simulations, such as: centre of mass coordinates, grid space and exhaustiveness of the search algorithm (grid box dimensions are provided in Table 1). Interaction analysis was performed using PoseView web31. The docked complexes were ranked according to their energy scores using the python script developed by the Scripps Research Institute vina_screen_get_top.py. All the docked poses were visually inspected in PyMol and were further analysed using PoseView31,49. PoseView determines four types of interactions namely; (i) hydrogen bonds, (ii) hydrophobic, (iii) metal interactions and (iv) π-stacking interactions.
SLC22A2 gene was characterised to contain five haplotypes in the South African Xhosa population. Haplotype structure and frequency within the study population is shown in Table 2. The five identified haplotypes are defined by three nsSNPs, i.e. rs316019 (A270S); rs8177516 (R400C) and rs8177517 (K432Q). Pairwise sequence alignments of SLC22A2 protein and the identified variants are shown in Fig. 1.
Pairwise protein sequence alignments of SLC22A2 variants: (a) rs316019 (A270S), (b) rs8177516 (R400C) and (c) rs8177517 (K432Q) (variant position indicated by red arrow).
The functional effects of the identified aa substitutions on protein structure and function were assessed using SIFT, PolyPhen-2 and I-mutant algorithms (Table 3). Both SIFT and PolyPhen-2 predicted that the function of SLC22A2 would be affected by the presence of the R400C variant, whilst the protein function would most likely not be affected by the presence of the K432Q variant. However, these algorithms showed contradictory results for the remaining A270S variant. The I-mutant algorithm which measures the degree of protein destabilization predicted that all of the identified aa substitutions would decrease protein stability.
A 3D structure was generated by AlphaFold 2 for the SLC22A2 protein (Fig S1). Figure 2 shows the outward facing conformation of the 3D structure of SLC22A2 and Fig. 3 shows the locations of the three variants relative to the drug binding site identified in the South African Xhosa population. The overall topology (consisting of alpha helices and beta sheets) of SLC22A2 is representative of a transmembrane protein involved in drug transport. The QMEANBrane predicted the generated template to be within the expected range (0–1) for a transmembrane structure. The Ramachandran plot indicated that 94.4% of residues had favourable dihedral angle conformations and the ERRAT overall quality factor score was 97.55% suggesting only ~ 2% error rate within the protein model.
3D structure of SLC22A2 protein embedded in the POPC lipid bilayer membrane: the SLC22A2 protein structure (coloured by secondary structure) in complex with the 5RE inhibitor (shown as spheres in magenta) phosphate lipid head groups shown as orange spheres.
3D structure of SLC22A2 in complex with 5RE generated by AlphaFold 2: the three mutated residues identified in the South African Xhosa population are shown as sticks (coloured deep-teal) and the 5RE inhibitor in magenta. Solid black arrows point to the locations of the mutations identified to characterise the haplotypes existing in the study population.
All five SLC22A2-5RE haplotype systems reached equilibrium after 30 ns based on the backbone RMSD values (Fig. 4A). The mean and standard deviation values for the change in protein RMSD backbone atoms for haplotype 1 (WT-5RE), haplotype 2, 3, 4 and 5 were 0.60 ± 0.12 nm, 0.54 ± 0.13 nm, 0.60 ± 0.14 nm, 0.50 ± 0.13 nm and 0.44 ± 0.06 nm, respectively (Fig. 4A). The RMSD values were the highest for the haplotype systems 1 and 3 compared to haplotype systems 2, 4 and 5. Furthermore, the RMSF fluctuation values were the lowest for haplotypes 5, and 3 each having 0.17 ± 0.12 nm, 0.22 ± 0.23 nm, while haplotypes 2 had an RMSF value of 0.24 ± 0.21 nm compared to haplotypes 1 and 4 each having the largest RMSF value of 0.29 ± 0.22 nm (Fig. 4B). The protein residues of haplotypes 1 and 4 showed four regions (residues R1:62-65, R2:80-85, R3:99-106 and R4:108-113), of high flexibility, while none of the other haplotypes showed any regions of high flexibility (Fig 4B). Interestingly, none of the mutated residues and the most important binding site residues showed high flexibility values, suggesting the active site is stable and not undergoing large conformational changes (Fig.4B).
Statistical analysis of haplotypes over 200 ns simulation. (A) The backbone RMSD deviation of SLC22A2-5RE in the five observed haplotypes [ordinate is RMSD (nm) and the abscissa is time (ns)]. (B) RMSF deviation for protein residues for the five haplotype systems over the last 100 ns. Regions of high flexibility are boxed in red and labelled as R1, 2, 3 and 4 on the figure.
The distance values between the protein residues and 5RE did not show any significant differences between the five systems. The average number of hydrogen bonds formed between SLC22A2-5RE was the highest for haplotype 2 with 2.29 contacts being formed compared to haplotypes 1, 4, 5 and 3 each having 1.53, 1.07, 0.19 and 0.15 bonds, respectively. The non-bonded interaction energy calculated between the protein and the drug was the highest for haplotype 1 (−261.67 ± 42.14 kJ/Mol) followed by haplotype 2, 4 and 3 each having −237.85 ± 30.62 kJ/Mol, −217.18 ± 35.93 kJ/Mol and −215.40 ± 39.80 kJ/Mol while haplotype 5 had the lowest non-bonded interaction energy with a score of −166.35 ± 36.17 kJ/Mol. However, contrasting results were observed for the binding free energy calculations for the different complexes. Here, the total free energy of binding calculated for haplotype 5 was the highest with −106.05 ± 14.32 kJ/Mol, followed by haplotype 4 having a total binding free energy score of −96.79 ± 14.86 kJ/Mol (Table 4). Haplotypes 3 and 2 demonstrated total binding free energies of −93.62 ± 1.74 kJ/Mol and −92.74 ± 14.96 kJ/Mol respectively, while haplotype 1 had the weakest total binding free energy with 5RE of −85.69 ± 14.05 kJ/Mol (Table 4). In summary, the inhibitor 5RE bound the strongest to haplotypes 5 and 2 compared to haplotype 1 (WT) (Table 4).
The average area per lipid reached convergence and fluctuated around a stable mean and standard deviation values for each system ranging between 0.79 ± 0.01 and 0.81 ± 0.01 nm (Fig. 5 and Table 5). The means and SDs for each of the various lipid bilayer properties are summarized in Table 5. The bilayer thickness values was found to shift towards smaller values for haplotype 4 (8.42 ± 4.91 nm) followed by haplotype 5 (8.89 ± 5.18 nm), then haplotype 3 had the largest values (9.02 ± 5.26 nm) and haplotype 2 (8.90 ± 5.18 nm), compared to haplotype 1 having (8.89 ± 5.18 nm). The lateral diffusion coefficient was the lowest for the haplotype 3 (0.86 ± 0.28 cm2/s) system, followed by haplotype 5 (1.01 ± 0.22 cm2/s), then haplotype 1 (1.16 ± 0.20 cm2/s) and haplotype 4 (1.31 ± 0.20 cm2/s), while haplotype 2 had the highest, reaching 1.33 ± 0.14 cm2/s.
Area per lipid fluctuation for the five haplotype systems over 200 ns.
The deuterium order parameters for the POPC carbon acyl chains sn-1 and sn-2 indicated similar results (Figs. S2–S6). For all the systems of chain 1, atoms 4 and 5 had the lowest deuterium order parameter values while atoms 2, 3, 6–8 had higher values, suggesting that atoms 2, 3, 6–8 are not undergoing any large structural changes. For chain 2, the atoms for haplotypes 2, 3, 4 and 5 showed similar results with terminal atoms (6 and 7) had lower deuterium order parameter values, while the central atoms 2–5 had higher values, suggesting that the central atoms are maintaining their structure, while haplotype 1 displayed low deuterium order parameter values for atoms 3, 4 and 6 and higher values for atoms 2, 5, 7 and 8. Calculating the average number of hydrogen bonds formed between the protein and POPC lipids indicated that haplotype 2 formed 27.90, haplotype 4 formed 27.81 and haplotype 3 formed 27.51. Haplotype 5 formed the least amount of hydrogen bonds compared to haplotype 1 each having 14.71 and 22.87, respectively. Interestingly, the non-bonded interaction energy results were similar to the hydrogen bond analysis with haplotype 5 showing the weakest total energy of −4727.42 ± 380.49 KJ/Mol and haplotypes 4 and 2 the largest values of −5031.12 ± 354.91 KJ/Mol and −4998.96 ± 391.41 KJ/Mol compared to haplotype 1 having −4982.29 ± 363.99 KJ/Mol. The final 200 ns snapshot structures for the five haplotype conformations were overlaid and showed similar outward facing conformations for the inner membrane structures with the largest difference being in the outer head region of the protein (Fig. S7).
The molecular docking of the 11 drugs to the equilibrated structures at 100 ns for the five structures were performed to determine any difference in drug binding affinity and number of interactions between the different haplotypes. Some drugs are known anti-diabetic compounds such as Metformin, Glibenclamide and Gliclazide. The remaining ligands, (i.e. Cimetidine, Creatinine, Dolutegravir, Isavuconazole, Ranitidine, Ranolazine, Trimethoprim and Vandetanib), can be classified as substrates &/or inhibitors of SLC22A2. In haplotypes 1 and 2, four drugs (Isavuconazole, Dolutegravir, Glibenclamide and Gliclazide), whilst in haplotypes 3 and 4 (Glibenclamide, Dolutegravir, Isavuconazole and Vandetanib), in haplotype 5 (Glibenclamide, Dolutegravir, Isavuconazole and Ranolazine) and 5RE consistently ranked among the top binders based on higher binding affinity and number of interactions compared to Metformin that ranked second lowest (Tables 6, 7, 8, 9 and 10). Creatinine demonstrated the lowest docking scores for all five haplotypes.
Important binding residues for the four top binding ligands and the 5RE inhibitor in haplotype 1 included Trp218, Tyr447 and Gln242 (Table 6 and Fig. 6). The residue Tyr447 was responsible for the strong binding of the four top binding ligands and the 5RE inhibitor in haplotypes 2 and 4 (Tables 7, 9 and Figs. 7, 9). In haplotype 3 aa residues Tyr245 and Ser471 featured as important for binding of the top four ligands and the 5RE inhibitor (Table 8 and Fig. 8). The residues Gln242 and Tyr447 were identified as essential in ligand binding in haplotype 5 (Table 10 and Fig. 10).
The study showed that Glibenclamide, Gliclazide, Dolutegravir, Isavuconazole, Vandetanib as well as Ranolazine to have stronger affinity for the SLC22A2 transporter than Metformin and the other six compounds (Tables 6, 7, 8, 9, 10). In addition to this, Metformin showed only hydrogen bond interactions in all haplotypes (Tables 6, 7, 8, 9, 10). In all five haplotypes, the top four ligands, Metformin and the 5RE inhibitor bound in the active site (Figs. 6, 7, 8, 9 and 10).
Ligand overlay for haplotype 1 SLC22A2 protein. (A) Image was generated using Pymol with the first line anti-diabetic drug, Metformin, shown in blue. The four top binding ligands Isavuconazole, Dolutegravir, Glibenclamide, and Gliclazide are shown in deep teal, grey, hot-pink and green respectively. The 5RE inhibitor is indicated in magenta. (B,C) 2D PoseView illustration, showing interacting residues for Metformin and the 5RE inhibitor respectively—dotted black lines denote hydrogen bonds, solid green lines denote hydrophobic interactions and dotted green lines denote π stacking &/or aromatic interactions.
Ligand overlay for haplotype 2 SLC22A2 protein. (A) Image was generated using Pymol with the first line anti-diabetic drug, Metformin, shown in blue. The four top binding ligands Isavuconazole, Dolutegravir, Glibenclamide, and Gliclazide are shown in deep teal, grey, hot-pink and green respectively. The 5RE inhibitor is indicated in magenta. (B,C) 2D PoseView illustration, showing interacting residues for Metformin and the 5RE inhibitor respectively—dotted black lines denote hydrogen bonds, solid green lines denote hydrophobic interactions and dotted green lines denote π stacking &/or aromatic interactions.
Ligand overlay for haplotype 3 SLC22A2 protein. (A) Image was generated using Pymol with the first line anti-diabetic drug, Metformin, shown in blue. The four top binding ligands Dolutegravir, Glibenclamide, Vandetanib and Isavuconazole are shown in grey, hot-pink, orange and deep teal respectively. The 5RE inhibitor is indicated in magenta. (B,C) 2D PoseView illustration, showing interacting residues for Metformin and the 5RE inhibitor respectively—dotted black lines denote hydrogen bonds, solid green lines denote hydrophobic interactions and dotted green lines denote π stacking &/or aromatic interactions.
Ligand overlay for haplotype 4 SLC22A2 protein. (A) Image was generated using Pymol with the first line anti-diabetic drug, Metformin, shown in blue. The four top binding ligands Isavuconazole, Dolutegravir, Glibenclamide, and Vandetanib are shown in deep teal, grey, hot-pink, and orange respectively. The 5RE inhibitor is indicated in magenta. (B,C) 2D PoseView illustration, showing interacting residues for Metformin and the 5RE inhibitor respectively—dotted black lines denote hydrogen bonds, solid green lines denote hydrophobic interactions and dotted green lines denote π stacking &/or aromatic interactions.
Ligand overlay for haplotype 5 SLC22A2 protein. (A) Image was generated using Pymol with the first line anti-diabetic drug, Metformin, shown in blue. The four top binding ligands Dolutegravir, Glibenclamide, Isavuconazole and Ranolazine are shown in grey, hot-pink, deep teal, and yellow respectively. The 5RE inhibitor is indicated in magenta. (B,C) 2D PoseView illustration, showing interacting residues for Metformin and the 5RE inhibitor respectively—dotted black lines denote hydrogen bonds, solid green lines denote hydrophobic interactions and dotted green lines denote π stacking &/or aromatic interactions.
3D molecular modelling is routinely used in the discovery, development and design of drugs. Statistical data obtained from multiple protein conformations are considered more reliable than data obtained from single static structures, since proteins can occur in different conformations50. The predicted 3D protein model of SLC22A2 was built using a template generated from multiple aa sequence alignments using AlphaFold 2 (Fig. S1). The constructed 3D protein model of SLC22A2 successfully satisfied several quality checks suggesting the protein model is reliable and accurate for use in molecular dynamic simulations and docking studies (Fig. 1). The predicted docking scores of ligands ranged between −10.4 and −3.8 kcal/Mol across haplotypes (Tables 6, 7, 8, 9 and 10). Larger ligands (such as: Dolutegravir, Isavuconazole, Ranolazine, Trimethoprim, Vandetanib, Glibenclamide and Gliclazide) tend to show higher docking scores in comparison to smaller ligands (such as: Metformin, Cimetidine, Ranitidine and Creatinine). This can be attributed to larger ligands having more residues available to interact with in the binding pocket. Thus, it was observed that Metformin which is a small ligand, tends to have a lower docking score in comparison to its counter-parts Glibenclamide and Gliclazide in all haplotypes (Tables 6, 7, 8, 9 and 10). Glibenclamide and Gliclazide are second-generation sulfonylureas used in the treatment of T2DM. Mutations in the CYP2C9 gene have been shown to influence the efficacy of their function45,46. They are commonly used in combination therapy with Metformin but can also be used as monotherapy. Combination therapy of Metformin and Glibenclamide as well as Metformin and Gliclazide has been shown to be more effective in the management of T2DM in patients who have not responded effectively to monotherapies available51,52,53.
The MD results showed that haplotypes 1, 2, 3 and 4 had higher RMSD and RMSF values suggesting they were less stable and more flexible compared to haplotype 5 being more stable based on lower RMSD and RMSF values. However, in contrast to the binding free energy haplotype 2 showed a higher average number of hydrogen bonds and non-bonded interaction energy with 5RE compared to haplotype 1, whilst haplotypes 3 and 5 only showed the least amount of average number of hydrogen bonds and non-bonded interaction energy compared to haplotype 1 and 4. The non-bonded interaction energy should be interpreted with caution as it excludes the polar solvation terms and solvent accessible surface area terms which is included in the free energy of binding calculation. Analysis of the POPC lipid bilayer indicated that the lipid structures for all five haplotype systems remained stable throughout the simulation, based on area per lipid and deuterium order parameters. Haplotypes 1, 2, 3 and 5 showed very little difference in membrane thickness with only haplotype 4 displaying smaller membrane thickness and larger lateral diffusion values suggesting higher mobility of lipids resulting in stronger POPC and protein interaction. Additionally, haplotypes 2, 3 and 4 showed higher binding affinity to the POPC lipid bilayer based on the number of hydrogen bonds and non-bonded interaction energy while haplotype 5 displayed weaker binding affinity to the POPC lipid bilayer.
Furthermore, haplotypes 4 and 5 also demonstrated the highest number of interactions with inhibitor 5RE and free binding energy compared to haplotype 1. For the docking of metformin to the five haplotypes again haplotype 4 and 5 higher docking scores compared to haplotype 1 and only haplotype 4 showed slightly higher number of interactions with haplotype 4 residues. Overall the docking results for the top four compounds to haplotypes 4 and 5 showed higher docking scores and in some instances higher number of interactions compared to haplotype 1. However, the binding interactions observed between proteins and ligands needs to be validated experimentally using X-ray or Cryo-EM methods. One of the most common SNPs of SLC22A2 is rs316019 (A270S variant) which displays a minor allele frequency (MAF) of 80% in all ethnicities and has been associated with varied transport activity54. Sajib et al (2018) demonstrated through in silico analysis that the 270A variant of SLC22A2 (i.e. WT) was better at binding Metformin and other drugs in comparison to the 270S variant (i.e. haplotype 2) since it had a more open and wider active site conformation20. Furthermore, three studies by Yoon et al. (2013), Wang et al. (2008) and Song et al. (2008) demonstrated individuals with the haplotype 2 genotype (270S) to produce a lowered Metformin clearance rate in comparison to individuals with the WT haplotype in Korean and Chinese subjects55,56,57. The above-mentioned studies could imply that the reduced Metformin clearance observed in the 270S genotype can be attributed to the smaller active site (as predicted by Sajib et al (2018)) making Metformin binding more efficient.
However, in contrast, Chen et al. demonstrated the 270S variant to produce a higher clearance of Metformin compared to the WT in Caucasian and African American ancestries58. In addition to this, Dujíc et al., Christensen et al. and Tzvetkov showed no statistically significant association between the 270S variant and glycaemic response to Metformin monotherapy59,60,61. The data generated herein demonstrates that haplotypes 4 and 5 could possibly have a stronger affinity for ligands as shown in our MD study with 5RE binding free energy, suggesting Metformin clearance can be reduced based on a higher number of interactions. Regardless of the contradiction in data, in population groups where polymorphism(s) may influence the response to Metformin therapy attention should be given to the phenotypic response of individual patients.
Haplotype 4 has the aa substitution (A270S and K432Q) and 5 (R400C) with substitutions A270S and R400C predicted to be deleterious and damaging to the protein function by SIFT and PolyPhen2 algorithms (Table 2). If we assume, that haplotypes 4 and 5 are most likely to cause a slower Metformin clearance rate because of their stronger affinity for ligand(s) in comparison to haplotype 1, individuals with these haplotypes are most likely to develop complications associated with Metformin toxicity. Drugs of the biguanide class, i.e. Metformin, if accumulated in the body cause an increase in plasma lactate level in a concentration-dependent manner by inhibiting mitochondrial respiration20,62,63. This could ultimately result in developing a life-threatening condition known as metformin-associated lactic acidosis (MALA)62,64. Therefore, dose adjustments of Metformin may prove beneficial for patients of haplotypes 4 and 5 in order to reduce the toxicity of Metformin and improve its efficacy. Trimethoprim has been shown to significantly reduce the clearance of both Metformin and Creatinine65. Creatinine is a by-product of muscle metabolism and creatinine serum levels are used as an indirect indicator of kidney function. Since Creatinine is a substrate of SLC22A2 secreted into the urine, the binding of inhibitors to SLC22A2 protein can cause a rise in serum creatinine levels that may be incorrectly interpreted as impaired kidney function54.
The findings in this study suggest haplotypes 4 and 5 have a higher affinity for drug binding based on binding free energy values and interaction analysis. However, a limitation in this study is the exclusion of entropy within binding free energy calculations using the g_mmpbsa tool and not including implicit membrane solute dielectrics for accurate polar solvation energy calculations. This will be considered in future studies. In some instances, such as with Metformin, the increased binding affinity can prove detrimental to the patient resulting in Metformin toxicity. Provided that organic cationic drugs are also used to treat various cancers and HIV, the results generated herein can be used to expand pharmacogenetic information available for African-specific populations, potentially impacting healthcare on the African continent. Future investigations will include cellular uptake analysis of Metformin to determine each haplotype’s transport activity identified for the SLC22A2 gene in the Xhosa population of South Africa.
In this study, the structural and functional effects of mutation(s) present within haplotypes were evaluated within the SLC22A2 protein. The protein models generated for each haplotype of the SLC22A2 gene indicated that haplotypes 4 and 5 of the Xhosa population of South Africa might impair the transport of Metformin via the SLC22A2 protein transporter. In addition to this, we identified other drugs which may outcompete or inhibit Metformin transport via SLC22A2. We also observed, that larger ligands tend to preferentially bind in the conserved active site location of SLC22A2. Based on these evaluations, it can be deduced that T2DM patients should be cautious with the usage of over-the-counter drugs since it may lead to Metformin accumulation in the blood, resulting in increased lactate concentration and the increased likelihood of the development of MALA. It is recommended that future studies should experimentally validate the binding affinity of each haplotype for Metformin. In addition, results may be used to develop and/or improve pharmacogenomic profiling systems used for individualized therapeutic approaches employed in the treatment and management of T2DM.
The protein sequence analysed during the current study is available in the UniProt repository with accession number O15244. The genetic polymorphisms analysed during the current study are available in the dbSNP respository with accession number rs316019, rs8177516 and rs8177517.
International Diabetes Federation (IDF). Diabetes Atlas 8th Edition 2017. https://www.idf.org/our-network/regions-members/africa/welcome.html. Accessed 15 July 2018 (2018).
Singh, S., Usman, K. & Banerjee, M. Pharmacogenetic studies update in type 2 diabetes mellitus. World J. Diabetes. 7, 302. https://doi.org/10.4239/wjd.v7.i15.302 (2016).
Article PubMed PubMed Central Google Scholar
Inzucchi, S. E. et al. Management of hyperglycaemia in type 2 diabetes: A patient-centered approach. Position statement of the American Diabetes Association (ADA) and the European Association for the Study of Diabetes (EASD). Diabetologia 55, 1577–1596. https://doi.org/10.1007/s00125-012-2534-0 (2012).
CAS Article PubMed Google Scholar
Topić, E. The role of pharmacogenetics in the treatment of diabetes mellitus. J. Med. Biochem. 33, 58–70. https://doi.org/10.2478/jomb-2013-0045 (2014).
Article Google Scholar
Chen, L. et al. Effect of lifestyle intervention in patients with type 2 diabetes: A meta-analysis. Metabolism 64, 338–347. https://doi.org/10.1016/j.metabol.2014.10.018 (2015).
CAS Article PubMed Google Scholar
Zia, A. Genetic susceptibility to type 2 diabetes and implications for therapy. J. Diabetes Metab. https://doi.org/10.4172/2155-6156.1000248 (2013).
Article Google Scholar
Mannino, G. C. & Sesti, G. Individualized therapy for type 2 diabetes: Clinical implications of pharmacogenetic data. Mol. Diagn. Ther. 16, 285–302. https://doi.org/10.1007/s40291-012-0002-7 (2012).
CAS Article PubMed Google Scholar
He, L., Vasiliou, K. & Nebert, D. W. Analysis and update of the human solute carrier (SLC) gene superfamily. Hum. Genomics. 3, 195–206. https://doi.org/10.1186/1479-7364-3-2-195 (2009).
CAS Article PubMed PubMed Central Google Scholar
Ciarimboli, G. & Schlatter, E. Regulation of organic cation transport. Pflugers Arch. Eur. J. Physiol. 449, 423–441. https://doi.org/10.1007/s00424-004-1355-5 (2005).
CAS Article Google Scholar
Saier, M. H. Families of transmembrane sugar transport proteins. Mol. Microbiol. 35, 699–710. https://doi.org/10.1046/j.1365-2958.2000.01759.x (2000).
CAS Article PubMed Google Scholar
Motohashi, H. & Inui, K. I. Organic cation transporter OCTs (SLC22) and MATEs (SLC47) in the human kidney. AAPS J. 15, 581–588. https://doi.org/10.1208/s12248-013-9465-7 (2013).
CAS Article PubMed PubMed Central Google Scholar
Koepsell, H., Lips, K. & Volk, C. Polyspecific organic cation transporters: Structure, function, physiological roles, and biopharmaceutical implications. Pharm. Res. 24, 1227–1251. https://doi.org/10.1007/s11095-007-9254-z (2007).
CAS Article PubMed Google Scholar
Brast, S. et al. The cysteines of the extracellular loop are crucial for trafficking of human organic cation transporter 2 to the plasma membrane and are involved in oligomerization. FASEB J. 26, 976–986. https://doi.org/10.1096/fj.11-180679 (2012).
CAS Article PubMed Google Scholar
Kimura, N. et al. Metformin is a superior substrate for renal organic cation transporter OCT2 rather than hepatic OCT1. Drug Metab. Pharmacokinet. 20, 379–386. https://doi.org/10.2133/dmpk.20.379 (2005).
CAS Article PubMed Google Scholar
Filipski, K. K., Loos, W. J., Verweij, J. & Sparreboom, A. Interaction of cisplatin with the human organic cation transporter 2. Clin. Cancer Res. 14, 3875–3880. https://doi.org/10.1158/1078-0432.CCR-07-4793 (2008).
CAS Article PubMed Google Scholar
Ciarimboli, G. et al. Organic cation transporter 2 mediates cisplatin-induced oto- and nephrotoxicity and is a target for protective interventions. Am. J. Pathol. 176, 1169–1180. https://doi.org/10.2353/ajpath.2010.090610 (2010).
CAS Article PubMed PubMed Central Google Scholar
Kim, S. et al. New data content and improved web interfaces. Nucleic Acids Res. 49(2021), D1388–D1395. https://doi.org/10.1093/nar/gkaa971 (2021).
CAS Article PubMed Google Scholar
Wang, Z. & Moult, J. SNPs, protein structure, and disease. Hum. Mutat. 17, 263–270. https://doi.org/10.1002/humu.22 (2001).
Article PubMed Google Scholar
Crawford, D. C. & Nickerson, D. A. Definition and clinical importance of haplotypes. Annu. Rev. Med. 56, 303–320. https://doi.org/10.1146/annurev.med.56.082103.104540 (2005).
CAS Article PubMed Google Scholar
Sajib, A. A., Islam, T., Paul, N. & Yeasmin, S. Interaction of rs316019 variants of SLC22A2 with metformin and other drugs- an in silico analysis. J. Genet. Eng. Biotechnol. 16, 769–775. https://doi.org/10.1016/j.jgeb.2018.01.003 (2018).
Article PubMed PubMed Central Google Scholar
Jacobs, C., Pearce, B., Du Plessis, M., Hoosain, N. & Benjeddou, M. Single nucleotide polymorphisms of the SLC22A2 gene within the Xhosa population of South Africa. Drug Metab. Pharmacokinet. 30, 457–460. https://doi.org/10.1016/j.dmpk.2015.11.002 (2015).
CAS Article PubMed Google Scholar
Ng, P. C. & Henikoff, S. Predicting deleterious amino acid substitutions. Genome Res. 11, 863–874. https://doi.org/10.1101/gr.176601 (2001).
CAS Article PubMed PubMed Central Google Scholar
Ng, P. C. & Henikoff, S. SIFT: Predicting amino acid changes that affect protein function. Nucleic Acids Res. 31, 3812–3814. https://doi.org/10.1093/nar/gkg509 (2003).
CAS Article PubMed PubMed Central Google Scholar
Vaser, R., Adusumalli, S., Leng, S. N., Sikic, M. & Ng, P. C. SIFT missense predictions for genomes. Nat. Protoc. 11, 1–9. https://doi.org/10.1038/nprot.2015.123 (2016).
CAS Article PubMed Google Scholar
Sunyaev, S. et al. Prediction of deleterious human alleles. Hum. Mol. Genet. 10, 591–597. https://doi.org/10.1093/hmg/10.6.591 (2001).
CAS Article PubMed Google Scholar
Adzhubei, I. A. et al. A method and server for predicting damaging missense mutations. Nat. Methods. 7, 248–249. https://doi.org/10.1038/nmeth0410-248 (2010).
CAS Article PubMed PubMed Central Google Scholar
Adzhubei, I., Jordan, D. M. & Sunyaev, S. R. Predicting functional effect of human missense mutations using PolyPhen-2. Curr. Protoc. Hum. Genet. https://doi.org/10.1002/0471142905.hg0720s76 (2013).
Article PubMed PubMed Central Google Scholar
Capriotti, E., Fariselli, P. & Casadio, R. I-Mutant2.0: Predicting stability changes upon mutation from the protein sequence or structure. Nucleic Acids Res. 33, 306–310. https://doi.org/10.1093/nar/gki375 (2005).
CAS Article Google Scholar
Capriotti, E., Fariselli, P., Rossi, I. & Casadio, R. A three-state prediction of single point mutations on protein stability changes. BMC Bioinform. 9, S6. https://doi.org/10.1186/1471-2105-9-S2-S6 (2008).
CAS Article Google Scholar
Jumper, J. et al. Highly accurate protein structure prediction with AlphaFold. Nature 596, 583–589. https://doi.org/10.1038/s41586-021-03819-2 (2021).
ADS CAS Article PubMed PubMed Central Google Scholar
Stierand, K. & Rarey, M. PoseView—Molecular interaction patterns at a glance. J. Cheminform. 2, 20146. https://doi.org/10.1186/1758-2946-2-s1-p50 (2010).
Article Google Scholar
Studer, G., Biasini, M. & Schwede, T. Assessing the local structural quality of transmembrane protein models using statistical potentials (QMEANBrane). Bioinformatics 30, 505–511. https://doi.org/10.1093/bioinformatics/btu457 (2014).
CAS Article Google Scholar
Colovos, C. & Yeates, T. O. Verification of protein structures: Patterns of nonbonded atomic interactions. Protein Sci. 2, 1511–1519. https://doi.org/10.1002/pro.5560020916 (1993).
CAS Article PubMed PubMed Central Google Scholar
Laskowski, R. A., MacArthur, M. W., Moss, D. S. & Thornton, J. M. PROCHECK: A program to check the stereochemical quality of protein structures. J. Appl. Crystallogr. 26, 283–291. https://doi.org/10.1107/s0021889892009944 (1993).
CAS Article Google Scholar
Laskowski, R. A., Rullmann, J. A. C., MacArthur, M. W., Kaptein, R. & Thornton, J. M. AQUA and PROCHECK-NMR: Programs for checking the quality of protein structures solved by NMR. J. Biomol. NMR. 8, 477–486. https://doi.org/10.1007/BF00228148 (1996).
CAS Article PubMed Google Scholar
Schrödinger, L., & D.W.P. [Internet]. http://www.pymol.org/pymo. (2020).
Lee, J. et al. CHARMM-GUI input generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM simulations using the CHARMM36 additive force field. J. Chem. Theory Comput. 12, 405–413. https://doi.org/10.1021/acs.jctc.5b00935 (2016).
CAS Article PubMed Google Scholar
Abraham, M. J. et al. GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX. 1–2, 19–25. https://doi.org/10.1016/J.SOFTX.2015.06.001 (2015).
ADS Article Google Scholar
Huang, J. et al. CHARMM36m: An improved force field for folded and intrinsically disordered proteins. Nat. Methods. 14, 71–73. https://doi.org/10.1038/nmeth.4067 (2017).
CAS Article PubMed Google Scholar
Vanommeslaeghe, K. et al. CHARMM general force field: A force field for drug-like molecules compatible with the CHARMM all-atom additive biological force fields. J. Comput. Chem. https://doi.org/10.1002/jcc.21367 (2009).
Article Google Scholar
Parrinello, M. & Rahman, A. Polymorphic transitions in single crystals: A new molecular dynamics method. J. Appl. Phys. 52, 7182–7190. https://doi.org/10.1063/1.328693 (1981).
ADS CAS Article Google Scholar
Essmann, U. et al. A smooth particle mesh Ewald method. J. Chem. Phys. 103, 8577–8593. https://doi.org/10.1063/1.470117 (1995).
ADS CAS Article Google Scholar
Kumari, R., Kumar, R. & Lynn, A. G-mmpbsa—A GROMACS tool for high-throughput MM-PBSA calculations. J. Chem. Inf. Model. 54, 1951–1962. https://doi.org/10.1021/ci500020m (2014).
CAS Article PubMed Google Scholar
Gapsys, V., de Groot, B. L. & Briones, R. Computational analysis of local membrane properties. J. Comput. Aided. Mol. Des. 27, 845–858. https://doi.org/10.1007/s10822-013-9684-0 (2013).
ADS CAS Article PubMed PubMed Central Google Scholar
Surendiran, A. et al. Influence of CYP2C9 gene polymorphisms on response to glibenclamide in type 2 diabetes mellitus patients. Eur. J. Clin. Pharmacol. 67, 797–801. https://doi.org/10.1007/s00228-011-1013-8 (2011).
CAS Article PubMed Google Scholar
Banu, H., Renuka, N. & Vasanthakumar, G. Reduced catalytic activity of human CYP2C9 natural alleles for gliclazide: Molecular dynamics simulation and docking studies. Biochimie 93, 1028–1036. https://doi.org/10.1016/j.biochi.2011.02.008 (2011).
CAS Article PubMed Google Scholar
Trott, O. & Olson, A. J. AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J. Comput. Chem. https://doi.org/10.1002/jcc.21334 (2009).
Article Google Scholar
Morris, G. M. et al. AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility. J. Comput. Chem. https://doi.org/10.1002/jcc.21256 (2009).
Article PubMed PubMed Central Google Scholar
Delano, W. The PyMOL Molecular Graphics System. (2002).
Hernández-santoyo, A., Tenorio-barajas, A. Y., Altuzar, V., Vivanco-cid, H. & Mendoza-barrera, C. Protein–protein and protein–ligand docking. Protein Eng. Technol. Appl. 29, 63–81 (2013).
Google Scholar
Tosi, F. et al. Combination treatment with metformin and glibenclamide versus single-drug therapies in type 2 diabetes mellitus: A randomized, double-blind, comparative study. Metabolism 52, 862–867. https://doi.org/10.1016/S0026-0495(03)00101-X (2003).
CAS Article PubMed Google Scholar
Marre, M., Howlett, H., Lehert, P. & Allavoine, T. Improved glycaemic control with metformin-glibenclamide combined tablet therapy (glucovance in type 2 diabetic patients inadequately controlled on metformin). Diabet. Med. 19, 673–680. https://doi.org/10.1046/j.1464-5491.2002.00774.x (2002).
CAS Article PubMed Google Scholar
Vilar, L. et al. Comparison of metformin, gliclazide MR and rosiglitazone in monotherapy and in combination for type 2 diabetes. Arq. Bras. Endocrinol. Metabol. 54, 311–318. https://doi.org/10.1590/s0004-27302010000300010 (2010).
Article PubMed Google Scholar
Zazuli, Z. et al. The impact of genetic polymorphisms in organic cation transporters on renal drug disposition. Int. J. Mol. Sci. 21, 1–23. https://doi.org/10.3390/ijms21186627 (2020).
CAS Article Google Scholar
Yoon, H., Cho, H. Y., Yoo, H. D., Kim, S. M. & Lee, Y. B. Influences of organic cation transporter polymorphisms on the population pharmacokinetics of metformin in healthy subjects. AAPS J. 15, 571–580. https://doi.org/10.1208/s12248-013-9460-z (2013).
CAS Article PubMed PubMed Central Google Scholar
Wang, Z. J., Yin, O. Q. P., Tomlinson, B. & Chow, M. S. S. OCT2 polymorphisms and in-vivo renal functional consequence: Studies with metformin and cimetidine. Pharmacogenet. Genomics 18, 637–645. https://doi.org/10.1097/FPC.0b013e328302cd41 (2008).
CAS Article PubMed Google Scholar
Song, I. S. et al. Genetic variants of the organic cation transporter 2 influence the disposition of metformin. Clin. Pharmacol. Ther. 84, 559–562. https://doi.org/10.1038/clpt.2008.61 (2008).
CAS Article PubMed Google Scholar
Chen, Y. et al. Effect of genetic variation in the organic cation transporter 2 on the renal elimination of metformin. Pharmacogenet. Genomics 19, 497–504. https://doi.org/10.1097/FPC.0b013e32832cc7e9 (2009).
CAS Article PubMed PubMed Central Google Scholar
Dujic, T. et al. Variants in pharmacokinetic transporters and glycemic response to metformin: A metgen meta-analysis. Clin. Pharmacol. Ther. 101, 763–772. https://doi.org/10.1002/cpt.567 (2017).
CAS Article PubMed PubMed Central Google Scholar
Christensen, M. M. H. et al. A gene-gene interaction between polymorphisms in the OCT2 and MATE1 genes influences the renal clearance of metformin. Pharmacogenet. Genomics 23, 526–534. https://doi.org/10.1097/FPC.0b013e328364a57d (2013).
CAS Article PubMed Google Scholar
Tzvetkov, M. V. et al. The effects of genetic polymorphisms in the organic cation transporters OCT1, OCT2, and OCT3 on the renal clearance of metformin. Clin. Pharmacol. Ther. 86, 299–306. https://doi.org/10.1038/clpt.2009.92 (2009).
CAS Article PubMed Google Scholar
Defronzo, R., Fleming, G. A., Chen, K. & Bicsak, T. A. Metformin-associated lactic acidosis: Current perspectives on causes and risk. Metabolism 65, 20–29. https://doi.org/10.1016/j.metabol.2015.10.014 (2016).
CAS Article PubMed Google Scholar
Piel, S., Ehinger, J. K., Elmér, E. & Hansson, M. J. Metformin induces lactate production in peripheral blood mononuclear cells and platelets through specific mitochondrial complex I inhibition. Acta Physiol. 213, 171–180. https://doi.org/10.1111/apha.12311 (2015).
CAS Article Google Scholar
Huang, W., Castelino, R. L. & Peterson, G. M. Adverse event notifications implicating metformin with lactic acidosis in Australia. J. Diabetes Complic. 29, 1261–1265. https://doi.org/10.1016/j.jdiacomp.2015.06.001 (2015).
Article Google Scholar
Grün, B. et al. Trimethoprim-metformin interaction and its genetic modulation by OCT2 and MATE1 transporters. Br. J. Clin. Pharmacol. 76, 787–796. https://doi.org/10.1111/bcp.12079 (2013).
CAS Article PubMed PubMed Central Google Scholar
Download references
The work reported herein was made possible through funding by the South African Medical Research Council through its Division of Research Capacity Development under funding received from the South African National Treasury. The content hereof is the sole responsibility of the authors and do not necessarily represent the official views of the SAMRC or the funders. Appreciation goes to Dr Ruben Cloete from the South African National Bioinformatics Institute (SANBI), University of the Western Cape, South Africa for performing the structural modelling and molecular dynamic simulations.
Precision Medicine Unit, Department of Biotechnology, Faculty of Natural Sciences, University of the Western Cape, Robert Sobukwe Road, Bellville, 7535, South Africa
Zainonesa Abrahams-October & Mongi Benjeddou
Biomedical Research and Innovation Platform, South African Medical Research Council, Tygerberg, Cape Town, 7505, South Africa
Rabia Johnson
Division of Medical Physiology, Faculty of Medicine and Health Sciences, Stellenbosch University, Tygerberg, 7505, South Africa
Rabia Johnson
South African Medical Research Council Bioinformatics Unit, South African National Bioinformatics Institute, University of the Western Cape, Private Bag X17, Bellville, Cape Town, 7535, South Africa
Ruben Cloete
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
Z.A.O. and R.C. conceptualized the project, wrote the main manuscript, performed analysis and prepared all figures. R.J., M.B. and R.C. provided supervision. All authors reviewed the manuscript.
Correspondence to Ruben Cloete.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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 http://creativecommons.org/licenses/by/4.0/.
Reprints and Permissions
Abrahams-October, Z., Johnson, R., Benjeddou, M. et al. The determination of the effect(s) of solute carrier family 22-member 2 (SLC22A2) haplotype variants on drug binding via molecular dynamic simulation systems. Sci Rep 12, 16936 (2022). https://doi.org/10.1038/s41598-022-21291-4
Download citation
Received: 02 February 2022
Accepted: 26 September 2022
Published: 08 October 2022
DOI: https://doi.org/10.1038/s41598-022-21291-4
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Advertisement
© 2022 Springer Nature Limited
Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.