Genetic interaction mapping informs integrative structure determination of protein complexes

Determining structures of protein complexes is crucial for understanding cellular functions. Here, we describe an integrative structure determination approach that relies on in vivo measurements of genetic interactions. We construct phenotypic profiles for point mutations crossed against gene deletions or exposed to environmental perturbations, followed by converting similarities between two profiles into an upper bound on the distance between the mutated residues. We determine the structure of the yeast histone H3-H4 complex based on similar to 500,000 genetic interactions of 350 mutants. We then apply the method to subunits Rpb1-Rpb2 of yeast RNA polymerase II and subunits RpoB-RpoC of bacterial RNA polymerase. The accuracy is comparable to that based on chemical cross-links; using restraints from both genetic interactions and cross-links further improves model accuracy and precision. The approach provides an efficient means to augment integrative structure determination with in vivo observations

Unveiling a novel transient druggable pocket in BACE-1 through molecular simulations: conformational analysis and binding mode of multisite inhibitors

The critical role of BACE-1 in the formation of neurotoxic ß-amyloid peptides in the brain makes it an attractive target for an efficacious treatment of Alzheimer’s disease. However, the development of clinically useful BACE-1 inhibitors has proven to be extremely challeng- ing. In this study we examine the binding mode of a novel potent inhibitor (compound 1, with IC50 80 nM) designed by synergistic combination of two fragments—huprine and rhein— that individually are endowed with very low activity against BACE-1. Examination of crystal structures reveals no appropriate binding site large enough to accommodate 1. Therefore we have examined the conformational flexibility of BACE-1 through extended molecular dynamics simulations, paying attention to the highly flexible region shaped by loops 8–14, 154–169 and 307–318. The analysis of the protein dynamics, together with studies of pocket druggability, has allowed us to detect the transient formation of a secondary binding site, which contains Arg307 as a key residue for the interaction with small molecules, at the edge of the catalytic cleft. The formation of this druggable “floppy” pocket would enable the bind- ing of multisite inhibitors targeting both catalytic and secondary sites. Molecular dynamics simulations of BACE-1 bound to huprine-rhein hybrid compounds support the feasibility of this hypothesis. The results provide a basis to explain the high inhibitory potency of the two enantiomeric forms of 1, together with the large dependence on the length of the oligo- methylenic linker. Furthermore, the multisite hypothesis has allowed us to rationalize the inhibitory potency of a series of tacrine-chromene hybrid compounds, specifically regarding the apparent lack of sensitivity of the inhibition constant to the chemical modifications intro- duced in the chromene unit. Overall, these findings pave the way for the exploration of novel functionalities in the design of optimized BACE-1 multisite inhibitors

In silico exploration of Red Sea Bacillus genomes for natural product biosynthetic gene clusters

Background: The increasing spectrum of multidrug-resistant bacteria is a major global public health concern, necessitating discovery of novel antimicrobial agents. Here, members of the genus Bacillus are investigated as a potentially attractive source of novel antibiotics due to their broad spectrum of antimicrobial activities. We specifically focus on a computational analysis of the distinctive biosynthetic potential of Bacillus paralicheniformis strains isolated from the Red Sea, an ecosystem exposed to adverse, highly saline and hot conditions. Results: We report the complete circular and annotated genomes of two Red Sea strains, B. paralicheniformis Bac48 isolated from mangrove mud and B. paralicheniformis Bac84 isolated from microbial mat collected from Rabigh Harbor Lagoon in Saudi Arabia. Comparing the genomes of B. paralicheniformis Bac48 and B. paralicheniformis Bac84 with nine publicly available complete genomes of B. licheniformis and three genomes of B. paralicheniformis, revealed that all of the B. paralicheniformis strains in this study are more enriched in nonribosomal peptides (NRPs). We further report the first computationally identified trans-acyltransferase (trans-AT) nonribosomal peptide synthetase/polyketide synthase (PKS/ NRPS) cluster in strains of this species. Conclusions:B. paralicheniformis species have more genes associated with biosynthesis of antimicrobial bioactive compounds than other previously characterized species of B. licheniformis, which suggests that these species are better potential sources for novel antibiotics. Moreover, the genome of the Red Sea strain B. paralicheniformis Bac48 is more enriched in modular PKS genes compared to B. licheniformis strains and other B. paralicheniformis strains. This may be linked to adaptations that strains surviving in the Red Sea underwent to survive in the relatively hot and saline ecosystems

Synthetic biology to access and expand nature's chemical diversity

Bacterial genomes encode the biosynthetic potential to produce hundreds of thousands of complex molecules with diverse applications, from medicine to agriculture and materials. Accessing these natural products promises to reinvigorate drug discovery pipelines and provide novel routes to synthesize complex chemicals. The pathways leading to the production of these molecules often comprise dozens of genes spanning large areas of the genome and are controlled by complex regulatory networks with some of the most interesting molecules being produced by non-model organisms. In this Review, we discuss how advances in synthetic biology — including novel DNA construction technologies, the use of genetic parts for the precise control of expression and for synthetic regulatory circuits — and multiplexed genome engineering can be used to optimize the design and synthesis of pathways that produce natural products

Comparative genomics reveals phylogenetic distribution patterns of secondary metabolites in Amycolatopsis species

