Membrane association of the PTEN tumor suppressor: neutron scattering and MD simulations reveal the structure of protein-membrane complexes

Cancer and Autism: How PTEN Mutations Degrade Function at the Membrane and Isoform Expression in the Human Brain

Mutations causing loss of PTEN lipid phosphatase activity can promote cancer, benign tumors (PHTS), and neurodevelopmental disorders (NDDs). Exactly how they preferentially trigger distinct phenotypic outcomes has been puzzling. Here, we demonstrate that PTEN mutations differentially allosterically bias P loop dynamics and its connection to the catalytic site, affecting catalytic activity. NDD-related mutations are likely to sample conformations of the functional wild-type state, while sampled conformations for the strong, cancer-related driver mutation hotspots favor catalysis-primed conformations, suggesting that NDD mutations are likely to be weaker, and our large-scale simulations show why. Prenatal PTEN isoform expression data suggest exons 5 and 7, which harbor NDD mutations, as cancer-risk carriers. Since cancer requires more than a single mutation, our conformational and genomic analysis helps discover how same protein mutations can foster different clinical manifestations, articulates a role for co-occurring background latent driver mutations, and uncovers relationships of splicing isoform expression to life expectancy.

How PTEN mutations degrade function at the membrane and life expectancy of carriers of mutations in the human brain

PTEN dysfunction, caused by loss of lipid phosphatase activity or deletion, promotes pathologies, cancer, benign tumors, and neurodevelopmental disorders (NDDs). Despite efforts, exactly how the mutations trigger distinct phenotypic outcomes, cancer or NDD, has been puzzling. It has also been unclear how to distinguish between mutations harbored by isoforms, are they cancer or NDDs-related. Here we address both. We demonstrate that PTEN mutations differentially allosterically bias P-loop dynamics and its connection to the catalytic site, affecting catalytic activity. NDD-related mutations are likely to sample conformations present in the wild-type, while sampled conformations sheltering cancer-related hotspots favor catalysis-prone conformations, suggesting that NDD mutations are weaker. Analysis of isoform expression data indicates that if the transcript has NDD-related mutations, alone or in combination with cancer hotspots, there is high prenatal expression. If no mutations within the measured days, low expression levels. Cancer mutations promote stronger signaling and cell proliferation; NDDs' are weaker, influencing brain cell differentiation. Further, exon 5 is impacted by NDD or non-NDD mutations, while exon 7 is exclusively impacted by NDD mutations. Our comprehensive conformational and genomic analysis helps discover how same allele mutations can foster different clinical manifestations and uncovers correlations of splicing isoform expression to life expectancy.

Comparative Protein Structural Network Analysis Reveals C-Terminal Tail Phosphorylation Structural Communication Fingerprint in PTEN-Associated Mutations in Autism and Cancer

PTEN (phosphatase and tensin homolog deleted on chromosome 10) is a tightly regulated dual-specificity phosphatase and key regulator of the PI3K/AKT/mTOR signaling pathway. PTEN phosphorylation at its carboxy-terminal tail (CTT) serine/threonine cluster negatively regulates its tumor suppressor function by inducing a stable, closed, and inactive conformation. Germline PTEN mutations predispose individuals to PTEN hamartoma tumor syndrome (PHTS), a rare inherited cancer syndrome and, intriguingly, one of the most common causes of autism spectrum disorder (ASD). However, the mechanistic details that govern phosphorylated CTT catalytic conformational dynamics in the context of PHTS-associated mutations are unknown. Here, we utilized a comparative protein structure network (PSN)-based approach to investigate PTEN CTT phosphorylation-induced conformational dynamics specific to PTEN-ASD compared to PTEN-cancer phenotypes. Results from our study show differences in structural flexibility, inter-residue contacts, and allosteric communication patterns mediated by CTT phosphorylation, differentiating PTEN-ASD and PTEN-cancer phenotypes. Further, we identified perturbations among global metapaths and community network connections within the active site and inter-domain regions, indicating the significance of these regions in transmitting information across the PSN. Together, our studies provide a mechanistic underpinning of allosteric regulation through the coupled interplay of CTT phosphorylation conformational dynamics in PTEN-ASD and PTEN-cancer mutations. Importantly, the detailed atomistic interactions and structural consequences of PTEN variants reveal potential allosteric druggable target sites as a viable and currently unexplored treatment approach for individuals with different PHTS-associated mutations.

