An atypical GABARAP binding module drives the pro-autophagic potential of the AML-associated NPM1c variant.

The nucleolar scaffold protein NPM1 is a multifunctional regulator of cellular homeostasis, genome integrity, and stress response. NPM1 mutations, known as NPM1c variants promoting its aberrant cytoplasmic localization, are the most frequent genetic alterations in acute myeloid leukemia (AML). A hallmark of AML cells is their dependency on elevated autophagic flux. Here, we show that NPM1 and NPM1c induce the autophagy-lysosome pathway by activating the master transcription factor TFEB, thereby coordinating the expression of lysosomal proteins and autophagy regulators. Importantly, both NPM1 and NPM1c bind to autophagy modifiers of the GABARAP subfamily through an atypical binding module preserved within its N terminus. The propensity of NPM1c to induce autophagy depends on this module, likely indicating that NPM1c exerts its pro-autophagic activity by direct engagement with GABARAPL1. Our data report a non-canonical binding mode of GABARAP family members that drives the pro-autophagic potential of NPM1c, potentially enabling therapeutic options.

Ubiquitination regulates ER-phagy and remodelling of endoplasmic reticulum.

The endoplasmic reticulum (ER) undergoes continuous remodelling via a selective autophagy pathway, known as ER-phagy1. ER-phagy receptors have a central role in this process2, but the regulatory mechanism remains largely unknown. Here we report that ubiquitination of the ER-phagy receptor FAM134B within its reticulon homology domain (RHD) promotes receptor clustering and binding to lipidated LC3B, thereby stimulating ER-phagy. Molecular dynamics (MD) simulations showed how ubiquitination perturbs the RHD structure in model bilayers and enhances membrane curvature induction. Ubiquitin molecules on RHDs mediate interactions between neighbouring RHDs to form dense receptor clusters that facilitate the large-scale remodelling of lipid bilayers. Membrane remodelling was reconstituted in vitro with liposomes and ubiquitinated FAM134B. Using super-resolution microscopy, we discovered FAM134B nanoclusters and microclusters in cells. Quantitative image analysis revealed a ubiquitin-mediated increase in FAM134B oligomerization and cluster size. We found that the E3 ligase AMFR, within multimeric ER-phagy receptor clusters, catalyses FAM134B ubiquitination and regulates the dynamic flux of ER-phagy. Our results show that ubiquitination enhances RHD functions via receptor clustering, facilitates ER-phagy and controls ER remodelling in response to cellular demands.

Heteromeric clusters of ubiquitinated ER-shaping proteins drive ER-phagy.

Membrane-shaping proteins characterized by reticulon homology domains play an important part in the dynamic remodelling of the endoplasmic reticulum (ER). An example of such a protein is FAM134B, which can bind LC3 proteins and mediate the degradation of ER sheets through selective autophagy (ER-phagy)1. Mutations in FAM134B result in a neurodegenerative disorder in humans that mainly affects sensory and autonomic neurons2. Here we report that ARL6IP1, another ER-shaping protein that contains a reticulon homology domain and is associated with sensory loss3, interacts with FAM134B and participates in the formation of heteromeric multi-protein clusters required for ER-phagy. Moreover, ubiquitination of ARL6IP1 promotes this process. Accordingly, disruption of Arl6ip1 in mice causes an expansion of ER sheets in sensory neurons that degenerate over time. Primary cells obtained from Arl6ip1-deficient mice or from patients display incomplete budding of ER membranes and severe impairment of ER-phagy flux. Therefore, we propose that the clustering of ubiquitinated ER-shaping proteins facilitates the dynamic remodelling of the ER during ER-phagy and is important for neuronal maintenance.

Identification of potential selective autophagy receptors from protein-content profiling of autophagosomes.

