结核杆菌的药物外排泵

外排泵基因过度表达促使细胞内的药物转移到细胞外,胞内药物浓度下降,致使药物临床疗效降低或无效。外排泵上有多个与识别底物结构无关的阴离子凹穴(底物结合位点),通过静电引力非特异性结合各种含阳离子的药物,在能量协助下将药物排到细胞外,进而导致细菌抗药性的产生。研究表明,结核杆菌基因组中存在外排泵基因,但仅分离到几个有功效的外排泵。本文就相关研究进展进行综述。     

Plasmid-based controls to detect rpoB mutations in Mycobacterium tuberculosis by quantitative polymerase chain reaction-high-resolution melting

Quantitative polymerase chain reaction-high-resolution melting (qPCR-HRM) analysis was used to screen for mutations related to drug resistance in Mycobacterium tuberculosis. We detected the C526T and C531T mutations in the rifampicin resistance-determining region (RRDR) of the rpoB gene with qPCR-HRM using plasmid-based controls. A segment of the RRDR region from M. tuberculosis H37Rv and from strains carrying C531T or C526T mutations in the rpoB were cloned into pGEM-T vector and these vectors were used as controls in the qPCR-HRM analysis of 54 M. tuberculosis strains. The results were confirmed by DNA sequencing and showed that recombinant plasmids can replace genomic DNA as controls in the qPCR-HRM assay. Plasmids can be handled outside of biosafety level 3 facilities, reducing the risk of contamination and the cost of the assay. Plasmids have a high stability, are normally maintained in Escherichia coli and can be extracted in large amounts.     

Virtual screening to identify novel potential inhibitors for Glutamine synthetase of Mycobacterium tuberculosis

Abstract

Glutamine synthetase (GS) of Mycobacterium tuberculosis (Mtb) is an essential enzyme which is involved in nitrogen metabolism and cell wall synthesis. It is involved in the inhibition of phagosome-lysosome fusion by preventing acidification. Targeting GS can be helpful to control the infection of Mtb. In order to identify potential inhibitors, we screened chemical libraries (56,400 compounds of ChEMBL anti-mycobacterial, 1596 FDA approved drugs, 419 Natural product and 916 phytochemical) against this target using the virtual screening approach. Screening by molecular docking identified ten top-ranked compounds as GSMtb inhibitors and they were compared with known inhibitors (as control). Since GS enzyme (GSHs) is also present in human. We have compared the protein sequence of GS from Mtb and human using the P-BLAST in NCBI. We found ～27% identity in between these two sequences, so we also compared the binding affinity of inhibitor between Mtb and human. Finally, we identified top two compounds namely CHEMBL387509, CHEMBL226198 from ChEMBL anti-mycobacterial dataset, and Eriocitrin and Malvidin from phytochemical dataset which showed lees binding affinity towards GSHs whereas Pamidronate, and Phentermine from FDA approved drugs and (-)-Quinic Acid, Hexopyranuronic acid, Quebrachit, and Castanospermine from natural product showed protein-ligand interaction with Mtb protein while no interaction with GSHs. The top two docked complexes were subjected to molecular dynamic simulation to understand the stability of the molecule. Further, we calculated the binding free energy of the docked complex and analyzed hydrogen bond, salt bridge, pie stacking, and hydrophobic interaction in the docking region. These ligands exhibited very good binding affinity GSMtb enzymes. Therefore, these ligands are novel and drug-likeness compounds, and they may be potential inhibitors of M tuberculosis.

Communicated by Ramaswamy H. Sarma     

Genetics-directed drug discovery for combating Mycobacterium tuberculosis infection

