Structural behavior of monomer of SARS-CoV-2 spike protein during initial stage of adsorption on graphene

Spike glycoprotein of the SARS-CoV-2 virus and its structure play a crucial role in the infections of cells containing angiotensin-converting enzyme 2 (ACE2) as well as in the interactions of this virus with surfaces. Protection against viruses and often even their deactivation is one of the great varieties of graphene applications. The structural changes of the non-glycosylated monomer of the spike glycoprotein trimer (denoted as S-protein in this work) triggered by its adsorption onto graphene at the initial stage are investigated by means of atomistic molecular dynamics simulations. The adsorption of the S-protein happens readily during the first 10 ns. The shape of the S-protein becomes more prolate during the adsorption, but this trend, albeit less pronounced, is observed also for the freely relaxing S-protein in water. The receptor-binding domain (RBD) of the free and adsorbed S-protein manifests itself as the most rigid fragment of the whole S-protein. The adsorption even enhances the rigidity of the whole S-protein as well as its subunits. Only one residue of the RBD involved in the specific interactions with ACE2 during the cell infection is involved in the direct contact of the adsorbed S-protein with the graphene. The new intramolecular hydrogen bonds formed during the S-protein adsorption replace the S-protein-water hydrogen bonds; this trend, although less apparent, is observed also during the relaxation of the free S-protein in water. In the initial phase, the secondary structure of the RBD fragment specifically interacting with ACE2 receptor is not affected during the S-protein adsorption onto the graphene.     

Computational Screening of Phenylamino-Phenoxy-Quinoline Derivatives against the Main Protease of SARS-CoV-2 Using Molecular Docking and the ONIOM Method

In the search for new anti-HIV-1 agents, two forms of phenylamino-phenoxy-quinoline derivatives have been synthesized, namely, 2-phenylamino-4-phenoxy-quinoline and 6-phenylamino-4-phenoxy-quinoline. In this study, the binding interactions of phenylamino-phenoxy-quinoline derivatives and six commercially available drugs (hydroxychloroquine, ritonavir, remdesivir, S-217622, N3, and PF-07321332) with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main protease (M^pro) were investigated using molecular docking and the ONIOM method. The molecular docking showed the hydrogen bonding and hydrophobic interactions of all the compounds in the pocket of SARS-CoV-2 main protease (M^pro), which plays an important role for the division and proliferation of the virus into the cell. The binding free energy values between the ligands and M^pro ranged from −7.06 to −10.61 kcal/mol. The molecular docking and ONIOM results suggested that 4-(2′,6′-dimethyl-4′-cyanophenoxy)-2-(4″-cyanophenyl)-aminoquinoline and 4-(4′-cyanophenoxy)-2-(4″-cyanophenyl)-aminoquinoline have low binding energy values and appropriate molecular properties; moreover, both compounds could bind to M^pro via hydrogen bonding and Pi-Pi stacking interactions with amino acid residues, namely, HIS41, GLU166, and GLN192. These amino acids are related to the proteolytic cleavage process of the catalytic triad mechanisms. Therefore, this study provides important information for further studies on synthetic quinoline derivatives as antiviral candidates in the treatment of SARS-CoV-2.     

Modeling the interactions of a peptide-major histocompatibility class I ligand with its receptors. II. Cross-reaction between a monoclonal antibody and two αβ T cell receptors

