GAG positioning on IL-1RI; A mechanism regulated by dual effect of glycosylation.

IL-1RI is the signaling receptor for the IL-1 family of cytokines that are involved in establishment of the innate and acquired immune systems. Glycosylated extracellular (EC) domain of the IL-1RI binds to agonist such as IL-1ß or antagonist ligands and the accessory protein to form the functional signaling complex. Dynamics and ligand binding of the IL-1RI is influenced by presence of the glycosaminoglycans (GAGs) of the EC matrix. Here a combination of molecular dockings and molecular dynamics simulations of the unglycosylated, partially N-glycosylated and fully N-glycosylated IL-1RI EC domain in the apo, GAG-bound and IL-1ß-bound states were carried out to explain the co-occurring dynamical effect of receptor's glycosylation and GAGs. It was shown that the IL-1RI adopts two types of "extended" and "locked" conformations in its dynamical pattern, and glycosylation maintains the receptor in the latter form. Maintaining the receptor in the locked conformation disfavors IL-1ß binding by burying its two binding site on the IL-1RI EC domain. Glycosylation disfavors GAG binding to the extended IL-1RI EC domain by sterically limiting the GAGs degrees of freedom in targeting its binding site, while it favors GAG binding to the locked IL-1RI by favorable packing interactions.

Role of N-glycosylation in EGFR ectodomain ligand binding.

The epidermal growth factor receptor (EGFR) is a tyrosine kinase protein, overexpressed in several cancers. The extracellular domain of EGFR is known to be heavily glycosylated. Growth factor (mostly epidermal growth factor or EGF) binding activates EGFR. This occurs by inducing the transition from the autoinhibited tethered conformation to an extended conformation of the monomeric form of EGFR and by stabilizing the flexible preformed dimer. Activated EGFR adopts a back-to-back dimeric conformation after binding of another homologous receptor to its extracellular domain as the dimeric partner. Several antibodies inhibit EGFR by targeting the growth factor binding site or the dimeric interfaces. Glycosylation has been shown to be important for modulating the stability and function of EGFR. Here, atomistic MD simulations show that N-glycosylation of the EGFR extracellular domain plays critical roles in the binding of growth factors, monoclonal antibodies, and the dimeric partners to the monomeric EGFR extracellular domain. N-glycosylation results in the formation of several noncovalent interactions between the glycans and EGFR extracellular domain near the EGF binding site. This stabilizes the growth factor binding site, resulting in stronger interactions (electrostatic) between the growth factor and EGFR. N-glycosylation also helps maintain the dimeric interface and plays distinct roles in binding of antibodies to spatially separated epitopes of the EGFR extracellular domain. Analysis of SNP data suggests the possibility of altered glycosylation with functional consequences. Proteins 2017; 85:1529-1549. © 2017 Wiley Periodicals, Inc.

Computational design and investigation of the monomeric spike SARS-CoV-2-ferritin nanocage vaccine stability and interactions.

Since the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) outbreak, several solutions have been proposed to manage the disease. The most viable option for controlling this virus is to produce effective vaccines. Most of the current SARS-CoV-2 vaccines have focused on the infusion spike protein. Spike exists as a trimer and plays a vital role in infecting host cells by binding to the Angiotensin-Converting Enzyme 2 (ACE2) receptor through its Receptor Binding Domain (RBD). Ferritin protein, a naturally occurring iron-storage protein, has gained attention for vaccine production due to its self-assembling property, non-toxic nature, and biocompatibility. Ferritin nanocages have recently been employed in the development of a SARS-CoV-2 vaccination eliciting not only long-term protective memory cells but also a sustained antibody response. In this study, a combination of in silico investigations including molecular docking, molecular dynamics simulations, and immune simulations were carried out to computationally model the monomeric spike protein on the ferritin nanocage as well as to evaluate its stability and interactions for the first time. The structural dynamics of the modeled complex demonstrated noticeable stability. In particular, the Receptor Binding Domain (RBD) and ferritin within the monomeric spike-ferritin complex illustrated significant stability. The lack of alterations in the secondary structure further supported the overall steadiness of the complex. The decline in the distance between ferritin and spike suggests a strong interaction over time. The cross-correlation matrices revealed that the monomeric spike and ferritin move towards each other supporting the stable interaction between spike and ferritin. Further, the orientation of monomeric spike protein within the ferritin unit facilitated the exposure of critical epitopes, specifically upward active Receptor Binding Domain (RBD), enabling effective interactions with the ACE2 receptor. The immune simulations of the model indicated high-level stimulations of both cellular and humoral immunity in the human body. It was also found that the employed model is effective regardless of the mutated spikes in different variants. These findings shed light on the current status of the SARS-CoV-2-ferritin nanoparticle vaccines and could be used as a framework for other similar vaccine designs.

