Role of Epistasis in Amikacin, Kanamycin, Bedaquiline, and Clofazimine Resistance in Mycobacterium tuberculosis Complex.

Antibiotic resistance among bacterial pathogens poses a major global health threat. Mycobacterium tuberculosis complex (MTBC) is estimated to have the highest resistance rates of any pathogen globally. Given the low growth rate and the need for a biosafety level 3 laboratory, the only realistic avenue to scale up drug susceptibility testing (DST) for this pathogen is to rely on genotypic techniques. This raises the fundamental question of whether a mutation is a reliable surrogate for phenotypic resistance or whether the presence of a second mutation can completely counteract its effect, resulting in major diagnostic errors (i.e., systematic false resistance results). To date, such epistatic interactions have only been reported for streptomycin that is now rarely used. By analyzing more than 31,000 MTBC genomes, we demonstrated that the eis C-14T promoter mutation, which is interrogated by several genotypic DST assays endorsed by the World Health Organization, cannot confer resistance to amikacin and kanamycin if it coincides with loss-of-function (LoF) mutations in the coding region of eis. To our knowledge, this represents the first definitive example of antibiotic reversion in MTBC. Moreover, we raise the possibility that mmpR (Rv0678) mutations are not valid markers of resistance to bedaquiline and clofazimine if these coincide with an LoF mutation in the efflux pump encoded by mmpS5 (Rv0677c) and mmpL5 (Rv0676c).

Application of Targeted Next-Generation Sequencing Assay on a Portable Sequencing Platform for Culture-Free Detection of Drug-Resistant Tuberculosis from Clinical Samples.

Targeted next-generation sequencing (tNGS) has emerged as a comprehensive alternative to existing methods for drug susceptibility testing (DST) of Mycobacterium tuberculosis from patient sputum samples for clinical diagnosis of drug-resistant tuberculosis (DR-TB). However, the complexity of sequencing platforms has limited their uptake in low-resource settings. The goal of this study was to evaluate the use of the tNGS-based DST solution Genoscreen Deeplex Myc-TB, for use on the compact, low-cost Oxford Nanopore Technologies MinION sequencer. One hundred four DNA samples extracted from smear-positive sputum sediments, previously sequenced using the Deeplex assay on an Illumina MiniSeq, were resequenced on MinION after applying a custom library preparation. MinION read quality, mapping statistics, and variant calling were computed using an in-house pipeline and compared to the reference MiniSeq data. The average percentage of MinION reads mapped to an H37RV reference genome was 90.8%, versus 99.5% on MiniSeq. The mean depths of coverage were 4,151× and 4,177× on MinION and MiniSeq, respectively, with heterogeneous distribution across targeted genes. Composite reference coverage breadth was >99% for both platforms. We observed full concordance between technologies in reporting the clinically relevant drug-resistant markers, including full gene deletions. In conclusion, we demonstrated that the workflow and sequencing data obtained from Deeplex on MinION are comparable to those for the MiniSeq, despite the higher raw error rates on MinION, with the added advantage of MinION's portability, versatility, and low capital costs. Targeted NGS on MinION is a promising DST solution for rapidly providing clinically relevant data to manage complex DR-TB cases.

NanoShaper-VMD interface: computing and visualizing surfaces, pockets and channels in molecular systems.

SUMMARY: NanoShaper is a program specifically aiming the construction and analysis of the molecular surface of nanoscopic systems. It uses ray-casting for parallelism and it performs analytical computations whenever possible to maximize robustness and accuracy of the approach. Among the other features, NanoShaper provides volume, surface area, including that of internal cavities, for any considered molecular system. It identifies pockets via a very intuitive definition based on the concept of probe radius, intrinsic to the definition of the solvent excluded surface. We show here that, with a suitable choice of the parameters, the same approach can also permit the visualisation of molecular channels. NanoShaper has now been interfaced with the widely used molecular visualization software VMD, further enriching its already well furnished toolset. AVAILABILITY AND IMPLEMENTATION: VMD is available at http://www.ks.uiuc.edu/Research/vmd/. NanoShaper, its documentation, tutorials and supporting programs are available at http://concept.iit.it/downloads. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

