Proteomics analysis of the p.G849D variant in neurexin 2 alpha may reveal insight into Parkinson's disease pathobiology.

Parkinson's disease (PD), the fastest-growing neurological disorder globally, has a complex etiology. A previous study by our group identified the p.G849D variant in neurexin 2 (NRXN2), encoding the synaptic protein, NRXN2α, as a possible causal variant of PD. Therefore, we aimed to perform functional studies using proteomics in an attempt to understand the biological pathways affected by the variant. We hypothesized that this may reveal insight into the pathobiology of PD. Wild-type and mutant NRXN2α plasmids were transfected into SH-SY5Y cells. Thereafter, total protein was extracted and prepared for mass spectrometry using a Thermo Scientific Fusion mass spectrometer equipped with a Nanospray Flex ionization source. The data were then interrogated against the UniProt H. sapiens database and afterward, pathway and enrichment analyses were performed using in silico tools. Overexpression of the wild-type protein led to the enrichment of proteins involved in neurodegenerative diseases, while overexpression of the mutant protein led to the decline of proteins involved in ribosomal functioning. Thus, we concluded that the wild-type NRXN2α may be involved in pathways related to the development of neurodegenerative disorders, and that biological processes related to the ribosome, transcription, and tRNA, specifically at the synapse, could be an important mechanism in PD. Future studies targeting translation at the synapse in PD could therefore provide further information on the pathobiology of the disease.

Neurexin 2 p.G849D variant, implicated in Parkinson's disease, increases reactive oxygen species, and reduces cell viability and mitochondrial membrane potential in SH-SY5Y cells.

Parkinson's disease (PD) is a neurodegenerative movement disorder, affecting 1-2% of the human population over 65. A previous study by our group identified a p.G849D variant in neurexin 2α (NRXN2) co-segregating with PD, prompting validation of its role using experimental methods. This novel variant had been found in a South African family with autosomal dominant PD. NRXN2α is an essential synaptic maintenance protein with multiple functional roles at the synaptic cleft. The aim of the present study was to investigate the potential role of the translated protein NRXN2α and the observed mutant in PD by performing functional studies in an in vitro model. Wild-type and mutant NRXN2α plasmids were transfected into SH-SY5Y cells to assess the effect of the mutant on cell viability and apoptosis [(3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) Assay; ApoTox-Glo™ Triplex Assay)], mitochondrial membrane potential (MMP; MitoProbe™ JC-1 Assay), mitochondrial network analysis (MitoTracker®) and reactive oxygen species (ROS; ROS-Glo™ H2O2 Assay). Cells transfected with the mutant NRXN2α plasmid showed decreased cell viability and MMP. They also exhibited increased ROS production. However, these cells showed no changes in mitochondrial fragmentation. Our findings led us to speculate that the p.G849D variant may be involved in a toxic feedback loop leading to neuronal death in PD. Mitochondrial dysfunction and synaptic dysfunction have been linked to PD. Therefore, findings from this exploratory study are in line with previous studies connecting these two processes and warrants further investigation into the role of this variant in other cellular and animal models.

The determination of the effect(s) of solute carrier family 22-member 2 (SLC22A2) haplotype variants on drug binding via molecular dynamic simulation systems.

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 Flexible, Extended Coil of the PDZ-Binding Motif of the Three Deadly Human Coronavirus E Proteins Plays a Role in Pathogenicity.

