Molecular Dynamics Simulations of Adsorption of SARS-CoV-2 Spike Protein on Polystyrene Surface

Molecular and Energetic Descriptions of the Plasma Protein Adsorption onto the PVC Surface: Implications for Biocompatibility in Medical Devices

Protein adsorption on material surfaces plays a key role in the biocompatibility of medical devices. Therefore, understanding the complex interplay of physicochemical factors driving this kind of biofouling is paramount for advancing biomaterial design. In this study, we investigated the interaction of the most prominent plasma proteins with polyvinyl chloride (PVC) as one of the ubiquitous materials in medical devices. Through molecular docking, we identified human serum albumin (HSA) as a plasma protein with the highest affinity for adsorption onto the PVC surface with the binding energy of -25.9 kJ mol-1. Subsequently, utilizing triplicate molecular dynamics (MD) simulations (0.5 μs each), we quantitatively analyzed the interactions between HSA and PVC, probing potential structural changes in the protein upon adsorption. Our findings revealed that water-mediated hydrogen bonds and van der Waals forces are key contributors in stabilizing HSA onto the surface of PVC without significant alteration to its secondary and tertiary structures. The observed distribution of water molecules further highlights the importance of the hydration layer in facilitating and modulating protein-polymer interactions. We further evaluated the thermodynamic properties governing the adsorption process by calculating the potential of mean force (PMF) along the direction normal to the surface. The computed Gibbs free energy of adsorption at 300 K (-507.4 kJ/mol) indicated a thermodynamically favored and spontaneous process. Moreover, our investigations across different temperatures (290 to 310 K) consistently showed an enthalpy-driven adsorption process.

Adsorption-Driven Deformation and Footprints of the RBD Proteins in SARS-CoV-2 Variants on Biological and Inanimate Surfaces

Respiratory viruses, carried through airborne microdroplets, frequently adhere to surfaces, including plastics and metals. However, our understanding of the interactions between viruses and materials remains limited, particularly in scenarios involving polarizable surfaces. Here, we investigate the role of the receptor-binding domain (RBD) of the spike protein mutations on the adsorption of SARS-CoV-2 to hydrophobic and hydrophilic surfaces employing molecular simulations. To contextualize our findings, we contrast the interactions on inanimate surfaces with those on native biological interfaces, specifically the angiotensin-converting enzyme 2. Notably, we identify a 2-fold increase in structural deformations for the protein's receptor binding motif (RBM) onto inanimate surfaces, indicative of enhanced shock-absorbing mechanisms. Furthermore, the distribution of adsorbed amino acids (landing footprints) on the inanimate surface reveals a distinct regional asymmetry relative to the biological interface, with roughly half of the adsorbed amino acids arranged in opposite sites. In spite of the H-bonds formed at the hydrophilic substrate, the simulations consistently show a higher number of contacts and interfacial area with the hydrophobic surface, where the wild-type RBD adsorbs more strongly than the Delta or Omicron RBDs. In contrast, the adsorption of Delta and Omicron to hydrophilic surfaces was characterized by a distinctive hopping-pattern. The novel shock-absorbing mechanisms identified in the virus adsorption on inanimate surfaces show the embedded high-deformation capacity of the RBD without losing its secondary structure, which could lead to current experimental strategies in the design of virucidal surfaces.

Probing Enzymatic PET Degradation: Molecular Dynamics Analysis of Cutinase Adsorption and Stability

Understanding the mechanisms influencing poly(ethylene terephthalate) (PET) biodegradation is crucial for developing innovative strategies to accelerate the breakdown of this persistent plastic. In this study, we employed all-atom molecular dynamics simulation to investigate the adsorption process of the LCC-ICCG cutinase enzyme onto the PET surface. Our results revealed that hydrophobic, π-π, and H bond interactions, specifically involving aliphatic, aromatic, and polar uncharged amino acids, were the primary driving forces for the adsorption of the cutinase enzyme onto PET. Additionally, we observed a negligible change in the enzyme's tertiary structure during the interaction with PET (RMSD = 1.35 Å), while its secondary structures remained remarkably stable. Quantitative analysis further demonstrated that there is about a 24% decrease in the number of enzyme-water hydrogen bonds upon adsorption onto the PET surface. The significance of this study lies in unraveling the molecular intricacies of the adsorption process, providing valuable insights into the initial steps of enzymatic PET degradation.

Intracellular Protein Adsorption Behavior and Biological Effects of Polystyrene Nanoplastics in THP-1 Cells