The recombinant antibody, pSAN13.4.1, has a unique T cell like specificity; it binds an Influenza Hemagglutinin octapeptide (Ha_255–262) in an MHC (H-2K^k)-restricted manner, and a detailed comparison of the fine specificity of pSAN13.4.1 with the fine specificity of two Ha_255–262-specific, H-2K^k-restricted T cell hybridomas has supported this contention. A three-dimensional model of pSAN13.4.1 has been derived by homology modeling techniques. Subsequently, the structure of the pSAN13.4.1 antibody in complex with the antigenic Ha-K^k ligand was derived after a flexible and automated docking of the MHC-peptide pair into the Fab combining site. Interestingly, the most energetically favored binding mode shows numerous analogies to the recently determined recognition of class I MHC-peptide complexes by T cell receptors (TCRs). The pSAN13.4.1 also binds diagonally across the MHC binding groove but is more deeply anchored to the peptide-MHC (pep/MHC) ligand than TCRs, notably through numerous interactions of its heavy chain. The present model accounts well for the experimentally determined binding affinity of a set of 144 single amino acid substituted Ha analogues and the observed shared specificity between the pSAN antibody and two different T cell receptors for the Ha-K^k antigenic ligand. Analogies and differences between Fab and TCR recognition are explained by dissecting the binding role of each chain of the immune receptors as well as the contribution of all peptide amino acids.     

Structural Insight into the Binding of Cyanovirin-N with the Spike Glycoprotein,Mpro and PLpro of SARS-CoV-2: Protein–Protein Interactions,Dynamics Simulations and Free Energy Calculations

The emergence of COVID-19 continues to pose severe threats to global public health. The pandemic has infected over 171 million people and claimed more than 3.5 million lives to date. We investigated the binding potential of antiviral cyanobacterial proteins including cyanovirin-N, scytovirin and phycocyanin with fundamental proteins involved in attachment and replication of SARS-CoV-2. Cyanovirin-N displayed the highest binding energy scores (−16.8 ± 0.02 kcal/mol, −12.3 ± 0.03 kcal/mol and −13.4 ± 0.02 kcal/mol, respectively) with the spike protein, the main protease (M^pro) and the papainlike protease (PL^pro) of SARS-CoV-2. Cyanovirin-N was observed to interact with the crucial residues involved in the attachment of the human ACE2 receptor. Analysis of the binding affinities calculated employing the molecular mechanics-Poisson–Boltzmann surface area (MM-PBSA) approach revealed that all forms of energy, except the polar solvation energy, favourably contributed to the interactions of cyanovirin-N with the viral proteins. With particular emphasis on cyanovirin-N, the current work presents evidence for the potential inhibition of SARS-CoV-2 by cyanobacterial proteins, and offers the opportunity for in vitro and in vivo experiments to deploy the cyanobacterial proteins as valuable therapeutics against COVID-19.     

Virtual Screening of Natural Chemical Databases to Search for Potential ACE2 Inhibitors

The angiotensin-converting enzyme II (ACE2) is a multifunctional protein in both health and disease conditions, which serves as a counterregulatory component of RAS function in a cardioprotective role. ACE2 modulation may also have relevance to ovarian cancer, diabetes, acute lung injury, fibrotic diseases, etc. Furthermore, since the outbreak of the coronavirus disease in 2019 (COVID-19), ACE2 has been recognized as the host receptor of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The receptor binding domain of the SARS-CoV-2 S-protein has a strong interaction with ACE2, so ACE2 may be a potent drug target to prevent the virus from invading host cells for anti-COVID-19 drug discovery. In this study, structure- and property-based virtual screening methods were combined to filter natural product databases from ChemDiv, TargetMol, and InterBioScreen to find potential ACE2 inhibitors. The binding affinity between protein and ligands was predicted using both Glide SP and XP scoring functions and the MM-GBSA method. ADME properties were also calculated to evaluate chemical drug-likeness. Then, molecular dynamics (MD) simulations were performed to further explore the binding modes between the highest-potential compounds and ACE2. Results showed that the compounds 154-23-4 and STOCK1N-07141 possess potential ACE2 inhibition activities and deserve further study.     

Aromatic Rings as Molecular Determinants for the Molecular Recognition of Protein Kinase Inhibitors

