Computer simulations on oxidative stress-induced reactions in SARS-CoV-2 spike glycoprotein: a multi-scale approach.

Oxidative stress, which occurs when an organism is exposed to an adverse stimulus that results in a misbalance of antioxidant and pro-oxidants species, is the common denominator of diseases considered as a risk factor for SARS-CoV-2 lethality. Indeed, reactive oxygen species caused by oxidative stress have been related to many virus pathogenicity. In this work, simulations have been performed on the receptor binding domain of SARS-CoV-2 spike glycoprotein to study what residues are more susceptible to be attacked by ·OH, which is one of the most reactive radicals associated to oxidative stress. The results indicate that isoleucine (ILE) probably plays a crucial role in modification processes driven by radicals. Accordingly, QM/MM-MD simulations have been conducted to study both the ·OH-mediated hydrogen abstraction of ILE residues and the induced modification of the resulting ILE radical through hydroxylation or nitrosylation reactions. All in all, in silico studies show the importance of the chemical environment triggered by oxidative stress on the modifications of the virus, which is expected to help for foreseeing the identification or development of antioxidants as therapeutic drugs.

In silico study of substrate chemistry effect on the tethering of engineered antibodies for SARS-CoV-2 detection: Amorphous silica vs gold.

The influence of the properties of different solid substrates on the tethering of two antibodies, IgG1-CR3022 and IgG1-S309, which were specifically engineered for the detection of SARS-CoV-2, has been examined at the molecular level using conventional and accelerated Molecular Dynamics (cMD and aMD, respectively). Two surfaces with very different properties and widely used in immunosensors for diagnosis, amorphous silica and the most stable facet of the face-centered cubic gold structure, have been considered. The effects of such surfaces on the structure and orientation of the immobilized antibodies have been determined by quantifying the tilt and hinge angles that describe the orientation and shape of the antibody, respectively, and the dihedrals that measure the relative position of the antibody arms with respect to the surface. Results show that the interactions with amorphous silica, which are mainly electrostatic due to the charged nature of the surface, help to preserve the orientation and structure of the antibodies, especially of the IgG1-CR3022, indicating that the primary sequence of those antibodies also plays some role. Instead, short-range van der Waals interactions with the inert gold surface cause a higher degree tilting and fraying of the antibodies with respect to amorphous silica. The interactions between the antibodies and the surface also affect the correlation among the different angles and dihedrals, which increases with their strength. Overall, results explain why amorphous silica substrates are frequently used to immobilize antibodies in immunosensors.

Unravelling the molecular interactions between the SARS-CoV-2 RBD spike protein and various specific monoclonal antibodies.

Vaccination against SARS-CoV-2 just started in most of the countries. However, the development of specific vaccines against SARS-CoV-2 is not the only approach to control the virus and monoclonal antibodies (mAbs) start to merit special attention as a therapeutic option to treat COVID-19 disease. Here, the main conformations and interactions between the receptor-binding domain (RBD) of spike glycoprotein of SARS-CoV-2 (S protein) with two mAbs (CR3022 and S309) and the ACE2 cell receptor are studied as the main representatives of three different epitopes on the RBD of S protein. The combined approach of 1 µs accelerated molecular dynamics (aMD) and ab-initio hybrid molecular dynamics is used to identify the most predominant interactions under physiological conditions. Results allow to determine the main receptor-binding mapping, hydrogen bonding network and salt bridges in the most populated antigen-antibody interface conformations. The deep knowledge on the protein-protein interactions involving mAbs and ACE2 receptor with the spike glycoprotein of SARS-CoV-2 increases background knowledge to speed up the development of new vaccines and therapeutic drugs.

IgG1-b12-HIV-gp120 Interface in Solution: A Computational Study.

The use of broadly neutralizing antibodies against human immunodeficiency virus type 1 (HIV-1) has been shown to be a promising therapeutic modality in the prevention of HIV infection. Understanding the b12-gp120 binding mechanism under physiological conditions may assist the development of more broadly effective antibodies. In this work, the main conformations and interactions between the receptor-binding domain (RBD) of spike glycoprotein gp120 of HIV-1 and the IgG1-b12 mAb are studied. Accelerated molecular dynamics (aMD) and ab initio hybrid molecular dynamics have been combined to determine the most persistent interactions between the most populated conformations of the antibody-antigen complex under physiological conditions. The results show the most persistent receptor-binding mapping in the conformations of the antibody-antigen interface in solution. The binding-free-energy decomposition reveals a small enhancement in the contribution played by the CDR-H3 region to the b12-gp120 interface compared to the crystal structure.

In silico antibody engineering for SARS-CoV-2 detection.

