Shedding Light on the Inhibitory Mechanisms of SARS-CoV-1/CoV-2 Spike Proteins by ACE2-Designed Peptides

Uncovering the Association Mechanism between Two Intrinsically Flexible Proteins

The understanding of protein-protein interaction mechanisms is key to the atomistic description of cell signaling pathways and for the development of new drugs. In this context, the mechanism of intrinsically disordered proteins folding upon binding has attracted attention. The VirB9 C-terminal domain (VirB9Ct) and the VirB7 N-terminal motif (VirB7Nt) associate with VirB10 to form the outer membrane core complex of the Type IV Secretion System injectisome. Despite forming a stable and rigid complex, VirB7Nt behaves as a random coil, while VirB9Ct is intrinsically dynamic in the free state. Here we combined NMR, stopped-flow fluorescence, and computer simulations using structure-based models to characterize the VirB9Ct-VirB7Nt coupled folding and binding mechanism. Qualitative data analysis suggested that VirB9Ct preferentially binds to VirB7Nt by way of a conformational selection mechanism at lower temperatures. However, at higher temperatures, energy barriers between different VirB9Ct conformations are more easily surpassed. Under these conditions the formation of non-native initial encounter complexes may provide alternative pathways toward the native complex conformation. These observations highlight the intimate relationship between folding and binding, calling attention to the fact that the two molecular partners must search for the most favored intramolecular and intermolecular interactions on a rugged and funnelled conformational energy landscape, along which multiple intermediates may lead to the final native state.

Coordinate-Dependent Drift-Diffusion Reveals the Kinetic Intermediate Traps of Top7-Based Proteins

The computer-designed Top7 served as a scaffold to produce immunoreactive proteins by grafting of the 2F5 HIV-1 antibody epitope (Top7-2F5) followed by biotinylation (Top7-2F5-biotin). The resulting nonimmunoglobulin affinity proteins were effective in inducing and detecting the HIV-1 antibody. However, the grafted Top7-2F5 design led to protein aggregation, as opposed to the soluble biotinylated Top7-2F5-biotin. The structure-based model predicted that the thermodynamic cooperativity of Top7 increases after grafting and biotin-labeling, reducing their intermediate state populations. In this work, the folding kinetic traps that might contribute to the aggregation propensity are investigated by the diffusion theory. Since the engineered proteins have similar sequence and structural homology, they served as protein models to study the kinetic intermediate traps that were uncovered by characterizing the position-dependent drift-velocity (v(Q)) and the diffusion (D(Q)) coefficients. These coordinate-dependent coefficients were taken into account to obtain the folding and transition path times over the free energy transition states containing the intermediate kinetic traps. This analysis may be useful to predict the aggregated kinetic traps of scaffold-epitope proteins that might compose novel diagnostic and therapeutic platforms.

Biotinylation Eliminates the Intermediate State of Top7 Designed with an HIV-1 Epitope

Broadly neutralizing antibodies against HIV-1 are rare with the 2F5 antibody being one of the most protective. Insertion of an antibody epitope into a stable and small protein scaffold overcomes many of the obstacles found to produce antibodies. However, the design leads to grafting of epitopes that may cause protein aggregation. Here, I investigated the 2F5 epitope grafted into the Top7 as the scaffold in which the resulting immunoreactive protein precipitates along the storage time, as opposed to its completely soluble biotinylated version. Molecular dynamics showed that biotinylation eliminates the intermediate state of the scaffold-epitope Top7-2F5 by switching a noncooperative to a cooperative folding. The aggregation propensity of the Top7-designed proteins is examined in light of thermodynamic cooperativity and kinetic traps along the decreasing depth of the intermediate ensemble in the free energy landscape. This protocol may predict stable and soluble scaffold-epitopes with the purpose of composing novel therapeutic and diagnostic platforms.

Computation of the Binding Energies between Human ACE2 and Spike RBDs of the Original Strain, Delta and Omicron Variants of the SARS-CoV-2: A DFT Simulation Approach