The less virulent human (h) coronaviruses (CoVs) 229E, NL63, OC43, and HKU1 cause mild, self-limiting respiratory tract infections, while the more virulent SARS-CoV-1, MERS-CoV, and SARS-CoV-2 have caused severe outbreaks. The CoV envelope (E) protein, an important contributor to the pathogenesis of severe hCoV infections, may provide insight into this disparate severity of the disease. We, therefore, generated full-length E protein models for SARS-CoV-1 and -2, MERS-CoV, HCoV-229E, and HCoV-NL63 and docked C-terminal peptides of each model to the PDZ domain of the human PALS1 protein. The PDZ-binding motif (PBM) of the SARS-CoV-1 and -2 and MERS-CoV models adopted a more flexible, extended coil, while the HCoV-229E and HCoV-NL63 models adopted a less flexible alpha helix. All the E peptides docked to PALS1 occupied the same binding site and the more virulent hCoV E peptides generally interacted more stably with PALS1 than the less virulent ones. We hypothesize that the increased flexibility of the PBM in the more virulent hCoVs facilitates more stable binding to various host proteins, thereby contributing to more severe disease. This is the first paper to model full-length 3D structures for both the more virulent and less virulent hCoV E proteins, providing novel insights for possible drug and/or vaccine development.

Molecular Modeling of Subtype-Specific Tat Protein Signatures to Predict Tat-TAR Interactions That May Be Involved in HIV-Associated Neurocognitive Disorders.

HIV-1 is responsible for a spectrum of neurocognitive deficits defined as HIV-associated neurocognitive disorders (HAND). The HIV transactivator of transcription (Tat) protein plays a key role in the neuropathophysiology of HAND. The Tat protein functions by transactivation of viral genes through its interaction with the transactivation response (TAR) RNA element. Subtype-specific Tat protein signatures including C31S, R57S and Q63E present in Tat subtype C has previously been linked to a lowered neuropathophysiology compared to Tat subtype B. In this study, we attempted to understand the molecular mechanism by which Tat subtype-specific variation, particularly, C31S, R57S, and Q63E influence the Tat-TAR interaction. We performed molecular modeling to generate accurate three-dimensional protein structures of the HIV-1 Tat subtypes C and B using the Swiss model webserver. Thereafter, we performed a molecular docking of the TAR RNA element to each of the Tat subtypes B and C protein structures using the HDOCK webserver. Our findings indicate that Tat subtype B had a higher affinity for the TAR RNA element compared to Tat subtype C based on a higher docking score of -187.37, a higher binding free energy value of -9834.63 ± 216.17 kJ/mol, and a higher number of protein-nucleotide interactions of 26. Furthermore, Tat subtype B displayed more flexible regions when bound to the TAR element and this flexibility could account for the stronger affinity of Tat subtype B to TAR. From the Tat signatures linked to neuropathogenesis, only R57/R57S are involved in Tat-TAR interaction. Due to the lack of electrostatic interactions observed between Tat subtype C and TAR, weaker affinity is observed, and this may contribute to a lower level of neuropathophysiology observed in subtype C infection.

In silico repurposing of a Novobiocin derivative for activity against latency associated Mycobacterium tuberculosis drug target nicotinate-nucleotide adenylyl transferase (Rv2421c).

Nicotinamide-nucleotide adenylyl transferase (Rv2421c) was selected as a potential drug target, because it has been shown, in vitro, to be essential for Mycobacterium tuberculosis growth. It is conserved between mycobacterium species, is up-regulated during dormancy, has a known 3D crystal structure and has no known human homologs. A model of Rv2421c in complex with nicotinic acid adenine dinucleotide and magnesium ion was constructed and subject tovirtual ligand screening against the Prestwick Chemical Library and the ZINC database, which yielded 155 potential hit molecules. Of the 155 compounds identified five were pursued further using an IC50 based 3D-QSAR study. The 3D-QSAR model validated the inhibition properties of the five compounds based on R2 value of 0.895 and Q2 value of 0.944 compared to known inhibitors of Rv2421c. Higher binding affinities was observed for the novel ZINC13544129 and two FDA approved compounds (Novobiocin sodium salt, Sulfasalazine). Similarly, the total interaction energy was found to be the highest for Cromolyn disodium system (-418.88 kJ/mol) followed by Novobiocin (-379.19 kJ/mol) and Sulfasalazine with (-330.13 kJ/mol) compared to substrate DND having (-185.52 kJ/mol). Subsequent in vitro testing of the five compounds identified Novobiocin sodium salt with activity against Mycobacterium tuberculosis at 50 µM, 25µM and weakly at 10µM concentrations. Novobiocin salt interacts with a MG ion and active site residues His20, Thr86, Gly107 and Leu164 similar to substrate DND of Mycobacterium tuberculosis Rv2421c. Additional in silico structural analysis of known Novobiocin sodium salt derivatives against Rv2421c suggest Coumermycin as a promising alternative for the treatment of Mycobacterium tuberculosis based on large number of hydrogen bond interactions with Rv2421c similar in comparison to Novobiocin salt and substrate DND.

