The DprE1 enzyme, one of the most vulnerable targets of Mycobacterium tuberculosis

The re-emergence of tuberculosis in recent years led the World Health Organization (WHO) to launch the Stop TB Strategy program. Beside repurposing the existing drugs and exploring novel molecular combinations, an essential step to face the burden of tuberculosis will be to develop new drugs by identifying vulnerable bacterial targets. Recent studies have focused on decaprenylphosphoryl-d-ribose oxidase (DprE1) of Mycobacterium tuberculosis, an essential enzyme involved in cell wall metabolism, for which new promising molecules have proved efficacy as antitubercular agents. This review summarizes the state of the art concerning DprE1 in terms of structure, enzymatic activity and inhibitors. This enzyme is emerging as one of the most vulnerable target in M. tuberculosis.     

Human CD8+ T-cells Recognizing Peptides from Mycobacterium tuberculosis (Mtb) Presented by HLA-E Have an Unorthodox Th2-like,Multifunctional, Mtb Inhibitory Phenotype and Represent a Novel Human T-cell Subset

Mycobacterial antigens are not exclusively presented to T-cells by classical HLA-class Ia and HLA-class II molecules, but also through alternative antigen presentation molecules such as CD1a/b/c, MR1 and HLA-E. We recently described mycobacterial peptides that are presented in HLA-E and recognized by CD8⁺ T-cells. Using T-cell cloning, phenotyping, microbiological, functional and RNA-expression analyses, we report here that these T-cells can exert cytolytic or suppressive functions, inhibit mycobacterial growth, yet express GATA3, produce Th2 cytokines (IL-4,-5,-10,-13) and activate B-cells via IL-4. In TB patients, Mtb specific cells were detectable by peptide-HLA-E tetramers, and IL-4 and IL-13 were produced following peptide stimulation. These results identify a novel human T-cell subset with an unorthodox, multifunctional Th2 like phenotype and cytolytic or regulatory capacities, which is involved in the human immune response to mycobacteria and demonstrable in active TB patients’ blood. The results challenge the current dogma that only Th1 cells are able to inhibit Mtb growth and clearly show that Th2 like cells can strongly inhibit outgrowth of Mtb from human macrophages. These insights significantly expand our understanding of the immune response in infectious disease.     

Carbohydrate Recognition by RpfB from Mycobacterium tuberculosis Unveiled by Crystallographic and Molecular Dynamics Analyses

Resuscitation of Mtb is crucial to the etiology of Tuberculosis, because latent tuberculosis is estimated to affect one-third of the world population. The resuscitation-promoting factor RpfB is mainly responsible for Mtb resuscitation from dormancy. Given the impact of latent Tuberculosis, RpfB represents an interesting target for tuberculosis drug discovery. However, no molecular models of substrate binding and catalysis are hitherto available for this enzyme. Here, we identified key interactions involved in substrate binding to RpfB by combining x-ray diffraction studies and computational approaches. The crystal structure of RpfB catalytic domain in complex with N,N′,N″-triacetyl-chitotriose, as described here, provides the first, to our knowledge, atomic representation of ligand recognition by RpfB and demonstrates that the strongest interactions are established by the N-acetylglucosamine moiety in the central region of the enzyme binding cleft. Molecular dynamics analyses provided information on the dynamic behavior of protein-substrate interactions and on the role played by the solvent in RpfB function. These data combined with sequence conservation analysis suggest that Glu-292 is the sole residue crucial for catalysis, implying that RpfB acts via the formation of an oxocarbenium ion rather than a covalent intermediate. Present data represent a solid base for the design of effective drug inhibitors of RpfB. Moreover, homology models were generated for the catalytic domains of all members of the Mtb Rpf family (RpfA-E). The analysis of these models unveiled analogies and differences among the different members of the Rpf protein family.Abbreviations used: Mtb, Mycobacterium tuberculosis; RpfB, Resuscitation promoting factor B; PDB, Protein Data Bank; RpfBc, catalytic domain of RpfB; NAG3, N,N'',N"-triacetyl-chitotriose; NAG6, hexa-N- acetylchitohexaose; MD, molecular dynamics; RMSF, root mean-square fluctuation     

Molecular Dynamics Simulation of CaF2 Structure in Crystalline,Superionic and Molten Phases