Structural and biophysical properties of farnesylated KRas interacting with the chaperone SmgGDS-558

KRas is a small GTPase and membrane-bound signaling protein. Newly synthesized KRas is post-translationally modified with a membrane-anchoring prenyl group. KRas chaperones are therapeutic targets in cancer due to their participation in trafficking oncogenic KRas to membranes. SmgGDS splice variants are chaperones for small GTPases with basic residues in their hypervariable domain (HVR), including KRas. SmgGDS-607 escorts pre-prenylated small GTPases, while SmgGDS-558 escorts prenylated small GTPases. We provide a structural description of farnesylated and fully processed KRas (KRas-FMe) in complex with SmgGDS-558 and define biophysical properties of this interaction. Surface plasmon resonance measurements on biomimetic model membranes quantified the thermodynamics of the interaction of SmgGDS with KRas, and small-angle x-ray scattering was used to characterize complexes of SmgGDS-558 and KRas-FMe structurally. Structural models were refined using Monte Carlo and molecular dynamics simulations. Our results indicate that SmgGDS-558 interacts with the HVR and the farnesylated C-terminus of KRas-FMe, but not its G-domain. Therefore, SmgGDS-558 interacts differently with prenylated KRas than prenylated RhoA, whose G-domain was found in close contact with SmgGDS-558 in a recent crystal structure. Using immunoprecipitation assays, we show that SmgGDS-558 binds the GTP-bound, GDP-bound, and nucleotide-free forms of farnesylated and fully processed KRas in cells, consistent with SmgGDS-558 not engaging the G-domain of KRas. We found that the dissociation constant, Kd, for KRas-FMe binding to SmgGDS-558 is comparable with that for the KRas complex with PDEδ, a well-characterized KRas chaperone that also does not interact with the KRas G-domain. These results suggest that KRas interacts in similar ways with the two chaperones SmgGDS-558 and PDEδ. Therapeutic targeting of the SmgGDS-558/KRas complex might prove as useful as targeting the PDEδ/KRas complex in KRas-driven cancers.

Structural and Dynamic Effects of PTEN C-Terminal Tail Phosphorylation

The phosphatase and tensin homologue deleted on chromosome 10 (PTEN) tumor suppressor gene encodes a tightly regulated dual-specificity phosphatase that serves as the master regulator of PI3K/AKT/mTOR signaling. The carboxy-terminal tail (CTT) is key to regulation and harbors multiple phosphorylation sites (Ser/Thr residues 380-385). CTT phosphorylation suppresses the phosphatase activity by inducing a stable, closed conformation. However, little is known about the mechanisms of phosphorylation-induced CTT-deactivation dynamics. Using explicit solvent microsecond molecular dynamics simulations, we show that CTT phosphorylation leads to a partially collapsed conformation, which alters the secondary structure of PTEN and induces long-range conformational rearrangements that encompass the active site. The active site rearrangements prevent localization of PTEN to the membrane, precluding lipid phosphatase activity. Notably, we have identified phosphorylation-induced allosteric coupling between the interdomain region and a hydrophobic site neighboring the active site in the phosphatase domain. Collectively, the results provide a mechanistic understanding of CTT phosphorylation dynamics and reveal potential druggable allosteric sites in a previously believed clinically undruggable protein.

Binding of Ca2+-independent C2 domains to lipid membranes: A multi-scale molecular dynamics study

