Network-based elucidation of colon cancer drug resistance mechanisms by phosphoproteomic time-series analysis.

Aberrant signaling pathway activity is a hallmark of tumorigenesis and progression, which has guided targeted inhibitor design for over 30 years. Yet, adaptive resistance mechanisms, induced by rapid, context-specific signaling network rewiring, continue to challenge therapeutic efficacy. Leveraging progress in proteomic technologies and network-based methodologies, we introduce Virtual Enrichment-based Signaling Protein-activity Analysis (VESPA)-an algorithm designed to elucidate mechanisms of cell response and adaptation to drug perturbations-and use it to analyze 7-point phosphoproteomic time series from colorectal cancer cells treated with clinically-relevant inhibitors and control media. Interrogating tumor-specific enzyme/substrate interactions accurately infers kinase and phosphatase activity, based on their substrate phosphorylation state, effectively accounting for signal crosstalk and sparse phosphoproteome coverage. The analysis elucidates time-dependent signaling pathway response to each drug perturbation and, more importantly, cell adaptive response and rewiring, experimentally confirmed by CRISPR knock-out assays, suggesting broad applicability to cancer and other diseases.

Network-based elucidation of colon cancer drug resistance by phosphoproteomic time-series analysis.

Aberrant signaling pathway activity is a hallmark of tumorigenesis and progression, which has guided targeted inhibitor design for over 30 years. Yet, adaptive resistance mechanisms, induced by rapid, context-specific signaling network rewiring, continue to challenge therapeutic efficacy. By leveraging progress in proteomic technologies and network-based methodologies, over the past decade, we developed VESPA-an algorithm designed to elucidate mechanisms of cell response and adaptation to drug perturbations-and used it to analyze 7-point phosphoproteomic time series from colorectal cancer cells treated with clinically-relevant inhibitors and control media. Interrogation of tumor-specific enzyme/substrate interactions accurately inferred kinase and phosphatase activity, based on their inferred substrate phosphorylation state, effectively accounting for signal cross-talk and sparse phosphoproteome coverage. The analysis elucidated time-dependent signaling pathway response to each drug perturbation and, more importantly, cell adaptive response and rewiring that was experimentally confirmed by CRISPRko assays, suggesting broad applicability to cancer and other diseases.

Integrating 3D structural information into systems biology.

Systems biology is a data-heavy field that focuses on systems-wide depictions of biological phenomena necessarily sacrificing a detailed characterization of individual components. As an example, genome-wide protein interaction networks are widely used in systems biology and continuously extended and refined as new sources of evidence become available. Despite the vast amount of information about individual protein structures and protein complexes that has accumulated in the past 50 years in the Protein Data Bank, the data, computational tools, and language of structural biology are not an integral part of systems biology. However, increasing effort has been devoted to this integration, and the related literature is reviewed here. Relationships between proteins that are detected via structural similarity offer a rich source of information not available from sequence similarity, and homology modeling can be used to leverage Protein Data Bank structures to produce 3D models for a significant fraction of many proteomes. A number of structure-informed genomic and cross-species (i.e., virus-host) interactomes will be described, and the unique information they provide will be illustrated with a number of examples. Tissue- and tumor-specific interactomes have also been developed through computational strategies that exploit patient information and through genetic interactions available from increasingly sensitive screens. Strategies to integrate structural information with these alternate data sources will be described. Finally, efforts to link protein structure space with chemical compound space offer novel sources of information in drug design, off-target identification, and the identification of targets for compounds found to be effective in phenotypic screens.

Education and Outreach in Physical Sciences in Oncology.

Physical sciences are often overlooked in the field of cancer research. The Physical Sciences in Oncology Initiative was launched to integrate physics, mathematics, chemistry, and engineering with cancer research and clinical oncology through education, outreach, and collaboration. Here, we provide a framework for education and outreach in emerging transdisciplinary fields.

Oncoprotein-specific molecular interaction maps (SigMaps) for cancer network analyses.

Tumor-specific elucidation of physical and functional oncoprotein interactions could improve tumorigenic mechanism characterization and therapeutic response prediction. Current interaction models and pathways, however, lack context specificity and are not oncoprotein specific. We introduce SigMaps as context-specific networks, comprising modulators, effectors and cognate binding-partners of a specific oncoprotein. SigMaps are reconstructed de novo by integrating diverse evidence sources-including protein structure, gene expression and mutational profiles-via the OncoSig machine learning framework. We first generated a KRAS-specific SigMap for lung adenocarcinoma, which recapitulated published KRAS biology, identified novel synthetic lethal proteins that were experimentally validated in three-dimensional spheroid models and established uncharacterized crosstalk with RAB/RHO. To show that OncoSig is generalizable, we first inferred SigMaps for the ten most mutated human oncoproteins and then for the full repertoire of 715 proteins in the COSMIC Cancer Gene Census. Taken together, these SigMaps show that the cell's regulatory and signaling architecture is highly tissue specific.