Protein kinases are key enzymes in many signal transduction pathways, and play a crucial role in cellular proliferation, differentiation, and various cell regulatory processes. However, aberrant function of kinases has been associated with cancers and many other diseases. Consequently, competitive inhibition of the ATP binding site of protein kinases has emerged as an effective means of curing these diseases. Over the past three decades, thousands of protein kinase inhibitors (PKIs) with varying molecular frames have been developed. Large-scale data mining of the Protein Data Bank resulted in a database of 2139 non-redundant high-resolution X-ray crystal structures of PKIs bound to protein kinases. This provided us with a unique opportunity to study molecular determinants for the molecular recognition of PKIs. A chemoinformatic analysis of 2139 PKIs resulted in findings that PKIs are “flat” molecules with high aromatic ring counts and low fractions of sp³ carbon. All but one PKI possessed one or more aromatic rings. More importantly, it was found that the average weighted hydrogen bond count is inversely proportional to the number of aromatic rings. Based on this linear relationship, we put forward the exchange rule of hydrogen bonding interactions and non-bonded π-interactions. Specifically, a loss of binding affinity caused by a decrease in hydrogen bonding interactions is compensated by a gain in binding affinity acquired by an increase in aromatic ring-originated non-bonded interactions (i.e., π–π stacking interactions, CH–π interactions, cation–π interactions, etc.), and vice versa. The very existence of this inverse relationship strongly suggests that both hydrogen bonding and aromatic ring-originated non-bonded interactions are responsible for the molecular recognition of PKIs. As an illustration, two representative PKI–kinase complexes were employed to examine the relative importance of different modes of non-bonded interactions for the molecular recognition of PKIs. For this purpose, two FDA-approved PKI drugs, ibrutinib and lenvatinib, were chosen. The binding pockets of both PKIs were thoroughly examined to identify all non-bonded intermolecular interactions. Subsequently, the strengths of interaction energies between ibrutinib and its interacting residues in tyrosine kinase BTK were quantified by means of the double hybrid DFT method B2PLYP. The resulting energetics for the binding of ibrutinib in tyrosine kinase BTK showed that CH–π interactions and π–π stacking interactions between aromatic rings of the drug and hydrophobic residues in its binding pocket dominate the binding interactions. Thus, this work establishes that, in addition to hydrogen bonding, aromatic rings function as important molecular determinants for the molecular recognition of PKIs. In conclusion, our findings support the following pharmacophore model for ATP-competitive kinase inhibitors: a small molecule features a scaffold of one or more aromatic rings which is linked with one or more hydrophilic functional groups. The former has the structural role of acting as a scaffold and the functional role of participating in aromatic ring-originated non-bonded interactions with multiple hydrophobic regions in the ATP binding pocket of kinases. The latter ensure water solubility and form hydrogen bonds with the hinge region and other hydrophilic residues of the ATP binding pocket.     

Binding of human ACE2 and RBD of Omicron enhanced by unique interaction patterns among SARS-CoV-2 variants of concern

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus causing COVID-19, has continued to mutate and spread worldwide despite global vaccination efforts. In particular, the Omicron variant, first identified in South Africa in late November 2021, has become the dominant strain worldwide. Compared to the original strain identified in Wuhan, Omicron features 50 genetic mutations, with 15 mutations in the receptor-binding domain (RBD) of the spike protein, which binds to the human angiotensin-converting enzyme 2 (ACE2) receptor for viral entry. However, it is not completely understood how these mutations alter the interaction and binding strength between the Omicron RBD and ACE2. In this study, we used a combined steered molecular dynamics (SMD) simulation and experimental microscale thermophoresis (MST) approach to quantify the interaction between Omicron RBD and ACE2. We report that the Omicron brings an enhanced RBD-ACE2 interface through N501Y, Q498R, and T478K mutations; the changes further lead to unique interaction patterns, reminiscing the features of previously dominated variants, Alpha (N501Y) and Delta (L452R and T478K). Among the Q493K and Q493R, we report that Q493R shows stronger binding to ACE2 than Q493K due to increased interactions. Our MST data confirmed that the Omicron mutations in RBD are associated with a five-fold higher binding affinity to ACE2 compared to the RBD of the original strain. In conclusion, our results could help explain the Omicron variant's prevalence in human populations, as higher interaction forces or affinity for ACE2 likely promote greater viral binding and internalization, leading to increased infectivity.     