In-silico structural analysis of Heterocephalus glaber amyloid beta: an anti-Alzheimer's peptide.

Heterocephalus glaber, known as the Naked mole-rat, has an extraordinary immunity to Alzheimer's disease. The pathological hallmark of Alzheimer's disease is cerebral accumulations of plaques, consisting of self-aggregated amyloid beta peptides. Homo sapiens and H. glaber amyloid beta peptides are different in only one amino acid. Herein, computational structural analyses were carried out to determine whether plaque development in H. glaber is prevented by the replacement of His13 with Arg13 in the amyloid beta peptide. AlphaFold2 was used to predict the structure of the H. glaber amyloid beta peptide. HADDOCK and Hex were used to self-dock the peptides and dock ions on peptides, respectively. Illustrations were made by PyMol and ChimeraX. Using VMD, we calculated the radius of gyration. The phylogenetic analysis was conducted by Mega. The results showed an accurate structure with two alpha helices separated by a short coil for H. glaber. Self-docking of the two amyloid beta peptides demonstrated a globular conformation in the H. glaber dimer, implying the unlikeliness of amyloid beta peptides' self-aggregation to form fibrillar structures. This conformational state resulted in lower electrostatic energy compared to H. sapiens, contributing to H. glaber's lower tendency for fibril and, ultimately, plaque formation. Phylogenetic analysis confirmed that amyloid precursor protein is highly conserved in each taxon of rodentia and primata. This study provides insight into the connection between the structure of H. glaber amyloid beta and its plaque formation properties, showing that the Arg13 in H. glaber leads to fibril instability, and might prevent senile plaque accumulation.

Oligomer formation of SARS-CoV-2 ORF8 through 73YIDI76 motifs regulates immune response and non-infusion antiviral interactions.

Introduction: Open Reading Frame 8 (ORF8) is a 121 amino acid length SARS-CoV-2 specific accessory protein that plays crucial roles in viral infectivity, and pathogenesis. Current SARS-CoV-2 treatments focus on spike or RNA-dependent RNA polymerase proteins. Hence, directing attention to ORF8 yields substantial benefits for innovative non-infusional therapeutics. Functional ORF8 is proposed to form oligomers via a crystallographic contact centered by 73YIDI76 motifs. Methods: Hence, the structure and atomistic interactions of trimeric and tetrameric ORF8 oligomeric forms were modeled by means of thorough molecular modeling and molecular dynamics simulations. Results: Results show that trimeric and tetrameric oligomers are stabilized by the interaction of ß4-ß5 (47-83) loops. 73YIDI76 motifs are involved in obtaining the oligomerization interfaces. It is shown that the tetramers which resemble a doughnut-like construction are the most stabilized oligomeric forms. Where four ß4-ß5 loops form the interfaces between two dimers. Each monomer links to two others through ß4-ß5 loops and a covalent Cys20-Cys20 bridge. Epitope mapping, binding site predictions, and solvent-accessible surface area analyses of different ORF8 forms show that the B-cell, MHC-I, and drug epitopes stay exposed in oligomeric forms. Discussion: Approving that the viral infectivity is expanded upon ORF8 oligomerization and the regions involved in oligomerization can be considered as therapeutic targets.

