GRB2 is a BECN1 interacting protein that regulates autophagy.

GRB2 is an adaptor protein of HER2 (and several other tyrosine kinases), which we identified as a novel BECN1 (Beclin 1) interacting partner. GRB2 co-immunoprecipitated with BECN1 in several breast cancer cell lines and regulates autophagy through a mechanism involving the modulation of the class III PI3Kinase VPS34 activity. In ovo studies in a CAM (Chicken Chorioallantoic Membrane) model indicated that GRB2 knockdown, as well as overexpression of GRB2 loss-of-function mutants (Y52A and S86A-R88A) compromised tumor growth. These differences in tumor growth correlated with differential autophagy activity, indicating that autophagy effects might be related to the effects on tumorigenesis. Our data highlight a novel function of GRB2 as a BECN1 binding protein and a regulator of autophagy.

Membrane association of the ATG8 conjugation machinery emerges as a key regulatory feature for autophagosome biogenesis.

Autophagy is a highly conserved intracellular pathway that is essential for survival in all eukaryotes. In healthy cells, autophagy is used to remove damaged intracellular components, which can be as simple as unfolded proteins or as complex as whole mitochondria. Once the damaged component is captured, the autophagosome engulfs it and closes, isolating the content from the cytoplasm. The autophagosome then fuses with the late endosome and/or lysosome to deliver its content to the lysosome for degradation. Formation of the autophagosome, sequestration or capture of content, and closure all require the ATG proteins, which constitute the essential core autophagy protein machinery. This brief 'nutshell' will highlight recent data revealing the importance of small membrane-associated domains in the ATG proteins. In particular, recent findings from two parallel studies reveal the unexpected key role of α-helical structures in the ATG8 conjugation machinery and ATG8s. These studies illustrate how unique membrane association modules can control the formation of autophagosomes.

ARL8B mediates lipid droplet contact and delivery to lysosomes for lipid remobilization.

Lipid droplets (LDs) play a crucial role in maintaining cellular lipid balance by storing and delivering lipids as needed. However, the intricate lipolytic pathways involved in LD turnover remain poorly described, hindering our comprehension of lipid catabolism and related disorders. Here, we show a function of the small GTPase ARL8B in mediating LD turnover in lysosomes. ARL8B-GDP localizes to LDs, while ARL8-GTP predominantly favors lysosomes. GDP binding induces a conformation with an exposed N-terminal amphipathic helix, enabling ARL8B to bind to LDs. By associating with LDs and lysosomes, and with its property to form a heterotypic complex, ARL8B mediates LD-lysosome contacts and efficient lipid transfer between these organelles. In human macrophages, this ARL8B-dependent LD turnover mechanism appears as the major lipolytic pathway. Our finding opens exciting possibilities for understanding the molecular mechanisms underlying LD degradation and its potential implications for inflammatory disorders.

AI-based AlphaFold2 significantly expands the structural space of the autophagy pathway.

ABBREVIATIONS: AF2: AlphaFold2; AF2-Mult: AlphaFold2 multimer; ATG: autophagy-related; CTD: C-terminal domain; ECTD: extreme C-terminal domain; FR: flexible region; MD: molecular dynamics; NTD: N-terminal domain; pLDDT: predicted local distance difference test; UBL: ubiquitin-like.

Autophagosome membrane expansion is mediated by the N-terminus and cis-membrane association of human ATG8s.

Autophagy is an essential catabolic pathway which sequesters and engulfs cytosolic substrates via autophagosomes, unique double-membraned structures. ATG8 proteins are ubiquitin-like proteins recruited to autophagosome membranes by lipidation at the C-terminus. ATG8s recruit substrates, such as p62, and play an important role in mediating autophagosome membrane expansion. However, the precise function of lipidated ATG8 in expansion remains obscure. Using a real-time in vitro lipidation assay, we revealed that the N-termini of lipidated human ATG8s (LC3B and GABARAP) are highly dynamic and interact with the membrane. Moreover, atomistic MD simulation and FRET assays indicate that N-termini of LC3B and GABARAP associate in cis on the membrane. By using non-tagged GABARAPs, we show that GABARAP N-terminus and its cis-membrane insertion are crucial to regulate the size of autophagosomes in cells irrespectively of p62 degradation. Our study provides fundamental molecular insights into autophagosome membrane expansion, revealing the critical and unique function of lipidated ATG8.