Abstract

We describe here a number of molecular dynamics simulations on calcium fluoride over a range of temperatures spanning the transitions to the superionic and molten state. The simulation temperatures are 1400, 1590, 1800, 2000, 2200 K. By using the bond spherical harmonics method with equal neighbor number, we have studied the structure and bond orientation of cation sublattice and anion sublattice in superionic conductor CaF₂. The bond order parameters Q₁ have been calculated both for the computer generated instataneous configurations from the simulation system and for the standard configurations from the normal distribution model of bond orientation. The comparison of Q₁ between the molecular dynamics simulation and the normal distribution model shows that not only the cation sublattice but also the anion sublattice can be described by the normal distribution model. The cations keep their original fcc frame, but in the anion case there is a great deal of random distortions from the original anion sublattice.     

The Role of Electrostatic Interaction in Triggering the Unraveling of Stable Helix 1 in Normal Prion Protein. A Molecular Dynamics Simulation Investigation

Abstract

The conversion of normal prion protein (PrP^C) into scrapie isoform (PrP^Sc) is a key event in the pathogenesis of prion diseases. However, the conversion mechanism has given rise to much controversy. For instance, there is much debate on the behavior of helix 1 (H1) in the conversion. A series of experiments demonstrated that H1 in isolated state was very stable under a variety of conditions. But, other experiments indicated that helices 2 and 3 rather than H1 were retained in PrP^Sc. In this paper, molecular dynamics (MD) simulation is employed to investigate the dynamic behavior of H1. It is revealed that although the helix 1 of Human PrP^C (HuPrP^C) is very stable in the isolated state, it becomes unstable when incorporated into native HuPrP^C, which likely results from the long-range electrostatic interaction between Asp147 and Arg208 located in the helices 1 and 3, respectively. This explanation is supported by experimental evaluation and MD simulation on D147N mutant of HuPrP^C that the mutant becomes a little more stable than the wild type HuPrP^C. This finding not only help to reconcile the existing debate on the role of helix 1 in the PrP^C→PrP^Sc transition, but also reveals a possible mechanism for triggering the PrP^C→PrP^Sc conversion.     

Structure and Stability of the Insulin Dimer Investigated by Molecular Dynamics Simulation

Abstract

Molecular dynamics simulation indicates that the dynamical behaviour of the insulin dimer is asymmetric. Atomic level knowledge of the interaction modes and protein conformation in the solvation state identifies dynamical structures, held by hydrogen bonds that stabilize, mainly in one monomer, the interaction between the chains. Dynamic cross-correlation analysis shows that the two insulin monomers behave asymmetrically and are almost independent. Solvation energy, calculated to evaluate the contribute of each interface residue to the dimer association pattern, well compares with the experimental association state found in protein mutants indicating that this parameter is an important factor to explain the association properties of mutated insulin dimers.     

Structure and Mechanical Characterization of DNA i-Motif Nanowires by Molecular Dynamics Simulation

We studied the structure and mechanical properties of DNA i-motif nanowires by means of molecular dynamics computer simulations. We built up to 230 nm-long nanowires, based on a repeated TC₅ sequence from crystallographic data, fully relaxed and equilibrated in water. The unusual C⋅C⁺ stacked structure, formed by four ssDNA strands arranged in an intercalated tetramer, is here fully characterized both statically and dynamically. By applying stretching, compression, and bending deformations with the steered molecular dynamics and umbrella sampling methods, we extract the apparent Young’s and bending moduli of the nanowire, as well as estimates for the tensile strength and persistence length. According to our results, the i-motif nanowire shares similarities with structural proteins, as far as its tensile stiffness, but is closer to nucleic acids and flexible proteins, as far as its bending rigidity is concerned. Furthermore, thanks to its very thin cross section, the apparent tensile toughness is close to that of a metal. Besides their yet to be clarified biological significance, i-motif nanowires may qualify as interesting candidates for nanotechnology templates, due to such outstanding mechanical properties.     

Structure and Mechanical Characterization of DNA i-Motif Nanowires by Molecular Dynamics Simulation