Structural evaluation of an alternative Protein A biomimetic ligand for antibody purification

Affinity chromatography is one of the most common techniques employed at the industrial-scale for antibody purification. In particular, the purification of human immunoglobulin G (hIgG) has gained relevance with the immobilization of its natural binding counterpart—Staphylococcus aureus Protein A (SpA) or with the recent development of biomimetic affinity ligands, namely triazine-based ligands. These ligands have been developed in order to overcome economic and leaching issues associated to SpA. The most recent triazine-based ligand—TPN-BM, came up as an analogue of 2-(3-amino-phenol)-6-(4-amino-1-naphthol)-4-chloro-sym-triazine ligand also known as ligand 22/8 with improved physico-chemical properties and a greener synthetic route. This work intends to evaluate the potential of TPN-BM as an alternative affinity ligand towards antibody recognition and binding, namely IgG, at an atomic level, since it has already been tested, after immobilization onto chitosan-based monoliths and demonstrated interesting affinity behaviour for this purpose. Herein, combining automated molecular docking and molecular dynamics simulations it was predicted that TPN-BM has high propensity to bind IgG through the same binding site found in the crystallographic structure of SpA_IgG complex, as well as theoretically predicted for ligand 22/8_IgG complex. Furthermore, it was found that TPN-BM established preferential interactions with aromatic residues at the Fab domain (Trp 50, Tyr 53, Tyr 98 and Trp 100), while in the Fc domain the main interactions are based on hydrogen bonds with pH sensitive residues at operational regime for binding and elution like histidines (His 460, His 464, His 466). Moreover, the pH dependence of TPN-BM_IgG complex formation was more evident for the Fc domain, where at pH 3 the protonation state and consequently the charge alteration of histidine residues located at the IgG binding site induced ligand detachment which explains the optimal elution condition at this pH observed experimentally.     

Discovery of SARS-CoV-2 Mpro peptide inhibitors from modelling substrate and ligand binding

The main protease (M^pro) of SARS-CoV-2 is central to viral maturation and is a promising drug target, but little is known about structural aspects of how it binds to its 11 natural cleavage sites. We used biophysical and crystallographic data and an array of biomolecular simulation techniques, including automated docking, molecular dynamics (MD) and interactive MD in virtual reality, QM/MM, and linear-scaling DFT, to investigate the molecular features underlying recognition of the natural M^pro substrates. We extensively analysed the subsite interactions of modelled 11-residue cleavage site peptides, crystallographic ligands, and docked COVID Moonshot-designed covalent inhibitors. Our modelling studies reveal remarkable consistency in the hydrogen bonding patterns of the natural M^pro substrates, particularly on the N-terminal side of the scissile bond. They highlight the critical role of interactions beyond the immediate active site in recognition and catalysis, in particular plasticity at the S2 site. Building on our initial M^pro-substrate models, we used predictive saturation variation scanning (PreSaVS) to design peptides with improved affinity. Non-denaturing mass spectrometry and other biophysical analyses confirm these new and effective ‘peptibitors’ inhibit M^pro competitively. Our combined results provide new insights and highlight opportunities for the development of M^pro inhibitors as anti-COVID-19 drugs.

The main protease (M^pro) of SARS-CoV-2 is central to viral maturation and is a promising drug target. In silico methods reveal structural aspects of how it binds to its 11 natural cleavage sites, the design of novel peptide inhibitors, and insights into drug design.     

Probing conformational hotspots for the recognition and intervention of protein complexes by lysine reactivity profiling

Probing the conformational and functional hotspot sites within aqueous native protein complexes is still a challenging task. Herein, a mass spectrometry (MS)-based two-step isotope labeling-lysine reactivity profiling (TILLRP) strategy is developed to quantify the reactivities of lysine residues and probe the molecular details of protein–protein interactions as well as evaluate the conformational interventions by small-molecule active compounds. The hotspot lysine sites that are crucial to the SARS-CoV-2 S1–ACE2 combination could be successfully probed, such as S1 Lys⁴¹⁷ and Lys⁴⁴⁴. Significant alteration of the reactivities of lysine residues at the interaction interface of S1-RBD Lys³⁸⁶–Lys⁴⁶² was observed during the formation of complexes, which might be utilized as indicators for investigating the S1-ACE2 dynamic recognition and intervention at the molecular level in high throughput.