Emerging evidence implicating a role for neurexins in neurodegenerative and neuropsychiatric disorders.

Synaptopathies are brain disorders characterized by dysfunctional synapses, which are specialized junctions between neurons that are essential for the transmission of information. Synaptic dysfunction can occur due to mutations that alter the structure and function of synaptic components or abnormal expression levels of a synaptic protein. One class of synaptic proteins that are essential to their biology are cell adhesion proteins that connect the pre- and post-synaptic compartments. Neurexins are one type of synaptic cell adhesion molecule that have, recently, gained more pathological interest. Variants in both neurexins and their common binding partners, neuroligins, have been associated with several neuropsychiatric disorders. In this review, we summarize some of the key physiological functions of the neurexin protein family and the protein networks they are involved in. Furthermore, examination of published literature has implicated neurexins in both neuropsychiatric and neurodegenerative disorders. There is a clear link between neurexins and neuropsychiatric disorders, such as autism spectrum disorder and schizophrenia. However, multiple expression studies have also shown changes in neurexin expression in several neurodegenerative disorders, including Alzheimer's disease and Parkinson's disease. Therefore, this review highlights the potential importance of neurexins in brain disorders and the importance of doing more targeted studies on these genes and proteins.

Aptamer-Based Diagnostic Systems for the Rapid Screening of TB at the Point-of-Care.

The transmission of Tuberculosis (TB) is very rapid and the burden it places on health care systems is felt globally. The effective management and prevention of this disease requires that it is detected early. Current TB diagnostic approaches, such as the culture, sputum smear, skin tuberculin, and molecular tests are time-consuming, and some are unaffordable for low-income countries. Rapid tests for disease biomarker detection are mostly based on immunological assays that use antibodies which are costly to produce, have low sensitivity and stability. Aptamers can replace antibodies in these diagnostic tests for the development of new rapid tests that are more cost effective; more stable at high temperatures and therefore have a better shelf life; do not have batch-to-batch variations, and thus more consistently bind to a specific target with similar or higher specificity and selectivity and are therefore more reliable. Advancements in TB research, in particular the application of proteomics to identify TB specific biomarkers, led to the identification of a number of biomarker proteins, that can be used to develop aptamer-based diagnostic assays able to screen individuals at the point-of-care (POC) more efficiently in resource-limited settings.

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

BACKGROUND: 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. METHODS: 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. RESULTS: 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. CONCLUSIONS: 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.

Prioritization of candidate genes for a South African family with Parkinson's disease using in-silico tools.

Parkinson's disease (PD) is a neurodegenerative disorder exhibiting Mendelian inheritance in some families. Next-generation sequencing approaches, including whole exome sequencing (WES), have revolutionized the field of Mendelian disorders and have identified a number of PD genes. We recruited a South African family with autosomal dominant PD and used WES to identify a possible pathogenic mutation. After filtration and prioritization, we found five potential causative variants in CFAP65, RTF1, NRXN2, TEP1 and CCNF. The variant in NRXN2 was selected for further analysis based on consistent prediction of deleteriousness across computational tools, not being present in unaffected family members, ethnic-matched controls or public databases, and its expression in the substantia nigra. A protein model for NRNX2 was created which provided a three-dimensional (3D) structure that satisfied qualitative mean and global model quality assessment scores. Trajectory analysis showed destabilizing effects of the variant on protein structure, indicated by high flexibility of the LNS-6 domain adopting an extended conformation. We also found that the known substrate N-acetyl-D-glucosamine (NAG) contributed to restoration of the structural stability of mutant NRXN2. If NRXN2 is indeed found to be the causal gene, this could reveal a new mechanism for the pathobiology of PD.