C2 domains facilitate protein interactions with lipid bilayers in either a Ca²⁺-dependent or -independent manner. We used molecular dynamics (MD) simulations to explore six Ca²⁺-independent C2 domains, from KIBRA, PI3KC2α, RIM2, PTEN, SHIP2, and Smurf2. In coarse-grained MD simulations these C2 domains formed transient interactions with zwitterionic bilayers, compared with longer-lived interactions with anionic bilayers containing phosphatidylinositol bisphosphate (PIP₂). Type I C2 domains bound non-canonically via the front, back, or side of the β sandwich, whereas type II C2 domains bound canonically, via the top loops. C2 domains interacted strongly with membranes containing PIP₂, causing bound anionic lipids to cluster around the protein. Binding modes were refined via atomistic simulations. For PTEN and SHIP2, CG simulations of their phosphatase plus C2 domains with PIP₂-containing bilayers were also performed, and the roles of the two domains in membrane localization compared. These studies establish a simulation protocol for membrane-recognition proteins.

•

Binding of Ca²⁺-independent C2 domains to membranes was explored by MD simulation

•

C2 domains from KIBRA, PI3KC2α, RIM2, PTEN, SHIP2, and Smurf2 were compared

•

C2 domains formed longer-lived interactions with lipid bilayers containing PIP₂

•

For PTEN and SHIP2, simulations of their phosphatase plus C2 domains were performed

The mechanism of full activation of tumor suppressor PTEN at the phosphoinositide-enriched membrane

Tumor suppressor PTEN, the second most highly mutated protein in cancer, dephosphorylates signaling lipid PIP₃ produced by PI3Ks. Excess PIP₃ promotes cell proliferation. The mechanism at the membrane of this pivotal phosphatase is unknown hindering drug discovery. Exploiting explicit solvent simulations, we tracked full-length PTEN trafficking from the cytosol to the membrane. We observed its interaction with membranes composed of zwitterionic phosphatidylcholine, anionic phosphatidylserine, and phosphoinositides, including signaling lipids PIP₂ and PIP₃. We tracked its moving away from the zwitterionic and getting absorbed onto anionic membrane that harbors PIP₃. We followed it localizing on microdomains enriched in signaling lipids, as PI3K does, and observed PIP₃ allosterically unfolding the N-terminal PIP₂ binding domain, positioning it favorably for the polybasic motif interaction with PIP₂. Finally, we determined PTEN catalytic action at the membrane, all in line with experimental observations, deciphering the mechanisms of how PTEN anchors to the membrane and restrains cancer.

•

PTEN localizes on membrane microdomains enriched in phosphoinositides, as PI3K does

•

Full PTEN activation requires both signaling lipids, PIP₂ and PIP₃

•

Strong salt bridge interactions sustain stable PTEN membrane localization

•

Substrate-induced P loop conformational change implicates PTEN catalytic activity

The MemMoRF database for recognizing disordered protein regions interacting with cellular membranes

Protein and lipid membrane interactions play fundamental roles in a large number of cellular processes (e.g. signalling, vesicle trafficking, or viral invasion). A growing number of examples indicate that such interactions can also rely on intrinsically disordered protein regions (IDRs), which can form specific reversible interactions not only with proteins but also with lipids. We named IDRs involved in such membrane lipid-induced disorder-to-order transition as MemMoRFs, in an analogy to IDRs exhibiting disorder-to-order transition upon interaction with protein partners termed Molecular Recognition Features (MoRFs). Currently, both the experimental detection and computational characterization of MemMoRFs are challenging, and information about these regions are scattered in the literature. To facilitate the related investigations we generated a comprehensive database of experimentally validated MemMoRFs based on manual curation of literature and structural data. To characterize the dynamics of MemMoRFs, secondary structure propensity and flexibility calculated from nuclear magnetic resonance chemical shifts were incorporated into the database. These data were supplemented by inclusion of sentences from papers, functional data and disease-related information. The MemMoRF database can be accessed via a user-friendly interface at https://memmorf.hegelab.org, potentially providing a central resource for the characterization of disordered regions in transmembrane and membrane-associated proteins.

Membrane surface recognition by the ASAP1 PH domain and consequences for interactions with the small GTPase Arf1

