Evolutionary sequence and structural basis for the distinct conformational landscapes of Tyr and Ser/Thr kinases.

Protein kinases are molecular machines with rich sequence variation that distinguishes the two main evolutionary branches - tyrosine kinases (TKs) from serine/threonine kinases (STKs). Using a sequence co-variation Potts statistical energy model we previously concluded that TK catalytic domains are more likely than STKs to adopt an inactive conformation with the activation loop in an autoinhibitory "folded" conformation, due to intrinsic sequence effects. Here we investigated the structural basis for this phenomenon by integrating the sequence-based model with structure-based molecular dynamics (MD) to determine the effects of mutations on the free energy difference between active and inactive conformations, using a novel thermodynamic cycle involving many (n=108) protein-mutation free energy perturbation (FEP) simulations in the active and inactive conformations. The sequence and structure-based results are consistent and support the hypothesis that the inactive conformation "DFG-out Activation Loop Folded", is a functional regulatory state that has been stabilized in TKs relative to STKs over the course of their evolution via the accumulation of residue substitutions in the activation loop and catalytic loop that facilitate distinct substrate binding modes in trans and additional modes of regulation in cis for TKs.

Kinetic coevolutionary models predict the temporal emergence of HIV-1 resistance mutations under drug selection pressure.

Drug resistance in HIV type 1 (HIV-1) is a pervasive problem that affects the lives of millions of people worldwide. Although records of drug-resistant mutations (DRMs) have been extensively tabulated within public repositories, our understanding of the evolutionary kinetics of DRMs and how they evolve together remains limited. Epistasis, the interaction between a DRM and other residues in HIV-1 protein sequences, is key to the temporal evolution of drug resistance. We use a Potts sequence-covariation statistical-energy model of HIV-1 protein fitness under drug selection pressure, which captures epistatic interactions between all positions, combined with kinetic Monte-Carlo simulations of sequence evolutionary trajectories, to explore the acquisition of DRMs as they arise in an ensemble of drug-naive patient protein sequences. We follow the time course of 52 DRMs in the enzymes protease, RT, and integrase, the primary targets of antiretroviral therapy. The rates at which DRMs emerge are highly correlated with their observed acquisition rates reported in the literature when drug pressure is applied. This result highlights the central role of epistasis in determining the kinetics governing DRM emergence. Whereas rapidly acquired DRMs begin to accumulate as soon as drug pressure is applied, slowly acquired DRMs are contingent on accessory mutations that appear only after prolonged drug pressure. We provide a foundation for using computational methods to determine the temporal evolution of drug resistance using Potts statistical potentials, which can be used to gain mechanistic insights into drug resistance pathways in HIV-1 and other infectious agents.

Evolutionary sequence and structural basis for the distinct conformational landscapes of Tyr and Ser/Thr kinases.

Protein kinases are molecular machines with rich sequence variation that distinguishes the two main evolutionary branches - tyrosine kinases (TKs) from serine/threonine kinases (STKs). Using a sequence co-variation Potts statistical energy model we previously concluded that TK catalytic domains are more likely than STKs to adopt an inactive conformation with the activation loop in an autoinhibitory "folded" conformation, due to intrinsic sequence effects. Here we investigated the structural basis for this phenomenon by integrating the sequence-based model with structure-based molecular dynamics (MD) to determine the effects of mutations on the free energy difference between active and inactive conformations, using a novel thermodynamic cycle involving many (n=108) protein-mutation free energy perturbation (FEP) simulations in the active and inactive conformations. The sequence and structure-based results are consistent and support the hypothesis that the inactive conformation "DFG-out Activation Loop Folded", is a functional regulatory state that has been stabilized in TKs relative to STKs over the course of their evolution via the accumulation of residue substitutions in the activation loop and catalytic loop that facilitate distinct substrate binding modes in trans and additional modes of regulation in cis for TKs.

Potts Hamiltonian Models and Molecular Dynamics Free Energy Simulations for Predicting the Impact of Mutations on Protein Kinase Stability.

