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.

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.

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.

A conserved degron containing an amphipathic helix regulates the cholesterol-mediated turnover of human squalene monooxygenase, a rate-limiting enzyme in cholesterol synthesis.

Cholesterol biosynthesis in the endoplasmic reticulum (ER) is tightly controlled by multiple mechanisms to regulate cellular cholesterol levels. Squalene monooxygenase (SM) is the second rate-limiting enzyme in cholesterol biosynthesis and is regulated both transcriptionally and post-translationally. SM undergoes cholesterol-dependent proteasomal degradation when cholesterol is in excess. The first 100 amino acids of SM (designated SM N100) are necessary for this degradative process and represent the shortest cholesterol-regulated degron identified to date. However, the fundamental intrinsic characteristics of this degron remain unknown. In this study, we performed a series of deletions, point mutations, and domain swaps to identify a 12-residue region (residues Gln-62-Leu-73), required for SM cholesterol-mediated turnover. Molecular dynamics and circular dichroism revealed an amphipathic helix within this 12-residue region. Moreover, 70% of the variation in cholesterol regulation was dependent on the hydrophobicity of this region. Of note, the earliest known Doa10 yeast degron, Deg1, also contains an amphipathic helix and exhibits 42% amino acid similarity with SM N100. Mutating SM residues Phe-35/Ser-37/Leu-65/Ile-69 into alanine, based on the key residues in Deg1, blunted SM cholesterol-mediated turnover. Taken together, our results support a model whereby the amphipathic helix in SM N100 attaches reversibly to the ER membrane depending on cholesterol levels; with excess, the helix is ejected and unravels, exposing a hydrophobic patch, which then serves as a degradation signal. Our findings shed new light on the regulation of a key cholesterol synthesis enzyme, highlighting the conservation of critical degron features from yeast to humans.

Probing molecular forces in multi-component physiological membranes.

Biological membranes are remarkably heterogeneous, composed of diverse lipid mixtures with distinct chemical structure and composition. By combining molecular dynamics simulations and the newly developed Lipid-Force Distribution Analysis (L-FDA), we explore force transmission in complex multi-component membrane models mimicking eukaryotic organelles. We found that the chemical-moiety based segmentation at membrane interfaces revealed a distinctive distribution of bonded and non-bonded forces in diverse membrane environment. Our molecular stress analysis could have far-reaching implications in describing the relationship between membrane mechanical properties and functional states of chemically distinct lipids.

Decoding Structural Properties of a Partially Unfolded Protein Substrate: En Route to Chaperone Binding.

Many proteins comprising of complex topologies require molecular chaperones to achieve their unique three-dimensional folded structure. The E.coli chaperone, GroEL binds with a large number of unfolded and partially folded proteins, to facilitate proper folding and prevent misfolding and aggregation. Although the major structural components of GroEL are well defined, scaffolds of the non-native substrates that determine chaperone-mediated folding have been difficult to recognize. Here we performed all-atomistic and replica-exchange molecular dynamics simulations to dissect non-native ensemble of an obligate GroEL folder, DapA. Thermodynamics analyses of unfolding simulations revealed populated intermediates with distinct structural characteristics. We found that surface exposed hydrophobic patches are significantly increased, primarily contributed from native and non-native ß-sheet elements. We validate the structural properties of these conformers using experimental data, including circular dichroism (CD), 1-anilinonaphthalene-8-sulfonic acid (ANS) binding measurements and previously reported hydrogen-deutrium exchange coupled to mass spectrometry (HDX-MS). Further, we constructed network graphs to elucidate long-range intra-protein connectivity of native and intermediate topologies, demonstrating regions that serve as central "hubs". Overall, our results implicate that genomic variations (or mutations) in the distinct regions of protein structures might disrupt these topological signatures disabling chaperone-mediated folding, leading to formation of aggregates.

