Benchmarking Refined and Unrefined AlphaFold2 Structures for Hit Discovery.

The recently developed AlphaFold2 (AF2) algorithm predicts proteins' 3D structures from amino acid sequences. The open AlphaFold protein structure database covers the complete human proteome. Using an industry-leading molecular docking method (Glide), we investigated the virtual screening performance of 37 common drug targets, each with an AF2 structure and known holo and apo structures from the DUD-E data set. In a subset of 27 targets where the AF2 structures are suitable for refinement, the AF2 structures show comparable early enrichment of known active compounds (avg. EF 1%: 13.0) to apo structures (avg. EF 1%: 11.4) while falling behind early enrichment of the holo structures (avg. EF 1%: 24.2). With an induced-fit protocol (IFD-MD), we can refine the AF2 structures using an aligned known binding ligand as the template to improve the performance in structure-based virtual screening (avg. EF 1%: 18.9). Glide-generated docking poses of known binding ligands can also be used as templates for IFD-MD, achieving similar improvements (avg. EF 1% 18.0). Thus, with proper preparation and refinement, AF2 structures show considerable promise for in silico hit identification.

The YΦ motif defines the structure-activity relationships of human 20S proteasome activators.

The 20S proteasome (20S) facilitates turnover of most eukaryotic proteins. Substrate entry into the 20S first requires opening of gating loops through binding of HbYX motifs that are present at the C-termini of certain proteasome activators (PAs). The HbYX motif has been predominantly characterized in the archaeal 20S, whereas little is known about the sequence preferences of the human 20S (h20S). Here, we synthesize and screen ~120 HbYX-like peptides, revealing unexpected differences from the archaeal system and defining the h20S recognition sequence as the Y-F/Y (YФ) motif. To gain further insight, we create a functional chimera of the optimized sequence, NLSYYT, fused to the model activator, PA26E102A. A cryo-EM structure of PA26E102A-h20S is used to identify key interactions, including non-canonical contacts and gate-opening mechanisms. Finally, we demonstrate that the YФ sequence preferences are tuned by valency, allowing multivalent PAs to sample greater sequence space. These results expand the model for termini-mediated gating and provide a template for the design of h20S activators.

Toward Lung Tumor Localization based on Strain Variability of Lung Surface during Video-Assisted Thoracoscopic Surgery.

Tumor localization, especially in case of minimally invasive lung tumor resection surgery, is extremely challenging due to the continuous motion of the organ. This motion can be troublesome as it results in spatial discrepancy corresponding to preoperative and intraoperative tumor location. In order to characterize lung tissue stiffness for the purpose of lung tumor localization, in this paper, we present a novel characterization approach based on variability in resistance of the healthy region vs. the tumorous region resulting from lung motion. The proposed approach is numerically validated on a Finite Element (FE) model of the lung with varying surface stiffnesses, where higher stiffness represents tumor and lower stiffness corresponds to healthy lung tissue. The numerical simulation validates the sensitivity of our mechanism for different grades of tumors by demonstrating that the strain on the healthy tissue is 31.8 and 67.1 times higher than that on the tumor surface for a selected relative stiffness variation of 3.6x and 24.4x respectively, at a pressure of 1.6 KPa. Additionally, a framework is developed to validate the proposed approach in a video of a video-assisted thoracoscopic surgery (VATS), where multiple landmarks on the lung surface are tracked. This enables us to quantify the motion of points residing on healthy surface and tumorous surface. The motion data is further analyzed to study the relative surface strain, and it is shown that the proposed approach differentiates a tumor from healthy surface.

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 interprotein and 26 intraprotein 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.

Dynamic Allostery Mediated by a Conserved Tryptophan in the Tec Family Kinases.

Bruton's tyrosine kinase (Btk) is a Tec family non-receptor tyrosine kinase that plays a critical role in immune signaling and is associated with the immunological disorder X-linked agammaglobulinemia (XLA). Our previous findings showed that the Tec kinases are allosterically activated by the adjacent N-terminal linker. A single tryptophan residue in the N-terminal 17-residue linker mediates allosteric activation, and its mutation to alanine leads to the complete loss of activity. Guided by hydrogen/deuterium exchange mass spectrometry results, we have employed Molecular Dynamics simulations, Principal Component Analysis, Community Analysis and measures of node centrality to understand the details of how a single tryptophan mediates allostery in Btk. A specific tryptophan side chain rotamer promotes the functional dynamic allostery by inducing coordinated motions that spread across the kinase domain. Either a shift in the rotamer population, or a loss of the tryptophan side chain by mutation, drastically changes the coordinated motions and dynamically isolates catalytically important regions of the kinase domain. This work also identifies a new set of residues in the Btk kinase domain with high node centrality values indicating their importance in transmission of dynamics essential for kinase activation. Structurally, these node residues appear in both lobes of the kinase domain. In the N-lobe, high centrality residues wrap around the ATP binding pocket connecting previously described Catalytic-spine residues. In the C-lobe, two high centrality node residues connect the base of the R- and C-spines on the αF-helix. We suggest that the bridging residues that connect the catalytic and regulatory architecture within the kinase domain may be a crucial element in transmitting information about regulatory spine assembly to the catalytic machinery of the catalytic spine and active site.

