N501Y mutation of spike protein in SARS-CoV-2 strengthens its binding to receptor ACE2

SARS-CoV-2 is spreading around the world for the past year. Recently, several variants such as B.1.1.7 (Alpha), B.1.351 (Beta), and P.1 (Gamma), sharing a key mutation N501Y on the RBD, appear to be more infectious to humans. To understand the underlying mechanism, we performed cell surface binding assay, kinetics study, single-molecule technique, and computational method to investigate the interaction between these RBD (mutations) and ACE2. Remarkably, RBD with the N501Y mutation exhibited a considerably stronger interaction, with a faster association rate and slower dissociation rate. Consistently, atomic force microscopy-based single-molecule force microscopy quantifies their strength showing a higher binding probability and unbinding force for the mutation. Molecular dynamics simulations of RBD-ACE2 complexes indicated that the N501Y introduced additional π-π and π-cation interaction for the higher force/interaction. Taken together, we suggested that the reinforced interaction from N501Y mutation in RBD should play an essential role in the higher transmission of SARS-CoV-2 variants and future mutations in the RBD of the virus should be under surveillance.

Molecular Mechanics Study of Flow and Surface Influence in Ligand–Protein Association

Ligand–protein association is the first and critical step for many biological and chemical processes. This study investigated the molecular association processes under different environments. In biology, cells have different compartments where ligand–protein binding may occur on a membrane. In experiments involving ligand–protein binding, such as the surface plasmon resonance and continuous flow biosynthesis, a substrate flow and surface are required in experimental settings. As compared with a simple binding condition, which includes only the ligand, protein, and solvent, the association rate and processes may be affected by additional ligand transporting forces and other intermolecular interactions between the ligand and environmental objects. We evaluated these environmental factors by using a ligand xk263 binding to HIV protease (HIVp) with atomistic details. Using Brownian dynamics simulations, we modeled xk263 and HIVp association time and probability when a system has xk263 diffusion flux and a non-polar self-assembled monolayer surface. We also examined different protein orientations and accessible surfaces for xk263. To allow xk263 to access to the dimer interface of immobilized HIVp, we simulated the system by placing the protein 20Å above the surface because immobilizing HIVp on a surface prevented xk263 from contacting with the interface. The non-specific interactions increased the binding probability while the association time remained unchanged. When the xk263 diffusion flux increased, the effective xk263 concentration around HIVp, xk263–HIVp association time and binding probability decreased non-linearly regardless of interacting with the self-assembled monolayer surface or not. The work sheds light on the effects of the solvent flow and surface environment on ligand–protein associations and provides a perspective on experimental design.

New insights into the microscopic interactions associated with the physical mechanism of action of highly diluted biologics

In this investigation, we report the effect on the microscopic dynamics and interactions of the cytokine interferon gamma (IFN-g) and antibodies to IFN-g (anti-IFN-g) and to the interferon gamma receptor 1 (anti-IFNGR1) prepared in exceptionally dilute solutions of initial proteins. Using both THz spectroscopy and molecular dynamics (MD) simulations we have uncovered that the high dilution method of sample preparation results in the reorganization of the sample surface residue dynamics at the solvent-protein interface that leads to both structural and kinetic heterogeneous dynamics that ultimately create interactions that enhance the binding probability of the antigen binding site. Our results indicate that the modified interfacial dynamics of anti-IFN-g and anti-IFGNR1 that we probe experimentally are directly associated with alterations in the complementarity regions of the distinct antibodies that designate both antigen-antibody affinity and recognition.

An evolutionary analysis of the SARS-CoV-2 genomes from the countries in the same meridian (Preprint)

UNSTRUCTURED In the current study we analyzed the genomes of SARS-CoV-2 strains isolated from Italy, Sweden, Congo (countries in the same meridian) and Brazil, as outgroup country. Evolutionary analysis revealed codon 9628 under episodic selective pressure for all four countries, suggesting it as a key site for the virus evolution. Belonging to the P0DTD3 (Y14_SARS2) uncharacterized protein 14, further investigation has been conducted showing the codon mutation as responsible for the helical modification in the secondary structure. According to the predictions done, the codon is placed into the more ordered region of the gene (41-59) and close the area acting as transmembrane (54-67), suggesting its involvement into the attachment phase of the virus. The predicted structures of P0DTD3 mutated and not confirmed the importance of the codon to define the protein structure and the ontological analysis of the protein emphasized that the mutation enhances the binding probability.