Cyclin-Dependent Kinase 5 Regulates cPLA2 Activity and Neuroinflammation in Parkinson's Disease.

Hyperactivation of cyclin-dependent kinase 5 (Cdk5) by p25, contributes to neuroinflammation causing neurodegeneration in Parkinson's disease (PD) and Alzheimer's disease. However, the mechanism by which Cdk5 induces neuroinflammation in the PD brain is largely unexplored. Here, we show that Cdk5 phosphorylates cytosolic phospholipase A2 (cPLA2) at Thr-268 and Ser-505 sites lead to its activation and generation of eicosanoid products. Mutational studies using site-directed mutagenesis and molecular simulations show that the architecture of the protein changes on each single-point mutation. Interestingly, double mutations also led to a severe decline in the activity of cPLA2 and to the disruption of its translocation to the plasma membrane. Further, the brain lysates of transgenic PD mouse models show hyperactivation of Cdk5, resulting in enhanced phosphorylation of Thr-268 and Ser-505 of cPLA2 and its heightened activity, confirming the findings observed in the cell culture model of PD. These phosphorylation sites of cPLA2 and Cdk5 could be explored as the future therapeutic targets against neuroinflammation in PD. Further, conjoint transcriptomic analysis of the publicly available human PD datasets strengthens the hypothesis that genes of the arachidonic acid, prostaglandin synthesis, and inflammatory pathways are significantly upregulated in the case of PD patients compared with that of healthy control subjects.

Altered TMPRSS2 usage by SARS-CoV-2 Omicron impacts infectivity and fusogenicity.

The SARS-CoV-2 Omicron BA.1 variant emerged in 20211 and has multiple mutations in its spike protein2. Here we show that the spike protein of Omicron has a higher affinity for ACE2 compared with Delta, and a marked change in its antigenicity increases Omicron's evasion of therapeutic monoclonal and vaccine-elicited polyclonal neutralizing antibodies after two doses. mRNA vaccination as a third vaccine dose rescues and broadens neutralization. Importantly, the antiviral drugs remdesivir and molnupiravir retain efficacy against Omicron BA.1. Replication was similar for Omicron and Delta virus isolates in human nasal epithelial cultures. However, in lung cells and gut cells, Omicron demonstrated lower replication. Omicron spike protein was less efficiently cleaved compared with Delta. The differences in replication were mapped to the entry efficiency of the virus on the basis of spike-pseudotyped virus assays. The defect in entry of Omicron pseudotyped virus to specific cell types effectively correlated with higher cellular RNA expression of TMPRSS2, and deletion of TMPRSS2 affected Delta entry to a greater extent than Omicron. Furthermore, drug inhibitors targeting specific entry pathways3 demonstrated that the Omicron spike inefficiently uses the cellular protease TMPRSS2, which promotes cell entry through plasma membrane fusion, with greater dependency on cell entry through the endocytic pathway. Consistent with suboptimal S1/S2 cleavage and inability to use TMPRSS2, syncytium formation by the Omicron spike was substantially impaired compared with the Delta spike. The less efficient spike cleavage of Omicron at S1/S2 is associated with a shift in cellular tropism away from TMPRSS2-expressing cells, with implications for altered pathogenesis.

Spatio-temporal dynamics of intra-host variability in SARS-CoV-2 genomes.

During the course of the COVID-19 pandemic, large-scale genome sequencing of SARS-CoV-2 has been useful in tracking its spread and in identifying variants of concern (VOC). Viral and host factors could contribute to variability within a host that can be captured in next-generation sequencing reads as intra-host single nucleotide variations (iSNVs). Analysing 1347 samples collected till June 2020, we recorded 16 410 iSNV sites throughout the SARS-CoV-2 genome. We found ∼42% of the iSNV sites to be reported as SNVs by 30 September 2020 in consensus sequences submitted to GISAID, which increased to ∼80% by 30th June 2021. Following this, analysis of another set of 1774 samples sequenced in India between November 2020 and May 2021 revealed that majority of the Delta (B.1.617.2) and Kappa (B.1.617.1) lineage-defining variations appeared as iSNVs before getting fixed in the population. Besides, mutations in RdRp as well as RNA-editing by APOBEC and ADAR deaminases seem to contribute to the differential prevalence of iSNVs in hosts. We also observe hyper-variability at functionally critical residues in Spike protein that could alter the antigenicity and may contribute to immune escape. Thus, tracking and functional annotation of iSNVs in ongoing genome surveillance programs could be important for early identification of potential variants of concern and actionable interventions.