The receptor binding domain (RBD) of SARS-CoV-2 binds to human ACE2 leading to infection. In this study, the complexes that are formed by the attachment of the SARS-CoV-2 spike RBDs of the original strain, delta and omicron variants to the human ACE2 are investigated via density functional theory (DFT) simulations to obtain binding energies. The DFT computations are performed without fragmenting the interfaces to involve longer-range interactions for improved accuracy, which is one of the primary features of the approach used in this study. Basis set superposition error corrections and van der Waals dispersions are also included in the DFT simulations. The binding energies of the SARS-CoV-2 spike RBDs of the original strain, delta and omicron variants to the human ACE2 are computed as -4.76, -6.68, and -11.77 eV, respectively. These binding energy values indicate that the binding of the omicron variant to the ACE2 is much more favorable than the binding of the original strain and the delta variant, which constitute a molecular reason for the takeover of the omicron variant. The binding energies and the decomposition of these energies found in this study are expected to aid in the development of neutralizing agents.

The Impact of Angiotensin-Converting Enzyme-2/Angiotensin 1-7 Axis in Establishing Severe COVID-19 Consequences

The COVID-19 pandemic has put a tremendous stress on the medical community over the last two years. Managing the infection proved a lot more difficult after several research communities started to recognize the long-term effects of this disease. The cellular receptor for the virus was identified as angiotensin-converting enzyme-2 (ACE2), a molecule responsible for a wide array of processes, broadly variable amongst different organs. Angiotensin (Ang) 1-7 is the product of Ang II, a decaying reaction catalysed by ACE2. The effects observed after altering the level of ACE2 are essentially related to the variation of Ang 1-7. The renin-angiotensin-aldosterone system (RAAS) is comprised of two main branches, with ACE2 representing a crucial component of the protective part of the complex. The ACE2/Ang (1-7) axis is well represented in the testis, heart, brain, kidney, and intestine. Infection with the novel SARS-CoV-2 virus determines downregulation of ACE2 and interrupts the equilibrium between ACE and ACE2 in these organs. In this review, we highlight the link between the local effects of RAAS and the consequences of COVID-19 infection as they arise from observational studies.

Characterization of binding interactions of SARS-CoV-2 spike protein and DNA-peptide nanostructures

Binding interactions of the spike proteins of the severe acute respiratory syndrome corona virus 2 (SARS-CoV-2) to a peptide fragment derived from the human angiotensin converting enzyme 2 (hACE2) receptor are investigated. The peptide is employed as capture moiety in enzyme linked immunosorbent assays (ELISA) and quantitative binding interaction measurements that are based on fluorescence proximity sensing (switchSENSE). In both techniques, the peptide is presented on an oligovalent DNA nanostructure, in order to assess the impact of mono- versus trivalent binding modes. As the analyte, the spike protein and several of its subunits are tested as well as inactivated SARS-CoV-2 and pseudo viruses. While binding of the peptide to the full-length spike protein can be observed, the subunits RBD and S1 do not exhibit binding in the employed concentrations. Variations of the amino acid sequence of the recombinant full-length spike proteins furthermore influence binding behavior. The peptide was coupled to DNA nanostructures that form a geometric complement to the trimeric structure of the spike protein binding sites. An increase in binding strength for trimeric peptide presentation compared to single peptide presentation could be generally observed in ELISA and was quantified in switchSENSE measurements. Binding to inactivated wild type viruses could be shown as well as qualitatively different binding behavior of the Alpha and Beta variants compared to the wild type virus strain in pseudo virus models.

Exploring the Spike-hACE 2 Residue-Residue Interaction in Human Coronaviruses SARS-CoV-2, SARS-CoV, and HCoV-NL63

Coronaviruses (CoVs) have been responsible for three major outbreaks since the beginning of the 21st century, and the emergence of the recent COVID-19 pandemic has resulted in considerable efforts to design new therapies against coronaviruses. Thus, it is crucial to understand the structural features of their major proteins related to the virus-host interaction. Several studies have shown that from the seven known CoV human pathogens, three of them use the human Angiotensin-Converting Enzyme 2 (hACE-2) to mediate their host's cell entry: SARS-CoV-2, SARS-CoV, and HCoV-NL63. Therefore, we employed quantum biochemistry techniques within the density function theory (DFT) framework and the molecular fragmentation with conjugate caps (MFCC) approach to analyze the interactions between the hACE-2 and the spike protein-RBD of the three CoVs in order to map the hot-spot residues that form the recognition surface for these complexes and define the similarities and differences in the interaction scenario. The total interaction energy evaluated showed a good agreement with the experimental binding affinity order: SARS-2 > SARS > NL63. A detailed investigation revealed the energetically most relevant regions of hACE-2 and the spike protein for each complex, as well as the key residue-residue interactions. Our results provide valuable information to deeply understand the structural behavior and binding site characteristics that could help to develop antiviral therapeutics that inhibit protein-protein interactions between CoVs S protein and hACE-2.