BiKi Life Sciences: A New Suite for Molecular Dynamics and Related Methods in Drug Discovery.

In this paper, we introduce the BiKi Life Sciences suite. This software makes it easy for computational medicinal chemists to run ad hoc molecular dynamics protocols in a novel and task-oriented environment; as a notebook, BiKi (acronym of Binding Kinetics) keeps memory of any activity together with dependencies among them. It offers unique accelerated protein-ligand binding/unbinding methods and other useful tools to gain actionable knowledge from molecular dynamics simulations and to simplify the drug discovery process.

Investigating drug-target association and dissociation mechanisms using metadynamics-based algorithms.

CONSPECTUS: This Account highlights recent advances and discusses major challenges in the field of drug-target recognition, binding, and unbinding studied using metadynamics-based approaches, with particular emphasis on their role in structure-based design. Computational chemistry has significantly contributed to drug design and optimization in an extremely broad range of areas, including prediction of target druggability and drug likeness, de novo design, fragment screening, ligand docking, estimation of binding affinity, and modulation of ADMET (absorption, distribution, metabolism, excretion, toxicity) properties. Computationally driven drug discovery must continuously adapt to keep pace with the evolving knowledge of the factors that modulate the pharmacological action of drugs. There is thus an urgent need for novel computational approaches that integrate the vast amount of complex information currently available for small (bio)organic compounds, biologically relevant targets and their complexes, while also accounting accurately for the thermodynamics and kinetics of drug-target association, the intrinsic dynamical behavior of biomolecular systems, and the complexity of protein-protein networks. Understanding the mechanism of drug binding to and unbinding from biological targets is fundamental for optimizing lead compounds and designing novel biologically active ones. One major challenge is the accurate description of the conformational complexity prior to and upon formation of drug-target complexes. Recently, enhanced sampling methods, including metadynamics and related approaches, have been successfully applied to investigate complex mechanisms of drugs binding to flexible targets. Metadynamics is a family of enhanced sampling techniques aimed at enhancing the rare events and reconstructing the underlying free energy landscape as a function of a set of order parameters, usually referred to as collective variables. Studies of drug binding mechanisms have predicted the most probable association and dissociation pathways and the related binding free energy profile. In addition, the availability of an efficient open-source implementation, running on cost-effective GPU (i.e., graphical processor unit) architectures, has considerably decreased the learning curve and the computational costs of the methods, and increased their adoption by the community. Here, we review the recent contributions of metadynamics and other enhanced sampling methods to the field of drug-target recognition and binding. We discuss how metadynamics has been used to search for transition states, to predict binding and unbinding paths, to treat conformational flexibility, and to compute free energy profiles. We highlight the importance of such predictions in drug discovery. Major challenges in the field and possible solutions will finally be discussed.

Oxidation-induced structural changes of ceruloplasmin foster NGR motif deamidation that promotes integrin binding and signaling.

Asparagine deamidation occurs spontaneously in proteins during aging; deamidation of Asn-Gly-Arg (NGR) sites can lead to the formation of isoAsp-Gly-Arg (isoDGR), a motif that can recognize the RGD-binding site of integrins. Ceruloplasmin (Cp), a ferroxidase present in the cerebrospinal fluid (CSF), contains two NGR sites in its sequence: one exposed on the protein surface ((568)NGR) and the other buried in the tertiary structure ((962)NGR). Considering that Cp can undergo oxidative modifications in the CSF of neurodegenerative diseases, we investigated the effect of oxidation on the deamidation of both NGR motifs and, consequently, on the acquisition of integrin binding properties. We observed that the exposed (568)NGR site can deamidate under conditions mimicking accelerated Asn aging. In contrast, the hidden (962)NGR site can deamidate exclusively when aging occurs under oxidative conditions, suggesting that oxidation-induced structural changes foster deamidation at this site. NGR deamidation in Cp was associated with gain of integrin-binding function, intracellular signaling, and cell pro-adhesive activity. Finally, Cp aging in the CSF from Alzheimer disease patients, but not in control CSF, causes Cp deamidation with gain of integrin-binding function, suggesting that this transition might also occur in pathological conditions. In conclusion, both Cp NGR sites can deamidate during aging under oxidative conditions, likely as a consequence of oxidative-induced structural changes, thereby promoting a gain of function in integrin binding, signaling, and cell adhesion.