Protein kinase Cθ C2 domain is a phosphotyrosine binding module that plays a key role in its activation.

Protein kinase Cθ (PKCθ) is a novel PKC that plays a key role in T lymphocyte activation. To understand how PKCθ is regulated in T cells, we investigated the properties of its N-terminal C2 domain that functions as an autoinhibitory domain. Our measurements show that a Tyr(P)-containing peptide derived from CDCP1 binds the C2 domain of PKCθ with high affinity and activates the enzyme activity of the intact protein. The Tyr(P) peptide also binds the C2 domain of PKCδ tightly, but no enzyme activation was observed with PKCδ. Mutations of PKCθ-C2 residues involved in Tyr(P) binding abrogated the enzyme activation and association of PKCθ with Tyr-phosphorylated full-length CDCP1 and severely inhibited the T cell receptor/CD28-mediated activation of a PKCθ-dependent reporter gene in T cells. Collectively, these studies establish the C2 domain of PKCθ as a Tyr(P)-binding domain and suggest that the domain may play a major role in PKCθ activation via its Tyr(P) binding.

Molecular dynamics simulations of PIP2 and PIP3 in lipid bilayers: determination of ring orientation, and the effects of surface roughness on a Poisson-Boltzmann description.

Molecular dynamics (MD) simulations of phosphatidylinositol (4,5)-bisphosphate (PIP2) and phosphatidylinositol (3,4,5)-trisphosphate (PIP3) in 1-palmitoyl 2-oleoyl phosphatidylcholine (POPC) bilayers indicate that the inositol rings are tilted approximately 40 degrees with respect to the bilayer surface, as compared with 17 degrees for the P-N vector of POPC. Multiple minima were obtained for the ring twist (analogous to roll for an airplane). The phosphates at position 1 of PIP2 and PIP3 are within an Angström of the plane formed by the phosphates of POPC; lipids in the surrounding shell are depressed by 0.5-0.8 A, but otherwise the phosphoinositides do not substantially perturb the bilayer. Finite size artifacts for ion distributions are apparent for systems of approximately 26 waters/lipid, but, based on simulations with a fourfold increase of the aqueous phase, the phosphoinositide positions and orientations do not show significant size effects. Electrostatic potentials evaluated from Poisson-Boltzmann (PB) calculations show a strong dependence of potential height and ring orientation, with the maxima on the -25 mV surfaces (17.1 +/- 0.1 A for PIP2 and 19.4 +/- 0.3 A for PIP3) occurring near the most populated orientations from MD. These surfaces are well above the background height of 10 A estimated for negatively charged cell membranes, as would be expected for lipids involved in cellular signaling. PB calculations on microscopically flat bilayers yield similar maxima as the MD-based (microscopically rough) systems, but show less fine structure and do not clearly indicate the most probable regions. Electrostatic free energies of interaction with pentalysine are also similar for the rough and flat systems. These results support the utility of a rigid/flat bilayer model for PB-based studies of PIP2 and PIP3 as long as the orientations are judiciously chosen.

Structure of an acetyl-CoA binding protein from Staphylococcus aureus representing a novel subfamily of GCN5-related N-acetyltransferase-like proteins.

We have determined the solution NMR structure of SACOL2532, a putative GCN5-like N-acetyltransferase (GNAT) from Staphylococcus aureus. SACOL2532 was shown to bind both CoA and acetyl-CoA, and structures with and without bound CoA were determined. Based on analysis of the structure and sequence, a subfamily of small GCN5-related N-acetyltransferase (GNAT)-like proteins can be defined. Proteins from this subfamily, which is largely congruent with COG2388, are characterized by a cysteine residue in the acetyl-CoA binding site near the acetyl group, by their small size in relation to other GNATs, by a lack of obvious substrate binding site, and by a distinct conformation of bound CoA in relation to other GNATs. Subfamily members are found in many bacterial and eukaryotic genomes, and in some archaeal genomes. Whereas other GNATs transfer the acetyl group of acetyl-CoA directly to an aliphatic amine, the presence of the conserved cysteine residue suggests that proteins in the COG2388 GNAT-subfamily transfer an acetyl group from acetyl-CoA to one or more presently unidentified aliphatic amines via an acetyl (cysteine) enzyme intermediate. The apparent absence of a substrate-binding region suggests that the substrate is a macromolecule, such as another protein, or that a second protein subunit providing a substrate-binding region must combine with SACOL2532 to make a fully functional N-acetyl transferase.

