A Computational Pipeline for Accurate Prioritization of Protein-Protein Binding Candidates in High-Throughput Protein Libraries.

Identifying the interactome for a protein of interest is challenging due to the large number of possible binders. High-throughput experimental approaches narrow down possible binding partners but often include false positives. Furthermore, they provide no information about what the binding region is (e.g., the binding epitope). We introduce a novel computational pipeline based on an AlphaFold2 (AF) Competitive Binding Assay (AF-CBA) to identify proteins that bind a target of interest from a pull-down experiment and the binding epitope. Our focus is on proteins that bind the Extraterminal (ET) domain of Bromo and Extraterminal domain (BET) proteins, but we also introduce nine additional systems to show transferability to other peptide-protein systems. We describe a series of limitations to the methodology based on intrinsic deficiencies of AF and AF-CBA to help users identify scenarios where the approach will be most useful. Given the method's speed and accuracy, we anticipate its broad applicability to identify binding epitope regions among potential partners, setting the stage for experimental verification.

Structure Determination of Challenging Protein-Peptide Complexes Combining NMR Chemical Shift Data and Molecular Dynamics Simulations.

Intrinsically disordered regions of proteins often mediate important protein-protein interactions. However, the folding-upon-binding nature of many polypeptide-protein interactions limits the ability of modeling tools to predict the three-dimensional structures of such complexes. To address this problem, we have taken a tandem approach combining NMR chemical shift data and molecular simulations to determine the structures of peptide-protein complexes. Here, we use the MELD (Modeling Employing Limited Data) technique applied to polypeptide complexes formed with the extraterminal domain (ET) of bromo and extraterminal domain (BET) proteins, which exhibit a high degree of binding plasticity. This system is particularly challenging as the binding process includes allosteric changes across the ET receptor upon binding, and the polypeptide binding partners can adopt different conformations (e.g., helices and hairpins) in the complex. In a blind study, the new approach successfully modeled bound-state conformations and binding poses, using only protein receptor backbone chemical shift data, in excellent agreement with experimentally determined structures for moderately tight (Kd ∼100 nM) binders. The hybrid MELD + NMR approach required additional peptide ligand chemical shift data for weaker (Kd ∼250 µM) peptide binding partners. AlphaFold also successfully predicts the structures of some of these peptide-protein complexes. However, whereas AlphaFold can provide qualitative peptide rankings, MELD can directly estimate relative binding affinities. The hybrid MELD + NMR approach offers a powerful new tool for structural analysis of protein-polypeptide complexes involving disorder-to-order transitions upon complex formation, which are not successfully modeled with most other complex prediction methods, providing both the 3D structures of peptide-protein complexes and their relative binding affinities.

Disrupting MLV integrase:BET protein interaction biases integration into quiescent chromatin and delays but does not eliminate tumor activation in a MYC/Runx2 mouse model.

Murine leukemia virus (MLV) integrase (IN) lacking the C-terminal tail peptide (TP) loses its interaction with the host bromodomain and extraterminal (BET) proteins and displays decreased integration at promoter/enhancers and transcriptional start sites/CpG islands. MLV lacking the IN TP via an altered open reading frame was used to infect tumorigenesis mouse model (MYC/Runx2) animals to observe integration patterns and phenotypic effects, but viral passage resulted in the restoration of the IN TP through small deletions. Mice subsequently infected with an MLV IN lacking the TP coding sequence (TP-) showed an improved median survival by 15 days compared to wild type (WT) MLV infection. Recombination with polytropic endogenous retrovirus (ERV), Pmv20, was identified in seven mice displaying both fast and slow tumorigenesis, highlighting the strong selection within the mouse to maintain the full-length IN protein. Mapping the genomic locations of MLV in tumors from an infected mouse with no observed recombination with ERVs, TP-16, showed fewer integrations at TSS and CpG islands, compared to integrations observed in WT tumors. However, this mouse succumbed to the tumor in relatively rapid fashion (34 days). Analysis of the top copy number integrants in the TP-16 tumor revealed their proximity to known MLV common insertion site genes while maintaining the MLV IN TP- genotype. Furthermore, integration mapping in K562 cells revealed an insertion preference of MLV IN TP- within chromatin profile states associated with weakly transcribed heterochromatin with fewer integrations at histone marks associated with BET proteins (H3K4me1/2/3, and H3K27Ac). While MLV IN TP- showed a decreased overall rate of tumorigenesis compared to WT virus in the MYC/Runx2 model, MLV integration still occurred at regions associated with oncogenic driver genes independently from the influence of BET proteins, either stochastically or through trans-complementation by functional endogenous Gag-Pol protein.