We studied the structure and mechanical properties of DNA i-motif nanowires by means of molecular dynamics computer simulations. We built up to 230 nm-long nanowires, based on a repeated TC₅ sequence from crystallographic data, fully relaxed and equilibrated in water. The unusual C⋅C⁺ stacked structure, formed by four ssDNA strands arranged in an intercalated tetramer, is here fully characterized both statically and dynamically. By applying stretching, compression, and bending deformations with the steered molecular dynamics and umbrella sampling methods, we extract the apparent Young’s and bending moduli of the nanowire, as well as estimates for the tensile strength and persistence length. According to our results, the i-motif nanowire shares similarities with structural proteins, as far as its tensile stiffness, but is closer to nucleic acids and flexible proteins, as far as its bending rigidity is concerned. Furthermore, thanks to its very thin cross section, the apparent tensile toughness is close to that of a metal. Besides their yet to be clarified biological significance, i-motif nanowires may qualify as interesting candidates for nanotechnology templates, due to such outstanding mechanical properties.     

Molecular Dynamics Simulation Studies of dTTP Binding and Catalysis Mediated by YhdE Dimerization

YhdE is a Maf-like (multicopy associated filamentation) protein that primarily acts as dTTPase to hydrolyze dTTP into dTMP and two phosphate molecules in cell metabolism pathway. Two crystal structures of YhdE have been previously determined, representing the open and closed active site conformations, respectively. Based on the structures, we have carried out molecular dynamics simulations and free energy calculations to investigate dTTP binding to and hydrolysis by YhdE. Our results suggest that YhdE closed state is structurally more compact than its open state at room temperature. YhdE open state is a favorable conformation for dTTP binding and closed state is a structurally favorable conformation for catalytic reaction. This observation is supported by the structure of YhdE homolog in complex with a nucleotide analog. Free energy calculations reveal that YhdE dimerization occurs preferentially in dTTP binding and is favorable for successive cooperative reaction. The key residues R11, R12 and K80, are found to contribute to the substrate stabilization. Further, YhdE dimerization and binding of dTTP induce the cooperative effect through a direct allosteric communication network in YhdE from the dTTP binding sites in the catalytic center to the intermolecular β-strand in YhdE dimer.     

Molecular structure and peptidoglycan recognition of Mycobacterium tuberculosis ArfA (Rv0899)

Mycobacterium tuberculosis ArfA (Rv0899) is a membrane protein encoded by an operon that is required for supporting bacterial growth in acidic environments. Its C-terminal domain (C domain) shares significant sequence homology with the OmpA-like family of peptidoglycan-binding domains, suggesting that its physiological function in acid stress protection may be related to its interaction with the mycobacterial cell wall. Previously, we showed that ArfA forms three independently structured modules, and we reported the structure of its central domain (B domain). Here, we describe the high-resolution structure and dynamics of the C domain, we identify ArfA as a peptidoglycan-binding protein and we elucidate the molecular basis for its specific recognition of diaminopimelate-type peptidoglycan. The C domain of ArfA adopts the characteristic fold of the OmpA-like family. It exhibits pH-dependent conformational dynamics (with significant heterogeneity at neutral pH and a more ordered structure at acidic pH), which could be related to its acid stress response. The C domain associates tightly with polymeric peptidoglycan isolated from M. tuberculosis and also associates with a soluble peptide intermediate of peptidoglycan biosynthesis. This enabled us to characterize the peptidoglycan binding site where five highly conserved ArfA residues, including two key arginines, establish the specificity for diaminopimelate- but not Lys-type peptidoglycan. ArfA is the first peptidoglycan-binding protein to be identified in M. tuberculosis. Its functions in acid stress protection and peptidoglycan binding suggest a link between the acid stress response and the physicochemical properties of the mycobacterial cell wall.     

ZNC(C)PR的构象及分子动力学模拟

本文用分子动力学方法和质子间距约束结合建立了ZNC(C)PR的构象。质子间距是由二维核磁共振的NOE效应得到的。结果表明ZNC(C)PR分子的主链是伸展的,但Asn2的侧链与Arg5的侧链相互靠近,使整个分子成紧密结构。在所建模型的基础上,对ZNC(C)PR进行了分子动力学摸拟。发现第一位的焦谷氨酸环很灵活,这可能意味着焦谷氨酸环对ZNC(C)PR的活性并不重要,这与实验结果是一致的。     