Mycobacterium tuberculosis (Mtb), the pathogen of tuberculosis (TB), is one of the most infectious bacteria in the world. The traditional strategy to combat TB involves targeting the pathogen directly; however, the rapid evolution of drug resistance lessens the efficiency of this anti-TB method. Therefore, in recent years, some researchers have turned to an alternative anti-TB strategy, which hinders Mtb infection through targeting host genes. In this work, using a theoretical genetic analysis, we identified 170 Mtb infection-associated genes from human genetic variations related to Mtb infection. Then, the agents targeting these genes were identified to have high potential as anti-TB drugs. In particular, the agents that can target multiple Mtb infection-associated genes are more druggable than the single-target counterparts. These potential anti-TB agents were further screened by gene expression data derived from connectivity map. As a result, some agents were revealed to have high interest for experimental evaluation. This study not only has important implications for anti-TB drug discovery, but also provides inspirations for streamlining the pipeline of modern drug discovery.     

Zika virus NS5 protein potential inhibitors: an enhanced in silico approach in drug discovery

The re-emerging Zika virus (ZIKV) is an arthropod-borne virus that has been described to have explosive potential as a worldwide pandemic. The initial transmission of the virus was through a mosquito vector, however, evolving modes of transmission has allowed the spread of the disease over continents. The virus has already been linked to irreversible chronic central nervous system conditions. The concerns of the scientific and clinical community are the consequences of Zika viral mutations, thus suggesting the urgent need for viral inhibitors. There have been large strides in vaccine development against the virus but there are still no FDA approved drugs available. Rapid rational drug design and discovery research is fundamental in the production of potent inhibitors against the virus that will not just mask the virus, but destroy it completely. In silico drug design allows for this prompt screening of potential leads, thus decreasing the consumption of precious time and resources. This study demonstrates an optimized and proven screening technique in the discovery of two potential small molecule inhibitors of ZIKV Methyltransferase and RNA dependent RNA polymerase. This in silico ‘per-residue energy decomposition pharmacophore’ virtual screening approach will be critical in aiding scientists in the discovery of not only effective inhibitors of Zika viral targets, but also a wide range of anti-viral agents.     

结核分枝杆菌三种耐药基因的检测方法

建立3种结核分枝杆菌耐药基因的检测方法,探讨耐药基因突变与耐药性的关系。将58株临床分离株均做聚合酶链反应-单链构象多态性分析(PCR—SSCP)和传统药物敏感试验。结果表明,结核分枝杆菌耐药基因突变与耐药水平有密切联系,绝大多数结核分枝杆菌耐药基因突变发生在高耐药株,少部分在低耐药株发生基因突变。     

Structural insight into Mycobacterium tuberculosis maltosyl transferase inhibitors: pharmacophore-based virtual screening,docking, and molecular dynamics simulations

Pharmacophore-based virtual screening, subsequent docking, and molecular dynamics (MD) simulations have been done to identify potential inhibitors of maltosyl transferase of Mycobacterium tuberculosis (mtb GlgE). Ligand and structure-based pharmacophore models representing its primary binding site (pbs) and unique secondary binding site 2 (sbs2), respectively, were constructed based on the three dimensional structure of mtb GlgE. These pharmacophore models were further used for screening of ZINC and antituberculosis compounds database (ATD). Virtually screened molecules satisfying Lipinski’s rule of five were then analyzed using docking studies and have identified 23 molecules with better binding affinity than its natural substrate, maltose. Four top scoring ligands from ZINC and ATD that either binds to pbs or sbs2 have been subjected to 10 ns each MD simulations and binding free energy calculations. Results of these studies have confirmed stable protein ligand binding. Results reported in the article are likely to be helpful in antitubercular therapeutic development research.     

Structure-based screening and molecular dynamics simulations offer novel natural compounds as potential inhibitors of Mycobacterium tuberculosis isocitrate lyase