Single-point mutations in kinase proteins can affect their stability and fitness, and computational analysis of these effects can provide insights into the relationships among protein sequence, structure, and function for this enzyme family. To assess the impact of mutations on protein stability, we used a sequence-based Potts Hamiltonian model trained on a kinase family multiple-sequence alignment (MSA) to calculate the statistical energy (fitness) effects of mutations and compared these against relative folding free energies (ΔΔGs) calculated from all-atom molecular dynamics free energy perturbation (FEP) simulations in explicit solvent. The fitness effects of mutations in the Potts model (ΔEs) showed good agreement with experimental thermostability data (Pearson r = 0.68), similar to the correlation we observed with ΔΔGs predicted from structure-based relative FEP simulations. Recognizing the possible advantages of using Potts models to rapidly estimate protein stability effects of kinase mutations seen in cancer genomics data, we used the Potts statistical energy model to estimate the stability effects of 65 conservative and nonconservative mutations across three distinct kinases (Wee1, Abl1, and Cdc7) with somatic mutations reported in the Genomic Data Commons (GDC) database. The ΔEs of these mutations calculated from the Potts model are consistent with the corresponding ΔΔGs from FEP simulations (Pearson ratio of 0.72). The agreement between these methods suggests that the Potts model may be used as a sequence-based tool for high-throughput screening of mutational effects as part of a computational pipeline for predicting the stability effects of mutations. We also demonstrate how the scalability of the fitness-based Potts model calculations permits analyses that are not easily accessed using FEP simulations. To this end, we employed site-saturation mutagenesis in the Potts model in order to investigate the relative stability effects of mutations seen in different cancer evolutionary scenarios. We used this approach to analyze the effects of drug pressure in Abl kinase by contrasting the relative fitness penalties of somatic mutations seen in miscellaneous cancer types with those calculated for mutations associated with cancer drug resistance. We observed that, in contrast to somatic mutations of Abl seen in various tumors that appear to have evolved neutrally, cancer mutations that evolved under drug pressure in Abl-targeted therapies tend to preserve enzyme stability.

Elucidating the Molecular Determinants of the Binding Modes of a Third-Generation HIV-1 Integrase Strand Transfer Inhibitor: The Importance of Side Chain and Solvent Reorganization.

The first- and second-generation clinically used HIV-1 integrase (IN) strand transfer inhibitors (INSTIs) are key components of antiretroviral therapy (ART), which work by blocking the integration step in the HIV-1 replication cycle that is catalyzed by a nucleoprotein assembly called an intasome. However, resistance to even the latest clinically used INSTIs is beginning to emerge. Developmental third-generation INSTIs, based on naphthyridine scaffolds, are promising candidates to combat drug-resistant viral variants. Among these novel INSTIs, compound 4f exhibits two distinct conformations when binding with intasomes from HIV-1 and the closely related prototype foamy virus (PFV) despite the high structural similarity of their INSTI binding pockets. The molecular mechanism and the key active site residues responsible for these differing binding modes in closely related intasomes remain elusive. To unravel the molecular determinants governing the two distinct binding modes, we applied a novel molecular dynamics-based free energy method that utilizes alchemical pathways to overcome the sampling challenges associated with transitioning between the two bound conformations of ligand 4f within the crowded environments of the INSTI binding pockets in these intasomes. The calculated conformational free energies successfully recapitulate the experimentally observed binding mode preferences in the two viral intasomes. Analysis of the simulated structures suggests that the observed binding mode preferences are caused by amino acid residue differences in both the front and the central catalytic sub-pocket of the INSTI binding site in HIV-1 and PFV. Additional free energy calculations on mutants of HIV-1 and PFV revealed that while both sub-pockets contribute to binding mode selection, the central sub-pocket plays a more important role. These results highlight the importance of both side chain and solvent reorganization, as well as the conformational entropy in determining the ligand binding mode, and will help inform the development of more effective INSTIs for combatting drug-resistant viral variants.

A clinical trial of therapeutic vaccination in lymphoma with serial tumor sampling and single-cell analysis.