Computational Analysis of the SARS-CoV-2 RBD-ACE2-Binding Process Based on MD and the 3D-RISM Theory

The binding process of angiotensin-converting enzyme 2 (ACE2) to the receptor-binding domain (RBD) of the severe acute respiratory syndrome-like coronavirus 2 spike protein was investigated using molecular dynamics simulation and the three-dimensional reference interaction-site model theory. The results suggested that the protein-binding process consists of a protein-protein approaching step, followed by a local structural rearrangement step. In the approaching step, the interprotein interaction energy decreased as the proteins approached each other, whereas the solvation free energy increased. As the proteins approached, the glycan of ACE2 first established a hydrogen bond with the RBD. Thereafter, the number of interprotein hydrogen bonds increased rapidly. The solvation free energy increased because of the desolvation of the protein as it approached its partner. The spatial distribution function of the solvent revealed the presence of hydrogen bonds bridged by water molecules on the RBD-ACE2 interface. Finally, principal component analysis revealed that ACE2 showed a pronounced conformational change, whereas there was no significant change in RBD.

Biotin-painted proteins have thermodynamic stability switched by kinetic folding routes

Biotin-labeled proteins are widely used as tools to study protein-protein interactions and proximity in living cells. Proteomic methods broadly employ proximity-labeling technologies based on protein biotinylation in order to investigate the transient encounters of biomolecules in subcellular compartments. Biotinylation is a post-translation modification in which the biotin molecule is attached to lysine or tyrosine residues. So far, biotin-based technologies proved to be effective instruments as affinity and proximity tags. However, the influence of biotinylation on aspects such as folding, binding, mobility, thermodynamic stability, and kinetics needs to be investigated. Here, we selected two proteins [biotin carboxyl carrier protein (BCCP) and FKBP3] to test the influence of biotinylation on thermodynamic and kinetic properties. Apo (without biotin) and holo (biotinylated) protein structures were used separately to generate all-atom structure-based model simulations in a wide range of temperatures. Holo BCCP contains one biotinylation site, and FKBP3 was modeled with up to 23 biotinylated lysines. The two proteins had their estimated thermodynamic stability changed by altering their energy landscape. In all cases, after comparison between the apo and holo simulations, differences were observed on the free-energy profiles and folding routes. Energetic barriers were altered with the density of states clearly showing changes in the transition state. This study suggests that analysis of large-scale datasets of biotinylation-based proximity experiments might consider possible alterations in thermostability and folding mechanisms imposed by the attached biotins.

Structure/Function Analysis of Truncated Amino-Terminal ACE2 Peptide Analogs That Bind to SARS-CoV-2 Spike Glycoprotein

The global burden of the SARS-CoV-2 pandemic is thought to result from a high viral transmission rate. Here, we consider mechanisms that influence host cell-virus binding between the SARS-CoV-2 spike glycoprotein (SPG) and the human angiotensin-converting enzyme 2 (ACE2) with a series of peptides designed to mimic key ACE2 hot spots through adopting a helical conformation analogous to the N-terminal α1 helix of ACE2, the region experimentally shown to bind to the SARS-CoV-2 receptor-binding domain (RBD). The approach examines putative structure/function relations by assessing SPG binding affinity with surface plasmon resonance (SPR). A cyclic peptide (c[KFNHEAEDLFEKLM]) was characterized in an α-helical conformation with micromolar affinity (KD = 500 µM) to the SPG. Thus, stabilizing the helical structure of the 14-mer through cyclization improves binding to SPG by an order of magnitude. In addition, end-group peptide analog modifications and residue substitutions mediate SPG binding, with net charge playing an apparent role. Therefore, we surveyed reported viral variants, and a correlation of increased positive charge with increased virulence lends support to the hypothesis that charge is relevant to enhanced viral fusion. Overall, the structure/function relationship informs the importance of conformation and charge for virus-binding analog design.