A molecular switch required for retrovirus assembly participates in the hexagonal immature lattice.

In the Rous sarcoma virus (RSV) Gag protein, the 25 amino-acid residues of the p10 domain immediately upstream of the CA domain are essential for immature particle formation. We performed systematic mutagenesis on this region and found excellent correlation between the amino-acid side chains required for in vitro assembly and those that participate in the p10-CA dimer interface in a previously described crystal structure. We introduced exogenous cysteine residues that were predicted to form disulphide bonds across the dimer interface. Upon oxidation of immature particles, a disulphide-linked Gag hexamer was formed, implying that p10 participates in and stabilizes the immature Gag hexamer. This is the first example of a critical interaction between two different Gag domains. Molecular modeling of the RSV immature hexamer indicates that the N-terminal domains of CA must expand relative to the murine leukaemia virus mature hexamer to accommodate the p10 contact; this expansion is strikingly similar to recent cryotomography results for immature human immunodeficiency virus particles.

Structural and membrane binding analysis of the Phox homology domain of Bem1p: basis of phosphatidylinositol 4-phosphate specificity.

Phox homology (PX) domains, which have been identified in a variety of proteins involved in cell signaling and membrane trafficking, have been shown to interact with phosphoinositides (PIs) with different affinities and specificities. To elucidate the structural origin of the diverse PI specificity of PX domains, we determined the crystal structure of the PX domain from Bem1p that has been reported to bind phosphatidylinositol 4-phosphate (PtdIns(4)P). We also measured the membrane binding properties of the PX domain and its mutants by surface plasmon resonance and monolayer techniques and calculated the electrostatic potentials for the PX domain in the absence and presence of bound PtdIns(4)P. The Bem1p PX domain contains a signature PI-binding site optimized for PtdIns(4)P binding and also harbors basic and hydrophobic residues on the membrane-binding surface. The membrane binding of the Bem1p PX domain is initiated by nonspecific electrostatic interactions between the cationic membrane-binding surface of the domain and anionic membrane surfaces, followed by the membrane penetration of hydrophobic residues. Unlike other PX domains, the Bem1p PX domain has high intrinsic membrane penetrating activity in the absence of PtdIns(4)P, suggesting that the partial membrane penetration may occur before specific PtdIns(4)P binding and last after the removal of PtdIns(4)P under certain conditions. This structural and functional study of the PtdIns(4)P-binding Bem1p PX domain provides new insight into the diverse PI specificities and membrane-binding mechanisms of PX domains.

Binding of phosphoinositide-specific phospholipase C-zeta (PLC-zeta) to phospholipid membranes: potential role of an unstructured cluster of basic residues.

Phospholipase C-zeta (PLC-zeta) is a sperm-specific enzyme that initiates the Ca2+ oscillations in mammalian eggs that activate embryo development. It shares considerable sequence homology with PLC-delta1, but lacks the PH domain that anchors PLC-delta1 to phosphatidylinositol 4,5-bisphosphate, PIP2. Thus it is unclear how PLC-zeta interacts with membranes. The linker region between the X and Y catalytic domains of PLC-zeta, however, contains a cluster of basic residues not present in PLC-delta1. Application of electrostatic theory to a homology model of PLC-zeta suggests this basic cluster could interact with acidic lipids. We measured the binding of catalytically competent mouse PLC-zeta to phospholipid vesicles: for 2:1 phosphatidylcholine/phosphatidylserine (PC/PS) vesicles, the molar partition coefficient, K, is too weak to be of physiological significance. Incorporating 1% PIP2 into the 2:1 PC/PS vesicles increases K about 10-fold, to 5x10(3) M-1, a biologically relevant value. Expressed fragments corresponding to the PLC-zeta X-Y linker region also bind with higher affinity to polyvalent than monovalent phosphoinositides on nitrocellulose filters. A peptide corresponding to the basic cluster (charge=+7) within the linker region, PLC-zeta-(374-385), binds to PC/PS vesicles with higher affinity than PLC-zeta, but its binding is less sensitive to incorporating PIP2. The acidic residues flanking this basic cluster in PLC-zeta may account for both these phenomena. FRET experiments suggest the basic cluster could not only anchor the protein to the membrane, but also enhance the local concentration of PIP2 adjacent to the catalytic domain.

Electrostatic interactions drive membrane association of the human immunodeficiency virus type 1 Gag MA domain.