A rigorous framework for detecting SARS-CoV-2 spike protein mutational ensemble from genomic and structural features.

The recent release of SARS-CoV-2 genomic data from several countries has provided clues into the potential antigenic drift of the coronavirus population. In particular, the genomic instability observed in the spike protein necessitates immediate action and further exploration in the context of viral-host interactions. By temporally tracking 527,988 SARS-CoV-2 genomes, we identified invariant and hypervariable regions within the spike protein. We evaluated combination of mutations from SARS-CoV-2 lineages and found that maximum number of lineage-defining mutations were present in the N-terminal domain (NTD). Based on mutant 3D-structural models of known Variants of Concern (VOCs), we found that structural properties such as accessibility, secondary structural type, and intra-protein interactions at local mutation sites are greatly altered. Further, we observed significant differences between intra-protein networks of wild-type and Delta mutant, with the latter showing dense intra-protein contacts. Extensive molecular dynamics simulations of D614G mutant spike structure with hACE2 further revealed dynamic features with 47.7% of mutations mapping on flexible regions of spike protein. Thus, we propose that significant changes within spike protein structure have occurred that may impact SARS-CoV-2 pathogenesis, and repositioning of vaccine candidates is required to contain the spread of COVID-19 pathogen.

Genomic characterization and epidemiology of an emerging SARS-CoV-2 variant in Delhi, India.

Delhi, the national capital of India, experienced multiple severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) outbreaks in 2020 and reached population seropositivity of >50% by 2021. During April 2021, the city became overwhelmed by COVID-19 cases and fatalities, as a new variant, B.1.617.2 (Delta), replaced B.1.1.7 (Alpha). A Bayesian model explains the growth advantage of Delta through a combination of increased transmissibility and reduced sensitivity to immune responses generated against earlier variants (median estimates: 1.5-fold greater transmissibility and 20% reduction in sensitivity). Seropositivity of an employee and family cohort increased from 42% to 87.5% between March and July 2021, with 27% reinfections, as judged by increased antibody concentration after a previous decline. The likely high transmissibility and partial evasion of immunity by the Delta variant contributed to an overwhelming surge in Delhi.

Clinical, Serological, Whole Genome Sequence Analyses to Confirm SARS-CoV-2 Reinfection in Patients From Mumbai, India.

Background: SARS-CoV-2 infection may not provide long lasting post-infection immunity. While hundreds of reinfections have reported only a few have been confirmed. Whole genome sequencing (WGS) of the viral isolates from the different episodes is mandatory to establish reinfection. Methods: Nasopharyngeal (NP), oropharyngeal (OP) and whole blood (WB) samples were collected from paired samples of four individuals who were suspected of SARS-CoV-2 reinfection based on distinct clinical episodes and RT-PCR tests. Details from their case record files and investigations were documented. RNA was extracted from the NP and OP samples and subjected to WGS, and the nucleotide and amino acid sequences were subjected to genome and protein-based functional annotation analyses. Serial serology was performed for Anti-N IgG, Anti- S1 RBD IgG, and sVNT (surrogate virus neutralizing test). Findings: Three patients were more symptomatic with lower Ct values and longer duration of illness. Seroconversion was detected soon after the second episode in three patients. WGS generated a genome coverage ranging from 80.07 to 99.7%. Phylogenetic analysis revealed sequences belonged to G, GR and "Other" clades. A total of 42mutations were identified in all the samples, consisting of 22 non-synonymous, 17 synonymous, two in upstream, and one in downstream regions of the SARS-CoV-2 genome. Comparative genomic and protein-based annotation analyses revealed differences in the presence and absence of specific mutations in the virus sequences from the two episodes in all four paired samples. Interpretation: Based on the criteria of genome variations identified by whole genome sequencing and supported by clinical presentation, molecular and serological tests, we were able to confirm reinfections in two patients, provide weak evidence of reinfection in the third patient and unable to rule out a prolonged infection in the fourth. This study emphasizes the importance of detailed analyses of clinical and serological information as well as the virus's genomic variations while assessing cases of SARS-CoV-2 reinfection.