Accelerating the discovery of the beyond rule of five compounds that have high affinities toward SARS-CoV-2 spike RBD

The battle against SARS-CoV-2 coronavirus is the focal point for the global pandemic that has affected millions of lives worldwide. The need for effective and selective therapeutics for the treatment of the disease caused by SARS-CoV-2 is critical. Herein, we performed a hierarchical computational approach incorporating molecular docking studies, molecular dynamics simulations, absolute binding energy calculations, and steered molecular dynamics simulations for the discovery of potential compounds with high affinity towards SARS-CoV-2 spike RBD. By leveraging ZINC15 database, a total of 1282 in-clinical and FDA approved drugs were filtered out from nearly 0.5 million protomers of relatively large compounds (MW > 500, and LogP ≤ 5). Our results depict plausible mechanistic aspects related to the blockage of SARS-CoV-2 spike RBD by the top hits discovered. We found that the most promising candidates, namely, ZINC95628821, ZINC95617623, ZINC3979524, and ZINC261494658, strongly bind to the spike RBD and interfere with the human ACE2 receptor. These findings accelerate the rational design of selective inhibitors targeting the spike RBD protein of SARS-CoV-2.Communicated by Ramaswamy H. Sarma.

Systematic virtual screening in search of SARS CoV-2 inhibitors against spike glycoprotein: pharmacophore screening, molecular docking, ADMET analysis and MD simulations

In the absence of efficient anti-viral medications, the coronavirus disease 2019 (COVID-19), stemming from severe acute respiratory syndrome coronavirus-2 (SARS CoV-2), has spawned a worldwide catastrophe and global emergency. Amidst several anti-viral targets of COVID-19, spike glycoprotein has been recognized as an essential target for the viral entry into the host cell. In the search of effective SARS CoV-2 inhibitors acting against spike glycoprotein, the virtual screening of 175,851 ligands from the 2020.1 Asinex BioDesign library has been performed using in silico tools like SiteMap analysis, pharmacophore-based screening, molecular docking using different levels of precision, such as high throughput virtual screening, standard precision and extra precision, followed by absorption, distribution, metabolism, excretion and toxicity analysis, and molecular dynamics (MD) simulation. Following a molecular docking study, seventeen molecules (with a docking score of less than - 6.0) were identified having the substantial interactions with the catalytic amino acid and nucleic acid residues of spike glycoprotein at the binding site. In investigations using MD simulations for 10 ns, the hit molecules (1 and 2) showed adequate compactness and uniqueness, as well as satisfactory stability. These computational research findings have offered a key starting point in the field of design and development of novel SARS CoV-2 entry inhibitors with appropriate drug likeliness.

Assessment of activity of chalcone compounds as inhibitors of 3-chymotrypsin like protease (3CLPro) of SARS-CoV-2: in silico study

The COVID-19 is still pandemic due to emerging of various variant of concern of SARS-CoV2. Hence, it is devastating the world, causing significant economic as well as social chaos. This needs great effort to search and develop effective alternatives along with vaccination. Therefore, to continue drug discovery endeavors, we used chalcone derivatives to find an effective drug candidate against SARS-CoV2. Chalcone is a common simple scaffold that exists in many diets as well as in traditional medicine. Natural as well as synthetic chalcones have shown numerous interesting biological activities and are also effective in fighting various diseases. Hence, various computational methods were applied to find out potential inhibitors of 3CL^Pro using a library of 3000 compounds of chalcones. Firstly, the screening by structure-based pharmacophore model yielded 84 hits that were subjected to molecular docking. The top 10 docked compounds were characterized for stability by using 100 ns molecular dynamic (MD) simulation approach. Further, the binding free energy calculation by MMPBSA showed that four compounds bind to 3CL^Pro enzyme with high affinity, i.e., − 87.962 (kJ/mol), − 66.125 (kJ/mol), − 59.589 (kJ/mol), and − 66.728 (kJ/mol), respectively. Since chalcone is a common simple scaffold that is present in many diets as well as in traditional medicine, we suggest that screened compounds may emerge as promising drug candidates for SARS-CoV-2. These compounds may be investigated in vitro to evaluate the efficacy against SARS-CoV-2.