ABSTRACT: In situ vaccination (ISV) triggers an immune response to tumor-associated antigens at 1 tumor site, which can then tackle the disease throughout the body. Here, we report clinical and biological results of a phase 1/2 ISV trial in patients with low-grade lymphoma, combining an intratumoral toll-like receptor 9 (TLR9) agonist with local low-dose radiation and ibrutinib (an inhibitor of B- and T-cell kinases). Adverse events were predominately low grade. The overall response rate was 50%, including 1 complete response. All patients experienced tumor reduction at distant sites. Single-cell analyses of serial fine needle aspirates from injected and uninjected tumors revealed correlates of clinical response, such as lower CD47 and higher major histocompatibility complex class II expression on tumor cells, enhanced T-cell and natural killer cell effector function, and reduced immune suppression from transforming growth factor ß and inhibitory T regulatory 1 cells. Although changes at the local injected site were more pronounced, changes at distant uninjected sites were more often associated with clinical responses. Functional immune response assays and tracking of T-cell receptor sequences provided evidence of treatment-induced tumor-specific T-cell responses. Induction of immune effectors and reversal of negative regulators were both important in producing clinically meaningful tumor responses. The trial was registered at www.clinicaltrials.gov as #NCT02927964.

Elucidating the molecular determinants for binding modes of a third-generation HIV-1 integrase strand transfer inhibitor: Importance of side chain and solvent reorganization.

The first and second-generation clinically used HIV-1 integrase (IN) strand transfer inhibitors (INSTIs) are key components of antiretroviral therapy (ART), which work by blocking the integration step in the HIV-1 replication cycle that is catalyzed by a nucleoprotein assembly called an intasome. However, resistance to even the latest clinically used INSTIs is beginning to emerge. Developmental third-generation INSTIs, based on naphthyridine scaffold, are promising candidates to combat drug-resistant viral variants. Among these novel INSTIs, compound 4f exhibits two distinct conformations when binding to intasomes from HIV-1 and the closely related prototype foamy virus (PFV), despite the high structural similarity of their INSTI binding pockets. The molecular mechanism and the key active site residues responsible for these differing binding modes in closely related intasomes remain elusive. To unravel the molecular determinants governing the two distinct binding modes, we employ a novel molecular dynamics-based free energy approach that utilizes alchemical pathways to overcome the sampling challenges associated with transitioning between two ligand conformations within crowded environments along physical pathways. The calculated conformational free energies successfully recapitulate the experimentally observed binding mode preferences in the two viral intasomes. Analysis of the simulated structures suggests that the observed binding mode preferences are caused by amino acid residue differences in both the front and the central catalytic sub-pocket of the INSTI binding site in HIV-1 and PFV. Additional free energy calculations on mutants of HIV-1 and PFV revealed that while both sub-pockets contribute to the binding mode selection, the central sub-pocket plays a more important role. These results highlight the importance of both side chain and solvent reorganization, as well as the conformational entropy in determining the ligand binding mode and will help inform the development of more effective INSTIs for combatting drug-resistant viral variants.

Follicular lymphoma evolves with a surmountable dependency on acquired glycosylation motifs in the B-cell receptor.

ABSTRACT: An early event in the genesis of follicular lymphoma (FL) is the acquisition of new glycosylation motifs in the B-cell receptor (BCR) due to gene rearrangement and/or somatic hypermutation. These N-linked glycosylation motifs (N-motifs) contain mannose-terminated glycans and can interact with lectins in the tumor microenvironment, activating the tumor BCR pathway. N-motifs are stable during FL evolution, suggesting that FL tumor cells are dependent on them for their survival. Here, we investigated the dynamics and potential impact of N-motif prevalence in FL at the single-cell level across distinct tumor sites and over time in 17 patients. Although most patients had acquired at least 1 N-motif as an early event, we also found (1) cases without N-motifs in the heavy or light chains at any tumor site or time point and (2) cases with discordant N-motif patterns across different tumor sites. Inferring phylogenetic trees of the patients with discordant patterns, we observed that both N-motif-positive and N-motif-negative tumor subclones could be selected and expanded during tumor evolution. Comparing N-motif-positive with N-motif-negative tumor cells within a patient revealed higher expression of genes involved in the BCR pathway and inflammatory response, whereas tumor cells without N-motifs had higher activity of pathways involved in energy metabolism. In conclusion, although acquired N-motifs likely support FL pathogenesis through antigen-independent BCR signaling in most patients with FL, N-motif-negative tumor cells can also be selected and expanded and may depend more heavily on altered metabolism for competitive survival.