A mass spectrometry-based two-step isotope labeling-lysine reactivity profiling strategy is developed to probe the molecular details of protein–protein interactions and evaluate the conformational interventions by small-molecule active compounds.     

Inhibitor binding influences the protonation states of histidines in SARS-CoV-2 main protease

The main protease (M^pro) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an attractive target for antiviral therapeutics. Recently, many high-resolution apo and inhibitor-bound structures of M^pro, a cysteine protease, have been determined, facilitating structure-based drug design. M^pro plays a central role in the viral life cycle by catalyzing the cleavage of SARS-CoV-2 polyproteins. In addition to the catalytic dyad His41–Cys145, M^pro contains multiple histidines including His163, His164, and His172. The protonation states of these histidines and the catalytic nucleophile Cys145 have been debated in previous studies of SARS-CoV M^pro, but have yet to be investigated for SARS-CoV-2. In this work we have used molecular dynamics simulations to determine the structural stability of SARS-CoV-2 M^pro as a function of the protonation assignments for these residues. We simulated both the apo and inhibitor-bound enzyme and found that the conformational stability of the binding site, bound inhibitors, and the hydrogen bond networks of M^pro are highly sensitive to these assignments. Additionally, the two inhibitors studied, the peptidomimetic N3 and an α-ketoamide, display distinct His41/His164 protonation-state-dependent stabilities. While the apo and the N3-bound systems favored N_δ (HD) and N_ϵ (HE) protonation of His41 and His164, respectively, the α-ketoamide was not stably bound in this state. Our results illustrate the importance of using appropriate histidine protonation states to accurately model the structure and dynamics of SARS-CoV-2 M^pro in both the apo and inhibitor-bound states, a necessary prerequisite for drug-design efforts.

The main protease (M^pro) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an attractive target for antiviral therapeutics.     

Understanding the recognition mechanism of protein-RNA complexes using energy based approach

Protein-RNA interactions perform diverse functions within the cell. Understanding the recognition mechanism of protein-RNA complexes is a challenging task in molecular and computational biology. In this work, we have developed an energy based approach for identifying the binding sites and important residues for binding in protein-RNA complexes. The new approach considers the repulsive interactions as well as the effect of distance between the atoms in protein and RNA in terms of interaction energy, which are not considered in traditional distance based methods to identify the binding sites. We found that the positively charged, polar and aromatic residues are important for binding. These residues influence to form electrostatic, hydrogen bonding and stacking interactions. Our observation has been verified with the experimental binding specificity of protein-RNA complexes and found good agreement with experiments. Further, the propensities of residues/nucleotides in the binding sites of proteins/RNA and their atomic contributions have been derived. Based on these results we have proposed a novel mechanism for the recognition of protein-RNA complexes: the charged and polar residues in proteins initiate recognition with RNA by making electrostatic and hydrogen bonding interactions between them; the aromatic side chains tend to form aromatic-aromatic interactions and the hydrophobic residues aid to stabilize the complex.     

Plant-Derived Antiviral Compounds as Potential Entry Inhibitors against Spike Protein of SARS-CoV-2 Wild-Type and Delta Variant: An Integrative in SilicoApproach

