Interferon-γ as a Potential Inhibitor of SARS-CoV-2 ORF6 Accessory Protein.

The ORF6 protein of the SARS-CoV-2 virus plays a crucial role in blocking the innate immune response of the infected cells by inhibiting interferon pathways. Additionally, it binds to and immobilises the RAE1 protein on the cytoplasmic membranes, thereby blocking mRNA transport from the nucleus to the cytoplasm. In all these cases, the host cell proteins are tethered by the flexible C-terminus of ORF6. A possible strategy to inhibit the biological activity of ORF6 is to bind its C-terminus with suitable ligands. Our in silico experiments suggest that hIFNγ binds the ORF6 protein with high affinity, thus impairing its interactions with RAE1 and, consequently, its activity in viral invasion. The in vitro studies reported here reveal a shift of the localisation of RAE1 in ORF6 overexpressing cells upon treatment with hIFNγ from predominantly cytoplasmic to mainly nuclear, resulting in the restoration of the export of mRNA from the nucleus. We also explored the expression of GFP in transfected-with-ORF6 cells by means of fluorescence microscopy and qRT-PCR, finding that treatment with hIFNγ unblocks the mRNA trafficking and reinstates the GFP expression level. The ability of the cytokine to block ORF6 is also reflected in minimising its negative effects on DNA replication by reducing accumulated RNA-DNA hybrids. Our results, therefore, suggest hIFNγ as a promising inhibitor of the most toxic SARS-CoV-2 protein.

Insights into the SARS-CoV-2 ORF6 Mechanism of Action.

ORF6 is responsible for suppressing the immune response of cells infected by the SARS-CoV-2 virus. It is also the most toxic protein of SARS-CoV-2, and its actions are associated with the viral pathogenicity. Here, we study in silico and in vitro the structure of the protein, its interaction with RAE1 and the mechanism of action behind its high toxicity. We show both computationally and experimentally that SARS-CoV-2 ORF6, embedded in the cytoplasmic membranes, binds to RAE1 and sequesters it in the cytoplasm, thus depleting its availability in the nucleus and impairing nucleocytoplasmic mRNA transport. This negatively affects the cellular genome stability by compromising the cell cycle progression into the S-phase and by promoting the accumulation of RNA-DNA hybrids. Understanding the multiple ways in which ORF6 affects DNA replication may also have important implications for elucidating the pathogenicity of SARS-CoV-2 and developing therapeutic strategies to mitigate its deleterious effects on host cells.

Immunohistochemical phenotype of colorectal carcinoma in patients with KRAS mutation and mismatch repair status.

INTRODUCTION: Aberrant expression of CK7/CK20/CDX2 is reported in percentage of colorectal carcinomas (CRC).

Heparan Sulfate Facilitates Binding of hIFNγ to Its Cell-Surface Receptor hIFNGR1.

Human interferon-gamma (hIFNγ) is a crucial signaling molecule with an important role in the initialization and development of the immune response of the host. However, its aberrant activity is also associated with the progression of a multitude of autoimmune and other diseases, which determines the need for effective inhibitors of its activity. The development of such treatments requires proper understanding of the interaction of hIFNγ to its cell-surface receptor hIFNGR1. Currently, there is no comprehensive model of the mechanism of this binding process. Here, we employ molecular dynamics simulations to study on a microscopic level the process of hIFNγ-hIFNGR1 complex formation in different scenarios. We find that the two molecules alone fail to form a stable complex, but the presence of heparan-sulfate-like oligosaccharides largely facilitates the process by both demobilizing the highly flexible C-termini of the cytokine and assisting in the proper positioning of its globule between the receptor subunits. An antiproliferative-activity assay on cells depleted from cell-surface heparan sulfate (HS) sulfation together with the phosphorylation levels of the signal transducer and activator of transcription STAT1 confirms qualitatively the simulation-based multistage complex-formation model. Our results reveal the key role of HS and its proteoglycans in all processes involving hIFNγ signalling.

Immunohistochemical expression of CK20, CK7, and CDX2 in colorectal carcinoma in correlation with pathomorphological characteristics.

INTRODUCTION: Colorectal carcinoma is the third most common cancer worldwide. The usual immunophenotype of colorectal adenocarcinoma is CDX2 positive, CK20 positive, and CK7 negative. Aberrant expression is reported in a variety of colorectal carcinomas but its relation to morphological variables and survival data is still unclear.

