High-throughput computational discovery of inhibitory protein fragments with AlphaFold.

Peptides can bind to specific sites on larger proteins and thereby function as inhibitors and regulatory elements. Peptide fragments of larger proteins are particularly attractive for achieving these functions due to their inherent potential to form native-like binding interactions. Recently developed experimental approaches allow for high-throughput measurement of protein fragment inhibitory activity in living cells. However, it has thus far not been possible to predict de novo which of the many possible protein fragments bind their protein targets, let alone act as inhibitors. We have developed a computational method, FragFold, that employs AlphaFold to predict protein fragment binding to full-length protein targets in a high-throughput manner. Applying FragFold to thousands of fragments tiling across diverse proteins revealed peaks of predicted binding along each protein sequence. These predictions were compared with experimentally measured peaks of inhibitory activity in E. coli. We establish that our approach is a sensitive predictor of protein fragment function: Evaluating inhibitory fragments derived from known protein-protein interaction interfaces, we found 87% were predicted by FragFold to bind in a native-like mode. Across full protein sequences, 68% of FragFold-predicted binding peaks match experimentally measured inhibitory peaks. This is true even when the underlying inhibitory mechanism is unclear from existing structural data, and we find FragFold is able to predict novel binding modes for inhibitory fragments of unknown structure, explaining previous genetic and biochemical data for these fragments. The success rate of FragFold demonstrates that this computational approach should be broadly applicable for discovering inhibitory protein fragments across proteomes.

Mapping functional regions of essential bacterial proteins with dominant-negative protein fragments.

Massively parallel measurements of dominant-negative inhibition by protein fragments have been used to map protein interaction sites and discover peptide inhibitors. However, the underlying principles governing fragment-based inhibition have thus far remained unclear. Here, we adapted a high-throughput inhibitory fragment assay for use in Escherichia coli, applying it to a set of 10 essential proteins. This approach yielded single amino acid resolution maps of inhibitory activity, with peaks localized to functionally important interaction sites, including oligomerization interfaces and folding contacts. Leveraging these data, we performed a systematic analysis to uncover principles of fragment-based inhibition. We determined a robust negative correlation between susceptibility to inhibition and cellular protein concentration, demonstrating that inhibitory fragments likely act primarily by titrating native protein interactions. We also characterized a series of trade-offs related to fragment length, showing that shorter peptides allow higher-resolution mapping but suffer from lower inhibitory activity. We employed an unsupervised statistical analysis to show that the inhibitory activities of protein fragments are largely driven not by generic properties such as charge, hydrophobicity, and secondary structure, but by the more specific characteristics of their bespoke macromolecular interactions. Overall, this work demonstrates fundamental characteristics of inhibitory protein fragment function and provides a foundation for understanding and controlling protein interactions in vivo.

Effects of sequence motifs in the yeast 3' untranslated region determined from massively parallel assays of random sequences.

BACKGROUND: The 3' untranslated region (UTR) plays critical roles in determining the level of gene expression through effects on activities such as mRNA stability and translation. Functional elements within this region have largely been identified through analyses of native genes, which contain multiple co-evolved sequence features. RESULTS: To explore the effects of 3' UTR sequence elements outside of native sequence contexts, we analyze hundreds of thousands of random 50-mers inserted into the 3' UTR of a reporter gene in the yeast Saccharomyces cerevisiae. We determine relative protein expression levels from the fitness of transformants in a growth selection. We find that the consensus 3' UTR efficiency element significantly boosts expression, independent of sequence context; on the other hand, the consensus positioning element has only a small effect on expression. Some sequence motifs that are binding sites for Puf proteins substantially increase expression in the library, despite these proteins generally being associated with post-transcriptional downregulation of native mRNAs. Our measurements also allow a systematic examination of the effects of point mutations within efficiency element motifs across diverse sequence backgrounds. These mutational scans reveal the relative in vivo importance of individual bases in the efficiency element, which likely reflects their roles in binding the Hrp1 protein involved in cleavage and polyadenylation. CONCLUSIONS: The regulatory effects of some 3' UTR sequence features, like the efficiency element, are consistent regardless of sequence context. In contrast, the consequences of other 3' UTR features appear to be strongly dependent on their evolved context within native genes.

Seeds of their own destruction: Dominant-negative peptide screening yields functional insight and therapeutic leads.

Systematic, high-throughput screening for "dominant-negative" protein fragments is an emerging method for mapping functional regions of the parental protein in vivo. In this issue of Cell Systems, Ford et al. apply this approach to 65 cancer drivers, providing functional insights and demonstrating therapeutic potential for several dominant-negative peptides.

Comprehensive sequence-to-function mapping of cofactor-dependent RNA catalysis in the glmS ribozyme.