Integrated genomic view of SARS-CoV-2 in India.

Background: India first detected SARS-CoV-2, causal agent of COVID-19 in late January 2020, imported from Wuhan, China. From March 2020 onwards, the importation of cases from countries in the rest of the world followed by seeding of local transmission triggered further outbreaks in India. Methods: We used ARTIC protocol-based tiling amplicon sequencing of SARS-CoV-2 (n=104) from different states of India using a combination of MinION and MinIT sequencing from Oxford Nanopore Technology to understand how introduction and local transmission occurred. Results: The analyses revealed multiple introductions of SARS-CoV-2 genomes, including the A2a cluster from Europe and the USA, A3 cluster from Middle East and A4 cluster (haplotype redefined) from Southeast Asia (Indonesia, Thailand and Malaysia) and Central Asia (Kyrgyzstan). The local transmission and persistence of genomes A4, A2a and A3 was also observed in the studied locations. The most prevalent genomes with patterns of variance (confined in a cluster) remain unclassified, and are here proposed as A4-clade based on its divergence within the A cluster. Conclusions: The viral haplotypes may link their persistence to geo-climatic conditions and host response. Multipronged strategies including molecular surveillance based on real-time viral genomic data is of paramount importance for a timely management of the pandemic.

Human LC3 and GABARAP subfamily members achieve functional specificity via specific structural modulations.

Autophagy is a conserved adaptive cellular pathway essential to maintain a variety of physiological functions. Core components of this machinery are the six human Atg8 orthologs that initiate formation of appropriate protein complexes. While these proteins are routinely used as indicators of autophagic flux, it is presently not possible to discern their individual biological functions due to our inability to predict specific binding partners. In our attempts towards determining downstream effector functions, we developed a computational pipeline to define structural determinants of human Atg8 family members that dictate functional diversity. We found a clear evolutionary separation between human LC3 and GABARAP subfamilies and also defined a novel sequence motif responsible for their specificity. By analyzing known protein structures, we observed that functional modules or microclusters reveal a pattern of intramolecular network, including distinct hydrogen bonding of key residues (F52/Y49; a subset of HP2) that may directly modulate their interaction preferences. Multiple molecular dynamics simulations were performed to characterize how these proteins interact with a common protein binding partner, PLEKHM1. Our analysis showed remarkable differences in binding modes via intrinsic protein dynamics, with PLEKHM1-bound GABARAP complexes showing less fluctuations and higher number of contacts. We further mapped 373 genomic variations and demonstrated that distinct cancer-related mutations are likely to lead to significant structural changes. Our findings present a quantitative framework to establish factors underlying exquisite specificity of human Atg8 proteins, and thus facilitate the design of precise modulators.Abbreviations: Atg: autophagy-related; ECs: evolutionary constraints; GABARAP: GABA type A receptor-associated protein; HsAtg8: human Atg8; HP: hydrophobic pocket; KBTBD6: kelch repeat and BTB domain containing 6; LIR: LC3-interacting region; MAP1LC3/LC3: microtubule associated protein 1 light chain 3; MD: molecular dynamics; HIV-1 Nef: human immunodeficiency virus type 1 negative regulatory factor; PLEKHM1: pleckstrin homology and RUN domain containing M1; RMSD: root mean square deviation; SQSTM1/p62: sequestosome 1; WDFY3/ALFY: WD repeat and FYVE domain containing 3.

Revealing Conformational Transitions in G-Protein-Coupled Receptor Rhodopsin upon Phosphorylation.