The wild-type SARS-CoV-2 has continuously evolved into several variants with increased transmissibility and virulence. The Delta variant which was initially identified in India created a devastating impact throughout the country during the second wave. While the efficacy of the existing vaccines against the latest SARS-CoV-2 variants remains unclear, extensive research is being carried out to develop potential antiviral drugs through approaches like in silico screening and drug-repurposing. This study aimed to conduct the docking-based virtual screening of 50 potential phytochemical compounds against a Spike glycoprotein of the wild-type and the Delta SARS-CoV-2 variant. Subsequently, molecular docking was performed for the five best compounds, such as Lupeol, Betulin, Hypericin, Corilagin, and Geraniin, along with synthetic controls. From the results obtained, it was evident that Lupeol exhibited a remarkable binding affinity towards the wild-type Spike protein (−8.54 kcal/mol), while Betulin showed significant binding interactions with the mutated Spike protein (−8.83 kcal/mol), respectively. The binding energy values of the selected plant compounds were slightly higher than that of the controls. Key hydrogen bonding and hydrophobic interactions of the resulting complexes were visualized, which explained their greater binding affinity against the target proteins—the Delta S protein of SARS-CoV-2, in particular. The lower RMSD, the RMSF values of the complexes and the ligands, Rg, H-bonds, and the binding free energies of the complexes together revealed the stability of the complexes and significant binding affinities of the ligands towards the target proteins. Our study suggests that Lupeol and Betulin could be considered as potential ligands for SARS-CoV-2 spike antagonists. Further experimental validations might provide new insights for the possible antiviral therapeutic interventions of the identified lead compounds and their analogs against COVID-19 infection.     

Does SARS-CoV-2 affect neurodegenerative disorders? TLR2, a potential receptor for SARS-CoV-2 in the CNS

The coronavirus (COVID-19) pandemic, caused by severe acute respiratory system coronavirus 2 (SARS-CoV-2), has created significant challenges for scientists seeking to understand the pathogenic mechanisms of SARS-CoV-2 infection and to identify the best therapies for infected patients. Although ACE2 is a known receptor for the virus and has been shown to mediate viral entry into the lungs, accumulating reports highlight the presence of neurological symptoms resulting from infection. As ACE2 expression is low in the central nervous system (CNS), these neurological symptoms are unlikely to be caused by ACE2-virus binding. In this review, we will discuss a proposed interaction between SARS-CoV-2 and Toll-like receptor 2 (TLR2) in the CNS. TLR2 is an innate immune receptor that recognizes exogenous microbial components but has also been shown to interact with multiple viral components, including the envelope (E) protein of SARS-CoV-2. In addition, TLR2 plays an important role in the pathogenesis of neurodegenerative diseases such as Alzheimer’s disease (AD) and Parkinson’s disease (PD). Based on these observations, we hypothesize that TLR2 may play a critical role in the response to SARS-CoV-2 infiltration in the CNS, thereby resulting in the induction or acceleration of AD and PD pathologies in patients.Subject terms: Neurodegeneration, Neurodegenerative diseases     

In silico study of novel niclosamide derivatives,SARS-CoV-2 nonstructural proteins catalytic residue-targeting small molecules drug candidates

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2)-mediated coronavirus disease 2019 (COVID-19) infection remains a global pandemic and health emergency with overwhelming social and economic impacts throughout the world. Therapeutics for COVID-19 are limited to only remdesivir; therefore, there is a need for combined, multidisciplinary efforts to develop new therapeutic molecules and explore the effectiveness of existing drugs against SARS-CoV-2. In the present study, we reported eight (SCOV-L-02, SCOV-L-09, SCOV-L-10, SCOV-L-11, SCOV-L-15, SCOV-L-18, SCOV-L-22, and SCOV-L-23) novel structurally related small-molecule derivatives of niclosamide (SCOV-L series) for their targeting potential against angiotensin-converting enzyme-2 (ACE2), type II transmembrane serine protease (TMPRSS2), and SARS-COV-2 nonstructural proteins (NSPs) including NSP5 (3CLpro), NSP3 (PLpro), and RdRp. Our correlation analysis suggested that ACE2 and TMPRSS2 modulate host immune response via regulation of immune-infiltrating cells at the site of tissue/organs entries. In addition, we identified some TMPRSS2 and ACE2 microRNAs target regulatory networks in SARS-CoV-2 infection and thus open up a new window for microRNAs-based therapy for the treatment of SARS-CoV-2 infection. Our in vitro study revealed that with the exception of SCOV-L-11 and SCOV-L-23 which were non-active, the SCOV-L series exhibited strict antiproliferative activities and non-cytotoxic effects against ACE2- and TMPRSS2-expressing cells. Our molecular docking for the analysis of receptor-ligand interactions revealed that SCOV-L series demonstrated high ligand binding efficacies (at higher levels than clinical drugs) against the ACE2, TMPRSS2, and SARS-COV-2 NSPs. SCOV-L-18, SCOV-L-15, and SCOV-L-09 were particularly found to exhibit strong binding affinities with three key SARS-CoV-2’s proteins: 3CLpro, PLpro, and RdRp. These compounds bind to the several catalytic residues of the proteins, and satisfied the criteria of drug-like candidates, having good adsorption, distribution, metabolism, excretion, and toxicity (ADMET) pharmacokinetic profile. Altogether, the present study suggests the therapeutic potential of SCOV-L series for preventing and managing SARs-COV-2 infection and are currently under detailed investigation in our lab.     

