The khmer package is a freely available software library for working efficiently with fixed length DNA words, or k-mers. khmer provides implementations of a probabilistic k-mer counting data structure, a compressible De Bruijn graph representation, De Bruijn graph partitioning, and digital normalization. khmer is implemented in C++ and Python, and is freely available under the BSD license at  https://github.com/dib-lab/khmer/.

In pretreatment isolates of M. tuberculosis with decrements of MIC values of isoniazid or rifampin below standard resistance breakpoints, higher MIC values were associated with a greater risk of relapse than lower MIC values. (Funded by the National Institute of Allergy and Infectious Diseases.).

Cell wall biosynthesis inhibitors have proven highly effective for treating tuberculosis (TB). We discovered and validated members of the indazole sulfonamide class of small molecules as inhibitors of Mycobacterium tuberculosis KasA—a key component for biosynthesis of the mycolic acid layer of the bacterium’s cell wall and the same pathway as that inhibited by the first-line antitubercular drug isoniazid (INH). One lead compound, DG167, demonstrated synergistic lethality in combination with INH and a transcriptional pattern consistent with bactericidality and loss of persisters. Our results also detail a novel dual-binding mechanism for this compound as well as substantial structure-activity relationships (SAR) that may help in lead optimization activities. Together, these results suggest that KasA inhibition, specifically, that shown by the DG167 series, may be developed into a potent therapy that can synergize with existing antituberculars.

Epistatic interactions between residues determine a protein’s adaptability and shape its evolutionary trajectory. When a protein experiences a changed environment, it is under strong selection to find a peak in the new fitness landscape. It has been shown that strong selection increases epistatic interactions as well as the ruggedness of the fitness landscape, but little is known about how the epistatic interactions change under selection in the long-term evolution of a protein. Here we analyze the evolution of epistasis in the protease of the human immunodeficiency virus type 1 (HIV-1) using protease sequences collected for almost a decade from both treated and untreated patients, to understand how epistasis changes and how those changes impact the long-term evolvability of a protein. We use an information-theoretic proxy for epistasis that quantifies the co-variation between sites, and show that positive information is a necessary (but not sufficient) condition that detects epistasis in most cases. We analyze the “fossils” of the evolutionary trajectories of the protein contained in the sequence data, and show that epistasis continues to enrich under strong selection, but not for proteins whose environment is unchanged. The increase in epistasis compensates for the information loss due to sequence variability brought about by treatment, and facilitates adaptation in the increasingly rugged fitness landscape of treatment. While epistasis is thought to enhance evolvability via valley-crossing early-on in adaptation, it can hinder adaptation later when the landscape has turned rugged. However, we find no evidence that the HIV-1 protease has reached its potential for evolution after 9 years of adapting to a drug environment that itself is constantly changing. We suggest that the mechanism of encoding new information into pairwise interactions is central to protein evolution not just in HIV-1 protease, but for any protein adapting to a changing environment.