Short Communication: HIV-1 Drug Resistance Mutation Analyses of Cameroon-Derived Integrase Sequences.

HIV-1 integrase (IN) is a primary target for combination antiretroviral therapy. Only a limited number of studies report on the emergence of resistance-associated mutations (RAMs) in Cameroon. We observed that 1.4% of sequence from treatment-naive patients had IN strand transfer inhibitor (INSTI) RAMs. These mutations confer resistance to raltegravir and elvitegravir. We also observed that 10.1% of the sequences have INSTI accessory RAMs. HIV-1 CRF02_AG was the predominant subtype (44.7%) in this study analyses. The occurrence of INSTI RAMs among the sequences at baseline needs to be monitored carefully.

Structural Comparison of Diverse HIV-1 Subtypes using Molecular Modelling and Docking Analyses of Integrase Inhibitors.

The process of viral integration into the host genome is an essential step of the HIV-1 life cycle. The viral integrase (IN) enzyme catalyzes integration. IN is an ideal therapeutic enzyme targeted by several drugs; raltegravir (RAL), elvitegravir (EVG), dolutegravir (DTG), and bictegravir (BIC) having been approved by the USA Food and Drug Administration (FDA). Due to high HIV-1 diversity, it is not well understood how specific naturally occurring polymorphisms (NOPs) in IN may affect the structure/function and binding affinity of integrase strand transfer inhibitors (INSTIs). We applied computational methods of molecular modelling and docking to analyze the effect of NOPs on the full-length IN structure and INSTI binding. We identified 13 NOPs within the Cameroonian-derived CRF02_AG IN sequences and further identified 17 NOPs within HIV-1C South African sequences. The NOPs in the IN structures did not show any differences in INSTI binding affinity. However, linear regression analysis revealed a positive correlation between the Ki and EC50 values for DTG and BIC as strong inhibitors of HIV-1 IN subtypes. All INSTIs are clinically effective against diverse HIV-1 strains from INSTI treatment-naïve populations. This study supports the use of second-generation INSTIs such as DTG and BIC as part of first-line combination antiretroviral therapy (cART) regimens, due to a stronger genetic barrier to the emergence of drug resistance.

Genome Sequencing of a Severe Acute Respiratory Syndrome Coronavirus 2 Isolate Obtained from a South African Patient with Coronavirus Disease 2019.

As a contribution to the global efforts to track and trace the ongoing coronavirus pandemic, here we present the sequence, phylogenetic analysis, and modeling of nonsynonymous mutations for a severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome that was detected in a South African patient with coronavirus disease 2019 (COVID-19).

Correction: Molecular dynamic simulations to investigate the structural impact of known drug resistance mutations on HIV-1C Integrase-Dolutegravir binding.

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

Molecular dynamic simulations to investigate the structural impact of known drug resistance mutations on HIV-1C Integrase-Dolutegravir binding.