Reactions of proton transfer and exchange with the participation of OH,SH, and CH acids and amines

A comparison of proton exchange reactions between OH, SH, and CH acids and the NH groups of trialkylammonium ions showed that regardless of the nature of the acid XH, the mechanism of exchange includes transfer of a proton in the ion pair N-H⁺ ... X^– as the slow step. At the fast steps of proton exchange XH- N⁺H, i.e., molecular exchange with breaking of a hydrogen bond X-H ... N and transfer of a proton along these bonds, differences appear in the properties of XH acids. In the sequence from OH to SH and CH acids, the hydrogen bonds X-H ... N are weakened. As a result of this, in the same sequence the kinetic acidity (k₂) decreases but the rate of molecular exchange (k_H) increases. The ratio between the values of k₂ and k_H is inverted when the strong bonds O-H ... N (k₂/k_H 1) are replaced by weak bonds C-H ... N (k₂/k_H 1). It was also established that the kinetic stability of the anions increases as the oxygen atoms are replaced by sulfur in the series RCOO^– < RCOS^– < R₂PSS^– as a result of the more effective delocalization of the negative charge on the diffuse orbitals of sulfur.Translated from Teoreticheskaya i Éksperimental'naya Khimiya, Vol. 23, No. 4, pp. 471–475, July–August 1987.     

In Silico Screening of Novel TMPRSS2 Inhibitors for Treatment of COVID-19

COVID-19, a pandemic caused by the virus SARS-CoV-2, has spread globally, necessitating the search for antiviral compounds. Transmembrane protease serine 2 (TMPRSS2) is a cell surface protease that plays an essential role in SARS-CoV-2 infection. Therefore, researchers are searching for TMPRSS2 inhibitors that can be used for the treatment of COVID-19. As such, in this study, based on the crystal structure, we targeted the active site of TMPRSS2 for virtual screening of compounds in the FDA database. Then, we screened lumacaftor and ergotamine, which showed strong binding ability, using 100 ns molecular dynamics simulations to study the stability of the protein–ligand binding process, the flexibility of amino acid residues, and the formation of hydrogen bonds. Subsequently, we calculated the binding free energy of the protein–ligand complex by the MM-PBSA method. The results show that lumacaftor and ergotamine interact with residues around the TMPRSS2 active site, and reached equilibrium in the 100 ns molecular dynamics simulations. We think that lumacaftor and ergotamine, which we screened through in silico studies, can effectively inhibit the activity of TMPRSS2. Our findings provide a basis for subsequent in vitro experiments, having important implications for the development of effective anti-COVID-19 drugs.     

Efficacy of Phytochemicals Derived from Avicennia officinalis for the Management of COVID-19: A Combined In Silico and Biochemical Study