Engineered immunoglobulin-G molecules (IgGs) are of wide interest for the development of detection elements in protein-based biosensors with clinical applications. The strategy usually employed for the de novo design of such engineered IgGs consists on merging fragments of the three-dimensional structure of a native IgG, which is immobilized on the biosensor surface, and of an antibody with an exquisite target specificity and affinity. In this work conventional and accelerated classical molecular dynamics (cMD and aMD, respectively) simulations have been used to propose two IgG-like antibodies for COVID-19 detection. More specifically, the crystal structure of the IgG1 B12 antibody, which inactivates the human immunodeficiency virus-1, has been merged with the structure of the antibody CR3022 Fab tightly bounded to SARS-CoV-2 receptor-binding domain (RBD) and the structure of the S309 antibody Fab fragment complexed with SARS-CoV-2 RBD. The two constructed antibodies, named IgG1-CR3022 and IgG1-S309, respectively, have been immobilized on a stable gold surface through a linker. Analyses of the influence of both the merging strategy and the substrate on the stability of the two constructs indicate that the IgG1-S309 antibody better preserves the neutralizing structure than the IgG1-CR3022 one. Overall, results indicate that the IgG1-S309 is appropriated for the generation of antibody based sensors for COVID-19 diagnosis.

Temperature effect on the SARS-CoV-2: A molecular dynamics study of the spike homotrimeric glycoprotein.

Rapid spread of SARS-CoV-2 virus have boosted the need of knowledge about inactivation mechanisms to minimize the impact of COVID-19 pandemic. Recent studies have shown that SARS-CoV-2 virus can be disabled by heating, the exposure time for total inactivation depending on the reached temperature (e.g. more than 45 min at 329 K or less than 5 min at 373 K. In spite of recent crystallographic structures, little is known about the molecular changes induced by the temperature. Here, we unravel the molecular basis of the effect of the temperature over the SARS-CoV-2 spike glycoprotein, which is a homotrimer with three identical monomers, by executing atomistic molecular dynamics (MD) simulations at 298, 310, 324, 338, 358 and 373 K. Furthermore, both the closed down and open up conformational states, which affect the accessibility of receptor binding domain, have been considered. Our results suggest that the spike homotrimer undergoes drastic changes in the topology of the hydrogen bonding interactions and important changes on the secondary structure of the receptor binding domain (RBD), while electrostatic interactions (i.e. salt bridges) are mainly preserved. The proposed inactivation mechanism has important implications for engineering new approaches to fight the SARS-CoV-2 coronavirus, as for example, cleaving or reorganizing the hydrogen bonds through chaotropic agents or nanoparticles with local surface resonant plasmon effect.

Tethering of the IgG1 Antibody to Amorphous Silica for Immunosensor Development: A Molecular Dynamics Study.

A key factor for improving the sensitivity and performance of immunosensors based on mechanical-plasmonic methods is the orientation of the antibody proteins immobilized on the inorganic surface. Although experimental techniques fail to determine surface phenomena at the molecular level, modern simulations open the possibility for improving our understanding of protein-surface interactions. In this work, replica exchange molecular dynamics (REMD) simulations have been used to model the IgG1 protein tethered onto the amorphous silica surface by considering a united-atom model and a relatively large system (2500 nm2 surface). Additional molecular dynamics (MD) simulations have been conducted to derive an atomistic model for the amorphous silica surface using the cristobalite crystal structure as a starting point and to examine the structure of the free IgG1 antibody in the solution for comparison when immobilized. Analyses of the trajectories obtained for the tethered IgG1, which was sampled considering 32 different temperatures, have been used to define the geometry of the protein with respect to the inorganic surface. The tilt angle of the protein with respect to the surface plane increases with temperature, the most populated values being 24, 66, and 87° at the lowest (250 K), room (298 K), and the highest (380 K) temperatures. This variation indicates that the importance of protein-surface interactions decreases with increasing temperature. The influence of the surface on the structure of the antibody is very significant in the constant region, which is directly involved in the tethering process, while it is relatively unimportant for the antigen-binding fragments, which are farthest from the surface. These results are expected to contribute to the development of improved mechanical-plasmonic sensor microarrays in the near future.

Amyloid-like Fibrils from a Diphenylalanine Capped with an Aromatic Fluorenyl.

The self-assembly behavior of a diphenylalanine amphiphile blocked at the C-terminus with a 9-fluorenylmethyl ester and stabilized at the N-terminus with a trifluoroacetate (TFA) anion, TFA·FF-OFm, has been examined. At low peptide concentration (0.5 mg/mL), long amyloid-like fibrils, which come from the fusion of two or more helical ribbons and/or thinner fibrils, organized in bundles or as individual entities are detected. Microbeam synchrotron radiation infrared spectroscopy has shown that TFA·FF-OFm molecules in amyloid-like fibrils arrange, forming antiparallel ß-sheets. Alteration of the experimental conditions to prioritize the thermodynamic contribution with respect to the kinetic one in the self-assembly process inhibits the organization of amyloid-like structures in favor of the formation of conventional fibrous structures. On the basis of experimental observations, a structural model where the individual antiparallel ß-sheets are oriented in parallel has been proposed for TFA·FF-OFm amyloid-like fibrils.