Resistance associated mutations (RAMs) threaten the long-term success of combination antiretroviral therapy (cART) outcomes for HIV-1 treatment. HIV-1 Integrase (IN) strand transfer inhibitors (INSTIs) have proven to be a viable option for highly specific HIV-1 therapy. The INSTI, Dolutegravir is recommended by the World Health Organization for use as first-line cART. This study aims to understand how RAMs affect the stability of IN, as well as the binding of the drug Dolutegravir to the catalytic pocket of the protein. A homology model of HIV-1 subtype C IN was successfully constructed and validated. The site directed mutator webserver was used to predict destabilizing and/or stabilizing effects of known RAMs while FoldX confirmed any changes in protein energy upon introduction of mutation. Also, interaction analysis was performed between neighbouring residues. Three mutations known to be associated with Raltegravir, Elvitegravir and Dolutegravir resistance were selected; E92Q, G140S and Y143R, for molecular dynamics simulations. The structural quality assessment indicated high reliability of the HIV-1C IN tetrameric structure, with more than 90% confidence in modelled regions. Change in free energy for the three mutants indicated different effects, while simulation analysis showed G140S to have the largest affect on protein stability and flexibility. This was further supported by weaker non-bonded pairwise interaction energy and binding free energy values between the drug DTG and E92Q, Y143R and G140S mutants suggesting reduced binding affinity, as indicated by interaction analysis in comparison to the WT. Our findings suggest the G140S mutant has the strongest effect on the HIV-1C IN protein structure and Dolutegravir binding. To the best of our knowledge, this is the first study that uses the consensus wild type HIV-1C IN sequence to build an accurate 3D model to understand the effect of three known mutations on DTG drug binding in a South Africa context.

Drug Resistance Mutations Against Protease, Reverse Transcriptase and Integrase Inhibitors in People Living With HIV-1 Receiving Boosted Protease Inhibitors in South Africa.

The South African national combination antiretroviral therapy (cART) roll-out program started in 2006, with over 4.4 million people accessing treatment since it was first introduced. HIV-1 drug resistance can hamper the success of cART. This study determined the patterns of HIV-1 drug-resistance associated mutations (RAMs) in People Living with HIV-1 (PLHIV-1). Receiving first (for children below 3 years of age) and second-line (for adults) cART regimens in South Africa. During 2017 and 2018, 110 patients plasma samples were selected, 96 samples including those of 17 children and infants were successfully analyzed. All patients were receiving a boosted protease inhibitor (bPI) as part of their cART regimen. The viral sequences were analyzed for RAMs through genotypic resistance testing. We performed genotypic resistance testing (GRT) for Protease inhibitors (PIs), Reverse transcriptase inhibitors (RTIs) and Integrase strand transfer inhibitors (InSTIs). Viral sequences were subtyped using REGAv3 and COMET. Based on the PR/RT sequences, HIV-1 subtypes were classified as 95 (99%) HIV-1 subtype C (HIV-1C) while one sample as 02_AG. Integrase sequencing was successful for 89 sequences, and all the sequences were classified as HIV-1C (99%, 88/89) except one sequence classified CRF02_AG, as observed in PR/RT. Of the 96 PR/RT sequences analyzed, M184V/I (52/96; 54%) had the most frequent RAM nucleoside reverse transcriptase inhibitor (NRTI). The most frequent non-nucleoside reverse transcriptase inhibitor (NNRTI) RAM was K103N/S (40/96, 42%). Protease inhibitor (PI) RAMs M46I and V82A were present in 12 (13%) of the sequences analyzed. Among the InSTI major RAM two (2.2%) sequences have Y143R and T97A mutations while one sample had T66I. The accessory RAM E157Q was identified in two (2.2%). The data indicates that the majority of the patients failed on bPIs didn't have any mutation; therefore adherence could be major issue in these groups of individuals. We propose continued viral load monitoring for better management of infected PLHIV.

HIV-1 Integrase Diversity and Resistance-Associated Mutations and Polymorphisms Among Integrase Strand Transfer Inhibitor-Naive HIV-1 Patients from Cameroon.

The World Health Organization (WHO) has put forth recommendations for the use of integrase (IN) strand transfer inhibitors (INSTIs) to be part of the first-line combination antiretroviral therapy regimen to treat HIV infections. The knowledge of pretreatment drug resistance against INSTIs is still scarce in resource-limited settings (RLS). We characterized the integrase gene to identify resistance-associated mutations (RAMs) in 56 INSTI-naive patient viral sequences from Cameroon. Study analysis used 37 sequences with fragment size ≥500 bp or of good quality .The majority of the sequences were identified as CRF02_AG 54.% (n = 20/37) and 45.9% (n = 17/37), other subtype viral sequences include (A, CRF36_cpx, F,G, and C). A total of 18.9% (n = 7/37) of the sequences had RAMs, with only 5.4% (n = 2/37) having major RAMs (Y143R/C/D/G and P145S), against INSTIs. Accessory RAMs were present in 8.1% (n = 3/37) of the sequences, of which one sequence contained solely E157Q, and another Q95K. One patient sequence had three accessory RAMs (G140E, E157Q, and G163R). We identified major RAMs to INSTIs, which might have a potential clinical impact to dolutegravir rollout in RLS, including Cameroon. This is the first study to describe RAMs among INSTI-naive people living with HIV-1 (PLHIV-1) infected with CRF02_AG and other subtypes in Cameroon.