The assembly of most retroviruses occurs at the plasma membrane. Membrane association is directed by MA, the N-terminal domain of the Gag structural protein. For human immunodeficiency virus type 1 (HIV-1), this association is mediated in part by a myristate fatty acid modification. Conflicting evidence has been presented on the relative importance of myristoylation, of ionic interactions between protein and membrane, and of Gag multimerization in membrane association in vivo. We addressed these questions biochemically by determining the affinity of purified myristoylated HIV-1 MA for liposomes of defined composition, both for monomeric and for dimeric forms of the protein. Myristoylation increases the barely detectable intrinsic affinity of the apo-protein for liposomes by only 10-fold, and the resulting affinity is still weak, similar to that of the naturally nonmyristoylated MA of Rous sarcoma virus. Membrane binding of HIV-1 MA is absolutely dependent on the presence of negatively charged lipid and is abrogated at high ionic strength. Forced dimerization of MA increases its membrane affinity by several orders of magnitude. When green fluorescent protein fusions of monomeric or dimeric MA are expressed in cells, the dimeric but not the monomeric protein becomes strongly membrane associated. Computational modeling supports these results and suggests a molecular mechanism for the modest effect of myristoylation on binding, wherein the membrane provides a hydrophobic environment for the myristate that is energetically similar to that provided by the protein. Overall, the results imply that the driving force for membrane association stems largely from ionic interactions between multimerized Gag and negatively charged phospholipids.

Measuring the quality of pediatric day surgery care.

This brief report describes the process of developing a valid and reliable questionnaire for quality measurement of pediatric day surgery care from the families' perspective. Questionnaire items were generated through a literature search and interviews with clinicians and parents. A computer-assisted telephone interview was used to administer the questionnaire to 448 parents within 72 hours of patient discharge. Tests of reliability and validity were administered, and questionnaire items were improved or omitted based on the results. The investigators plan to readminister the improved questionnaire to confirm its validity and reliability.

The role of electrostatics in protein-membrane interactions.

Many experimental, structural and computational studies have established the importance of nonspecific electrostatics as a driving force for peripheral membrane association. Here we focus on this component of protein/membrane interactions by using examples ranging from phosphoinositide signaling to retroviral assembly. We stress the utility of the collaboration of experiment and theory in identifying and quantifying the role of electrostatics not only in contributing to membrane association, but also in affecting subcellular targeting, in the control of membrane binding, and in the organization of proteins and lipids at membrane surfaces.

Retroviral matrix domains share electrostatic homology: models for membrane binding function throughout the viral life cycle.

The matrix domain (MA) of Gag polyproteins performs multiple functions throughout the retroviral life cycle. MA structures have an electropositive surface patch that is implicated in membrane association. Here, we use computational methods to demonstrate that electrostatic control of membrane binding is a central characteristic of all retroviruses. We are able to explain a wide range of experimental observations and provide a level of quantitative and molecular detail that has been inaccessible to experiment. We further predict that MA may exist in a variety of oligomerization states and propose mechanistic models for the effects of phosphoinositides and phosphorylation. The calculations provide a conceptual model for how non-myristoylated and myristoylated MAs behave similarly in assembly and disassembly. Hence, they provide a unified quantitative picture of the structural and energetic origins of the entire range of MA function and thus enhance, extend, and integrate previous observations on individual stages of the process.

NMR solution structure of Thermotoga maritima protein TM1509 reveals a Zn-metalloprotease-like tertiary structure.

The 150-residue protein TM1509 is encoded in gene YF09_THEMA of Thermotoga maritima. TM1509 has so far no functional annotation and belongs to protein family UPF0054 (PFAM accession number: PF02130) which contains at least 146 members. The NMR structure of TM1509 reveals an alpha+beta fold comprising a four stranded beta-sheet with topology A( upward arrow), B( upward arrow), D( upward arrow), C( downward arrow) as well as five alpha-helices I-V. The structures of most members of family PF02130 can be reliably constructed using the TM1509 NMR structure, demonstrating high leverage for exploration of fold space. A multiple sequence alignment of TM1509 with homologues of family UPF0054 shows that three polypeptide segments, as well as a putative zinc-binding consensus motif HGXLHLXGYDH located at the C-terminal end of alpha-helix IV, are highly conserved. The spatial arrangement of the three His residues of this UPF0054 consensus motif is similar to the arrangement found for the His residues in the HEXXHXXGXXH zinc-binding consensus motif of matrix metallo-proteases (MMPs). Moreover, the other conserved polypeptide segments form a large cavity which encloses the putative Zn-binding pocket and might confer specificity during catalysis. However, TM1509 and the other members of the UPF0054 family do not have the crucial Glu residue in position 2 of the MMP consensus motif. Intriguingly, the TM1509 structure indicates that the Asp in the UPF0054 consensus motif (Asp 111 in TM1509) may overtake the catalytic role of the Glu. This suggests that protein family UPF0054 might contain members of a hitherto uncharacterized class of metalloproteases.