G-protein-coupled receptors (GPCRs) have evolved as highly specialized cellular machinery that can dictate biological outcomes in response to diverse stimuli. Specifically, they induce multiple pathway responses upon structural perturbations induced at local protein sites. GPCRs utilize a concurrent strategy involving a central transmembrane topology and biochemical modifications for precise functional implementation. However, the specific role of the latter is not known due to the lack of precise probing techniques that can characterize receptor dynamics upon biochemical modifications. Phosphorylation is known to be one of the critical biochemical modifications in GPCRs that aids in receptor desensitization via arrestin binding. Here, we carry out all-atom molecular dynamics simulations of rhodopsin in a membrane environment to study its conformational dynamics induced upon phosphorylation. Interestingly, our comparative analysis of non-phosphorylated and phosphorylated rhodopsin structure demonstrated enhanced receptor stability upon phosphorylation at the C-terminal region that leads to the opening of the extracellular part of the transmembrane helices. In addition, monitoring the distinct number of phosphorylation states showed that having fewer phosphorylated residues does not bring about appropriate conformational changes in the extracellular region. Since phosphorylation results in receptor desensitization and recycling of the ligand, our findings provide significant insights into the conformational dynamics of the mechanism of ligand exit from the receptor.

Myg1 exonuclease couples the nuclear and mitochondrial translational programs through RNA processing.

Semi-autonomous functioning of mitochondria in eukaryotic cell necessitates coordination with nucleus. Several RNA species fine-tune mitochondrial processes by synchronizing with the nuclear program, however the involved components remain enigmatic. In this study, we identify a widely conserved dually localized protein Myg1, and establish its role as a 3'-5' RNA exonuclease. We employ mouse melanoma cells, and knockout of the Myg1 ortholog in Saccharomyces cerevisiae with complementation using human Myg1 to decipher the conserved role of Myg1 in selective RNA processing. Localization of Myg1 to nucleolus and mitochondrial matrix was studied through imaging and confirmed by sub-cellular fractionation studies. We developed Silexoseqencing, a methodology to map the RNAse trail at single-nucleotide resolution, and identified in situ cleavage by Myg1 on specific transcripts in the two organelles. In nucleolus, Myg1 processes pre-ribosomal RNA involved in ribosome assembly and alters cytoplasmic translation. In mitochondrial matrix, Myg1 processes 3'-termini of the mito-ribosomal and messenger RNAs and controls translation of mitochondrial proteins. We provide a molecular link to the possible involvement of Myg1 in chronic depigmenting disorder vitiligo. Our study identifies a key component involved in regulating spatially segregated organellar RNA processing and establishes the evolutionarily conserved ribonuclease as a coordinator of nucleo-mitochondrial crosstalk.

Distinct Intramolecular Hydrogen Bonding Dictates Antimicrobial Action of Membrane-Targeting Amphiphiles.

As mechanisms underpinning the molecular interactions between membrane-targeting antimicrobials and Gram-negative bacterial membranes at atomistic scale remain elusive, we used cholic acid (CA)-derived amphiphiles with different hydrophobicities as model antimicrobials and assessed the effect of their conformational flexibility on antimicrobial activity. Relative to other hydrophobic counterparts, a compound with a hexyl chain (6) showed the strongest binding with the lipopolysaccharide (LPS) of Gram-negative bacterial membranes and acted as an effective antimicrobial. Biomolecular simulations, validated by complementary approaches, revealed that specific intramolecular hydrogen bonding imparts conformationally rigid character to compound 6. This conformational stability of compound 6 allows minimum but specific interactions of the amphiphile with LPS that are a sum of exothermic processes like electrostatic interactions, membrane insertion, and endothermic contributions from disaggregation of LPS. Therefore, our study reveals that a membrane-targeting mechanism with the help of conformationally selective molecules offers a roadmap for developing future therapeutics against bacterial infections.

Deciphering the Role of Intramolecular Networking in Cholic Acid-Peptide Conjugates on the Lipopolysaccharide Surface in Combating Gram-Negative Bacterial Infections.