Selective autophagy receptors (SARs) are central to cellular homeostatic and organellar recycling pathways. Over the last two decades, more than 30 SARs have been discovered and validated using a variety of experimental approaches ranging from cell biology to biochemistry, including high-throughput imaging and screening methods. Yet, the extent of selective autophagy pathways operating under various cellular contexts, for example, under basal and starvation conditions, remains unresolved. Currently, our knowledge of all known SARs and their associated cargo components is fragmentary and limited by experimental data with varying degrees of resolution. Here, we use classical predictive and modeling approaches to integrate high-quality autophagosome content profiling data with disparate datasets. We identify a global set of potential SARs and their associated cargo components active under basal autophagy, starvation-induced, and proteasome-inhibition conditions. We provide a detailed account of cellular components, biochemical pathways, and molecular processes that are degraded via autophagy. Our analysis yields a catalog of new potential SARs that satisfy the characteristics of bonafide, well-characterized SARs. We categorize them by the subcellular compartments they emerge from and classify them based on their likely mode of action. Our structural modeling validates a large subset of predicted interactions with the human ATG8 family of proteins and shows characteristic, conserved LC3-interacting region (LIR)-LIR docking site (LDS) and ubiquitin-interacting motif (UIM)-UIM docking site (UDS) binding modes. Our analysis also revealed the most abundant cargo molecules targeted by these new SARs. Our findings expand the repertoire of SARs and provide unprecedented details into the global autophagic state of HeLa cells. Taken together, our findings provide motivation for the design of new experiments, testing the role of these novel factors in selective autophagy.

Role of FAM134 paralogues in endoplasmic reticulum remodeling, ER-phagy, and Collagen quality control.

Degradation of the endoplasmic reticulum (ER) via selective autophagy (ER-phagy) is vital for cellular homeostasis. We identify FAM134A/RETREG2 and FAM134C/RETREG3 as ER-phagy receptors, which predominantly exist in an inactive state under basal conditions. Upon autophagy induction and ER stress signal, they can induce significant ER fragmentation and subsequent lysosomal degradation. FAM134A, FAM134B/RETREG1, and FAM134C are essential for maintaining ER morphology in a LC3-interacting region (LIR)-dependent manner. Overexpression of any FAM134 paralogue has the capacity to significantly augment the general ER-phagy flux upon starvation or ER-stress. Global proteomic analysis of FAM134 overexpressing and knockout cell lines reveals several protein clusters that are distinctly regulated by each of the FAM134 paralogues as well as a cluster of commonly regulated ER-resident proteins. Utilizing pro-Collagen I, as a shared ER-phagy substrate, we observe that FAM134A acts in a LIR-independent manner and compensates for the loss of FAM134B and FAM134C, respectively. FAM134C instead is unable to compensate for the loss of its paralogues. Taken together, our data show that FAM134 paralogues contribute to common and unique ER-phagy pathways.

Membrane fusion and drug delivery with carbon nanotube porins.

Drug delivery mitigates toxic side effects and poor pharmacokinetics of life-saving therapeutics and enhances treatment efficacy. However, direct cytoplasmic delivery of drugs and vaccines into cells has remained out of reach. We find that liposomes studded with 0.8-nm-wide carbon nanotube porins (CNTPs) function as efficient vehicles for direct cytoplasmic drug delivery by facilitating fusion of lipid membranes and complete mixing of the membrane material and vesicle interior content. Fusion kinetics data and coarse-grained molecular dynamics simulations reveal an unusual mechanism where CNTP dimers tether the vesicles, pull the membranes into proximity, and then fuse their outer and inner leaflets. Liposomes containing CNTPs in their membranes and loaded with an anticancer drug, doxorubicin, were effective in delivering the drug to cancer cells, killing up to 90% of them. Our results open an avenue for designing efficient drug delivery carriers compatible with a wide range of therapeutics.

FAM134B-RHD Protein Clustering Drives Spontaneous Budding of Asymmetric Membranes.

Living cells constantly remodel the shape of their lipid membranes. In the endoplasmic reticulum (ER), the reticulon homology domain (RHD) of the reticulophagy regulator 1 (RETR1/FAM134B) forms dense autophagic puncta that are associated with membrane removal by ER-phagy. In molecular dynamics (MD) simulations, we find that FAM134B-RHD spontaneously forms clusters, driven in part by curvature-mediated attractions. At a critical size, as in a nucleation process, the FAM134B-RHD clusters induce the formation of membrane buds. The kinetics of budding depends sensitively on protein concentration and bilayer asymmetry. Our MD simulations shed light on the role of FAM134B-RHD in ER-phagy and show that membrane asymmetry can be used to modulate the kinetic barrier for membrane remodeling.

TBK1-mediated phosphorylation of LC3C and GABARAP-L2 controls autophagosome shedding by ATG4 protease.