Micro(nano)plastics (MNPs) are emerging pollutants that can adsorb pollutants in the environment and biological molecules and ultimately affect human health. However, the aspects of adsorption of intracellular proteins onto MNPs and its biological effects in cells have not been investigated to date. The present study revealed that 100 nm polystyrene nanoplastics (NPs) could be internalized by THP-1 cells and specifically adsorbed intracellular proteins. In total, 773 proteins adsorbed onto NPs with high reliability were identified using the proteomics approach and analyzed via bioinformatics to predict the route and distribution of NPs following cellular internalization. The representative proteins identified via the Kyoto Encyclopedia of Genes and Genomes pathway analysis were further investigated to characterize protein adsorption onto NPs and its biological effects. The analysis revealed that NPs affect glycolysis through pyruvate kinase M (PKM) adsorption, trigger the unfolded protein response through the adsorption of ribophorin 1 (RPN1) and heat shock 70 protein 8 (HSPA8), and are chiefly internalized into cells through clathrin-mediated endocytosis with concomitant clathrin heavy chain (CLTC) adsorption. Therefore, this work provides new insights and research strategies for the study of the biological effects caused by NPs.

Identification and mechanistic exploration of structural and conformational dynamics of NF-kB inhibitors: rationale insights from in silico and in vitro studies

Increased expression of target genes that code for proinflammatory chemical mediators results from a series of intracellular cascades triggered by activation of dysregulated NF-κB signaling pathway. Dysfunctional NF-kB signaling amplifies and perpetuates autoimmune responses in inflammatory diseases, including psoriasis. This study aimed to identify therapeutically relevant NF-kB inhibitors and elucidate the mechanistic aspects behind NF-kB inhibition. After virtual screening and molecular docking, five hit NF-kB inhibitors opted, and their therapeutic efficacy was examined using cell-based assays in TNF-α stimulated human keratinocyte cells. To investigate the conformational changes of target protein and inhibitor-protein interaction mechanisms, molecular dynamics (MD) simulations, binding free energy calculations together with principal component (PC) analysis, dynamics cross-correlation matrix analysis (DCCM), free energy landscape (FEL) analysis and quantum mechanical calculations were carried out. Among identified NF-kB inhibitors, myricetin and hesperidin significantly scavenged intracellular ROS and inhibited NF-kB activation. Analysis of the MD simulation trajectories of ligand-protein complexes revealed that myricetin and hesperidin formed energetically stabilized complexes with the target protein and were able to lock NF-kB in a closed conformation. Myricetin and hesperidin binding to the target protein significantly impacted conformational changes and internal dynamics of amino acid residues in protein domains. Tyr57, Glu60, Lys144 and Asp239 residues majorly contributed to locking the NF-kB in a closed conformation. The combinatorial approach employing in silico tools integrated with cell-based approaches substantiated the binding mechanism and NF-kB active site inhibition by the lead molecule myricetin, which can be explored as a viable antipsoriatic drug candidate associated with dysregulated NF-kB.Communicated by Ramaswamy H. Sarma.

Experimental and in silico evaluations of the possible molecular interaction between airborne particulate matter and SARS-CoV-2

During the early stage of the COVID-19 pandemic (winter 2020), the northern part of Italy has been significantly affected by viral infection compared to the rest of the country leading the scientific community to hypothesize that airborne particulate matter (PM) could act as a carrier for the SARS-CoV-2. To address this controversial issue, we first verified and demonstrated the presence of SARS-CoV-2 RNA genome on PM2.5 samples, collected in the city of Bologna (Northern Italy) in winter 2021. Then, we employed classical molecular dynamics (MD) simulations to investigate the possible recognition mechanism(s) between a newly modelled PM2.5 fragment and the SARS-CoV-2 Spike protein. The potential molecular interaction highlighted by MD simulations suggests that the glycans covering the upper Spike protein regions would mediate the direct contact with the PM2.5 carbon core surface, while a cloud of organic and inorganic PM2.5 components surround the glycoprotein with a network of non-bonded interactions resulting in up to 4769 total contacts. Moreover, a binding free energy of -207.2 ± 3.9 kcal/mol was calculated for the PM-Spike interface through the MM/GBSA method, and structural analyses also suggested that PM attachment does not alter the protein conformational dynamics. Although the association between the PM and SARS-CoV-2 appears plausible, this simulation does not assess whether these established interactions are sufficiently stable to carry the virus in the atmosphere, or whether the virion retains its infectiousness after the transport. While these key aspects should be verified by further experimental analyses, for the first time, this pioneering study gains insights into the molecular interactions between PM and SARS-CoV-2 Spike protein and will support further research aiming at clarifying the possible relationship between PM abundance and the airborne diffusion of viruses.