Massively parallel, quantitative measurements of biomolecular activity across sequence space can greatly expand our understanding of RNA sequence-function relationships. We report the development of an RNA-array assay to perform such measurements and its application to a model RNA: the core glmS ribozyme riboswitch, which performs a ligand-dependent self-cleavage reaction. We measure the cleavage rates for all possible single and double mutants of this ribozyme across a series of ligand concentrations, determining kcat and KM values for active variants. These systematic measurements suggest that evolutionary conservation in the consensus sequence is driven by maintenance of the cleavage rate. Analysis of double-mutant rates and associated mutational interactions produces a structural and functional mapping of the ribozyme sequence, revealing the catalytic consequences of specific tertiary interactions, and allowing us to infer structural rearrangements that permit certain sequence variants to maintain activity.

Neutrophil-like HL-60 cells expressing only GFP-tagged ß-actin exhibit nearly normal motility.

Observations of actin dynamics in living cells using fluorescence microscopy have been foundational in the exploration of the mechanisms underlying cell migration. We used CRISPR/Cas9 gene editing to generate neutrophil-like HL-60 cell lines expressing GFP-ß-actin from the endogenous locus (ACTB). In light of many previous reports outlining functional deficiencies of labeled actin, we anticipated that HL-60 cells would only tolerate a monoallelic edit, as biallelic edited cells would produce no normal ß-actin. Surprisingly, we recovered viable monoallelic GFP-ß-actin cells as well as biallelic edited GFP-ß-actin cells, in which one copy of the ACTB gene is silenced and the other contains the GFP tag. Furthermore, the edited cells migrate with similar speeds and persistence as unmodified cells in a variety of motility assays, and have nearly normal cell shapes. These results might partially be explained by our observation that GFP-ß-actin incorporates into the F-actin network in biallelic edited cells at similar efficiencies as normal ß-actin in unedited cells. Additionally, the edited cells significantly upregulate γ-actin, perhaps helping to compensate for the loss of normal ß-actin. Interestingly, biallelic edited cells have only modest changes in global gene expression relative to the monoallelic line, as measured by RNA sequencing. While monoallelic edited cells downregulate expression of the tagged allele and are thus only weakly fluorescent, biallelic edited cells are quite bright and well-suited for live cell microscopy. The nondisruptive phenotype and direct interpretability of this fluorescent tagging approach make it a promising tool for studying actin dynamics in these rapidly migrating and highly phagocytic cells.

Self-cleavage of the glmS ribozyme core is controlled by a fragile folding element.

Riboswitches modulate gene expression in response to small-molecule ligands. Switching is generally thought to occur via the stabilization of a specific RNA structure conferred by binding the cognate ligand. However, it is unclear whether any such stabilization occurs for riboswitches whose ligands also play functional roles, such as the glmS ribozyme riboswitch, which undergoes self-cleavage using its regulatory ligand, glucosamine 6-phosphate, as a catalytic cofactor. To address this question, it is necessary to determine both the conformational ensemble and its ligand dependence. We used optical tweezers to measure folding dynamics and cleavage rates for the core glmS ribozyme over a range of forces and ligand conditions. We found that the folding of a specific structural element, the P2.2 duplex, controls active-site formation and catalysis. However, the folded state is only weakly stable, regardless of cofactor concentration, supplying a clear exception to the ligand-based stabilization model of riboswitch function.

Single-molecule studies of riboswitch folding.

The folding dynamics of riboswitches are central to their ability to modulate gene expression in response to environmental cues. In most cases, a structural competition between the formation of a ligand-binding aptamer and an expression platform (or some other competing off-state) determines the regulatory outcome. Here, we review single-molecule studies of riboswitch folding and function, predominantly carried out using single-molecule FRET or optical trapping approaches. Recent results have supplied new insights into riboswitch folding energy landscapes, the mechanisms of ligand binding, the roles played by divalent ions, the applicability of hierarchical folding models, and kinetic vs. thermodynamic control schemes. We anticipate that future work, based on improved data sets and potentially combining multiple experimental techniques, will enable the development of more complete models for complex RNA folding processes. This article is part of a Special Issue entitled: Riboswitches.

The Bxb1 gp47 recombination directionality factor is required not only for prophage excision, but also for phage DNA replication.

Mycobacteriophage Bxb1 encodes a serine-integrase that catalyzes both integrative and excisive site-specific recombination. However, excision requires a second phage-encoded protein, gp47, which serves as a recombination directionality factor (RDF). The viability of a Bxb1 mutant containing an S153A substitution in gp47 that eliminates the RDF activity of Bxb1 gp47 shows that excision is not required for Bxb1 lytic growth. However, the inability to construct a Δ47 deletion mutant of Bxb1 suggests that gp47 provides a second function that is required for lytic growth, although the possibility of an essential cis-acting site cannot be excluded. Characterization of a mutant prophage of mycobacteriophage L5 in which gene 54 - a homologue of Bxb1 gene 47 - is deleted shows that it also is defective in induced lytic growth, and exhibits a strong defect in DNA replication. Bxb1 gp47 and its relatives are also unusual in containing conserved motifs associated with a phosphoesterase function, although we have not been able to show robust phosphoesterase activity of the proteins, and amino acid substitutions with the conserved motifs do not interfere with RDF activity. We therefore propose that Bxb1 gp47 and its relatives provide an important function in phage DNA replication that has been co-opted by the integration machinery of the serine-integrases to control the directionality of recombination.