Autophagy is a highly conserved catabolic process through which defective or otherwise harmful cellular components are targeted for degradation via the lysosomal route. Regulatory pathways, involving post-translational modifications such as phosphorylation, play a critical role in controlling this tightly orchestrated process. Here, we demonstrate that TBK1 regulates autophagy by phosphorylating autophagy modifiers LC3C and GABARAP-L2 on surface-exposed serine residues (LC3C S93 and S96; GABARAP-L2 S87 and S88). This phosphorylation event impedes their binding to the processing enzyme ATG4 by destabilizing the complex. Phosphorylated LC3C/GABARAP-L2 cannot be removed from liposomes by ATG4 and are thus protected from ATG4-mediated premature removal from nascent autophagosomes. This ensures a steady coat of lipidated LC3C/GABARAP-L2 throughout the early steps in autophagosome formation and aids in maintaining a unidirectional flow of the autophagosome to the lysosome. Taken together, we present a new regulatory mechanism of autophagy, which influences the conjugation and de-conjugation of LC3C and GABARAP-L2 to autophagosomes by TBK1-mediated phosphorylation.

Phospholipid Scramblases Remodel the Shape of Asymmetric Membranes.

The cell membrane and many organellar membranes are asymmetric and highly curved. In experiments, it is challenging to reconstitute and characterize membranes that differ in the lipid composition of their leaflets. Here we use molecular dynamics simulations to study the large-scale membrane shape changes associated with lipid shuttling between asymmetric leaflets. We exploit leaflet asymmetry to create a stable, near-spherical vesicle bud connected to a flat bilayer under periodic boundary conditions. Then we demonstrate how the lipid scramblase nhTMEM16 relaxes the lipid-number asymmetry. By mediating the flipping of lipids, this transmembrane protein dissipates the mechanochemical gradient between the leaflets and drives a large-scale membrane reorganization, converting the vesicle bud into a flat membrane. Our procedure to exploit bilayer asymmetry for simulations of highly curved membranes can be used to study the function of other lipid transporters and membrane-shaping proteins.

Membrane perforation by the pore-forming toxin pneumolysin.

Pneumolysin (PLY), a major virulence factor of Streptococcus pneumoniae, perforates cholesterol-rich lipid membranes. PLY protomers oligomerize as rings on the membrane and then undergo a structural transition that triggers the formation of membrane pores. Structures of PLY rings in prepore and pore conformations define the beginning and end of this transition, but the detailed mechanism of pore formation remains unclear. With atomistic and coarse-grained molecular dynamics simulations, we resolve key steps during PLY pore formation. Our simulations confirm critical PLY membrane-binding sites identified previously by mutagenesis. The transmembrane ß-hairpins of the PLY pore conformation are stable only for oligomers, forming a curtain-like membrane-spanning ß-sheet. Its hydrophilic inner face draws water into the protein-lipid interface, forcing lipids to recede. For PLY rings, this zone of lipid clearance expands into a cylindrical membrane pore. The lipid plug caught inside the PLY ring can escape by lipid efflux via the lower leaflet. If this path is too slow or blocked, the pore opens by membrane buckling, driven by the line tension acting on the detached rim of the lipid plug. Interestingly, PLY rings are just wide enough for the plug to buckle spontaneously in mammalian membranes. In a survey of electron cryo-microscopy (cryo-EM) and atomic force microscopy images, we identify key intermediates along both the efflux and buckling pathways to pore formation, as seen in the simulations.

Curvature induction and membrane remodeling by FAM134B reticulon homology domain assist selective ER-phagy.

FAM134B/RETREG1 is a selective ER-phagy receptor that regulates the size and shape of the endoplasmic reticulum. The structure of its reticulon-homology domain (RHD), an element shared with other ER-shaping proteins, and the mechanism of membrane shaping remain poorly understood. Using molecular modeling and molecular dynamics (MD) simulations, we assemble a structural model for the RHD of FAM134B. Through MD simulations of FAM134B in flat and curved membranes, we relate the dynamic RHD structure with its two wedge-shaped transmembrane helical hairpins and two amphipathic helices to FAM134B functions in membrane-curvature induction and curvature-mediated protein sorting. FAM134B clustering, as expected to occur in autophagic puncta, amplifies the membrane-shaping effects. Electron microscopy of in vitro liposome remodeling experiments support the membrane remodeling functions of the different RHD structural elements. Disruption of the RHD structure affects selective autophagy flux and leads to disease states.

Cryo-EM structure of the bifunctional secretin complex of Thermus thermophilus.