Single-cell multi-gene identification of somatic mutations and gene rearrangements in cancer.

In this proof-of-concept study, we developed a single-cell method that provides genotypes of somatic alterations found in coding regions of messenger RNAs and integrates these transcript-based variants with their matching cell transcriptomes. We used nanopore adaptive sampling on single-cell complementary DNA libraries to validate coding variants in target gene transcripts, and short-read sequencing to characterize cell types harboring the mutations. CRISPR edits for 16 targets were identified using a cancer cell line, and known variants in the cell line were validated using a 352-gene panel. Variants in primary cancer samples were validated using target gene panels ranging from 161 to 529 genes. A gene rearrangement was also identified in one patient, with the rearrangement occurring in two distinct tumor sites.

Mechanisms of HIV-1 integrase resistance to dolutegravir and potent inhibition of drug-resistant variants.

HIV-1 infection depends on the integration of viral DNA into host chromatin. Integration is mediated by the viral enzyme integrase and is blocked by integrase strand transfer inhibitors (INSTIs), first-line antiretroviral therapeutics widely used in the clinic. Resistance to even the best INSTIs is a problem, and the mechanisms of resistance are poorly understood. Here, we analyze combinations of the mutations E138K, G140A/S, and Q148H/K/R, which confer resistance to INSTIs. The investigational drug 4d more effectively inhibited the mutants compared with the approved drug Dolutegravir (DTG). We present 11 new cryo-EM structures of drug-resistant HIV-1 intasomes bound to DTG or 4d, with better than 3-Å resolution. These structures, complemented with free energy simulations, virology, and enzymology, explain the mechanisms of DTG resistance involving E138K + G140A/S + Q148H/K/R and show why 4d maintains potency better than DTG. These data establish a foundation for further development of INSTIs that potently inhibit resistant forms in integrase.

The molecular mechanism of CD81 antibody inhibition of metastasis.

Metastases are reduced in CD81KO mice. In addition, a unique anti-CD81 antibody, 5A6, inhibits metastasis in vivo and invasion and migration in vitro. Here, we probed the structural components of CD81 required for the antimetastatic activity induced by 5A6. We found that the removal of either cholesterol or the intracellular domains of CD81 did not affect inhibition by the antibody. We show that the uniqueness of 5A6 is due not to increased affinity but rather to its recognition of a specific epitope on the large extracellular loop of CD81. Finally, we present a number of CD81 membrane-associated partners that may play a role in mediating the 5A6 antimetastatic attributes, including integrins and transferrin receptors.

Targeted TLR9 Agonist Elicits Effective Antitumor Immunity against Spontaneously Arising Breast Tumors.

Spontaneous tumors that arise in genetically engineered mice recapitulate the natural tumor microenvironment and tumor-immune coevolution observed in human cancers, providing a more physiologically relevant preclinical model relative to implanted tumors. Similar to many cancer patients, oncogene-driven spontaneous tumors are often resistant to immunotherapy, and thus novel agents that can effectively promote antitumor immunity against these aggressive cancers show considerable promise for clinical translation, and their mechanistic assessment can broaden our understanding of tumor immunology. In this study, we performed extensive immune profiling experiments to investigate how tumor-targeted TLR9 stimulation remodels the microenvironment of spontaneously arising tumors during an effective antitumor immune response. To model the clinical scenario of multiple tumor sites, we used MMTV-PyMT transgenic mice, which spontaneously develop heterogeneous breast tumors throughout their 10 mammary glands. We found that i.v. administration of a tumor-targeting TLR9 agonist, referred to as PIP-CpG, induced a systemic T cell-mediated immune response that not only promoted regression of existing mammary tumors, but also elicited immune memory capable of delaying growth of independent newly arising tumors. Within the tumor microenvironment, PIP-CpG therapy initiated an inflammatory cascade that dramatically amplified chemokine and cytokine production, prompted robust infiltration and expansion of innate and adaptive immune cells, and led to diverse and unexpected changes in immune phenotypes. This study demonstrates that effective systemic treatment of an autochthonous multisite tumor model can be achieved using a tumor-targeted immunostimulant and provides immunological insights that will inform future therapeutic strategies.