Adenosine diphosphate-ribosylation factor (Arf) guanosine triphosphatase-activating proteins (GAPs) are enzymes that need to bind to membranes to catalyze the hydrolysis of guanosine triphosphate (GTP) bound to the small GTP-binding protein Arf. Binding of the pleckstrin homology (PH) domain of the ArfGAP With SH3 domain, ankyrin repeat and PH domain 1 (ASAP1) to membranes containing phosphatidylinositol 4,5-bisphosphate [PI(4,5)P₂] is key for maximum GTP hydrolysis but not fully understood. By combining nuclear magnetic resonance, neutron reflectometry, and molecular dynamics simulation, we show that binding of multiple PI(4,5)P₂ molecules to the ASAP1 PH domain (i) triggers a functionally relevant allosteric conformational switch and (ii) maintains the PH domain in a well-defined orientation, allowing critical contacts with an Arf1 mimic to occur. Our model provides a framework to understand how binding of the ASAP1 PH domain to PI(4,5)P₂ at the membrane may play a role in the regulation of ASAP1.

Laterally Resolved Small-Angle Scattering Intensity from Lipid Bilayer Simulations: An Exact and a Limited-Range Treatment

When combined, molecular simulations and small-angle scattering experiments are able to provide molecular-scale resolution of structure. Separately, scattering experiments provide only intermingled pair correlations between atoms, while molecular simulations are limited by model quality and the relatively short time scales that they can access. Their combined strength relies on agreement between the experimental spectra and those computed by simulation. To date, computing the neutron spectra from a molecular simulation of a lipid bilayer is straightforward only if the structure is approximated by laterally averaging the in-plane bilayer structure. However, this neglects all information about lateral heterogeneity, e.g., clustering of components in a lipid mixture. This paper presents two methods for computing the scattering intensity of simulated bilayers with in-plane heterogeneity, enabling a full treatment of both the transverse and lateral bilayer structure for the first time. The first method, termed the Dirac Brush, computes the exact spectra including spurious artifacts resulting from using information from neighboring periodic cells to account for the long-range structure of the bilayer. The second method, termed PFFT, applies a mean-field treatment in the field far from a scattering element, resulting in a correlation range that can be tuned (eliminating correlations with neighboring periodic images), but with computational cost that prohibits obtaining the exact (Dirac Brush) spectra. Following their derivation, the two methods are applied to a coarse-grained molecular simulation of a bilayer inhomogeneity, demonstrating the contributions of lateral correlations to the resulting spectra.

Structural Mechanisms of PTEN Regulation

The tumor suppressor phosphatase and tensin homolog on chromosome 10 (PTEN) is a tightly regulated enzyme responsible for dephosphorylating the progrowth lipid messenger molecule phosphatidylinositol 3,4,5-trisphosphate (PIP₃) on the plasma membrane. The carboxy-terminal tail (CTT) of PTEN is key for regulation of the enzyme. When phosphorylated, the unstructured CTT interacts with the phosphatase-C2 superdomain to inactivate the enzyme by preventing membrane association. PTEN mutations associated with cancer also inactivate the enzyme. Alternate translation-initiation sites generate extended isoforms of PTEN, such as PTEN-L that has multiple roles in cells. The extended amino-terminal region bears a signal sequence and a polyarginine sequence to facilitate exit from and entry into cells, respectively, and a membrane-binding helix that activates the enzyme. This amino-terminal region also facilitates mitochondrial and nucleolar localization. This review explores PTEN structure and its impact on localization and regulation.

Membrane Anchoring of Hck Kinase via the Intrinsically Disordered SH4-U and Length Scale Associated with Subcellular Localization

Src family kinases (SFKs) are a group of nonreceptor tyrosine kinases that are characterized by their involvement in critical signal transduction pathways. SFKs are often found attached to membranes, but little is known about the conformation of the protein in this environment. Here, solution nuclear magnetic resonance (NMR), neutron reflectometry (NR), and molecular dynamics (MD) simulations were employed to study the membrane interactions of the intrinsically disordered SH4 and Unique domains of the Src family kinase Hck. Through development of a procedure to combine the information from the different techniques, we were able produce a first-of-its-kind atomically detailed structural ensemble of a membrane-bound intrinsically disordered protein. Evaluation of the model demonstrated its consistency with previous work and provided insight into how SFK Unique domains act to differentiate the family members from one another. Fortuitously, the position of the ensemble on the membrane allowed the model to be combined with configurations of the multidomain Hck kinase previously determined from small-angle solution X-ray scattering to produce full-length models of membrane-anchored Hck. The resulting models allowed us to estimate that the kinase active site is positioned about 65 ± 35 Å away from the membrane surface, offering the first estimations of the length scale associated with the concept of SFK subcellular localization.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