ACE2 Peptide Fragment Interaction with Different S1 Protein Sites

We study the effect of the peptide QAKTFLDKFNHEAEDLFYQ on the kinetics of the SARS-CoV-2 spike protein S1 binding to angiotensin-converting enzyme 2 (ACE2), with the aim to characterize the interaction mechanism of the SARS-CoV2 virus with its host cell. This peptide corresponds to the sequence 24–42 of the ACE2 α1 domain, which marks the binding site for the S1 protein. The kinetics of S1-ACE2 complex formation was measured in the presence of various concentrations of the peptide using bio-layer interferometry. Formation of the S1-ACE2 complex was inhibited by the peptide in cases where it was preincubated with S1 protein before the binding experiment. The kinetic analysis of S1-ACE2 complex dissociation revealed that preincubation stabilized this complex, and this effect was dependent on the peptide concentration as well as the preincubation time. The results point to the formation of the ternary complex of S1 with ACE2 and the peptide. This is possible in the presence of another binding site for the S1 protein beside the receptor-binding domain for ACE2, which binds the peptide QAKTFLDKFNHEAEDLFYQ. Therefore, we conducted computational mapping of the S1 protein surface, revealing two additional binding sites located at some distance from the main receptor-binding domain on S1. We suggest the possibility to predict and test the short protein derived peptides for development of novel strategies in inhibiting virus infections.

Ninhydrin as a novel DNA hybridization indicator applied to a highly reusable electrochemical genosensor for Candida auris

This work reports a simple strategy for Candida auris genomic DNA (gDNA) detection, a multi-resistant fungus associated with nosocomial outbreaks in healthcare settings, presenting high mortality and morbidity rates. The platform was developed using gold electrode sensitized with specific DNA capture probe and ninhydrin as a novel DNA hybridization indicator. The genosensor was able to detect C. auris in urine sample by differential pulse voltammetry and electrochemical impedance spectroscopy. The biosensor's analytical performance was evaluated by differential pulse voltammetry, detecting up to 4.5 pg μL^-1 of C. auris gDNA in urine (1:10, V/V). Moreover, the genosensor was reused eight times with no loss in the current signal response. The genosensor showed selectivity and stability, maintaining 100% of its response up to 80 days of storage. In order to analyze interactions of single and double-stranded DNA with ninhydrin, SEM, AFM and molecular dynamics studies followed by docking simulations were performed. Theoretical calculations showed ninhydrin interactions more favorably with dsDNA in an A-T rich binding pocket rather than with the ssDNA. Therefore, the proposed system is a promising electrochemical detection device towards a more accurate detection of C. auris gDNA in biological samples.

Small Neutral Crowding Solute Effects on Protein Folding Thermodynamic Stability and Kinetics

Molecular crowding is a ubiquitous phenomenon in biological systems, with significant consequences on protein folding and stability. Small compounds, such as the osmolyte trimethylamine N-oxide (TMAO), can also present similar effects. To analyze the effects arising from these solute-like molecules, we performed a series of crowded coarse-grained folding simulations. Two well-known proteins were chosen, CI2 and SH3, modeled by the alpha-carbon-structure-based model. In the simulations, the crowding molecules were represented by low-sized neutral atom beads in different concentrations. The results show that a low level of the volume fraction occupied by neutral agents can change protein stability and folding kinetics for the two systems. However, the kinetics were shown to be unaffected in their respective folding temperatures. The faster kinetics correlates with changes in the folding route for systems at the same temperature, where non-native contacts were shown to be relevant. The transition states of the two systems with and without crowders are similar. In the case of SH3, there are differences in the structuring of two strands, which may be associated with the increase in its folding rate, in addition to the destabilization of the denatured ensemble. The present study also detected a crossover in the thermodynamic stability behavior, previously observed experimentally and theoretically. As the temperature increases, crowders change from destabilizing to stabilizing agents.