X-ray crystal structure of the N-terminal region of Moloney murine leukemia virus integrase and its implications for viral DNA recognition.

The retroviral integrase (IN) carries out the integration of a dsDNA copy of the viral genome into the host DNA, an essential step for viral replication. All IN proteins have three general domains, the N-terminal domain (NTD), the catalytic core domain, and the C-terminal domain. The NTD includes an HHCC zinc finger-like motif, which is conserved in all retroviral IN proteins. Two crystal structures of Moloney murine leukemia virus (M-MuLV) IN N-terminal region (NTR) constructs that both include an N-terminal extension domain (NED, residues 1-44) and an HHCC zinc-finger NTD (residues 45-105), in two crystal forms are reported. The structures of IN NTR constructs encoding residues 1-105 (NTR1-105 ) and 8-105 (NTR8-105 ) were determined at 2.7 and 2.15 Å resolution, respectively and belong to different space groups. While both crystal forms have similar protomer structures, NTR1-105 packs as a dimer and NTR8-105 packs as a tetramer in the asymmetric unit. The structure of the NED consists of three anti-parallel ß-strands and an α-helix, similar to the NED of prototype foamy virus (PFV) IN. These three ß-strands form an extended ß-sheet with another ß-strand in the HHCC Zn2+ binding domain, which is a unique structural feature for the M-MuLV IN. The HHCC Zn2+ binding domain structure is similar to that in HIV and PFV INs, with variations within the loop regions. Differences between the PFV and MLV IN NEDs localize at regions identified to interact with the PFV LTR and are compared with established biochemical and virological data for M-MuLV. Proteins 2017; 85:647-656. © 2016 Wiley Periodicals, Inc.

Repression of the Chromatin-Tethering Domain of Murine Leukemia Virus p12.

Murine leukemia virus (MLV) p12, encoded within Gag, binds the viral preintegration complex (PIC) to the mitotic chromatin. This acts to anchor the viral PIC in the nucleus as the nuclear envelope re-forms postmitosis. Mutations within the p12 C terminus (p12 PM13 to PM15) block early stages in viral replication. Within the p12 PM13 region (p12 60PSPMA65), our studies indicated that chromatin tethering was not detected when the wild-type (WT) p12 protein (M63) was expressed as a green fluorescent protein (GFP) fusion; however, constructs bearing p12-I63 were tethered. N-terminal truncations of the activated p12-I63-GFP indicated that tethering increased further upon deletion of p12 25DLLTEDPPPY34, which includes the late domain required for viral assembly. The p12 PM15 sequence (p12 70RREPP74) is critical for wild-type viral viability; however, virions bearing the PM15 mutation (p12 70AAAAA74) with a second M63I mutant were viable, with a titer 18-fold lower than that of the WT. The p12 M63I mutation amplified chromatin tethering and compensated for the loss of chromatin binding of p12 PM15. Rescue of the p12-M63-PM15 nonviable mutant with prototype foamy virus (PFV) and Kaposi's sarcoma herpesvirus (KSHV) tethering sequences confirmed the function of p1270-74 in chromatin binding. Minimally, full-strength tethering was seen with only p12 61SPIASRLRGRR71 fused to GFP. These results indicate that the p12 C terminus alone is sufficient for chromatin binding and that the presence of the p12 25DLLTEDPPPY34 motif in the N terminus suppresses the ability to tether. IMPORTANCE: This study defines a regulatory mechanism controlling the differential roles of the MLV p12 protein in early and late replication. During viral assembly and egress, the late domain within the p12 N terminus functions to bind host vesicle release factors. During viral entry, the C terminus of p12 is required for tethering to host mitotic chromosomes. Our studies indicate that the p12 domain including the PPPY late sequence temporally represses the p12 chromatin tethering motif. Maximal p12 tethering was identified with only an 11-amino-acid minimal chromatin tethering motif encoded at p1261-71 Within this region, the p12-M63I substitution switches p12 into a tethering-competent state, partially rescuing the p12-PM15 tethering mutant. A model for how this conformational change regulates early versus late functions is presented.