Probing nanomechanical interactions of SARS-CoV-2 variants Omicron and XBB with common surfaces

The emergence of SARS-CoV-2 variants has further raised concerns about viral transmission. A fundamental understanding of the intermolecular interactions between the coronavirus and different surfaces is needed to address the transmission of SARS-CoV-2 through respiratory droplet-contaminated surfaces or fomites. The receptor-binding domain (RBD) of the spike protein is a key target for the adhesion of SARS-CoV-2 on the surface. To understand the effect of mutations on adhesion, atomic force microscopy-based single-molecule force spectroscopy (AFM-SMFS) was used to quantify the interactions between wild-type, Omicron, and XBB with several surfaces. The measurement revealed that RBD exhibits relatively higher forces on paper and gold surfaces, with the average force being 1.5 times greater compared to that on plastic surface. In addition, the force elevation on paper and gold surfaces for the variants can reach ∼28% relative to the wild type. These findings enhance our understanding of the nanomechanical interactions of the virus on common surfaces.

Impact of E484Q and L452R Mutations on Structure and Binding Behavior of SARS-CoV-2 B.1.617.1 Using Deep Learning AlphaFold2, Molecular Docking and Dynamics Simulation

During the outbreak of COVID-19, many SARS-CoV-2 variants presented key amino acid mutations that influenced their binding abilities with angiotensin-converting enzyme 2 (hACE2) and neutralizing antibodies. For the B.1.617 lineage, there had been fears that two key mutations, i.e., L452R and E484Q, would have additive effects on the evasion of neutralizing antibodies. In this paper, we systematically investigated the impact of the L452R and E484Q mutations on the structure and binding behavior of B.1.617.1 using deep learning AlphaFold2, molecular docking and dynamics simulation. We firstly predicted and verified the structure of the S protein containing L452R and E484Q mutations via the AlphaFold2-calculated pLDDT value and compared it with the experimental structure. Next, a molecular simulation was performed to reveal the structural and interaction stabilities of the S protein of the double mutant variant with hACE2. We found that the double mutations, L452R and E484Q, could lead to a decrease in hydrogen bonds and higher interaction energy between the S protein and hACE2, demonstrating the lower structural stability and the worse binding affinity in the long dynamic evolutional process, even though the molecular docking showed the lower binding energy score of the S1 RBD of the double mutant variant with hACE2 than that of the wild type (WT) with hACE2. In addition, docking to three approved neutralizing monoclonal antibodies (mAbs) showed a reduced binding affinity of the double mutant variant, suggesting a lower neutralization ability of the mAbs against the double mutant variant. Our study helps lay the foundation for further SARS-CoV-2 studies and provides bioinformatics and computational insights into how the double mutations lead to immune evasion, which could offer guidance for subsequent biomedical studies.

Investigation of the effects of N-Acetylglucosamine on the stability of the spike protein in SARS-CoV-2 by molecular dynamics simulations

A lot of effort has been made in developing vaccine and therapeutic agents against the SARS-CoV-2, concentrating on the Spike protein that binds angiotensin-converting enzyme 2 on human cells. Nowadays, some researches study the role of the N-linked glycans as potential targets for vaccines and new agents. Due to the flexibility and diversity of the N-linked glycans, in this work, we focus on the N-Acetylglucosamine moiety, which is the precursor of nearly all eukaryotic glycans. We performed molecular dynamics simulations to study the effects of the N-Acetylglucosamine on the stability of the spike glycoprotein in SARS-CoV-2. After a 100 ns of simulation on the spike proteins without and with the N-Acetylglucosamine molecules, we found that the presence of N-Acetylglucosamine increases the local stability in their vicinity; even though their effect on the full structure is negligible. Thus; it can be inferred that the N-Acetylglucosamine moieties can potentially affect the interaction of the S protein with the ACE2 receptor. We also found that the S1 domain is more flexible than the S2 domain. We propose which of the experimentally observed glycans found on the spike may be more functional than the others. Detailed understanding of glycans is key for the development of new therapeutic strategies.

The influence of single-point mutation D614G on the binding process between human angiotensin-converting enzyme 2 and the SARS-CoV-2 spike protein-an atomistic simulation study