The recent coronavirus disease 2019 (COVID-19) pandemic is a global threat for healthcare management and the economic system, and effective treatments against the pathogenic severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus responsible for this disease have not yet progressed beyond the developmental phases. As drug refinement and vaccine progression require enormously broad investments of time, alternative strategies are urgently needed. In this study, we examined phytochemicals extracted from Avicennia officinalis and evaluated their potential effects against the main protease of SARS-CoV-2. The antioxidant activities of A. officinalis leaf and fruit extracts at 150 µg/mL were 95.97% and 92.48%, respectively. Furthermore, both extracts displayed low cytotoxicity levels against Artemia salina. The gas chromatography–mass spectroscopy analysis confirmed the identifies of 75 phytochemicals from both extracts, and four potent compounds, triacontane, hexacosane, methyl linoleate, and methyl palminoleate, had binding free energy values of −6.75, −6.7, −6.3, and −6.3 Kcal/mol, respectively, in complexes with the SARS-CoV-2 main protease. The active residues Cys145, Met165, Glu166, Gln189, and Arg188 in the main protease formed non-bonded interactions with the screened compounds. The root-mean-square difference (RMSD), root-mean-square fluctuations (RMSF), radius of gyration (Rg), solvent-accessible surface area (SASA), and hydrogen bond data from a molecular dynamics simulation study confirmed the docked complexes′ binding rigidity in the atomistic simulated environment. However, this study′s findings require in vitro and in vivo validation to ensure the possible inhibitory effects and pharmacological efficacy of the identified compounds.     

Synergistic inhibition of SARS-CoV-2 cell entry by otamixaban and covalent protease inhibitors: pre-clinical assessment of pharmacological and molecular properties

SARS-CoV-2, the cause of the COVID-19 pandemic, exploits host cell proteins for viral entry into human lung cells. One of them, the protease TMPRSS2, is required to activate the viral spike protein (S). Even though two inhibitors, camostat and nafamostat, are known to inhibit TMPRSS2 and block cell entry of SARS-CoV-2, finding further potent therapeutic options is still an important task. In this study, we report that a late-stage drug candidate, otamixaban, inhibits SARS-CoV-2 cell entry. We show that otamixaban suppresses TMPRSS2 activity and SARS-CoV-2 infection of a human lung cell line, although with lower potency than camostat or nafamostat. In contrast, otamixaban inhibits SARS-CoV-2 infection of precision cut lung slices with the same potency as camostat. Furthermore, we report that otamixaban''s potency can be significantly enhanced by (sub-) nanomolar nafamostat or camostat supplementation. Dominant molecular TMPRSS2-otamixaban interactions are assessed by extensive 109 μs of atomistic molecular dynamics simulations. Our findings suggest that combinations of otamixaban with supplemental camostat or nafamostat are a promising option for the treatment of COVID-19.

SARS-CoV-2, the cause of the COVID-19 pandemic, exploits host proteins for viral entry into human lung cells and is blocked by otamixaban in combination with a covalent protease inhibitor.     

Identification of SARS-CoV-2 Receptor Binding Inhibitors by In Vitro Screening of Drug Libraries

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the coronavirus disease 2019 (COVID-19) global pandemic. The first step of viral infection is cell attachment, which is mediated by the binding of the SARS-CoV-2 receptor binding domain (RBD), part of the virus spike protein, to human angiotensin-converting enzyme 2 (ACE2). Therefore, drug repurposing to discover RBD-ACE2 binding inhibitors may provide a rapid and safe approach for COVID-19 therapy. Here, we describe the development of an in vitro RBD-ACE2 binding assay and its application to identify inhibitors of the interaction of the SARS-CoV-2 RBD to ACE2 by the high-throughput screening of two compound libraries (LOPAC^®1280 and DiscoveryProbeTM). Three compounds, heparin sodium, aurintricarboxylic acid (ATA), and ellagic acid, were found to exert an effective binding inhibition, with IC50 values ranging from 0.6 to 5.5 µg/mL. A plaque reduction assay in Vero E6 cells infected with a SARS-CoV-2 surrogate virus confirmed the inhibition efficacy of heparin sodium and ATA. Molecular docking analysis located potential binding sites of these compounds in the RBD. In light of these findings, the screening system described herein can be applied to other drug libraries to discover potent SARS-CoV-2 inhibitors.