Conformational Dynamics and Allosteric Regulation Landscapes of Germline PTEN Mutations Associated with Autism Compared to Those Associated with Cancer

Individuals with germline PTEN tumor-suppressor variants have PTEN hamartoma tumor syndrome (PHTS). Clinically, PHTS has variable presentations; there are distinct subsets of PHTS-affected individuals, such as those diagnosed with autism spectrum disorder (ASD) or cancer. It remains unclear why mutations in one gene can lead to such seemingly disparate phenotypes. Therefore, we sought to determine whether it is possible to predict a given PHTS-affected individual's a priori risk of ASD, cancer, or the co-occurrence of both phenotypes. By integrating network proximity analysis performed on the human interactome, molecular simulations, and residue-interaction networks, we demonstrate the role of conformational dynamics in the structural communication and long-range allosteric regulation of germline PTEN variants associated with ASD or cancer. We show that the PTEN interactome shares significant overlap with the ASD and cancer interactomes, providing network-based evidence that PTEN is a crucial player in the biology of both disorders. Importantly, this finding suggests that a germline PTEN variant might perturb the ASD or cancer networks differently, thus favoring one disease outcome at any one time. Furthermore, protein-dynamic structural-network analysis reveals small-world structural communication mediated by highly conserved functional residues and potential allosteric regulation of PTEN. We identified a salient structural-communication pathway that extends across the inter-domain interface for cancer-only mutations. In contrast, the structural-communication pathway is predominantly restricted to the phosphatase domain for ASD-only mutations. Our integrative approach supports the prediction and potential modulation of the relevant conformational states that influence structural communication and long-range perturbations associated with mutational effects that lead to PTEN-ASD or PTEN-cancer phenotypes.

Neutron scattering in the biological sciences: progress and prospects

The scattering of neutrons can be used to provide information on the structure and dynamics of biological systems on multiple length and time scales. Pursuant to a National Science Foundation-funded workshop in February 2018, recent developments in this field are reviewed here, as well as future prospects that can be expected given recent advances in sources, instrumentation and computational power and methods. Crystallography, solution scattering, dynamics, membranes, labeling and imaging are examined. For the extraction of maximum information, the incorporation of judicious specific deuterium labeling, the integration of several types of experiment, and interpretation using high-performance computer simulation models are often found to be particularly powerful.

Sensitivity of peripheral membrane proteins to the membrane context: A case study of phosphatidylserine and the TIM proteins

There is a diverse class of peripheral membrane-binding proteins that specifically bind phosphatidylserine (PS), a lipid that signals apoptosis or cell fusion depending on the membrane context of its presentation. PS-receptors are specialized for particular PS-presenting pathways, indicating that they might be sensitive to the membrane context. In this review, we describe a combination of thermodynamic, structural, and computational techniques that can be used to investigate the mechanisms underlying this sensitivity. As an example, we focus on three PS-receptors of the T-cell Immunoglobulin and Mucin containing (TIM) protein family, which we have previously shown to differ in their sensitivity to PS surface density.

Complex biomembrane mimetics on the sub-nanometer scale

Biomimetic lipid vesicles are indispensable tools for gaining insight into the biophysics of cell physiology on the molecular level. The level of complexity of these model systems has steadily increased, and now spans from domain-forming lipid mixtures to asymmetric lipid bilayers. Here, we review recent progress in the development and application of elastic neutron and X-ray scattering techniques for studying these systems in situ and under physiologically relevant conditions on the nanometer to sub-nanometer length scales. In particular, we focus on: (1) structural details of coexisting liquid-ordered and liquid-disordered domains, including their thickness and lipid packing mismatch as a function of a size transition from nanoscopic to macroscopic domains; (2) membrane-mediated protein partitioning into lipid domains; (3) the role of the aqueous medium in tuning interactions between membranes and domains; and (4) leaflet-specific structure in asymmetric bilayers and passive lipid flip-flop.