Molecular Mechanism of the Anti-Inflammatory Action of Heparin.

Our objective is to reveal the molecular mechanism of the anti-inflammatory action of low-molecular-weight heparin (LMWH) based on its influence on the activity of two key cytokines, IFNγ and IL-6. The mechanism of heparin binding to IFNγ and IL-6 and the resulting inhibition of their activity were studied by means of extensive molecular-dynamics simulations. The effect of LMWH on IFNγ signalling inside stimulated WISH cells was investigated by measuring its antiproliferative activity and the translocation of phosphorylated STAT1 in the nucleus. We found that LMWH binds with high affinity to IFNγ and is able to fully inhibit the interaction with its cellular receptor. It also influences the biological activity of IL-6 by binding to either IL-6 or IL-6/IL-6Rα, thus preventing the formation of the IL-6/IL-6Rα/gp130 signalling complex. These findings shed light on the molecular mechanism of the anti-inflammatory action of LMWH and underpin its ability to influence favourably conditions characterised by overexpression of these two cytokines. Such conditions are not only associated with autoimmune diseases, but also with inflammatory processes, in particular with COVID-19. Our results put forward heparin as a promising means for the prevention and suppression of severe CRS and encourage further investigations on its applicability as an anti-inflammatory agent.

Self-Association of Antimicrobial Peptides: A Molecular Dynamics Simulation Study on Bombinin.

Antimicrobial peptides (AMPs) are a diverse group of membrane-active peptides which play a crucial role as mediators of the primary host defense against microbial invasion. Many AMPs are found to be fully or partially disordered in solution and to acquire secondary structure upon interaction with a lipid membrane. Here, we report molecular dynamics simulations studies on the solution behaviour of a specific AMP, bombinin H2. We show that in monomeric form in water solution the peptide is somewhat disordered and preferably adopts a helix-loop-helix conformation. However, when more than a single monomer is placed in the solution, the peptides self-associate in aggregates. Within the aggregate, the peptides provide each other with an amphipathic environment that mimics the water-membrane interface, which allows them to adopt a single-helix structure. We hypothesise that this is the mechanism by which bombinin H2 and, possibly, other small linear AMPs reach the target membrane in a functional folded state and are able to effectively exert their antimicrobial action on it.

Intramolecular Domain Movements of Free and Bound pMHC and TCR Proteins: A Molecular Dynamics Simulation Study.

The interaction of antigenic peptides (p) and major histocompatibility complexes (pMHC) with T-cell receptors (TCR) is one of the most important steps during the immune response. Here we present a molecular dynamics simulation study of bound and unbound TCR and pMHC proteins of the LC13-HLA-B*44:05-pEEYLQAFTY complex to monitor differences in relative orientations and movements of domains between bound and unbound states of TCR-pMHC. We generated local coordinate systems for MHC α1- and MHC α2-helices and the variable T-cell receptor regions TCR Vα and TCR Vß and monitored changes in the distances and mutual orientations of these domains. In comparison to unbound states, we found decreased inter-domain movements in the simulations of bound states. Moreover, increased conformational flexibility was observed for the MHC α2-helix, the peptide, and for the complementary determining regions of the TCR in TCR-unbound states as compared to TCR-bound states.

Molecular modeling of the effects of glycosylation on the structure and dynamics of human interferon-gamma.

Natural hIFNγ is a glycoprotein with two N-glycosylation sites in each monomer chain, which are independently and differentially glycosylated. Although glycosylation is not necessary for the activity of the cytokine, it was proposed that it protects the cytokine from proteolytic degradation and thus extends its circulatory half-life. Here, we report the development of model structures of glycosylated full-length native hIFNγ homodimers. Our aim is to shed light on the mechanism through which glycosylation preserves the integrity of the cytokine molecule. To this end, we employ molecular dynamics simulations to study the interaction of the carbohydrate chains with the receptor-binding sites in the cytokine and with its flexible highly positively charged C-termini. The glycans interact primarily with the globular part of the protein, but also occasionally form contacts with the solvent-exposed and sensitive to proteases C-terminal tails. We show that the glycans restrict the C-termini wagging motion into the solvent, limit their flexibility and keep them closer to the α-helical globule of hIFNγ, thus possibly protecting them from proteolytic processing.