Metadynamics Simulations Rationalise the Conformational Effects Induced by N-Methylation of RGD Cyclic Hexapeptides.

We combined metadynamics, docking and molecular mechanics/generalised born surface area (MM/GBSA) re-scoring methods to investigate the impact of single and multiple N-methylation on a set of RGD cyclopeptides displaying different affinity for integrin αIIbß3. We rationalised the conformational effects induced by N-methylation and its interplay with receptor affinity, obtaining good agreement with experimental data. This approach can be exploited before entering time-consuming and expensive synthesis and binding experiments.

Identification of a novel de novo deletion in RAF1 associated with biventricular hypertrophy in Noonan syndrome.

Biventricular hypertrophy (BVH) is a disease state characterized by the thickening of the ventricle walls. The differential diagnosis of BVH with other congenital and familial diseases in which increased ventricle wall thickness is a prominent clinical feature is fundamental due to its therapeutic and prognostic value, mainly during infancy. We describe a 2-month-old infant presenting BVH. Using exome sequencing, we identified a novel de novo 3-bp deletion in the RAF1 gene that is located in the binding active site for the 14-3-3 peptide. Based on docking calculations, we demonstrate that this novel mutation impairs protein/target binding, thus constitutively activating Ras signaling, which is a dysregulation associated with Noonan syndrome. Finally, our study underlines the importance of molecular modeling to understand the roles of novel mutations in pathogenesis.

NMR and computational methods in the structural and dynamic characterization of ligand-receptor interactions.

The recurrent failures in drug discovery campaigns, the asymmetry between the enormous financial investments and the relatively scarce results have fostered the development of strategies based on complementary methods. In this context in recent years the rigid lock-and-key binding concept had to be revisited in favour of a dynamic model of molecular recognition accounting for conformational changes of both the ligand and the receptor. The high level of complexity required by a dynamic description of the processes underlying molecular recognition requires a multidisciplinary investigation approach. In this perspective, the combination of nuclear magnetic resonance spectroscopy with molecular docking, conformational searches along with molecular dynamics simulations has given new insights into the dynamic mechanisms governing ligand receptor interactions, thus giving an enormous contribution to the identification and design of new and effective drugs. Herein a succinct overview on the applications of both NMR and computational methods to the structural and dynamic characterization of ligand-receptor interactions will be presented.

Single-cell atlas of ABCA7 loss-of-function reveals impaired neuronal respiration via choline-dependent lipid imbalances.

Loss-of-function (LoF) variants in the lipid transporter ABCA7 significantly increase the risk of Alzheimer's disease (odds ratio ∼2), yet the pathogenic mechanisms and the neural cell types affected by these variants remain largely unknown. Here, we performed single-nuclear RNA sequencing of 36 human post-mortem samples from the prefrontal cortex of 12 ABCA7 LoF carriers and 24 matched non-carrier control individuals. ABCA7 LoF was associated with gene expression changes in all major cell types. Excitatory neurons, which expressed the highest levels of ABCA7, showed transcriptional changes related to lipid metabolism, mitochondrial function, cell cycle-related pathways, and synaptic signaling. ABCA7 LoF-associated transcriptional changes in neurons were similarly perturbed in carriers of the common AD missense variant ABCA7 p.Ala1527Gly (n = 240 controls, 135 carriers), indicating that findings from our study may extend to large portions of the at-risk population. Consistent with ABCA7's function as a lipid exporter, lipidomic analysis of isogenic iPSC-derived neurons (iNs) revealed profound intracellular triglyceride accumulation in ABCA7 LoF, which was accompanied by a relative decrease in phosphatidylcholine abundance. Metabolomic and biochemical analyses of iNs further indicated that ABCA7 LoF was associated with disrupted mitochondrial bioenergetics that suggested impaired lipid breakdown by uncoupled respiration. Treatment of ABCA7 LoF iNs with CDP-choline (a rate-limiting precursor of phosphatidylcholine synthesis) reduced triglyceride accumulation and restored mitochondrial function, indicating that ABCA7 LoF-induced phosphatidylcholine dyshomeostasis may directly disrupt mitochondrial metabolism of lipids. Treatment with CDP-choline also rescued intracellular amyloid ß -42 levels in ABCA7 LoF iNs, further suggesting a link between ABCA7 LoF metabolic disruptions in neurons and AD pathology. This study provides a detailed transcriptomic atlas of ABCA7 LoF in the human brain and mechanistically links ABCA7 LoF-induced lipid perturbations to neuronal energy dyshomeostasis. In line with a growing body of evidence, our study highlights the central role of lipid metabolism in the etiology of Alzheimer's disease.