SARS-CoV-2 has continuously evolved as changes in the genetic code occur during replication of the genome, with some of the mutations leading to higher transmission among human beings. The spike aspartic acid-614 to glycine (D614G) substitution in the spike represents a "more transmissible form of SARS-CoV-2" and occurs in all SARS-CoV-2 mutants. However, the underlying mechanism of the D614G substitution in virus infectivity has remained unclear. In this paper, we adopt molecular simulations to study the contact processes of the D614G mutant and wild-type (WT) spikes with hACE2. The interaction areas with hACE2 for the two spikes are completely different by visualizing the whole binding processes. The D614G mutant spike moves towards the hACE2 faster than the WT spike. We have also found that the receptor-binding domain (RBD) and N-terminal domain (NTD) of the D614G mutant extend more outwards than those of the WT spike. By analyzing the distances between the spikes and hACE2, the changes of number of hydrogen bonds and interaction energy, we suggest that the increased infectivity of the D614G mutant is not possibly related to the binding strength, but to the binding velocity and conformational change of the mutant spike. This work reveals the impact of D614G substitution on the infectivity of the SARS-CoV-2, and hopefully could provide a rational explanation of interaction mechanisms for all the SARS-CoV-2 mutants.

Effect of surfactants on SARS-CoV-2: Molecular dynamics simulations

Surfactants are commonly used as disinfection agents in personal care products against bacteria and viruses, including SARS-CoV-2. However, there is a lack of understanding of the molecular mechanisms of the inactivation of viruses by surfactants. Here, we employ coarse grain (CG) and all-atom (AA) molecular dynamics simulations to investigate the interaction between general families of surfactants and the SARS-CoV-2 virus. To this end, we considered a CG model of a full virion. Overall, we found that surfactants have only a small impact on the virus envelope, being inserted into the envelope without dissolving it or generating pores, at the conditions considered here. However, we found that surfactants may induce a deep impact on the spike protein of the virus (responsible for its infectivity), easily covering it and inducing its collapse over the envelope surface of the virus. AA simulations confirmed that both negatively and positively charged surfactants are able to extensively adsorb over the spike protein and get inserted into the virus envelope. Our results suggest that the best strategy for the design of surfactants as virucidal agents will be to focus on those strongly interacting with the spike protein.

On the Sensitivity and Affinity of Gold, Silver, and Platinum Surfaces against the SARS-CoV-2 Virus: A Comparative Computational Study

The novel coronavirus disease and its complications have motivated the design of new sensors with the highest sensitivity, and affinity for the detection of the SARS-CoV-2 virus is considered in many research studies. In this research article, we employ full atomistic molecular dynamics (MD) models to study the interactions between the receptor binding domain (RBD) and spike protein of the coronavirus and different metals such as gold (Au), platinum (Pt), and silver (Ag) to analyze their sensitivity against this virus. The comparison between the RBD interactions with ACE2 (angiotensin-converting enzyme 2) and different metals indicates that metals have remarkable effects on the structural features and dynamical properties of the RBD. The binding site of the RBD has more affinity to the surfaces of gold, platinum, and silver than to the other parts of the protein. Moreover, the initial configuration of the RBD relative to the metal surface plays an important role in the stability of metal complexes with the RBD. The binding face of the protein to the metal surface has been changed in the presence of different metals. In other words, the residues of the RBD that participate in RBD interactions with the metals are different irrespective of the initial configurations in which the [Asn, Thr, Tyr], [Ser, Thr, Tyr], and [Asn, Asp, Tyr] residues of the protein have a greater affinity to Ag, Au, and Pt, respectively. The corresponding metals have a considerable affinity to the RBD, which due to strong interactions with the protein can change the secondary structure and structural features. Based on the obtained results during the complexation process between the protein and metals, the helical structure of the protein changes to the bend and antiparallel β-sheets. The calculated binding energies for the RBD complexes with silver, gold, and platinum are -95.03, -138.03, and -133.96 kcal·mol-1, respectively. The adsorption process of the spike protein on the surfaces of different metals represents similar results and indicates that the entire spike protein of the coronavirus forms a more stable complex with the gold surface compared with other metals. Moreover, the RBD of the spike protein has more interactions with the surfaces than with the other parts of the protein. Therefore, it is possible to predict the properties of the coronavirus on the metal surface based on the dynamical behavior of the RBD. Overall, our computational results confirm that the gold surface can be considered as an outstanding substrate for developing new sensors with the highest sensitivity against SARS-CoV-2.