Phosphorylation Requirement of Murine Leukemia Virus p12.

The p12 protein of murine leukemia virus (MLV) Gag is associated with the preintegration complex (PIC), and mutants of p12 (PM14) exhibit defects in nuclear entry/retention. Mutants of the phosphorylated serine 61 also have been reported to have defects in the early life cycle. Here we show that a phosphorylated peptide motif derived from human papillomavirus 8 (HPV-8), the E2 hinge region including residues 240 to 255, can functionally replace the main phosphorylated motif of MLV p12 and can rescue the viral titer of a strain with the lethal p12-PM14 mutation. Complementation with the HPV-8 E2 hinge motif generated multiple second-site mutations in live viral passage assays. Additional p12 phosphorylation sites were detected, including the late domain of p12 (PPPY) as well as the late domain/protease cleavage site of matrix (LYPAL), by mass spectrometry and Western blotting. Chromatin binding of p12-green fluorescent protein (GFP) fusion protein and functional complementation of p12-PM14 occurred in a manner independent of the E2 hinge region phosphorylation. Replacement of serine 61 by alanine within the minimal tethering domain (61SPMASRLRGRR71) maintained tethering, but in the context of the full-length p12, mutants with substitutions in S61 remained untethered and lost infectivity, indicating phosphorylation of p12 serine 61 functions to temporally regulate early and late p12 functions. IMPORTANCE: The p12 protein, required for both early and late viral functions, is the predominant phosphorylated viral protein of Moloney MLV and is required for virus viability. Our studies indicate that the N terminus of p12 represses the early function of the chromatin binding domain and that deletion of the N terminus activates chromatin binding in the wild-type Moloney MLV p12 protein. Mass spectrometry and mutagenesis studies suggest that phosphorylation of both the repression domain and the chromatin binding domain acts to temporally regulate this process at the appropriate stages during infection.

Structural and sequencing analysis of local target DNA recognition by MLV integrase.

Target-site selection by retroviral integrase (IN) proteins profoundly affects viral pathogenesis. We describe the solution nuclear magnetic resonance structure of the Moloney murine leukemia virus IN (M-MLV) C-terminal domain (CTD) and a structural homology model of the catalytic core domain (CCD). In solution, the isolated MLV IN CTD adopts an SH3 domain fold flanked by a C-terminal unstructured tail. We generated a concordant MLV IN CCD structural model using SWISS-MODEL, MMM-tree and I-TASSER. Using the X-ray crystal structure of the prototype foamy virus IN target capture complex together with our MLV domain structures, residues within the CCD α2 helical region and the CTD ß1-ß2 loop were predicted to bind target DNA. The role of these residues was analyzed in vivo through point mutants and motif interchanges. Viable viruses with substitutions at the IN CCD α2 helical region and the CTD ß1-ß2 loop were tested for effects on integration target site selection. Next-generation sequencing and analysis of integration target sequences indicate that the CCD α2 helical region, in particular P187, interacts with the sequences distal to the scissile bonds whereas the CTD ß1-ß2 loop binds to residues proximal to it. These findings validate our structural model and disclose IN-DNA interactions relevant to target site selection.

Role of Cysteines in Stabilizing the Randomized Receptor Binding Domains within Feline Leukemia Virus Envelope Proteins.