Structure based identification of novel inhibitors against ATP synthase of Mycobacterium tuberculosis: A combined in silico and in vitro study.

The shortcomings of conventional tuberculosis treatments resulting from the development of drug resistance in Mycobacterium tuberculosis drive a need for the formulation of novel therapeutic agents. The diarylquinoline class of drugs such as bedaquiline was recently approved for the treatment of multidrug-resistant strains of tuberculosis, primarily targeting c and ε subunits of the ATP synthases. Yet resistance to bedaquiline has already been reported. Therefore, Rv1311 was used as the target for the identification of possible inhibitors against the M. tuberculosis. The structure of Rv1311 was predicted and common feature pharmacophore models were generated which facilitated the identification of potential inhibitors in the ZINC database. The activities of the selected molecules were compared with known inhibitors of the ATP synthase using quantitative structure-activity relationship. The ZINC classified inhibitors showed comparable predicted activities with that of known inhibitors. Furthermore, the inhibitory behavior of the studied drug molecules was experimentally determined using in vitro techniques and showed the minimum inhibitory concentration as low as 25 µM. The resulted outcomes provide a deeper insight into the structural basis of Rv1311 inhibitions and can facilitate the process of drug design against tuberculosis.

Molecular modelling and simulation studies of the Mycobacterium tuberculosis multidrug efflux pump protein Rv1258c.

Mycobacterial efflux pumps play a major role in the emergence of antimycobacterial drug resistance. Of particular interest is the proteinaceous multi-drug efflux pump protein Rv1258c that encodes a tetracycline/aminoglycoside resistance (TAP-2)-like efflux pump which is active in susceptible and drug resistant Mycobacterium tuberculosis. Rv1258c is implicated in drug resistance to numerous antimycobacterials including first line drugs rifampicin and isoniazid as well as fluoroquinolone and aminoglycoside antibiotic classes. To date, compounds like verapamil and piperine have been shown to inhibit Rv1258c but no direct evidence for binding or mode of action exist. Therefore in the present study we generated an accurate 3D model of Rv1258c using MODELLER and validated its structure using molecular dynamic simulation studies with GROMACS software. The 3D-structures of Rv1258c and the homologous template 1pw4 were simulated within a POPE/POPG lipid bilayer and found to behave similar. Another important finding was the identification of one local energy minima state of the apo protein, which speaks to the flexibility of the protein and will be investigated further. Extraction of one of the open channel conformations of Rv1258c and blind docking of various structurally diverse putative inhibitors and substrates, allowed for the identification of a probable binding site. Spectinamide was found to bind to a different location on the outside surface of the protein suggesting its ability to avoid the efflux channel. We further identified 246 putative compounds that showed higher binding affinity values to Rv1258c compared to piperine and verapamil. Interaction analysis of the top 20 purchasable compounds identified crucial hydrogen bond interactions with Ser26, Ser45 and Glu243 as well as a pi-pi stacking interaction with Trp32 that accounted for the strong affinity of these compounds for Rv1258c. Future studies will entail purchasing a number of compounds for in vitro activity testing against Mycobacterium tuberculosis.

Author Correction: Analyses of HIV-1 integrase sequences prior to South African national HIV-treatment program and availability of integrase inhibitors in Cape Town, South Africa.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has been fixed in the paper.