Safety and Tolerability of Whole Soybean Products: A Dose-Escalating Clinical Trial in Older Adults with Obesity.

Soybean products have nutrients, dietary fiber, and phytoalexins beneficial for cardiovascular and overall health. Despite their high consumption in Asian populations, their safety in Western diets is debated. We conducted a dose-escalating clinical trial of the safety and tolerability of soybean products in eight older adults (70-85 years) with obesity. Whole green soybean pods grown under controlled conditions were processed to flour (WGS) at the United States Department of Agriculture using common cooking techniques such as slicing and heat treatment. WGS incorporated into food products was consumed at 10 g, 20 g, and 30 g/day for one week at each dose. The gastrointestinal outcomes, clinical biomarkers, and adverse events were evaluated. We explored the stimulation of phytoalexin (glyceollin) production in live viable soybean seeds (LSS-G). We compared the compositions of WGS and LSS-G with commercial soybean flour and its fermented and enzymatically hydrolyzed forms. We found that although 30 g WSG was well-tolerated, and it made participants feel full. Our processing produced glyceollins (267 µg/g) in LSS-G. Processing soybean flour decreased the iron content, but reduced the oligosaccharides, which could attenuate flatulence. Providing soybean flour at <30 g/day may be prudent for overall health and to prevent the exclusion of other food groups and nutrients in older adults with obesity.

Lysine-Derived Charge-Altering Releasable Transporters: Targeted Delivery of mRNA and siRNA to the Lungs.

Targeted delivery of nucleic acid therapeutics to the lungs could transform treatment options for pulmonary disease. We have previously developed oligomeric charge-altering releasable transporters (CARTs) for in vivo mRNA transfection and demonstrated their efficacy for use in mRNA-based cancer vaccination and local immunomodulatory therapies against murine tumors. While our previously reported glycine-based CART-mRNA complexes (G-CARTs/mRNA) show selective protein expression in the spleen (mouse, >99%), here, we report a new lysine-derived CART-mRNA complex (K-CART/mRNA) that, without additives or targeting ligands, shows selective protein expression in the lungs (mouse, >90%) following systemic IV administration. We further show that by delivering siRNA using the K-CART, we can significantly decrease expression of a lung-localized reporter protein. Blood chemistry and organ pathology studies demonstrate that K-CARTs are safe and well-tolerated. We report on the new step economical, organocatalytic synthesis (two steps) of functionalized polyesters and oligo-carbonate-co-α-aminoester K-CARTs from simple amino acid and lipid-based monomers. The ability to direct protein expression selectively in the spleen or lungs by simple, modular changes to the CART structure opens fundamentally new opportunities in research and gene therapy.

Evolutionary divergence in the conformational landscapes of tyrosine vs serine/threonine kinases.