A conserved isoleucine maintains the inactive state of Bruton's tyrosine kinase.

Despite high level of homology among non-receptor tyrosine kinases, different kinase families employ a diverse array of regulatory mechanisms. For example, the catalytic kinase domains of the Tec family kinases are inactive without assembly of the adjacent regulatory domains, whereas the Src kinase domains are autoinhibited by the assembly of similar adjacent regulatory domains. Using molecular dynamics simulations, biochemical assays, and biophysical approaches, we have uncovered an isoleucine residue in the kinase domain of the Tec family member Btk that, when mutated to the closely related leucine, leads to a shift in the conformational equilibrium of the kinase domain toward the active state. The single amino acid mutation results in measureable catalytic activity for the Btk kinase domain in the absence of the regulatory domains. We suggest that this isoleucine side chain in the Tec family kinases acts as a "wedge" that restricts the conformational space available to key regions in the kinase domain, preventing activation until the kinase domain associates with its regulatory subunits and overcomes the energetic barrier to activation imposed by the isoleucine side chain.

Removal of acutely hazardous pharmaceuticals from water using multi-template imprinted polymer adsorbent.

Molecularly imprinted polymer adsorbent has been prepared to remove a group of recalcitrant and acutely hazardous (p-type) chemicals from water and wastewaters. The polymer adsorbent exhibited twofold higher adsorption capacity than the commercially used polystyrene divinylbenzene resin (XAD) and powdered activated carbon adsorbents. Higher adsorption capacity of the polymer adsorbent was explained on the basis of high specific surface area formed during molecular imprinting process. Freundlich isotherms drawn showed that the adsorption of p-type chemicals onto polymer adsorbent was kinetically faster than the other reference adsorbents. Matrix effect on adsorption of p-type chemicals was minimal, and also polymer adsorbent was amenable to regeneration by washing with water/methanol (3:1, v/v) solution. The polymer adsorbent was unaltered in its adsorption capacity up to 10 cycles of adsorption and desorption, which will be more desirable in cost reduction of treatment compared with single-time-use activated carbon.

Linkage, mobility, and selfishness in the MazF family of bacterial toxins: a snapshot of bacterial evolution.

Prokaryotic MazF family toxins cooccur with cognate antitoxins having divergent DNA-binding folds and can be of chromosomal or plasmid origin. Sequence similarity search was carried out to identify the Toxin-Antitoxin (TA) operons of MazF family followed by sequence analysis and phylogenetic studies. The genomic DNA upstream of the TA operons was searched for the presence of regulatory motifs. The MazF family toxins showed a conserved hydrophobic pocket in a multibinding site and are present in pathogenic bacteria. The toxins of the MazF family are associated with four main types of cognate antitoxin partners and cluster as a subfamily on the branches of the phylogenetic tree. This indicates that transmission of the entire operon is the dominant mode of inheritance. The plasmid borne TA modules were interspersed between the chromosomal TA modules of the same subfamily, compatible with a frequent interchange of TA genes between the chromosome and the plasmid akin to that observed for antibiotic resistance gens. The split network of the MazF family toxins showed the AbrB-linked toxins as a hub of horizontal gene transfer. Distinct motifs are present in the upstream region of each subfamily. The presence of MazF family TA modules in pathogenic bacteria and identification of a conserved binding pocket are significant for the development of novel antibacterials to disrupt the TA interaction. However, the role of TAs in stress resistance needs to be established. Phylogenetic studies provide insight into the evolution of MazF family TAs and effect on the bacterial genome.

Modeling of the structure and interactions of the B. anthracis antitoxin, MoxX: deletion mutant studies highlight its modular structure and repressor function.

Our previous report on Bacillus anthracis toxin-antitoxin module (MoxXT) identified it to be a two component system wherein, PemK-like toxin (MoxT) functions as a ribonuclease (Agarwal S et al. JBC 285:7254-7270, 2010). The labile antitoxin (MoxX) can bind to/neutralize the action of the toxin and is also a DNA-binding protein mediating autoregulation. In this study, molecular modeling of MoxX in its biologically active dimeric form was done. It was found that it contains a conserved Ribbon-Helix-Helix (RHH) motif, consistent with its DNA-binding function. The modeled MoxX monomers dimerize to form a two-stranded antiparallel ribbon, while the C-terminal region adopts an extended conformation. Knowledge guided protein-protein docking, molecular dynamics simulation, and energy minimization was performed to obtain the structure of the MoxXT complex, which was exploited for the de novo design of a peptide capable of binding to MoxT. It was found that the designed peptide caused a decrease in MoxX binding to MoxT by 42% at a concentration of 2 µM in vitro. We also show that MoxX mediates negative transcriptional autoregulation by binding to its own upstream DNA. The interacting regions of both MoxX and DNA were identified in order to model their complex. The repressor activity of MoxX was found to be mediated by the 16 N-terminal residues that contains the ribbon of the RHH motif. Based on homology with other RHH proteins and deletion mutant studies, we propose a model of the MoxX-DNA interaction, with the antiparallel ß-sheet of the MoxX dimer inserted into the major groove of its cognate DNA. The structure of the complex of MoxX with MoxT and its own upstream regulatory region will facilitate design of molecules that can disrupt these interactions, a strategy for development of novel antibacterials.