Genetic barriers more than environmental associations explain Serratia marcescens population structure.

Bacterial species often comprise well-separated lineages, likely emerged and maintained by genetic isolation and/or ecological divergence. How these two evolutionary actors interact in the shaping of bacterial population structure is currently not fully understood. In this study, we investigate the genetic and ecological drivers underlying the evolution of Serratia marcescens, an opportunistic pathogen with high genomic flexibility and able to colonise diverse environments. Comparative genomic analyses reveal a population structure composed of five deeply-demarcated genetic clusters with open pan-genome but limited inter-cluster gene flow, partially explained by Restriction-Modification (R-M) systems incompatibility. Furthermore, a large-scale research on hundred-thousands metagenomic datasets reveals only a partial habitat separation of the clusters. Globally, two clusters only show a separate gene composition coherent with ecological adaptations. These results suggest that genetic isolation has preceded ecological adaptations in the shaping of the species diversity, an evolutionary scenario coherent with the Evolutionary Extended Synthesis.

Advantages of long- and short-reads sequencing for the hybrid investigation of the Mycobacterium tuberculosis genome.

Introduction: In the fight to limit the global spread of antibiotic resistance, computational challenges associated with sequencing technology can impact the accuracy of downstream analysis, including drug resistance identification, transmission, and genome resolution. About 10% of Mycobacterium tuberculosis (MTB) genome is constituted by the PE/PPE family, a GC-rich repetitive genome region. Although sequencing using short read technology is widely used, it is well recognized its limit in the PE/PPE regions due to the unambiguously mapping process onto the reference genome. The aim of this study was to compare the performances of short-reads (SRS), long-reads (LRS) and hybrid-reads (HYBR) based analysis over different common investigative tasks: genome coverage estimation, variant calling and cluster analysis, drug resistance detection and de novo assembly. Methods: For the study 13 model MTB clinical isolates were sequenced with both SRS and LRS. HYBR were produced correcting the long reads with the short reads. The fastq from the three approaches were then processed using a customized version of MTBseq for genome coverage estimation and variant calling and using two different assemblers for de novo assembly evaluation. Results: Estimation of genome coverage performances showed lower 8X breadth coverage for SRS respect to LRS and HYBR: considering the PE/PPE genes, SRS showed low results for the PE_PGRS family, while obtained acceptable coverage in PE and PPE genes; LRS and HYBR reached optimal coverages in PE/PPE genes. For variant calling HYBR showed the highest resolution, detecting the highest percentage of uniquely identified mutations compared to LRS and SRS. All three approaches agreed on the identification of two major clusters, with HYBR identifying an higher number of SNPs between the two clusters. Comparing the quality of the assemblies, HYBR and LRS obtained better results than SRS. Discussion: In conclusion, depending on the aim of the investigation, both SRS and LRS present complementary advantages and limitations implying that for a full resolution of MTB genomes, where all the mentioned analyses and both technologies are needed, the use of the HYBR approach represents a valid option and a well-rounded strategy.

Comprehensive Drug Resistance Characterization of Pulmonary Tuberculosis in Algeria: Insights on Mycobacterium tuberculosis Strains by Whole-Genome Sequencing.

In this study, we aimed to characterize drug-resistant strains by whole-genome sequencing (WGS), to describe the spreading lineages and the history of transmission. Drug susceptibility testing was performed by 96-well broth microdilution plates. The genomic DNA was extracted and purified; libraries were prepared and run on the Illumina NextSeq500 System. Among 82 isolates, 21 tuberculosis (TB) isolates (25.6%) were drug resistant, including 10 MDR and 4 pre-extensively drug-resistant (XDR)-TB. The mutation Ser315Thr in the katG gene was confirmed in 15 isolates. In rpoB, Ser450Leu and His445Asp mutations were the most common. Asp94Asn and Ala90Val mutations were reported in gyrA. The LAM family, the most TB drug resistant, was widely predominant in the north and the T sublineage in the south of the country. This study provides the first insight on TB drug resistance using WGS in Algeria and clearly describes the first pre-XDR-TB cases and lineage distribution across the country.