An evolutionary analysis of the SARS-CoV-2 genomes from the countries in the same meridian

In the current study we analyzed the genomes of SARS-CoV-2 strains isolated from Italy, Sweden, Congo (countries in the same meridian) and Brazil, as outgroup country. Evolutionary analysis revealed codon 9628 under episodic selective pressure for all four countries, suggesting it as a key site for the virus evolution. Belonging to the P0DTD3 (Y14_SARS2) uncharacterized protein 14, further investigation has been conducted showing the codon mutation as responsible for the helical modification in the secondary structure. According to the predictions done, the codon is placed into the more ordered region of the gene (41-59) and close the area acting as transmembrane (54-67), suggesting its involvement into the attachment phase of the virus. The predicted structures of P0DTD3 mutated and not confirmed the importance of the codon to define the protein structure and the ontological analysis of the protein emphasized that the mutation enhances the binding probability.

Nucleosome allostery in pioneer transcription factor binding

While recent experiments revealed that some pioneer transcription factors (TFs) can bind to their target DNA sequences inside a nucleosome, the binding dynamics of their target recognitions are poorly understood. Here we used the latest coarse-grained models and molecular dynamics simulations to study the nucleosome-binding procedure of the two pioneer TFs, Sox2 and Oct4. In the simulations for a strongly positioning nucleosome, Sox2 selected its target DNA sequence only when the target was exposed. Otherwise, Sox2 entropically bound to the dyad region nonspecifically. In contrast, Oct4 plastically bound on the nucleosome mainly in two ways. First, the two POU domains of Oct4 separately bound to the two parallel gyres of the nucleosomal DNA, supporting the previous experimental results of the partial motif recognition. Second, the POUSdomain of Oct4 favored binding on the acidic patch of histones. Then, simulating the TFs binding to a genomic nucleosome, theLIN28Bnucleosome, we found that the recognition of a pseudo motif by Sox2 induced the local DNA bending and shifted the population of the rotational position of the nucleosomal DNA. The redistributed DNA phase, in turn, changed the accessibility of a distant TF binding site, which consequently affected the binding probability of a second Sox2 or Oct4. These results revealed a nucleosomal DNA-mediated allosteric mechanism, through which one TF binding event can change the global conformation, and effectively regulate the binding of another TF at distant sites. Our simulations provide insights into the binding mechanism of single and multiple TFs on the nucleosome.

Length of Mucin-Like Domains Enhance Cell-Ebola Virus Adhesion by Increasing Binding Probability

The Ebola virus (EBOV) hijacks normal physiological processes by apoptotic mimicry in order to be taken up by the cell it infects. The initial adhesion of the virus to the cell is based on the interaction between T-cell immunoglobulin and mucin domain protein, TIM, on the cell-surface and phosphatidylserine (PS) on the viral outer surface. Therefore, it is important to understand the interaction between EBOV/PS and TIM, with selective blocking of the interaction as a potential therapy. Recent experimental studies have shown that for TIM-dependent EBOV entry, a Mucin-like Domain (MLD) with a length of at least 120 amino acids is required, possibly due to the increase of area of the PS-coated surface sampled. We examine this hypothesis by modeling the process of TIM-PS adhesion using a coarse-grained molecular model. We find that the strength of bound PS−TIM pairs is essentially independent of TIM length. TIMs with longer MLDs have higher average binding strengths because of an increase in the probability of binding between EBOV and TIM proteins. Similarly, we find that for larger persistence length (less flexible) the average binding force decreases, again because of a reduction in the probability of binding.Statement of SignificanceThis work studies the mechanism of TIM-dependent adhesion of the Ebola virus to a cell. Through coarse grained modeling we show that longer TIM stalks adhere more easily as they can sample a larger area, thus offering a mechanistic interpretation of an experimental finding. Better mechanistic understanding can lead to therapeutic ideas for blocking adhesion.