Inactive conformations of protein kinase catalytic domains where the DFG motif has a "DFG-out" orientation and the activation loop is folded present a druggable binding pocket that is targeted by FDA-approved 'type-II inhibitors' in the treatment of cancers. Tyrosine kinases (TKs) typically show strong binding affinity with a wide spectrum of type-II inhibitors while serine/threonine kinases (STKs) usually bind more weakly which we suggest here is due to differences in the folded to extended conformational equilibrium of the activation loop between TKs vs. STKs. To investigate this, we use sequence covariation analysis with a Potts Hamiltonian statistical energy model to guide absolute binding free-energy molecular dynamics simulations of 74 protein-ligand complexes. Using the calculated binding free energies together with experimental values, we estimated free-energy costs for the large-scale (~17-20 Å) conformational change of the activation loop by an indirect approach, circumventing the very challenging problem of simulating the conformational change directly. We also used the Potts statistical potential to thread large sequence ensembles over active and inactive kinase states. The structure-based and sequence-based analyses are consistent; together they suggest TKs evolved to have free-energy penalties for the classical 'folded activation loop' DFG-out conformation relative to the active conformation, that is, on average, 4-6 kcal/mol smaller than the corresponding values for STKs. Potts statistical energy analysis suggests a molecular basis for this observation, wherein the activation loops of TKs are more weakly 'anchored' against the catalytic loop motif in the active conformation and form more stable substrate-mimicking interactions in the inactive conformation. These results provide insights into the molecular basis for the divergent functional properties of TKs and STKs, and have pharmacological implications for the target selectivity of type-II inhibitors.

Contingency and Entrenchment of Drug-Resistance Mutations in HIV Viral Proteins.

The ability of HIV-1 to rapidly mutate leads to antiretroviral therapy (ART) failure among infected patients. Drug-resistance mutations (DRMs), which cause a fitness penalty to intrinsic viral fitness, are compensated by accessory mutations with favorable epistatic interactions which cause an evolutionary trapping effect, but the kinetics of this overall process has not been well characterized. Here, using a Potts Hamiltonian model describing epistasis combined with kinetic Monte Carlo simulations of evolutionary trajectories, we explore how epistasis modulates the evolutionary dynamics of HIV DRMs. We show how the occurrence of a drug-resistance mutation is contingent on favorable epistatic interactions with many other residues of the sequence background and that subsequent mutations entrench DRMs. We measure the time-autocorrelation of fluctuations in the likelihood of DRMs due to epistatic coupling with the sequence background, which reveals the presence of two evolutionary processes controlling DRM kinetics with two distinct time scales. Further analysis of waiting times for the evolutionary trapping effect to reverse reveals that the sequences which entrench (trap) a DRM are responsible for the slower time scale. We also quantify the overall strength of epistatic effects on the evolutionary kinetics for different mutations and show these are much larger for DRM positions than polymorphic positions, and we also show that trapping of a DRM is often caused by the collective effect of many accessory mutations, rather than a few strongly coupled ones, suggesting the importance of multiresidue sequence variations in HIV evolution. The analysis presented here provides a framework to explore the kinetic pathways through which viral proteins like HIV evolve under drug-selection pressure.

PARP14 is a novel target in STAT6 mutant follicular lymphoma.

The variable clinical course of follicular lymphoma (FL) is determined by the molecular heterogeneity of tumor cells and complex interactions within the tumor microenvironment (TME). IL-4 producing follicular helper T cells (TFH) are critical components of the FL TME. Binding of IL-4 to IL-4R on FL cells activates JAK/STAT signaling. We identified STAT6 mutations (STAT6MUT) in 13% of FL (N = 33/258), all clustered within the DNA binding domain. Gene expression data and immunohistochemistry showed upregulation of IL-4/STAT6 target genes in STAT6MUT FL, including CCL17, CCL22, and FCER2 (CD23). Functionally, STAT6MUT was gain-of-function by serial replating phenotype in pre-B CFU assays. Expression of STAT6MUT enhanced IL-4 induced FCER2/CD23, CCL17 and CCL22 expression and was associated with nuclear accumulation of pSTAT6. RNA sequencing identified PARP14 -a transcriptional switch and co-activator of STAT6- among the top differentially upregulated genes in IL-4 stimulated STAT6MUT lymphoma cells and in STAT6MUT primary FL cells. Quantitative chromatin immunoprecipitation (qChIP) demonstrated binding of STAT6MUT but not STAT6WT to the PARP14 promotor. Reporter assays showed increased IL-4 induced transactivation activity of STAT6MUT at the PARP14 promotor, suggesting a self-reinforcing regulatory circuit. Knock-down of PARP14 or PARP-inhibition abrogated the STAT6MUT gain-of-function phenotype. Thus, our results identify PARP14 as a novel therapeutic target in STAT6MUT FL.