Molecular Modeling of Mycobacterium Tuberculosis dUTpase: Docking and Catalytic Mechanism Studies

Abstract

Mycobacterium tuberculosis is a leading cause of infectious disease in the world today. This outlook is aggravated by a growing number of M. tuberculosis infections in individuals who are immunocompromised as a result of HIV infections. Thus, new and more potent anti-TB agents are necessary. Therefore, dUTpase was selected as a target enzyme to combat M. tuberculosis. In this work, molecular modeling methods involving docking and QM/MM calculations were carried out to investigate the binding orientation and predict binding affinities of some potential dUTpase inhibitors. Our results suggest that the best potential inhibitor investigated, among the compounds studied in this work, is the compound dUPNPP. Regarding the reaction mechanism, we concluded that the decisive stage for the reaction is the stage 1. Furthermore, it was also observed that the compounds with a ?1 electrostatic charge presented lower activation energy in relation to the compounds with a ?2 charge.     

Crystal Structure of AhpE from Mycobacterium tuberculosis, a 1-Cys peroxiredoxin

All living systems require protection against the damaging effects of reactive oxygen species. The genome of Mycobacterium tuberculosis, the cause of TB, encodes a number of peroxidases that are thought to be active against organic and inorganic peroxides, and are likely to play a key role in the ability of this organism to survive within the phagosomes of macrophages. The open reading frame Rv2238c in M.tuberculosis encodes a 153-residue protein AhpE, which is a peroxidase of the 1-Cys peroxiredoxin (Prx) family. The crystal structure of AhpE, determined at 1.87 A resolution (R(cryst)=0.179, R(free)=0.210), reveals a compact single-domain protein with a thioredoxin fold. AhpE forms both dimers and octamers; a tightly-associated dimer and a ring-like octamer, generated by crystallographic 4-fold symmetry. In this native structure, the active site Cys45 is in its oxidized, sulfenic acid (S-O-H) state. A second crystal form of AhpE, obtained after soaking in sodium bromide and refined at 1.90 A resolution (R(cryst)=0.242, R(free)=0.286), reveals the reduced structure. In this structure, a conformational change in an external loop, in two of the four molecules in the asymmetric unit, allows Arg116 to stabilise the Cys45 thiolate ion, and concomitantly closes a surface channel. This channel is identified as the likely binding site for a physiological reductant, and the conformational change is inferred to be important for the reaction cycle of AhpE.     

Molecular Interaction of Novel Compound 2-Methylheptyl Isonicotinate Produced by Streptomyces sp. 201 with Dihydrodipicolinate Synthase (DHDPS) Enzyme of Mycobacterium tuberculosis for its Antibacterial Activity

Antibiotic resistance is a growing problem in multi-drug-resistant tuberculosis which is caused by Mycobacterium tuberculosis (MTB). Hence there is an urgent need for designing or developing a novel or potent anti-tubercular agent. The Lysine/DAP biosynthetic pathway is a promising target because of its role in cell wall and amino acid biosynthesis. In our study we performed a molecular docking analysis of a novel antibacterial isolated from Streptomyces sp. 201 at three different binding site of dihydrodipicolinate synthase (DHDPS) enzyme of MTB. The molecular docking studies suggest that the novel molecule shows favourable interaction at the three different binding sites as compared to five experimentally known inhibitors of DHDPS.     

抗菌肽BLFcin6与不同磷脂膜之间相互作用的分子动力学模拟研究

抗菌肽具有广谱抗菌特性,有望成为抗生素较好的替代产品.研究抗菌肽的抗菌机制,可以为新型抗菌肽的设计提供指导.无论抗菌肽采用哪种抗菌机制,其首先要稳定地吸附到细胞膜之上.因此,本文利用分子动力学模拟方法比较了抗菌小肽BLFcin6与5种不同细胞膜之间的相互作用.对这5种细胞膜而言,小肽会很快结合在POPG膜和DPPC-CHOL膜的表面,倾向于进入DPPC膜的疏水内部,与POPC膜和POPC-CHOL膜的接触很少.考察相互作用能,小肽与POPG膜的相互作用最强,主要是小肽与细胞膜亲水头部存在静电相互作用;小肽与DPPC膜的疏水尾部的相互作用较强,但受胆固醇影响,小肽只结合在DPPC-CHOL膜表面.在结合过程中,小肽N端的Arg会先结合到细胞膜上,静电相互作用在小肽锚定细胞膜的过程中起关键作用.以上研究从原子水平上解释了为什么BLFcin6小肽具有抗菌作用,哪些残基起关键作用,也为进一步开展BLFcin6小肽及其衍生小肽的研究奠定基础.     