The Molecular Mechanism Underlying Recruitment and Insertion of Lipid-Anchored LC3 Protein into Membranes.

Lipid modification of cytoplasmic proteins initiates membrane engagement that triggers diverse cellular processes. Despite the abundance of lipidated proteins in the human proteome, the key determinants underlying membrane recognition and insertion are poorly understood. Here, we define the course of spontaneous membrane insertion of LC3 protein modified with phosphatidylethanolamine using multiple coarse-grain simulations. The partitioning of the lipid anchor chains proceeds through a concerted process, with its two acyl chains inserting one after the other. Concurrently, a conformational rearrangement involving the α-helix III of LC3, especially in the three basic residues Lys65, Arg68, and Arg69, ensures stable insertion of the phosphatidylethanolamine anchor into membranes. Mutational studies validate the crucial role of these residues, and further live-cell imaging analysis shows a substantial reduction in the formation of autophagic vesicles for the mutant proteins. Our study captures the process of water-favored LC3 protein recruitment to the membrane and thus opens, to our knowledge, new avenues to explore the cellular dynamics underlying vesicular trafficking.

Structure and dynamics of DRD4 bound to an agonist and an antagonist using in silico approaches.

Human dopamine receptor D4 (DRD4), a member of G-protein coupled receptor (GPCR) family, plays a central role in cell signaling and trafficking. Dysfunctional activity of DRD4 can lead to several psychiatric conditions and, therefore, represents target for many neurological disorders. However, lack of atomic structure impairs our understanding of the mechanism regulating its activity. Here, we report the modeled structure of DRD4 alone and in complex with dopamine and spiperone, its natural agonist and antagonist, respectively. To assess the conformational dynamics induced upon ligand binding, all-atom explicit solvent molecular dynamics simulations in membrane environment were performed. Comprehensive analyses of simulations reveal that agonist binding triggers a series of conformational changes in the transmembrane region, including rearrangement of residues, characteristic of transmission and tyrosine toggle molecular switches. Further, the trajectories indicate that a loop region in the intracellular region--ICL3, is significantly dynamic in nature, mainly due to the side-chain movements of conserved proline residues involved in SH3 binding domains. Interestingly, in dopamine-bound receptor simulation, ICL3 represents an open conformation ideal for G protein binding. The structural and dynamical information presented here suggest a mode of activation of DRD4, upon ligand binding. Our study will help in further understanding of receptor activation, as acquiring structural information is crucial for the design of highly selective DRD4 ligands.

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.

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.

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.

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.

Structured pathway across the transition state for peptide folding revealed by molecular dynamics simulations.

Small globular proteins and peptides commonly exhibit two-state folding kinetics in which the rate limiting step of folding is the surmounting of a single free energy barrier at the transition state (TS) separating the folded and the unfolded states. An intriguing question is whether the polypeptide chain reaches, and leaves, the TS by completely random fluctuations, or whether there is a directed, stepwise process. Here, the folding TS of a 15-residue ß-hairpin peptide, Peptide 1, is characterized using independent 2.5 µs-long unbiased atomistic molecular dynamics (MD) simulations (a total of 15 µs). The trajectories were started from fully unfolded structures. Multiple (spontaneous) folding events to the NMR-derived conformation are observed, allowing both structural and dynamical characterization of the folding TS. A common loop-like topology is observed in all the TS structures with native end-to-end and turn contacts, while the central segments of the strands are not in contact. Non-native sidechain contacts are present in the TS between the only tryptophan (W11) and the turn region (P7-G9). Prior to the TS the turn is found to be already locked by the W11 sidechain, while the ends are apart. Once the ends have also come into contact, the TS is reached. Finally, along the reactive folding paths the cooperative loss of the W11 non-native contacts and the formation of the central inter-strand native contacts lead to the peptide rapidly proceeding from the TS to the native state. The present results indicate a directed stepwise process to folding the peptide.

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.

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.

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.

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.

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.