Mycobacterium tuberculosis is the etiological agent of tuberculosis in humans and is responsible for more than two million deaths annually. M. tuberculosis isocitrate lyase (MtbICL) catalyzes the first step in the glyoxylate cycle, plays a pivotal role in the persistence of M. tuberculosis, which acts as a potential target for an anti-tubercular drug. To identify the potential anti-tuberculosis compound, we conducted a structure-based virtual screening of natural compounds from the ZINC database (n = 1,67,748) against the MtbICL structure. The ligands were docked against MtbICL in three sequential docking modes that resulted in 340 ligands having better docking score. These compounds were evaluated for Lipinski and ADMET prediction, and 27 compounds were found to fit well with re-docking studies. After refinement by molecular docking and drug-likeness analyses, three potential inhibitors (ZINC1306071, ZINC2111081, and ZINC2134917) were identified. These three ligands and the reference compounds were further subjected to molecular dynamics simulation and binding energy analyses to compare the dynamic structure of protein after ligand binding and the stability of the MtbICL and bound complexes. The binding free energy analyses were calculated to validate and capture the intermolecular interactions. The results suggested that the three compounds had a negative binding energy with ?96.462, ?143.549, and ?122.526 kJ mol^?1 for compounds with IDs ZINC1306071, ZINC2111081, and ZINC2134917, respectively. These lead compounds displayed substantial pharmacological and structural properties to be drug candidates. We concluded that ZINC2111081 has a great potential to inhibit MtbICL and would add to the drug discovery process against tuberculosis.     

Crystallographic and pre-steady-state kinetics studies on binding of NADH to wild-type and isoniazid-resistant enoyl-ACP(CoA) reductase enzymes from Mycobacterium tuberculosis

An understanding of isoniazid (INH) drug resistance mechanism in Mycobacterium tuberculosis should provide significant insight for the development of newer anti-tubercular agents able to control INH-resistant tuberculosis (TB). The inhA-encoded 2-trans enoyl-acyl carrier protein reductase enzyme (InhA) has been shown through biochemical and genetic studies to be the primary target for INH. In agreement with these results, mutations in the inhA structural gene have been found in INH-resistant clinical isolates of M.tuberculosis, the causative agent of TB. In addition, the InhA mutants were shown to have higher dissociation constant values for NADH and lower values for the apparent first-order rate constant for INH inactivation as compared to wild-type InhA. Here, in trying to identify structural changes between wild-type and INH-resistant InhA enzymes, we have solved the crystal structures of wild-type and of S94A, I47T and I21V InhA proteins in complex with NADH to resolutions of, respectively, 2.3A, 2.2A, 2.0 A, and 1.9A. The more prominent structural differences are located in, and appear to indirectly affect, the dinucleotide binding loop structure. Moreover, studies on pre-steady-state kinetics of NADH binding have been carried out. The results showed that the limiting rate constant values for NADH dissociation from the InhA-NADH binary complexes (k(off)) were eleven, five, and tenfold higher for, respectively, I21V, I47T, and S94A INH-resistant mutants of InhA as compared to INH-sensitive wild-type InhA. Accordingly, these results are proposed to be able to account for the reduction in affinity for NADH for the INH-resistant InhA enzymes.     

Dynamics of fluoroquinolones induced resistance in DNA gyrase of Mycobacterium tuberculosis

DNA gyrase is a validated target of fluoroquinolones which are key components of multidrug resistance tuberculosis (TB) treatment. Most frequent occurring mutations associated with high level of resistance to fluoroquinolone in clinical isolates of TB patients are A90V, D94G, and A90V–D94G (double mutant [DM]), present in the larger subunit of DNA Gyrase. In order to explicate the molecular mechanism of drug resistance corresponding to these mutations, molecular dynamics (MD) and mechanics approach was applied. Structure-based molecular docking of complex comprised of DNA bound with Gyrase A (large subunit) and Gyrase C (small subunit) with moxifloxacin (MFX) revealed high binding affinity to wild type with considerably high Glide XP docking score of ?7.88 kcal/mol. MFX affinity decreases toward single mutants and was minimum toward the DM with a docking score of ?3.82 kcal/mol. Docking studies were also performed against 8-Methyl-moxifloxacin which exhibited higher binding affinity against wild and mutants DNA gyrase when compared to MFX. Molecular Mechanics/Generalized Born Surface Area method predicted the binding free energy of the wild, A90V, D94G, and DM complexes to be ?55.81, ?25.87, ?20.45, and ?12.29 kcal/mol, respectively. These complexes were further subjected to 30 ns long MD simulations to examine significant interactions and conformational flexibilities in terms of root mean square deviation, root mean square fluctuation, and strength of hydrogen bond formed. This comparative drug interaction analysis provides systematic insights into the mechanism behind drug resistance and also paves way toward identifying potent lead compounds that could combat drug resistance of DNA gyrase due to mutations.     