UNLABELLED: Retargeting of gammaretroviral envelope proteins has shown promising results in the isolation of novel isolates with therapeutic potential. However, the optimal conditions required to obtain high-affinity retargeted envelope proteins with narrow tropism are not understood. This study highlights the advantage of constrained peptides within receptor binding domains and validates the random library screening technique of obtaining novel retargeted Env proteins. Using a modified vector backbone to screen the envelope libraries on 143B osteosarcoma cells, three novel and unique retargeted envelopes were isolated. The use of complex disulfide bonds within variable regions required for receptor binding is found within natural gammaretroviral envelope isolates. Interestingly, two of the isolates, named AII and BV2, have a pair of cysteines located within the randomized region of 11 amino acids similar to that identified within the CP Env, an isolate identified in a previous Env library screen on the human renal carcinoma Caki-1 cell line. The amino acids within the randomized region of AII and BV2 envelopes that are essential for viral infection have been identified in this study and include these cysteine residues. Through mutagenesis studies, the putative disulfide bond pairs including and beyond the randomized region were examined. In parallel, the disulfide bonds of CP Env were identified using mass spectrometry. The results indicate that this pair of cysteines creates the structural context to position key hydrophobic (F and W) and basic (K and H) residues critical for viral titer and suggest that AII, BV2, and CP internal cysteines bond together in distinct ways. IMPORTANCE: Retargeted gammaretroviral particles have broad applications for therapeutic use. Although great advances have been achieved in identifying new Env-host cell receptor pairs, the rules for designing optimal Env libraries are still unclear. We have found that isolates with an additional pair of cysteines within the randomized region have the highest transduction efficiencies. This emphasizes the importance of considering cysteine pairs in the design of new libraries. Furthermore, our data clearly indicate that these cysteines are essential for viral infectivity by presenting essential residues to the host cell receptor. These studies facilitate the screening of Env libraries for functional entry into target cells, allowing the identification of novel gammaretroviral Envs targeting alternative host cell receptors for gene and protein delivery.

Bimodal high-affinity association of Brd4 with murine leukemia virus integrase and mononucleosomes.

The importance of understanding the molecular mechanisms of murine leukemia virus (MLV) integration into host chromatin is highlighted by the development of MLV-based vectors for human gene-therapy. We have recently identified BET proteins (Brd2, 3 and 4) as the main cellular binding partners of MLV integrase (IN) and demonstrated their significance for effective MLV integration at transcription start sites. Here we show that recombinant Brd4, a representative of the three BET proteins, establishes complementary high-affinity interactions with MLV IN and mononucleosomes (MNs). Brd4(1-720) but not its N- or C-terminal fragments effectively stimulate MLV IN strand transfer activities in vitro. Mass spectrometry- and NMR-based approaches have enabled us to map key interacting interfaces between the C-terminal domain of BRD4 and the C-terminal tail of MLV IN. Additionally, the N-terminal fragment of Brd4 binds to both DNA and acetylated histone peptides, allowing it to bind tightly to MNs. Comparative analyses of the distributions of various histone marks along chromatin revealed significant positive correlations between H3- and H4-acetylated histones, BET protein-binding sites and MLV-integration sites. Our findings reveal a bimodal mechanism for BET protein-mediated MLV integration into select chromatin locations.

Altering murine leukemia virus integration through disruption of the integrase and BET protein family interaction.

We report alterations to the murine leukemia virus (MLV) integrase (IN) protein that successfully result in decreasing its integration frequency at transcription start sites and CpG islands, thereby reducing the potential for insertional activation. The host bromo and extraterminal (BET) proteins Brd2, 3 and 4 interact with the MLV IN protein primarily through the BET protein ET domain. Using solution NMR, protein interaction studies, and next generation sequencing, we show that the C-terminal tail peptide region of MLV IN is important for the interaction with BET proteins and that disruption of this interaction through truncation mutations affects the global targeting profile of MLV vectors. The use of the unstructured tails of gammaretroviral INs to direct association with complexes at active promoters parallels that used by histones and RNA polymerase II. Viruses bearing MLV IN C-terminal truncations can provide new avenues to improve the safety profile of gammaretroviral vectors for human gene therapy.