Background Genome mining tools have enabled us to predict biosynthetic gene clusters that might encode compounds with valuable functions for industrial and medical applications. With the continuously increasing number of genomes sequenced, we are confronted with an overwhelming number of predicted clusters. In order to guide the effective prioritization of biosynthetic gene clusters towards finding the most promising compounds, knowledge about diversity, phylogenetic relationships and distribution patterns of biosynthetic gene clusters is necessary. Results Here, we provide a comprehensive analysis of the model actinobacterial genus Amycolatopsis and its potential for the production of secondary metabolites. A phylogenetic characterization, together with a pan-genome analysis showed that within this highly diverse genus, four major lineages could be distinguished which differed in their potential to produce secondary metabolites. Furthermore, we were able to distinguish gene cluster families whose distribution correlated with phylogeny, indicating that vertical gene transfer plays a major role in the evolution of secondary metabolite gene clusters. Still, the vast majority of the diverse biosynthetic gene clusters were derived from clusters unique to the genus, and also unique in comparison to a database of known compounds. Our study on the locations of biosynthetic gene clusters in the genomes of Amycolatopsis’ strains showed that clusters acquired by horizontal gene transfer tend to be incorporated into non-conserved regions of the genome thereby allowing us to distinguish core and hypervariable regions in Amycolatopsis genomes. Conclusions Using a comparative genomics approach, it was possible to determine the potential of the genus Amycolatopsis to produce a huge diversity of secondary metabolites. Furthermore, the analysis demonstrates that horizontal and vertical gene transfer play an important role in the acquisition and maintenance of valuable secondary metabolites. Our results cast light on the interconnections between secondary metabolite gene clusters and provide a way to prioritize biosynthetic pathways in the search and discovery of novel compounds

High-quality MS/MS spectrum prediction for data-dependent and data-independent acquisition data analysis

Peptide fragmentation spectra are routinely predicted in the interpretation of mass-spectrometry-based proteomics data. However, the generation of fragment ions has not been understood well enough for scientists to estimate fragment ion intensities accurately. Here, we demonstrate that machine learning can predict peptide fragmentation patterns in mass spectrometers with accuracy within the uncertainty of measurement. Moreover, analysis of our models reveals that peptide fragmentation depends on long-range interactions within a peptide sequence. We illustrate the utility of our models by applying them to the analysis of both data-dependent and data-independent acquisition datasets. In the former case, we observe a q-value-dependent increase in the total number of peptide identifications. In the latter case, we confirm that the use of predicted tandem mass spectrometry spectra is nearly equivalent to the use of spectra from experimental libraries

Modeling of the TFIIH complex using chemical cross-links and electron microscopy (EM) density maps

TFIIH is essential for both RNA polymerase II transcription and DNA repair, and mutations in TFIIH can result in human disease. Here, we determine the molecular architecture of human and yeast TFIIH by an integrative approach using chemical crosslinking/mass spectrometry (CXMS) data, biochemical analyses, and previously published electron microscopy maps. We identified four new conserved "topological regions" that function as hubs for TFIIH assembly and more than 35 conserved topological features within TFIIH, illuminating a network of interactions involved in TFIIH assembly and regulation of its activities. We show that one of these conserved regions, the p62/Tfb1 Anchor region, directly interacts with the DNA helicase subunit XPD/Rad3 in native TFIIH and is required for the integrity and function of TFIIH. We also reveal the structural basis for defects in patients with xeroderma pigmentosum and trichothiodystrophy, with mutations found at the interface between the p62 Anchor region and the XPD subunit. For more information about how to reproduce this modeling, see https://salilab.org/tfiih or the README file

Molecular Details Underlying Dynamic Structures and Regulation of the Human 26S Proteasome

The 26S proteasome is the macromolecular machine responsible for ATP/ubiquitin dependent degradation. As aberration in proteasomal degradation has been implicated in many human diseases, structural analysis of the human 26S proteasome complex is essential to advance our understanding of its action and regulation mechanisms. In recent years, cross-linking mass spectrometry (XL-MS) has emerged as a powerful tool for elucidating structural topologies of large protein assemblies, with its unique capability of studying protein complexes in cells. To facilitate the identification of cross-linked peptides, we have previously developed a robust amine reactive sulfoxide-containing MS-cleavable cross-linker, disuccinimidyl sulfoxide (DSSO). To better understand the structure and regulation of the human 26S proteasome, we have established new DSSO-based in vivo and in vitro XL-MS workflows by coupling with HB-tag based affinity purification to comprehensively examine protein-protein interactions within the 26S proteasome. In total, we have identified 447 unique lysine-to-lysine linkages delineating 67 inter-protein and 26 intra-protein interactions, representing the largest cross-link dataset for proteasome complexes. In combination with EM maps and computational modeling, the architecture of the 26S proteasome was determined to infer its structural dynamics. In particular, three proteasome subunits Rpn1, Rpn6 and Rpt6 displayed multiple conformations that have not been previously reported. Additionally, cross-links between proteasome subunits and 15 proteasome interacting proteins including 9 known and 6 novel ones have been determined to demonstrate their physical interactions at the amino-acid level. Our results have provided new insights on the dynamics of the 26S human proteasome and the methodologies presented here can be applied to study other protein complexes. For more information about how to reproduce this modeling, see https://salilab.org/26S-PIPs or the README file