Specific mutations in the Mycobacterium tuberculosis rpoB gene are associated with increased dnaE2 expression

In Mycobacterium tuberculosis (MTB), rifampicin resistance is almost invariably due to mutations in the rpoB gene, whose function is critical for cell viability. Most of these mutations, at least initially, impair the fitness of the bacteria but confer a selective advantage when antibiotic pressure is exerted. Subsequent adaptation may be critical to restore fitness. The possibility was considered that MTB with mutations in the rpoB gene elicits a constitutive stress response, increasing the probability of subsequent adaptation. In order to test this hypothesis, the expression of recA and dnaE2, an inducible putative error-prone DNA polymerase, was determined in six different isogenic laboratory-generated rpoB-mutants of MTB. Expression levels were determined with real-time PCR and the data obtained were compared with those of the wild-type parent. In four of the six rpoB mutants, a two- to fivefold induction of dnaE2 was detected (P<0.05). Thus, the presence of specific mutations in rpoB is not only associated with impaired fitness but also results in a detectable, moderate yet persistent increase in the expression of dnaE2 but not recA.     

Identification of novel leads applying in silico studies for Mycobacterium multidrug resistant (MMR) protein

Multidrug efflux mechanism is the main cause of intrinsic drug resistance in bacteria. Mycobacterium multidrug resistant (MMR) protein belongs to small multidrug resistant family proteins (SMR), causing multidrug resistance to proton (H⁺)-linked lipophilic cationic drug efflux across the cell membrane. In the present work, MMR is treated as a novel target to identify new molecular entities as inhibitors for drug resistance in Mycobacterium tuberculosis. In silico techniques are applied to evaluate the 3D structure of MMR protein. The putative amino acid residues present in the active site of MMR protein are predicted. Protein–ligand interactions are studied by docking cationic ligands transported by MMR protein. Virtual screening is carried out with an in-house library of small molecules against the grid created at the predicted active site residues in the MMR protein. Absorption distribution metabolism and elimination (ADME) properties of the molecules with best docking scores are predicted. The studies with cationic ligands and those of virtual screening are analysed for identification of new lead molecules as inhibitors for drug resistance caused by the MMR protein.     

结核耐药基因的基因芯片检测及临床应用

目的:采用基因芯片技术对结核分枝杆菌中常见耐药基因rpoB、katG及inhA进行检测,以了解结核分枝杆菌的耐药情况,及基因芯片技术检测结核菌耐药基因的临床应用价值。方法:收集40例涂片抗酸染色阳性并经分枝杆菌菌种鉴定芯片鉴定为结核的样本进行结核耐药基因检测。结果:40例样本中,14例无法判读结果,占35%,检出26例,检出率为65%。其中,无突变的野生型21例,占52.5%;突变型5例,总突变率为12.5%;3例rpoB基因的531点单独突变(TCG→TTG),突变率为7.5%;2例katG基因的315点单独突变(AGC→ACC),突变率为5%。结论:结核耐药基因芯片试剂盒检测结核菌耐药基因时针对单个菌落,用痰样本直接检测耐药基因虽能简便快速地了解结核分枝杆菌的耐药情况,但会出现一些无法判读的结果,原因须进一步探讨。     

Impact of drug resistance on fitness of Mycobacterium tuberculosis strains of the W-Beijing genotype