Structural characterization of membrane-bound human immunodeficiency virus-1 Gag matrix with neutron reflectometry

The structural characterization of peripheral membrane proteins represents a tremendous challenge in structural biology due to their transient interaction with the membrane and the potential multitude of protein conformations during this interaction. Neutron reflectometry is uniquely suited to address this problem because of its ability to structurally characterize biological model systems nondestructively and under biomimetic conditions that retain full protein functionality. Being sensitive to only the membrane-bound fraction of a water-soluble peripheral protein, neutron reflectometry obtains a low-resolution average structure of the protein-membrane complex that is further refined using integrative modeling strategies. Here, the authors review the current technological state of biological neutron reflectometry exemplified by a detailed report on the structure determination of the myristoylated human immunodeficiency virus-1 (HIV-1) Gag matrix associated with phosphoserine-containing model membranes. The authors found that the HIV-1 Gag matrix is able to adopt different configurations at the membrane in a pH-dependent manner and that the myristate group orients the protein in a way that is conducive to PIP2-binding.

Deuterium Labeling Strategies for Creating Contrast in Structure-Function Studies of Model Bacterial Outer Membranes Using Neutron Reflectometry

Studying the outer membrane of Gram-negative bacteria is challenging due to the complex nature of its structure. Therefore, simplified models are required to undertake structure-function studies of processes that occur at the outer membrane/fluid interface. Model membranes can be created by immobilizing bilayers to solid supports such as gold or silicon surfaces, or as monolayers on a liquid support where the surface pressure and fluidity of the lipids can be controlled. Both model systems are amenable to having their structure probed by neutron reflectometry, a technique that provides a one-dimensional depth profile through a membrane detailing its thickness and composition. One of the strengths of neutron scattering is the ability to use contrast matching, allowing molecules containing hydrogen and those enriched with deuterium to be highlighted or matched out against the bulk isotopic composition of the solvent. Lipopolysaccharides, a major component of the outer membrane, can be isolated for incorporation into model membranes. Here, we describe the deuteration of lipopolysaccharides from rough strains of Escherichia coli for incorporation into model outer membranes, and how the use of deuterated materials enhances structural analysis of model membranes by neutron reflectometry.

Voltage-sensing phosphatase modulation by a C2 domain

The voltage-sensing phosphatase (VSP) is the first example of an enzyme controlled by changes in membrane potential. VSP has four distinct regions: the transmembrane voltage-sensing domain (VSD), the inter-domain linker, the cytosolic catalytic domain, and the C2 domain. The VSD transmits the changes in membrane potential through the inter-domain linker activating the catalytic domain which then dephosphorylates phosphatidylinositol phosphate (PIP) lipids. The role of the C2, however, has not been established. In this study, we explore two possible roles for the C2: catalysis and membrane-binding. The Ci-VSP crystal structures show that the C2 residue Y522 lines the active site suggesting a contribution to catalysis. When we mutated Y522 to phenylalanine, we found a shift in the voltage dependence of activity. This suggests hydrogen bonding as a mechanism of action. Going one step further, when we deleted the entire C2 domain, we found voltage-dependent enzyme activity was no longer detectable. This result clearly indicates the entire C2 is necessary for catalysis as well as for modulating activity. As C2s are known membrane-binding domains, we tested whether the VSP C2 interacts with the membrane. We probed a cluster of four positively charged residues lining the top of the C2 and suggested by previous studies to interact with phosphatidylinositol 4,5-bisphosphate [PI(4,5)P₂] (Kalli et al., 2014). Neutralizing those positive charges significantly shifted the voltage dependence of activity to higher voltages. We tested membrane binding by depleting PI(4,5)P₂ from the membrane using the 5HT2C receptor and found that the VSD motions as measured by voltage clamp fluorometry (VCF) were not changed. These results suggest that if the C2 domain interacts with the membrane to influence VSP function it may not occur exclusively through PI(4,5)P₂. Together, this data advances our understanding of the VSP C2 by demonstrating a necessary and critical role for the C2 domain in VSP function.