An electrostatic engine model for autoinhibition and activation of the epidermal growth factor receptor (EGFR/ErbB) family.

We propose a new mechanism to explain autoinhibition of the epidermal growth factor receptor (EGFR/ErbB) family of receptor tyrosine kinases based on a structural model that postulates both their juxtamembrane and protein tyrosine kinase domains bind electrostatically to acidic lipids in the plasma membrane, restricting access of the kinase domain to substrate tyrosines. Ligand-induced dimerization promotes partial trans autophosphorylation of ErbB1, leading to a rapid rise in intracellular [Ca(2+)] that can activate calmodulin. We postulate the Ca(2+)/calmodulin complex binds rapidly to residues 645--660 of the juxtamembrane domain, reversing its net charge from +8 to -8 and repelling it from the negatively charged inner leaflet of the membrane. The repulsion has two consequences: it releases electrostatically sequestered phosphatidylinositol 4,5-bisphosphate (PIP(2)), and it disengages the kinase domain from the membrane, allowing it to become fully active and phosphorylate an adjacent ErbB molecule or other substrate. We tested various aspects of the model by measuring ErbB juxtamembrane peptide binding to phospholipid vesicles using both a centrifugation assay and fluorescence correlation spectroscopy; analyzing the kinetics of interactions between ErbB peptides, membranes, and Ca(2+)/calmodulin using fluorescence stop flow; assessing ErbB1 activation in Cos1 cells; measuring fluorescence resonance energy transfer between ErbB peptides and PIP(2); and making theoretical electrostatic calculations on atomic models of membranes and ErbB juxtamembrane and kinase domains.

High-quality homology models derived from NMR and X-ray structures of E. coli proteins YgdK and Suf E suggest that all members of the YgdK/Suf E protein family are enhancers of cysteine desulfurases.

The structural biology of proteins mediating iron-sulfur (Fe-S) cluster assembly is central for understanding several important biological processes. Here we present the NMR structure of the 16-kDa protein YgdK from Escherichia coli, which shares 35% sequence identity with the E. coli protein SufE. The SufE X-ray crystal structure was solved in parallel with the YdgK NMR structure in the Northeast Structural Genomics (NESG) consortium. Both proteins are (1) key components for Fe-S metabolism, (2) exhibit the same distinct fold, and (3) belong to a family of at least 70 prokaryotic and eukaryotic sequence homologs. Accurate homology models were calculated for the YgdK/SufE family based on YgdK NMR and SufE crystal structure. Both structural templates contributed equally, exemplifying synergy of NMR and X-ray crystallography. SufE acts as an enhancer of the cysteine desulfurase activity of SufS by SufE-SufS complex formation. A homology model of CsdA, a desulfurase encoded in the same operon as YgdK, was modeled using the X-ray structure of SufS as a template. Protein surface and electrostatic complementarities strongly suggest that YgdK and CsdA likewise form a functional two-component desulfurase complex. Moreover, structural features of YgdK and SufS, which can be linked to their interaction with desulfurases, are conserved in all homology models. It thus appears very likely that all members of the YgdK/SufE family act as enhancers of Suf-S-like desulfurases. The present study exemplifies that "refined" selection of two (or more) targets enables high-quality homology modeling of large protein families.

A glutamate-alanine-leucine (EAL) domain protein of Salmonella controls bacterial survival in mice, antioxidant defence and killing of macrophages: role of cyclic diGMP.

Signature-tagged transposon mutagenesis of Salmonella with differential recovery from wild-type and immunodeficient mice revealed that the gene here named cdgR[for c-diguanylate (c-diGMP) regulator] is required for the bacterium to resist host phagocyte oxidase in vivo. CdgR consists solely of a glutamate-alanine-leucine (EAL) domain, a predicted cyclic diGMP (c-diGMP) phosphodiesterase. Disruption of cdgR decreased bacterial resistance to hydrogen peroxide and accelerated bacterial killing of macrophages. An ultrasensitive assay revealed c-diGMP in wild-type Salmonella with increased levels in the CdgR-deficient mutant. Thus, besides its known role in regulating cellulose synthesis and biofilm formation, bacterial c-diGMP also regulates host-pathogen interactions involving antioxidant defence and cytotoxicity.