Mycobacterium tuberculosis strains of the W-Beijing genotype became a common cause of tuberculosis during the past years and they are often associated with drug resistance. The biological factors facilitating the selection and wide dissemination of these strains are not known. To determine how acquisition of drug resistance affected growth of strains of the W-Beijing genotype, the growth of 55 M. tuberculosis isolates were studied using the BBL MGIT Mycobacteria Growth Indicator Tube and the BACTEC MGIT 960 System. Susceptible strains of non-Beijing genotypes were found to be the most fit strains. Drug-resistant strains of non-Beijing genotypes were more likely to grow slower than susceptible strains (P=0.001). Drug-resistant strains of the W-Beijing genotype had two tendencies of growth: some of them showed reduced growth compared to susceptible strains, while others did not show loss of fitness measured as growth.     

Horizontal transfer of a virulence operon to the ancestor of Mycobacterium tuberculosis

The contribution of interspecies horizontal gene transfer (HGT) to the evolution and virulence of Mycobacterium tuberculosis, the agent of tuberculosis in humans, has been barely investigated. Here we have studied the evolutionary history of the M. tuberculosis Rv0986-8 virulence operon recently identified, through functional genomics approaches, as playing an important role in parasitism of host phagocytic cells. We showed that among actinobacteria, this operon is specific to the M. tuberculosis complex and to ancestral Mycobacterium prototuberculosis species. These data, together with phylogenetic reconstruction and other in silico analyses, provided strong evidence that this operon has been acquired horizontally by the ancestor of M. tuberculosis, before the recent evolutionary bottleneck that preceded the clonal-like evolution of the M. tuberculosis complex. Genomic signature profiling further suggested that the transfer was plasmid mediated and that the operon originated from a gamma-proteobacterium donor species. Our study points out for the first time the contribution of HGT to the emergence of M. tuberculosis and close relatives as major pathogens. In addition, our data underline the importance of deciphering gene transfer networks in M. tuberculosis in order to better understand the evolutionary mechanisms involved in mycobacterial virulence.     

结核分枝杆菌对异烟肼耐药的分子机制

抗结核一线药物异烟肼是应用最广泛的抗结核药物之一,自1952年应用于临床以来,异烟肼就成了治疗结核和潜在感染的基础药物.有报道,我国异烟肼耐药已排在首位.结核分枝杆菌对异烟肼耐药的分子机制十分复杂,涉及katG、inhA、kasA、ndh、axyR等多种基因,本研究仅就此方面的研究作一综述.     

A population-based study of first and second-line drug-resistant tuberculosis in a high-burden area of the Mexico/United States border

The resistance of 139 Mycobacterium tuberculosis (MTB) isolates from the city of Monterrey, Northeast Mexico, to first and second-line anti-TB drugs was analysed. A total of 73 isolates were susceptible and 66 were resistant to anti-TB drugs. Monoresistance to streptomycin, isoniazid (INH) and ethambutol was observed in 29 cases. Resistance to INH was found in 52 cases and in 29 cases INH resistance was combined with resistance to two or three drugs. A total of 24 isolates were multidrug-resistant (MDR) resistant to at least INH and rifampicin and 11 MDR cases were resistant to five drugs. The proportion of MDR-TB among new TB cases in our target population was 0.72% (1/139 cases). The proportion of MDR-TB among previously treated cases was 25.18% (35/139 cases). The 13 polyresistant and 24 MDR isolates were assayed against the following seven second-line drugs: amikacin (AMK), kanamycin (KAN), capreomycin (CAP), clofazimine (CLF), ethionamide (ETH), ofloxacin (OFL) and cycloserine (CLS). Resistance to CLF, OFL or CLS was not observed. Resistance was detected to ETH (10.80%) and to AMK (2.70%), KAN (2.70%) and CAP (2.70%). One isolate of MDR with primary resistance was also resistant to three second-line drugs. Monterrey has a high prevalence of MDR-TB among previously treated cases and extensively drug-resistant-MTB strains may soon appear.     

In silico insights into the identification of potential novel angiogenic inhibitors against human VEGFR-2: a new SAR-based hierarchical clustering approach

Abstract