Predicting the Structure of the Flavodoxin from Eschericia coli by Homology Modeling,Distance Geometry and Molecular Dynamics

Abstract

As part of our on-going development of a method, based upon distance geometry calculations, for predicting the structures of proteins from the known structures of their homologues, we have predicted the structure of the 176 residue Flavodoxin from Escherichia coli. This prediction was based upon the crystal structures of the homologous Flavodoxins from Anacystis nidulans, Chondrus crispus, Desulfovibrio vulgaris and Clostridium beijerinckii, whose sequence identities with Escherichia coli were 44%, 33%, 23% and 16%, respectively. A total of 13,043 distance constraints among the alpha-carbons of the Escherichia coli structure were derived from the sequence alignments with the known structures, together with 8,893 distance constraints among backbone and sidechain atoms of adjacent residues, 978 between the alpha-carbons and selected atoms of the flavin mononucleotide cofactor, 116 constraints to enforce conserved hydrogen bonds, and 452 constraints on the torsion angles in conserved residues. An ensemble of ten random Escherichia coli structures was computed from these constraints, with an average root mean square coordinate deviation (RMSD) among the alpha carbons of 0.85 Ångstroms (excluding the first 1 and last 6 residues, which have no corresponding residues in any of the homologues and hence were unconstrained); the corresponding average heavy-atom RMSD was 1.60 Å.

Since the distance geometry calculations were performed without hydrogen atoms, protons were added to the resulting structures and these structures embedded in a 50 × 50 × 40 Å solvent box with periodic boundary conditions. They were then subjected to a 20 picosecond dynamical simulated annealing procedure, starting at 300 K and gradually reduced to 10K, in which all the distance and torsion angle constraints were maintained by means of harmonic restraint functions. This was followed up by 1000 iterations of unrestrained conjugate gradients minimization. The goal of this energy refinement procedure was not to drastically modify the structures in an attempt at a priori prediction, but merely to improve upon the predictions obtained from the geometric constraints, particularly with regard to their local backbone and sidechain conformations and their hydrogen bonds. The resulting structures differed from the respective starting structures by an average of 1.52 Å in their heavy atom RMSD's, while the average RMSD among the heavy atoms of residues 2-170 increased slightly to 1.66 Å. We hope these structures will be good enough to enable the phase problem to be solved for the crystallographic data that is now being collected on this protein.     

Direct and simultaneous identification of Mycobacterium tuberculosis complex (MTBC) and Mycobacterium tuberculosis (MTB) by rapid multiplex nested PCR-ICT assay

The Mycobacterium tuberculosis (MTB) shows different virulence and host infection range from other members of the M. tuberculosis complex (MTBC). Differential identification of MTB from MTBC is thus important in certain occasions. The currently commercially available molecular assays which use either IS6110 or 16S rDNA fragment as identification targets are mainly designed for identifying MTBC but not for MTB. Comparative genomic DNA analysis has provided valuable information on regions of difference (RD) present in MTB but not in other members of the MTBC. RD9 region is further suggested to be a potential target for differential identification of MTB from MTBC. In this study, using IS6110 and Rv3618 (belong to RD9) as the specific identification targets for MTBC and MTB, respectively, we developed and tested a multiplex nested PCR-ICT (immuno-chromatography test) assay for simultaneously and directly detecting not only MTBC but also MTB from 1500 clinical sputum specimens. The results were compared with traditional culture and biochemical identification results together with patients' clinical assessments. This assay showed a 95.5% sensitivity, 97.9% specificity, 2.1% false positive rate and 4.5% false negative rate towards detection of MTBC, and a 93.0% sensitivity, 99.8% specificity, 0.2% false positive rate and 7.0% false negative rate for detection of MTB. This detection system shows great potential in clinical application.