His-FLAG Tag as a Fusion Partner of Glycosylated Human Interferon-Gamma and Its Mutant: Gain or Loss?

In order to obtain glycosylated human interferon-gamma (hIFNγ) and its highly prone to aggregation mutant K88Q, a secretory expression in insect cells was employed. To facilitate recombinant proteins purification, detection, and stability the baculovirus expression vectors were constructed to bear N-terminal His6-FLAG tag. Although the obtained proteins were glycosylated, we found that their biological activity was 100 times lower than expected. Our attempts to recover the biological properties of both proteins by tag removal failed due to enterokinase resistance of the tag. Surprisingly, the tag was easily cleaved when the proteins were expressed in E. coli cells and the tag-free proteins showed fully restored activity. To shed light on this phenomenon we performed molecular dynamics simulations. The latter showed that the tags interact with the receptor binding domains and the flexible C-termini of the fusion proteins thus suppressing their complex formation with the hIFNγ receptor. We hypothesize that in the case of glycosylated proteins the tag/C-terminal interaction positions the FLAG peptide in close proximity to the glycans thus sterically impeding the enterokinase access to its recognition site.

Multistage modeling of protein dynamics with monomeric Myc oncoprotein as an example.

We propose to combine a mean-field approach with all-atom molecular dynamics (MD) into a multistage algorithm that can model protein folding and dynamics over very long time periods yet with atomic-level precision. As an example, we investigate an isolated monomeric Myc oncoprotein that has been implicated in carcinomas including those in colon, breast, and lungs. Under physiological conditions a monomeric Myc is presumed to be an example of intrinsically disordered proteins that pose a serious challenge to existing modeling techniques. We argue that a room-temperature monomeric Myc is in a dynamical state, it oscillates between different conformations that we identify. For this we adopt the Cα backbone of Myc in a crystallographic heteromer as an initial ansatz for the monomeric structure. We construct a multisoliton of the pertinent Landau free energy to describe the Cα profile with ultrahigh precision. We use Glauber dynamics to resolve how the multisoliton responds to repeated increases and decreases in ambient temperature. We confirm that the initial structure is unstable in isolation. We reveal a highly degenerate ground-state landscape, an attractive set towards which Glauber dynamics converges in the limit of vanishing ambient temperature. We analyze the thermal stability of this Glauber attractor using room-temperature molecular dynamics. We identify and scrutinize a particularly stable subset in which the two helical segments of the original multisoliton align in parallel next to each other. During the MD time evolution of a representative structure from this subset, we observe intermittent quasiparticle oscillations along the C-terminal α helix, some of which resemble a translating Davydov's Amide-I soliton. We propose that the presence of oscillatory motion is in line with the expected intrinsically disordered character of Myc.

Bloch spin waves and emergent structure in protein folding with HIV envelope glycoprotein as an example.

We inquire how structure emerges during the process of protein folding. For this we scrutinize collective many-atom motions during all-atom molecular dynamics simulations. We introduce, develop, and employ various topological techniques, in combination with analytic tools that we deduce from the concept of integrable models and structure of discrete nonlinear Schrödinger equation. The example we consider is an α-helical subunit of the HIV envelope glycoprotein gp41. The helical structure is stable when the subunit is part of the biological oligomer. But in isolation, the helix becomes unstable, and the monomer starts deforming. We follow the process computationally. We interpret the evolving structure both in terms of a backbone based Heisenberg spin chain and in terms of a side chain based XY spin chain. We find that in both cases the formation of protein supersecondary structure is akin the formation of a topological Bloch domain wall along a spin chain. During the process we identify three individual Bloch walls and we show that each of them can be modelled with a precision of tenths to several angstroms in terms of a soliton solution to a discrete nonlinear Schrödinger equation.

Spatiotemporal multistage consensus clustering in molecular dynamics studies of large proteins.