In this study, binding efficiency of new pyrrolopyrimidine structural analogs against human vascular endothelial growth factor receptor-2 (VEGFR-2) were elucidated using integrated in silico methods. Optimized high-resolution model of VEGFR-2 was generated and adopted for structure-based virtual screening approaches. Pyrrolopyrimidine inhibitor (CP15) associated compounds were screened from PubChem database and subjected to virtual screening and comparative docking methods against the receptor ligand-binding domain. Accordingly, high efficient compounds were clustered with similarity indices through PubChem structure cluster module using single-linkage algorithm. Moreover, pharmacokinetics including drug metabolism activities of high-binding leads under investigation was portrayed using ADMET and similarity ensemble analysis. Optimal energy orientations of the selected protein model have been shown to be reliable, and highly recommended for screening and docking studies. Docking and clustering strategies were shown that nineteen candidates as most effective binders for VEGFR-2 than CP15, and are grouped as three classes. Lys⁸⁶⁸, Glu⁸⁸⁵, Cys⁹¹⁹, His¹⁰²⁶, Arg¹⁰²⁷, Asp¹⁰⁴⁶, and Gly¹⁰⁴⁸ residues were predicted as novel hotspot residues, and participate in H-bonds, π-cation, π-stacking, halogen bonds, and salt-bridges formation with ligands. These additional bonds are contributing extent stability that holds the receptor structure at flexible state, this make difficult to any further conformational changes for evoking angiogenic signals. The ADMET and similarity ensemble analysis results were strongly indicated that thirteen candidates as best ligands for angiogenesis targets. Altogether, these findings indicate potential angiogenic templates and their binding levels with VEGFR-2; sorted viewpoints could be useful as a promising way to describe potential angiogenesis inhibitors with related molecular targets.     

A novel genotypic test for rapid detection of multidrug-resistant Mycobacterium tuberculosis isolates by a multiplex probe array

AIMS: To develop and evaluate a novel genotypic test for rapid detection of rifampicin and isoniazid resistance of multidrug-resistant (MDR) Mycobacterium tuberculosis isolates by a multiplex probe array. METHODS AND RESULTS: A multiplex probe array was designed for genotypic test to simultaneously screen the mutations of rpoB, katG, inhA and ahpC genes, associated with rifampin and isoniazid resistance in M. tuberculosis, with a probe detecting one of the recently confirmed genetic markers of isoniazid resistance ahpC-6 and -9 locus added. By using the genotypic test developed, 52 MDR isolates were identified, among which 46 isolates had mutations in rpoB (88.5%) and 45 at codon 315 of katG, regulatory region of inhA and oxyR-ahpC intergenic region (86.5%), whereas all 35 susceptible isolates identified showed a wild-type hybridization pattern. The sensitivity and specificity were 88.5% and 100% for rifampicin resistance, and 86.5% and 100% for isoniazid resistance, respectively. CONCLUSION: A rapid and simultaneous detection of rifampicin and isoniazid resistance caused by the mutations of rpoB, katG, inhA and ahpC genes in M. tuberculosis isolates could be achieved by a multiplex probe array developed. SIGNIFICANCE AND IMPACT OF THE STUDY: This genotypic test protocol has the potential to be developed on clinical application for the rapid detection of drug resistant M. tuberculosis isolates before an efficient chemotherapy is initiated.     

Combined pharmacophore modeling, 3D-QSAR and docking studies to identify novel HDAC inhibitors using drug repurposing

Abstract

Histone deacetylases (HDACs), a critical family of epigenetic enzymes, has emerged as a promising target for antitumor drugs. Here, we describe our protocol of virtual screening in identification of novel potential HDAC inhibitors through pharmacophore modeling, 3D-QSAR, molecular docking and molecular dynamics (MD) simulation. Considering the limitation of current virtual screening works, drug repurposing strategy was applied to discover druggable HDAC inhibitor. The ligand-based pharmacophore and 3D-QSAR models were established, and their reliability was validated by different methods. Then, the DrugBank database was screened, followed by molecular docking. MD simulation (100?ns) was performed to further study the stability of ligand binding modes. Finally, results indicated the hit DB03889 with high in silico inhibitory potency was suitable for further experimental analysis.

Communicated by Ramaswamy H. Sarma