Secretins form multimeric channels across the outer membrane of Gram-negative bacteria that mediate the import or export of substrates and/or extrusion of type IV pili. The secretin complex of Thermus thermophilus is an oligomer of the 757-residue PilQ protein, essential for DNA uptake and pilus extrusion. Here, we present the cryo-EM structure of this bifunctional complex at a resolution of ~7 Å using a new reconstruction protocol. Thirteen protomers form a large periplasmic domain of six stacked rings and a secretin domain in the outer membrane. A homology model of the PilQ protein was fitted into the cryo-EM map. A crown-like structure outside the outer membrane capping the secretin was found not to be part of PilQ. Mutations in the secretin domain disrupted the crown and abolished DNA uptake, suggesting a central role of the crown in natural transformation.

Carbon Nanotubes Mediate Fusion of Lipid Vesicles.

The fusion of lipid membranes is opposed by high energetic barriers. In living organisms, complex protein machineries carry out this biologically essential process. Here we show that membrane-spanning carbon nanotubes (CNTs) can trigger spontaneous fusion of small lipid vesicles. In coarse-grained molecular dynamics simulations, we find that a CNT bridging between two vesicles locally perturbs their lipid structure. Their outer leaflets merge as the CNT pulls lipids out of the membranes, creating an hourglass-shaped fusion intermediate with still intact inner leaflets. As the CNT moves away from the symmetry axis connecting the vesicle centers, the inner leaflets merge, forming a pore that completes fusion. The distinct mechanism of CNT-mediated membrane fusion may be transferable, providing guidance in the development of fusion agents, e.g., for the targeted delivery of drugs or nucleic acids.

Structural and mechanistic insights into human splicing factor SF3b complex derived using an integrated approach guided by the cryo-EM density maps.

Pre-mRNA splicing in eukaryotes is performed by the spliceosome, a highly complex macromolecular machine. SF3b is a multi-protein complex which recognizes the branch point adenosine of pre-mRNA as part of a larger U2 snRNP or U11/U12 di-snRNP in the dynamic spliceosome machinery. Although a cryo-EM map is available for human SF3b complex, the structure and relative spatial arrangement of all components in the complex are not yet known. We have recognized folds of domains in various proteins in the assembly and generated comparative models. Using an integrative approach involving structural and other experimental data, guided by the available cryo-EM density map, we deciphered a pseudo-atomic model of the closed form of SF3b which is found to be a "fuzzy complex" with highly flexible components and multiplicity of folds. Further, the model provides structural information for 5 proteins (SF3b10, SF3b155, SF3b145, SF3b130 and SF3b14b) and localization information for 4 proteins (SF3b10, SF3b145, SF3b130 and SF3b14b) in the assembly for the first time. Integration of this model with the available U11/U12 di-snRNP cryo-EM map enabled elucidation of an open form. This now provides new insights on the mechanistic features involved in the transition between closed and open forms pivoted by a hinge region in the SF3b155 protein that also harbors cancer causing mutations. Moreover, the open form guided model of the 5' end of U12 snRNA, which includes the branch point duplex, shows that the architecture of SF3b acts as a scaffold for U12 snRNA: pre-mRNA branch point duplex formation with potential implications for branch point adenosine recognition fidelity.

Protein flexibility in the light of structural alphabets.

Protein structures are valuable tools to understand protein function. Nonetheless, proteins are often considered as rigid macromolecules while their structures exhibit specific flexibility, which is essential to complete their functions. Analyses of protein structures and dynamics are often performed with a simplified three-state description, i.e., the classical secondary structures. More precise and complete description of protein backbone conformation can be obtained using libraries of small protein fragments that are able to approximate every part of protein structures. These libraries, called structural alphabets (SAs), have been widely used in structure analysis field, from definition of ligand binding sites to superimposition of protein structures. SAs are also well suited to analyze the dynamics of protein structures. Here, we review innovative approaches that investigate protein flexibility based on SAs description. Coupled to various sources of experimental data (e.g., B-factor) and computational methodology (e.g., Molecular Dynamic simulation), SAs turn out to be powerful tools to analyze protein dynamics, e.g., to examine allosteric mechanisms in large set of structures in complexes, to identify order/disorder transition. SAs were also shown to be quite efficient to predict protein flexibility from amino-acid sequence. Finally, in this review, we exemplify the interest of SAs for studying flexibility with different cases of proteins implicated in pathologies and diseases.

CLAP: a web-server for automatic classification of proteins with special reference to multi-domain proteins.