Fingolimod-Conjugated Charge-Altering Releasable Transporters Efficiently and Specifically Deliver mRNA to Lymphocytes In Vivo and In Vitro.

Charge-altering releasable transporters (CARTs) are a class of oligonucleotide delivery vehicles shown to be effective for delivery of messenger RNA (mRNA) both in vitro and in vivo. Here, we exploited the chemical versatility of the CART synthesis to generate CARTs containing the small-molecule drug fingolimod (FTY720) as a strategy to increase mRNA delivery and expression in lymphocytes through a specific ligand-receptor interaction. Fingolimod is an FDA-approved small-molecule drug that, upon in vivo phosphorylation, binds to the sphingosine-1-phosphate receptor 1 (S1P1), which is highly expressed on lymphocytes. Compared to its non-fingolimod-conjugated analogue, the fingolimod-conjugated CART achieved superior transfection of activated human and murine T and B lymphocytes in vitro. The higher transfection of the fingolimod-conjugated CARTs was lost when cells were exposed to a free fingolimod before transfection. In vivo, the fingolimod-conjugated CART showed increased mRNA delivery to marginal zone B cells and NK cells in the spleen, relative to CARTs lacking fingolimod. Moreover, fingolimod-CART-mediated mRNA delivery induces peripheral blood T-cell depletion similar to free fingolimod. Thus, we show that functionalization of CARTs with a pharmacologically validated small molecule can increase transfection of a cellular population of interest while conferring some of the targeting properties of the conjugated small molecule to the CARTs.

Intratumoral immunotherapy relies on B and T cell collaboration.

Antitumor T cell responses are the primary mediators of cancer immunotherapy. However, many other components of the immune system are needed for efficient T cell responses to be generated. Here, we developed a combinatorial approach where a Toll-like receptor 9 agonist (CpG) and Fc-fused IL-12 protein were injected together into just one of several tumor sites in a mouse. This combination led to body-wide (abscopal) therapeutic responses in multiple cancer models. These systemic responses were dependent not only on T cells but also on B cells. B cells were activated by the treatment and were required for optimal T cell activation. This cross-talk was dependent on MHC and was tumor antigen specific. The addition of an agonistic antibody against OX40 further enhanced T cell activation and therapeutic responses. Our data suggest that the combination of CpG, anti-OX40, and IL-12Fc may have success in patients with cancer and that B and T cell collaboration is crucial for the efficacy of this combination immunotherapy.

Structure-based virtual screening workflow to identify antivirals targeting HIV-1 capsid.

We have identified novel HIV-1 capsid inhibitors targeting the PF74 binding site. Acting as the building block of the HIV-1 capsid core, the HIV-1 capsid protein plays an important role in the viral life cycle and is an attractive target for antiviral development. A structure-based virtual screening workflow for hit identification was employed, which includes docking 1.6 million commercially-available drug-like compounds from the ZINC database to the capsid dimer, followed by applying two absolute binding free energy (ABFE) filters on the 500 top-ranked molecules from docking. The first employs the Binding Energy Distribution Analysis Method (BEDAM) in implicit solvent. The top-ranked compounds are then refined using the Double Decoupling method in explicit solvent. Both docking and BEDAM refinement were carried out on the IBM World Community Grid as part of the FightAIDS@Home project. Using this virtual screening workflow, we identified 24 molecules with calculated binding free energies between - 6 and - 12 kcal/mol. We performed thermal shift assays on these molecules to examine their potential effects on the stability of HIV-1 capsid hexamer and found that two compounds, ZINC520357473 and ZINC4119064 increased the melting point of the latter by 14.8 °C and 33 °C, respectively. These results support the conclusion that the two ZINC compounds are primary hits targeting the capsid dimer interface. Our simulations also suggest that the two hit molecules may bind at the capsid dimer interface by occupying a new sub-pocket that has not been exploited by existing CA inhibitors. The possible causes for why other top-scored compounds suggested by ABFE filters failed to show measurable activity are discussed.