The presence of lipopolysaccharide and emergence of drug resistance make the treatment of Gram-negative bacterial infections highly challenging. Herein, we present the synthesis and antibacterial activities of cholic acid-peptide conjugates (CAPs), demonstrating that valine-glycine dipeptide-derived CAP 3 is the most effective antimicrobial. Molecular dynamics simulations and structural analysis revealed that a precise intramolecular network of CAP 3 is maintained in the form of evolving edges, suggesting intramolecular connectivity. Further, we found high conformational rigidity in CAP 3 that confers maximum perturbations in bacterial membranes relative to other small molecules. Interestingly, CAP 3-coated catheters did not allow the formation of biofilms in mice, and treatment of wound infections with CAP 3 was able to clear the bacterial infection. Our results demonstrate that molecular conformation and internal connectivity are critical parameters to describe the antimicrobial nature of compounds, and the analysis presented here may serve as a general principle for the design of future antimicrobials.

Cholic Acid-Derived Amphiphile which Combats Gram-Positive Bacteria-Mediated Infections via Disintegration of Lipid Clusters.

Inappropriate and uncontrolled use of antibiotics results in the emergence of antibiotic resistance, thereby threatening the present clinical regimens to treat infectious diseases. Therefore, new antimicrobial agents that can prevent bacteria from developing drug resistance are urgently needed. Selective disruption of bacterial membranes is the most effective strategy for combating microbial infections as accumulation of genetic mutations will not allow for the emergence of drug resistance against these antimicrobials. In this work, we tested cholic acid (CA) derived amphiphiles tethered with different alkyl chains for their ability to combat Gram-positive bacterial infections. In-depth biophysical and biomolecular simulation studies suggested that the amphiphile with a hexyl chain (6) executes more effective interactions with Gram-positive bacterial membranes as compared to other hydrophobic counterparts. Amphiphile 6 is effective against multidrug resistant Gram-positive bacterial strains as well and does not allow the adherence of S. aureus on amphiphile 6 coated catheters implanted in mice. Further, treatment of wound infections with amphiphile 6 clears the bacterial infections. Therefore, the current study presents strategic guidelines in design and development of CA-derived membrane-targeting antimicrobials for Gram-positive bacterial infections.

A novel mutation alters the stability of PapA2 resulting in the complete abrogation of sulfolipids in clinical mycobacterial strains.

The analysis of whole genomes has revealed specific geographical distribution of Mycobacterium tuberculosis (Mtb) strains across the globe suggestive of unique niche dependent adaptive mechanisms. We provide an important correlation of a genome-based mutation to a molecular phenotype across two predominant clinical Mtb lineages of the Indian subcontinent. We have identified a distinct lineage specific mutation-G247C, translating into an alanine-proline conversion in the papA2 gene of Indo-oceanic lineage 1 (L1) Mtb strains, and restoration of cell wall sulfolipids by simple genetic complementation of papA2 from lineage 3 (L3) or from H37Rv (lineage 4-L4) attributed the loss of this glycolipid to this specific mutation in Indo-Oceanic L1 Mtb. The investigation of structure of Mtb PapA2 revealed a distinct nonribosomal peptide synthetase (NRPS) C domain conformation with an unconventional presence of a zinc binding motif. Surprisingly, the A83P mutation did not map to either the catalytic center in the N-terminal subdomain or any of the substrate-binding region of the protein. On the contrary, the inherent ability of mutant PapA2 to form insoluble aggregates and molecular simulations with the wild-type/mutant (Wt/mut) PapA2 purports an important role for the surface associated 83rd residue in protein conformation. This study demonstrates the importance of a critical structural residue in the papA2 protein of Mtb and helps establish a link between observed genomic alteration and its molecular consequence in the successful human pathogen Mtb. Significance We demonstrate the effect of a unique SNP in PapA2 gene of Indo-oceanic Mycobacterium tuberculosis (Mtb) strains leading to the loss of sulfolipid from these strains. By X-ray crystallographic analysis and molecular dynamics (MD) simulations, we show the importance of this residue in the global PapA2 structure. The presence of a Zn atom has not been reported before for this class of proteins. Here, we provide an important link between genomic alteration and its molecular consequence in Mtb highlighting one of the many adaptive mechanisms that have contributed to its success as a human pathogen. A high degree of identity with PapA1, 3, or 4 would help in interpreting the structure of these PapA proteins and other acyl transferases of other biological systems.