Structural arrangement of the active back-to-back dimer in N-glycosylated ErbB receptors is regulated by heterodimerization.

The human epidermal growth factor receptor (EGFR/ErbB) family consists of four members (ErbB1-4) and belongs to the superfamily of receptor tyrosine kinases (RTKs). The ErbB family members participate in multiple cellular pathways and are the key players in several cancers (brain, breast, lung etc.). Activation of these family members depends on their extracellular domains forming back-to-back hetero/homo dimers. Moreover, dimers are glycosylated, which is a crucial post-translational modification that affects the conformation and function of the protein. Here, molecular modeling and molecular docking are used to comprehensively investigate the dimerization mechanism in glycosylated back-to-back active dimer formation in the entire ErbB receptors for the first time. Results showed that 21 out of 37 clusters of active back-to-back dimers formed by all family members are through heterodimerization. Including; ErbB1-ErbB3/ErbB4, ErbB2-ErbB3/ErbB4 and ErbB3-ErbB4. Ranking ErbB2-ErbB3 as the most stabilized back-to-back dimeric construct. While glycan arrangements favor both homo/hetero dimerization at the dimeric interfaces, it promotes heterodimerization by stabilizing and packing the ligand binding sites of EGFR and ErbB2 respectively. These findings pave the path to future heterodimeric interface/glycan targeting rational anti-cancer drug designs for ErbB receptors.

Glycan-mediated functional assembly of IL-1RI: structural insights into completion of the current description for immune response.

Interleukin 1 Receptor type I (IL-1RI) is a multi-domain transmembrane receptor that triggers the inflammatory response. Understanding its detailed mechanism of action is crucial for treating immune disorders. IL-1RI is activated upon formation of its functional assembly that occurs by binding of the IL-1 cytokine and the accessory protein (Il-1RAcP) to it. X-ray crystallography, small-Angle X-ray Scattering and molecular dynamics simulation studies showed that IL-1RI adopts two types of 'compact' and 'extended' conformational states in its dynamical pattern. Furthermore, glycosylation has shown to play a critical role in its activation process. Here, classical and accelerated atomistic molecular dynamics were carried out to examine the role of full glycosylation of IL-1RI and IL-1RAcP in arrangement of the functional assembly. Simulations showed that the 'compact' and 'extended' IL-1RI form two types of 'cytokine-inaccessible-non-signaling' and 'cytokine-accessible-signaling' assemblies with the IL-1RacP, respectively that are both abiding in the presence of glycans. Suggesting that the cytokine binding to IL-1RI is not required for the formation of IL-1RI-IL-1RAcP complex and the 'compact' complex could act as a down-regulatory mechanism. The 'extended' complex is maintained by formation of several persistent hydrogen bonds between the IL-1RI-IL-1RAcP inter-connected glycans. Taken together, it was shown that full glycosylation regulates formation of the IL-1RI functional assembly and play critical role in cytokine biding and triggering the IL-1RI involved downstream pathways in the cell.Communicated by Ramaswamy H. Sarma.

In silico design of refined ferritin-SARS-CoV-2 glyco-RBD nanoparticle vaccine.

With the onset of Coronavirus disease 2019 (COVID-19) pandemic, all attention was drawn to finding solutions to cure the coronavirus disease. Among all vaccination strategies, the nanoparticle vaccine has been shown to stimulate the immune system and provide optimal immunity to the virus in a single dose. Ferritin is a reliable self-assembled nanoparticle platform for vaccine production that has already been used in experimental studies. Furthermore, glycosylation plays a crucial role in the design of antibodies and vaccines and is an essential element in developing effective subunit vaccines. In this computational study, ferritin nanoparticles and glycosylation, which are two unique facets of vaccine design, were used to model improved nanoparticle vaccines for the first time. In this regard, molecular modeling and molecular dynamics simulation were carried out to construct three atomistic models of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) receptor binding domain (RBD)-ferritin nanoparticle vaccine, including unglycosylated, glycosylated, and modified with additional O-glycans at the ferritin-RBD interface. It was shown that the ferritin-RBD complex becomes more stable when glycans are added to the ferritin-RBD interface and optimal performance of this nanoparticle can be achieved. If validated experimentally, these findings could improve the design of nanoparticles against all microbial infections.