Bedaquiline and clofazimine resistance in Mycobacterium tuberculosis: an in-vitro and in-silico data analysis.

BACKGROUND: Bedaquiline is a core drug for the treatment of multidrug-resistant tuberculosis; however, the understanding of resistance mechanisms is poor, which is hampering rapid molecular diagnostics. Some bedaquiline-resistant mutants are also cross-resistant to clofazimine. To decipher bedaquiline and clofazimine resistance determinants, we combined experimental evolution, protein modelling, genome sequencing, and phenotypic data. METHODS: For this in-vitro and in-silico data analysis, we used a novel in-vitro evolutionary model using subinhibitory drug concentrations to select bedaquiline-resistant and clofazimine-resistant mutants. We determined bedaquiline and clofazimine minimum inhibitory concentrations and did Illumina and PacBio sequencing to characterise selected mutants and establish a mutation catalogue. This catalogue also includes phenotypic and genotypic data of a global collection of more than 14 000 clinical Mycobacterium tuberculosis complex isolates, and publicly available data. We investigated variants implicated in bedaquiline resistance by protein modelling and dynamic simulations. FINDINGS: We discerned 265 genomic variants implicated in bedaquiline resistance, with 250 (94%) variants affecting the transcriptional repressor (Rv0678) of the MmpS5-MmpL5 efflux system. We identified 40 new variants in vitro, and a new bedaquiline resistance mechanism caused by a large-scale genomic rearrangement. Additionally, we identified in vitro 15 (7%) of 208 mutations found in clinical bedaquiline-resistant isolates. From our in-vitro work, we detected 14 (16%) of 88 mutations so far identified as being associated with clofazimine resistance and also seen in clinically resistant strains, and catalogued 35 new mutations. Structural modelling of Rv0678 showed four major mechanisms of bedaquiline resistance: impaired DNA binding, reduction in protein stability, disruption of protein dimerisation, and alteration in affinity for its fatty acid ligand. INTERPRETATION: Our findings advance the understanding of drug resistance mechanisms in M tuberculosis complex strains. We have established an extended mutation catalogue, comprising variants implicated in resistance and susceptibility to bedaquiline and clofazimine. Our data emphasise that genotypic testing can delineate clinical isolates with borderline phenotypes, which is essential for the design of effective treatments. FUNDING: Leibniz ScienceCampus Evolutionary Medicine of the Lung, Deutsche Forschungsgemeinschaft, Research Training Group 2501 TransEvo, Rhodes Trust, Stanford University Medical Scientist Training Program, National Institute for Health and Care Research Oxford Biomedical Research Centre, Oxford University Hospitals NHS Foundation Trust, Bill & Melinda Gates Foundation, Wellcome Trust, and Marie Sklodowska-Curie Actions.

Resistance patterns and transmission of mono- and polyresistant TB: clinical impact of WGS.

Objectives: Rapidly diagnosing drug-resistant TB is crucial for improving treatment and transmission control. WGS is becoming increasingly accessible and has added value to the diagnosis and treatment of TB. The aim of the study was to perform WGS to determine the rate of false-positive results of phenotypic drug susceptibility testing (pDST) and characterize the molecular mechanisms of resistance and transmission of mono- and polyresistant Mycobacterium (M.) tuberculosis. Methods: WGS was performed on 53 monoresistant and 25 polyresistant M. tuberculosis isolates characterized by pDST. Sequencing data were bioinformatically processed to infer mutations encoding resistance and determine the origin of resistance and phylogenetic relationship between isolates studied. Results: The data showed the variable sensitivity and specificity of WGS in comparison with pDST as the gold standard: isoniazid 92.7% and 92.3%; streptomycin 41.9% and 100.0%; pyrazinamide 15% and 94.8%; and ethambutol 75.0% and 98.6%, respectively. We found novel mutations encoding resistance to streptomycin (in gidB) and pyrazinamide (in kefB). Most isolates belonged to lineage 4 (80.1%) and the overall clustering rate was 11.5%. We observed lineage-specific gene variations encoding resistance to streptomycin and pyrazinamide. Conclusions: This study highlights the clinical potential of WGS in ruling out false-positive drug resistance following phenotypic or genetic drug testing, and recommend this technology together with the WHO catalogue in designing an optimal individualized treatment regimen and preventing the development of MDR TB. Our results suggest that resistance is primarily developed through spontaneous mutations or selective pressure.