The aim of this work is to find semi-rigid domains within large proteins as reference structures for fitting molecular dynamics trajectories. We propose an algorithm, multistage consensus clustering, MCC, based on minimum variation of distances between pairs of Cα-atoms as target function. The whole dataset (trajectory) is split into sub-segments. For a given sub-segment, spatial clustering is repeatedly started from different random seeds, and we adopt the specific spatial clustering with minimum target function: the process described so far is stage 1 of MCC. Then, in stage 2, the results of spatial clustering are consolidated, to arrive at domains stable over the whole dataset. We found that MCC is robust regarding the choice of parameters and yields relevant information on functional domains of the major histocompatibility complex (MHC) studied in this paper: the α-helices and ß-floor of the protein (MHC) proved to be most flexible and did not contribute to clusters of significant size. Three alleles of the MHC, each in complex with ABCD3 peptide and LC13 T-cell receptor (TCR), yielded different patterns of motion. Those alleles causing immunological allo-reactions showed distinct correlations of motion between parts of the peptide, the binding cleft and the complementary determining regions (CDR)-loops of the TCR. Multistage consensus clustering reflected functional differences between MHC alleles and yields a methodological basis to increase sensitivity of functional analyses of bio-molecules. Due to the generality of approach, MCC is prone to lend itself as a potent tool also for the analysis of other kinds of big data.

Geometry Dynamics of α -Helices in Different Class I Major Histocompatibility Complexes.

MHC α-helices form the antigen-binding cleft and are of particular interest for immunological reactions. To monitor these helices in molecular dynamics simulations, we applied a parsimonious fragment-fitting method to trace the axes of the α-helices. Each resulting axis was fitted by polynomials in a least-squares sense and the curvature integral was computed. To find the appropriate polynomial degree, the method was tested on two artificially modelled helices, one performing a bending movement and another a hinge movement. We found that second-order polynomials retrieve predefined parameters of helical motion with minimal relative error. From MD simulations we selected those parts of α-helices that were stable and also close to the TCR/MHC interface. We monitored the curvature integral, generated a ruled surface between the two MHC α-helices, and computed interhelical area and surface torsion, as they changed over time. We found that MHC α-helices undergo rapid but small changes in conformation. The curvature integral of helices proved to be a sensitive measure, which was closely related to changes in shape over time as confirmed by RMSD analysis. We speculate that small changes in the conformation of individual MHC α-helices are part of the intrinsic dynamics induced by engagement with the TCR.

Relative Movements of Domains in Large Molecules of the Immune System.

Molecular dynamics was used to simulate large molecules of the immune system (major histocompatibility complex class I, presented epitope, T-cell receptor, and a CD8 coreceptor.) To characterize the relative orientation and movements of domains local coordinate systems (based on principal component analysis) were generated and directional cosines and Euler angles were computed. As a most interesting result, we found that the presence of the coreceptor seems to influence the dynamics within the protein complex, in particular the relative movements of the two α-helices, Gα1 and Gα2.

Geometric analysis of alloreactive HLA α-helices.

Molecular dynamics (MD) is a valuable tool for the investigation of functional elements in biomolecules, providing information on dynamic properties and processes. Previous work by our group has characterized static geometric properties of the two MHC α-helices comprising the peptide binding region recognized by T cells. We build upon this work and used several spline models to approximate the overall shape of MHC α-helices. We applied this technique to a series of MD simulations of alloreactive MHC molecules that allowed us to capture the dynamics of MHC α-helices' steric configurations. Here, we discuss the variability of spline models underlying the geometric analysis with varying polynomial degrees of the splines.

Finding semirigid domains in biomolecules by clustering pair-distance variations.

Dynamic variations in the distances between pairs of atoms are used for clustering subdomains of biomolecules. We draw on a well-known target function for clustering and first show mathematically that the assignment of atoms to clusters has to be crisp, not fuzzy, as hitherto assumed. This reduces the computational load of clustering drastically, and we demonstrate results for several biomolecules relevant in immunoinformatics. Results are evaluated regarding the number of clusters, cluster size, cluster stability, and the evolution of clusters over time. Crisp clustering lends itself as an efficient tool to locate semirigid domains in the simulation of biomolecules. Such domains seem crucial for an optimum performance of subsequent statistical analyses, aiming at detecting minute motional patterns related to antigen recognition and signal transduction.

Relaxation estimation of RMSD in molecular dynamics immunosimulations.

Molecular dynamics simulations have to be sufficiently long to draw reliable conclusions. However, no method exists to prove that a simulation has converged. We suggest the method of "lagged RMSD-analysis" as a tool to judge if an MD simulation has not yet run long enough. The analysis is based on RMSD values between pairs of configurations separated by variable time intervals Δt. Unless RMSD(Δt) has reached a stationary shape, the simulation has not yet converged.