Cancer regulator EGFR-ErbB4 heterodimer is stabilized through glycans at the dimeric interface.

EGFR and ErbB4 are the only two members of cancer-regulating ErbB RTKs that maintain the activation of their extracellular ligand-binding domain and intracellular tyrosine kinase domain. EGFR and ErbB4 could form homo and heterodimers upon their activation. Heterodimerization triggers more diverse intracellular pathways compared to homodimerization. Moreover, it is known that N-glycosylation is crucial for the stabilization and activation of EGFR and ErbB4 receptors. Herein, atomistic molecular dynamics were simulated to study the EGFR-ErbB4 heterodimer in the glycosylated and unglycosylated states. It was shown that the EGFR-ErbB4 heterodimer is highly stabilized by glycosylation. The increased stability is most significant at the dimeric interfaces, regulated by packing of three glycans attached to EGFR (Asn337) and ErbB4 (Asn333, Asn523) at the dimeric interface. Finally, it is proposed that heterodimerization is the persistent key player in the EGFR and ErbB4 activation. Thus, targeting the heterodimers in future therapeutic designs could be a promising approach against drug resistance to ErbB-positive cancers.

Metabolomics analysis of milk thistle lipids to identify drought-tolerant genes.

Milk thistle is an oil and medicinal crop known as an alternative oil crop with a high level of unsaturated fatty acids, which makes it a favorable edible oil for use in food production. To evaluate the importance of Milk thistle lipids in drought tolerance, an experiment was performed in field conditions under three different water deficit levels (Field capacity (FC), 70% FC and 40% FC). After harvesting seeds of the plant, their oily and methanolic extracts were isolated, and subsequently, types and amounts of lipids were measured using GC-MS. Genes and enzymes engaged in biosynthesizing of these lipids were identified and their expression in Arabidopsis was investigated under similar conditions. The results showed that content of almost all measured lipids of milk thistle decreased under severe drought stress, but genes (belonged to Arabidopsis), which were involved in their biosynthetic pathway showed different expression patterns. Genes biosynthesizing lipids, which had significant amounts were selected and their gene and metabolic network were established. Two networks were correlated, and for each pathway, their lipids and respective biosynthesizing genes were grouped together. Four up-regulated genes including PXG3, LOX2, CYP710A1, PAL and 4 down-regulated genes including FATA2, CYP86A1, LACS3, PLA2-ALPHA were selected. The expression of these eight genes in milk thistle was similar to Arabidopsis under drought stress. Thus, PXG3, PAL, LOX2 and CYP86A1 genes that increased expression were selected for protein analysis. Due to the lack of protein structure of these genes in the milk thistle, modeling homology was performed for them. The results of molecular docking showed that the four proteins CYP86A1, LOX2, PAL and PXG3 bind to ligands HEM, 11O, ACT and LIG, respectively. HEM ligand was involved in production of secondary metabolites and dehydration tolerance, and HEM binding site remained conserved in various plants. CA ligands were involved in synthesis of cuticles and waxes. Overall, this study confirmed the importance of lipids in drought stress tolerance in milk thistle.

Identification of the key functional genes in salt-stress tolerance of Cyanobacterium Phormidium tenue using in silico analysis.