Activation Pathways and Free Energy Landscapes of the SARS-CoV-2 Spike Protein

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) uses a spike protein (S-protein) to recognize the receptor protein ACE2 of human cells and initiate infection, during which the conformational transition of the S-protein from inactive (down) state to active (up) state is one of the key molecular events determining the infectivity but the underlying mechanism remains poorly understood. In this work, we investigated the activation pathways and free energy landscape of the S-protein of SARS-CoV-2 and compared with those of the closely related counterpart SARS-CoV using molecular dynamics simulations. Our results revealed a large difference between the activation pathways of the two S-proteins. The transition from inactive to an active state for the S-protein of SARS-CoV-2 is more cooperative, involving simultaneous disruptions of several key interfacial hydrogen bonds, and the transition encounters a much higher free energy barrier. In addition, the conformational equilibrium of the SARS-CoV-2 S-protein is more biased to the inactive state compared to that of the SARS-CoV S-protein, suggesting the transient feature of the active state before binding to the receptor protein of the host cell. The key interactions contributing to the difference of the activation pathways and free energy landscapes were discussed. The results provide insights into the molecular mechanism involved in the initial stage of the SARS-CoV-2 infection.

Effective ACE2 peptide-nanoparticle conjugation and its binding with the SARS-Cov-2 RBD quantified by dynamic light scattering

The infection of coronavirus initiates with the binding between its spike protein receptor binding domain (RBD) and a human cellular receptor called angiotensin-converting enzyme 2 (ACE2). Here, we construct truncated ACE2 peptide-conjugated gold nanoparticles as antiviral scaffolds and study their binding with the SARS-CoV-2 RBD using dynamic light scattering (DLS). Systematic DLS analysis identifies the effective peptide-nanoparticle conjugation and its efficient, specific, and long-lasting multivalent binding towards the RBD with a binding affinity of 41 nM, indicating the potential of this antiviral platform to compete with natural ACE2-RBD interactions for viral blocking and showcasing an accessible approach to measure the binding constants and kinetics.

Engineering mesenchymal stromal cells with neutralizing and anti-inflammatory capability against SARS-CoV-2 infection

The emergence of the novel human severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has led to the pandemic of coronavirus disease 2019 (COVID-19), which has markedly affected global health and the economy. Both uncontrolled viral replication and a proinflammatory cytokine storm can cause severe tissue damage in patients with COVID-19. SARS-CoV-2 utilizes angiotensin-converting enzyme 2 (ACE2) as its entry receptor. In this study, we generated ACE2 extracellular domain-Fc and single-chain variable fragment-interleukin 6 (IL-6) single-chain variable fragment against IL-6 receptor (scFv-IL6R)-Fc fusion proteins to differentially neutralize viruses and ameliorate the cytokine storm. The human ACE2 (hACE2)_1-740-Fc fusion protein showed a potent inhibitory effect on pseudo-typed SARS-CoV-2 entry and a good safety profile in mice. In addition, scFv-IL6R-Fc strongly blocked IL-6 signal activation. We also established a mesenchymal stromal cell (MSC)-based hACE2_1-740-Fc and scFv-IL6R-Fc delivery system, which could serve as a potential therapy strategy for urgent clinical needs of patients with COVID-19.

Computational optimization of angiotensin-converting enzyme 2 for SARS-CoV-2 Spike molecular recognition

Since the beginning of the Covid19 pandemic, many efforts have been devoted to identifying approaches to neutralize SARS-CoV-2 replication within the host cell. A promising strategy to block the infection consists of using a mutant of the human receptor angiotensin-converting enzyme 2 (ACE2) as a decoy to compete with endogenous ACE2 for the binding to the SARS-CoV-2 Spike protein, which decreases the ability of the virus to enter the host cell. Here, using a computational framework based on the 2D Zernike formalism we investigate details of the molecular binding and evaluate the changes in ACE2-Spike binding compatibility upon mutations occurring in the ACE2 side of the molecular interface. We demonstrate the efficacy of our method by comparing our results with experimental binding affinities changes upon ACE2 mutations, separating ones that increase or decrease binding affinity with an Area Under the ROC curve ranging from 0.66 to 0.93, depending on the magnitude of the effects analyzed. Importantly, the iteration of our approach leads to the identification of a set of ACE2 mutants characterized by an increased shape complementarity with Spike. We investigated the physico-chemical properties of these ACE2 mutants and propose them as bona fide candidates for Spike recognition.