Molecular dynamics reveal that isoDGR-containing cyclopeptides are true αvß3 antagonists unable to promote integrin allostery and activation.

Ain't got that swing(-out): The cyclopeptide isoDGR is emerging as a new αvß3 integrin binding motif. Agreement between the results of computational and biochemical studies reveals that isoDGR-containing cyclopeptides are true αvß3 integrin antagonists that block αvß3 in its inactive conformation (see scheme). isoDGR-based ligands may give αvß3 antagonists without paradoxical effects.

Whole genome sequencing of Bulgarian rifampicin resistant Mycobacterium tuberculosis strains.

INTRODUCTION: The transmission of drug-resistant tuberculosis is one of the greatest challenges facing the global tuberculosis control. AIM: The aim of the study was to investigate the resent transmission of rifampicin resistant tuberculosis in Bulgaria and to describe the mutations related to the antimicrobials' resistance using whole genome sequencing. MATERIALS AND METHODS: As part of an ECDC funded pilot study for evaluation of the systematic use of whole genome sequencing (WGS) of Mycobacterium tuberculosis (MTB) surveillance (EUSeqMyTB), Bulgaria provided 65 rifampicin resistant isolates over a three years' timeframe (2017-2019) representing 87.5% of the notified rifampicin resistant cases. Drug resistance prediction and relatedness analysis of the resistant isolates was performed in collaboration with San Raffaele Scientific Institute, Milan, Italy. RESULTS: Almost all of the isolates were identified as Euro-American lineage (96.9%); 18.5% of the isolates were found to be resistant to fluoroquinolones, but no mutations conferring resistance to bedaquiline or linezolid could be identified. Less than half (43.3%) of the isolates were clustered (<5 SNPs distance) into a total of seven national SNP-based clusters, while a total of six isolates were found to be part of different cross-border clusters. All clustered cases originated from Bulgaria. CONCLUSIONS: WGS has proven to be a reliable tool for surveillance and tracing of recent transmission of tuberculosis and has the potential for resistance prediction for most of the antituberculosis drugs.

Whole genome sequencing of multidrug-resistant Mycobacterium tuberculosis isolates collected in the Czech Republic, 2005-2020.

The emergence and spread of resistant tuberculosis (TB) pose a threat to public health, so it is necessary to diagnose the drug-resistant forms in a clinically short time frame and closely monitor their transmission. In this study, we carried out a first whole genome sequencing (WGS)-based analysis of multidrug resistant (MDR) M. tuberculosis strains to explore the phylogenetic lineages diversity, drug resistance mechanisms, and ongoing transmission chains within the country. In total, 65 isolates phenotypically resistant to at least rifampicin and isoniazid collected in the Czech Republic in 2005-2020 were enrolled for further analysis. The agreement of the results obtained by WGS with phenotypic drug susceptibility testing (pDST) in the determination of resistance to isoniazid, rifampicin, pyrazinamide, streptomycin, second-line injectables and fluoroquinolones was more than 80%. Phylogenetic analysis of WGS data revealed that the majority of MDR M. tuberculosis isolates were the Beijing lineage 2.2.1 (n = 46/65; 70.8%), while the remaining strains belonged to Euro-American lineage. Cluster analysis with a predefined cut-off distance of less than 12 single nucleotide polymorphisms between isolates showed 19 isolates in 6 clusters (clustering rate 29.2%), located mainly in the region of the capital city of Prague. This study highlights the utility of WGS as a high-resolution approach in the diagnosis, characterization of resistance patterns, and molecular-epidemiological analysis of resistant TB in the country.