The development of artificial biocrust using cyanobacterium Phormidium tenue has been suggested as an effective strategy to prevent soil degradation. Here, a combination of in silico approaches with growth rate, photosynthetic pigment, morphology, and transcript analysis was used to identify specific genes and their protein products in response to 500 mM NaCl in P. tenue. The results show that 500 mM NaCl induces the expression of genes encoding glycerol-3-phosphate dehydrogenase (glpD) as a Flavoprotein, ribosomal protein S12 methylthiotransferase (rimO), and a hypothetical protein (sll0939). The constructed co-expression network revealed a group of abiotic stress-responsive genes. Using the Basic Local Alignment Search Tool (BLAST), the homologous proteins of rimO, glpD, and sll0939 were identified in the P. tenue genome. Encoded proteins of glpD, rimO, and DUF1622 genes, respectively, contain (DAO and DAO C), (UPF0004, Radical SAM and TRAM 2), and (DUF1622) domains. The predicted ligand included 22B and MG for DUF1622, FS5 for rimO, and FAD for glpD protein. There was no direct disruption in ligand-binding sites of these proteins by Na+, Cl-, or NaCl. The growth rate, photosynthetic pigment, and morphology of P. tenue were investigated, and the result showed an acceptable tolerance rate of this microorganism under salt stress. The quantitative real-time polymerase chain reaction (qRT-PCR) results revealed the up-regulation of glpD, rimO, and DUF1622 genes under salt stress. This is the first report on computational and experimental analyses of the glpD, rimO, and DUF1622 genes in P. tenue under salt stress to the best of our knowledge. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s13205-021-03050-w.

Glycosylation promotes the cancer regulator EGFR-ErbB2 heterodimer formation - molecular dynamics study.

ErbB family of receptor tyrosine kinases play significant roles in cellular differentiation and proliferation. Mutation or overexpression of these receptors leads to several cancers in humans. The family has four homologous members including EGFR, ErbB2, ErbB3, and ErbB4. From which all except the ErbB2 bind to growth factors via the extracellular domain to send signals to the cell. However, dimerization of the ErbB receptor occurs in extracellular, transmembrane, and intracellular domains. The ErbB receptors are known to form homodimers and heterodimers in the active form. Heterodimerization increases the variety of identified ligands and signaling pathways that can be activated by these receptors. Furthermore, glycosylation of the ErbB receptors has shown to be critical for their stability, ligand binding, and dimerization. Here, atomistic molecular dynamics simulations on the glycosylated and unglycosylated heterodimer showed that the EGFR-ErbB2 heterodimer is more stable in its dynamical pattern compared to the EGFR-EGFR homodimer. This increased stability is regulated by maintaining the dimeric interface by the attached glycans. It was also shown that the presence of various glycosylation sites within the ErbB2 growth factor binding site leads to occlusion of this site by the glycans that inhibit ligand binding to ErbB2 and participate in further stabilization of the heterodimer construct. Putting together, glycosylation seems to promote the heterodimer formation within the ErbB family members as the dominant molecular mechanism of activation for these receptors.

S494 O-glycosylation site on the SARS-CoV-2 RBD affects the virus affinity to ACE2 and its infectivity; a molecular dynamics study.

SARS-CoV-2 is a strain of Coronavirus family that caused the ongoing pandemic of COVID-19. Several studies showed that the glycosylation of virus spike (S) protein and the Angiotensin-Converting Enzyme 2 (ACE2) receptor on the host cell is critical for the virus infectivity. Molecular Dynamics (MD) simulations were used to explore the role of a novel mutated O-glycosylation site (D494S) on the Receptor Binding Domain (RBD) of S protein. This site was suggested as a key mediator of virus-host interaction. By exploring the dynamics of three O-glycosylated models and the control systems of unglcosylated S4944 and S494D complexes, it was shown that the decoration of S494 with elongated O-glycans results in stabilized interactions on the direct RBD-ACE2. Calculation of the distances between RBD and two major H1, H2 helices of ACE2 and the interacting pairs of amino acids in the interface showed that the elongated O-glycan maintains these interactions by forming several polar contacts with the neighbouring residues while it would not interfere in the direct binding interface. Relative binding free energy of RBD-ACE2 is also more favorable in the O-glycosylated models with longer glycans. The increase of RBD binding affinity to ACE2 depends on the size of attached O-glycan. By increasing the size of O-glycan, the RBD-ACE2 binding affinity will increase. Hence, this crucial factor must be taken into account for any further inhibitory approaches towards RBD-ACE2 interaction.