BACKGROUND: The function of a protein can be deciphered with higher accuracy from its structure than from its amino acid sequence. Due to the huge gap in the available protein sequence and structural space, tools that can generate functionally homogeneous clusters using only the sequence information, hold great importance. For this, traditional alignment-based tools work well in most cases and clustering is performed on the basis of sequence similarity. But, in the case of multi-domain proteins, the alignment quality might be poor due to varied lengths of the proteins, domain shuffling or circular permutations. Multi-domain proteins are ubiquitous in nature, hence alignment-free tools, which overcome the shortcomings of alignment-based protein comparison methods, are required. Further, existing tools classify proteins using only domain-level information and hence miss out on the information encoded in the tethered regions or accessory domains. Our method, on the other hand, takes into account the full-length sequence of a protein, consolidating the complete sequence information to understand a given protein better. RESULTS: Our web-server, CLAP (Classification of Proteins), is one such alignment-free software for automatic classification of protein sequences. It utilizes a pattern-matching algorithm that assigns local matching scores (LMS) to residues that are a part of the matched patterns between two sequences being compared. CLAP works on full-length sequences and does not require prior domain definitions.Pilot studies undertaken previously on protein kinases and immunoglobulins have shown that CLAP yields clusters, which have high functional and domain architectural similarity. Moreover, parsing at a statistically determined cut-off resulted in clusters that corroborated with the sub-family level classification of that particular domain family. CONCLUSIONS: CLAP is a useful protein-clustering tool, independent of domain assignment, domain order, sequence length and domain diversity. Our method can be used for any set of protein sequences, yielding functionally relevant clusters with high domain architectural homogeneity. The CLAP web server is freely available for academic use at http://nslab.mbu.iisc.ernet.in/clap/.

The relationship between classification of multi-domain proteins using an alignment-free approach and their functions: a case study with immunoglobulins.

Establishing functional relationships between multi-domain protein sequences is a non-trivial task. Traditionally, delineating functional assignment and relationships of proteins requires domain assignments as a prerequisite. This process is sensitive to alignment quality and domain definitions. In multi-domain proteins due to multiple reasons, the quality of alignments is poor. We report the correspondence between the classification of proteins represented as full-length gene products and their functions. Our approach differs fundamentally from traditional methods in not performing the classification at the level of domains. Our method is based on an alignment free local matching scores (LMS) computation at the amino-acid sequence level followed by hierarchical clustering. As there are no gold standards for full-length protein sequence classification, we resorted to Gene Ontology and domain-architecture based similarity measures to assess our classification. The final clusters obtained using LMS show high functional and domain architectural similarities. Comparison of the current method with alignment based approaches at both domain and full-length protein showed superiority of the LMS scores. Using this method we have recreated objective relationships among different protein kinase sub-families and also classified immunoglobulin containing proteins where sub-family definitions do not exist currently. This method can be applied to any set of protein sequences and hence will be instrumental in analysis of large numbers of full-length protein sequences.

Accurate prediction of interfacial residues in two-domain proteins using evolutionary information: implications for three-dimensional modeling.

With the preponderance of multidomain proteins in eukaryotic genomes, it is essential to recognize the constituent domains and their functions. Often function involves communications across the domain interfaces, and the knowledge of the interacting sites is essential to our understanding of the structure-function relationship. Using evolutionary information extracted from homologous domains in at least two diverse domain architectures (single and multidomain), we predict the interface residues corresponding to domains from the two-domain proteins. We also use information from the three-dimensional structures of individual domains of two-domain proteins to train naïve Bayes classifier model to predict the interfacial residues. Our predictions are highly accurate (∼85%) and specific (∼95%) to the domain-domain interfaces. This method is specific to multidomain proteins which contain domains in at least more than one protein architectural context. Using predicted residues to constrain domain-domain interaction, rigid-body docking was able to provide us with accurate full-length protein structures with correct orientation of domains. We believe that these results can be of considerable interest toward rational protein and interaction design, apart from providing us with valuable information on the nature of interactions.

Understanding the role of domain-domain linkers in the spatial orientation of domains in multi-domain proteins.

Inter-domain linkers (IDLs)' bridge flanking domains and support inter-domain communication in multi-domain proteins. Their sequence and conformational preferences enable them to carry out varied functions. They also provide sufficient flexibility to facilitate domain motions and, in conjunction with the interacting interfaces, they also regulate the inter-domain geometry (IDG). In spite of the basic intuitive understanding of the inter-domain orientations with respect to linker conformations and interfaces, we still do not entirely understand the precise relationship among the three. We show that IDG is evolutionarily well conserved and is constrained by the domain-domain interface interactions. The IDLs modulate the interactions by varying their lengths, conformations and local structure, thereby affecting the overall IDG. Results of our analysis provide guidelines in modelling of multi-domain proteins from the tertiary structures of constituent domain components.