Molecular Dynamics and MM-PBSA Analysis of the SARS-CoV-2 Gamma Variant in Complex with the hACE-2 Receptor

The SARS-CoV-2 virus, since its appearance in 2019, has caused millions of cases and deaths. To date, there is no effective treatment or a vaccine that is fully protective. Despite the efforts made by governments and health institutions around the globe to control its propagation, the evolution of the virus has accelerated, diverging into hundreds of variants. However, not all of them are variants of concern (VoC’s). VoC’s have appeared in different regions and throughout the two years of the pandemic they have spread around the world. Specifically, in South America, the gamma variant (previously known as P.1) appeared in early 2021, bringing with it a second wave of infections. This variant contains the N501Y, E484K and K417T mutations in the receptor binding domain (RBD) of the spike protein. Although these mutations have been described experimentally, there is still no clarity regarding their role in the stabilization of the complex with the human angiotensin converting enzyme 2 (hACE-2) receptor. In this article we dissect the influence of mutations on the interaction with the hACE-2 receptor using molecular dynamics and estimations of binding affinity through a screened version of the molecular mechanics Poisson Boltzmann surface area (MM-PBSA) and interaction entropy. Our results indicate that mutations E484K and K417T compensate each other in terms of binding affinity, while the mutation N501Y promotes a more convoluted effect. This effect consists in the adoption of a cis configuration in the backbone of residue Y495 within the RBD, which in turn promotes polar interactions with the hACE-2 receptor. These results not only correlate with experimental observations and complement previous knowledge, but also expose new features associated with the specific contribution of concerned mutations. Additionally, we propose a recipe to assess the residue-specific contribution to the interaction entropy.     

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.     

Fruit Bromelain-Derived Peptide Potentially Restrains the Attachment of SARS-CoV-2 Variants to hACE2: A Pharmacoinformatics Approach

Before entering the cell, the SARS-CoV-2 spike glycoprotein receptor-binding domain (RBD) binds to the human angiotensin-converting enzyme 2 (hACE2) receptor. Hence, this RBD is a critical target for the development of antiviral agents. Recent studies have discovered that SARS-CoV-2 variants with mutations in the RBD have spread globally. The purpose of this in silico study was to determine the potential of a fruit bromelain-derived peptide. DYGAVNEVK. to inhibit the entry of various SARS-CoV-2 variants into human cells by targeting the hACE binding site within the RBD. Molecular docking analysis revealed that DYGAVNEVK interacts with several critical RBD binding residues responsible for the adhesion of the RBD to hACE2. Moreover, 100 ns MD simulations revealed stable interactions between DYGAVNEVK and RBD variants derived from the trajectory of root-mean-square deviation (RMSD), radius of gyration (Rg), and root-mean-square fluctuation (RMSF) analysis, as well as free binding energy calculations. Overall, our computational results indicate that DYGAVNEVK warrants further investigation as a candidate for preventing SARS-CoV-2 due to its interaction with the RBD of SARS-CoV-2 variants.     

Molecular recognition of SARS-CoV-2 spike glycoprotein: quantum chemical hot spot and epitope analyses

Due to the COVID-19 pandemic, researchers have attempted to identify complex structures of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike glycoprotein (S-protein) with angiotensin-converting enzyme 2 (ACE2) or a blocking antibody. However, the molecular recognition mechanism—critical information for drug and antibody design—has not been fully clarified at the amino acid residue level. Elucidating such a microscopic mechanism in detail requires a more accurate molecular interpretation that includes quantum mechanics to quantitatively evaluate hydrogen bonds, XH/π interactions (X = N, O, and C), and salt bridges. In this study, we applied the fragment molecular orbital (FMO) method to characterize the SARS-CoV-2 S-protein binding interactions with not only ACE2 but also the B38 Fab antibody involved in ACE2-inhibitory binding. By analyzing FMO-based interaction energies along a wide range of binding interfaces carefully, we identified amino acid residues critical for molecular recognition between S-protein and ACE2 or B38 Fab antibody. Importantly, hydrophobic residues that are involved in weak interactions such as CH–O hydrogen bond and XH/π interactions, as well as polar residues that construct conspicuous hydrogen bonds, play important roles in molecular recognition and binding ability. Moreover, through these FMO-based analyses, we also clarified novel hot spots and epitopes that had been overlooked in previous studies by structural and molecular mechanical approaches. Altogether, these hot spots/epitopes identified between S-protein and ACE2/B38 Fab antibody may provide useful information for future antibody design, evaluation of the binding property of the SARS-CoV-2 variants including its N501Y, and small or medium drug design against the SARS-CoV-2.

Quantum chemical calculations investigated molecular recognition of SARS-CoV-2 spike glycoproteins including its N501Y variant for ACE2 and antibody. Hot spot and epitope analyses revealed key residues to design drugs and antibodies against COVID-19.     

Binding of Glycans to the SARS CoV-2 Spike Protein,an Open Question: NMR Data on Binding Site Localization,Affinity, and Selectivity

We have used NMR experiments to explore the binding of selected glycans and glycomimetics to the SARS CoV-2 spike glycoprotein (S-protein) and to its receptor binding domain (RBD). STD NMR experiments confirm the binding of sialoglycans to the S-protein of the prototypic Wuhan strain virus and yield dissociation constants in the millimolar range. The absence of STD effects for sialoglycans in the presence of the Omicron/BA.1 S-protein reflects a loss of binding as a result of S-protein evolution. Likewise, no STD effects are observed for the deletion mutant Δ143–145 of the Wuhan S-protein, thus supporting localization of the binding site in the N-terminal domain (NTD). The glycomimetics Oseltamivir and Zanamivir bind weakly to the S-protein of both virus strains. Binding of blood group antigens to the Wuhan S-protein cannot be confirmed by STD NMR. Using ¹H,¹⁵N TROSY HSQC-based chemical shift perturbation (CSP) experiments, we excluded binding of any of the ligands studied to the RBD of the Wuhan S-protein. Our results put reported data on glycan binding into perspective and shed new light on the potential role of glycan-binding to the S-protein.     

Understanding of the drug resistance mechanism of hepatitis C virus NS3/4A to paritaprevir due to D168N/Y mutations: A molecular dynamics simulation perspective

Hepatitis C virus (HCV) NS3/4A protease is an attractive target for the development of antiviral therapy. However, the evolution of antiviral drug resistance is a major problem for treatment of HCV infected patients. Understanding of drug-resistance mechanisms at molecular level is therefore very important for the guidance of further design of antiviral drugs with high efficiency and specificity. Paritaprevir is a potent inhibitor against HCV NS3/4A protease genotype 1a. Unfortunately, this compound is highly susceptible to the substitution at D168 in the protease. In this work, molecular dynamics simulations of paritaprevir complexed with wild-type (WT) and two mutated strains (D168 N and D168Y) were carried out. Due to such mutations, the inhibitor-protein hydrogen bonding between them was weakened and the salt-bridge network among residues R123, R155 and D168 responsible for inhibitor binding was disrupted. Moreover, the per-residue free energy decomposition suggested that the main contributions from key residues such as Q80, V132, K136, G137 and R155 were lost in the D168 N/Y mutations. These lead to a lower binding affinity of paritaprevir for D168 N/Y variants of the HCV NS3/4A protease, consistent with the experimental data. This detailed information could be useful for further design of high potency anti-HCV NS3/4A inhibitors.     

Affinity capillary electrophoresis and density functional theory employed for the characterization of hexaarylbenzene-based receptor complexation with alkali metal ions

In this study, affinity capillary electrophoresis (ACE) and quantum mechanical density functional theory (DFT) calculations were combined to investigate non-covalent binding interactions between the hexaarylbenzene-based receptor (R) and alkali metal ions, Rb(+) and Cs(+) , in methanol. The apparent binding (stability) constants (K(b) ) of the complexes of receptor R with alkali metal ions in the methanolic medium were determined by ACE from the dependence of effective electrophoretic mobility of the receptor R on the concentration of Rb(+) and Cs(+) ions in the BGE using a non-linear regression analysis. The receptor R formed relatively strong complexes both with rubidium (log K(b) =4.04±0.21) and cesium ions (log K(b) =3.72±0.22). The structural characteristics of the above alkali metal ion complexes with the receptor R were described by ab initio density functional theory calculations. These calculations have shown that the studied cations bind to the receptor R because they synergistically interact with the polar ethereal fence and with the central benzene ring via cation-π interaction.     

Prediction and mitigation of mutation threats to COVID-19 vaccines and antibody therapies

Antibody therapeutics and vaccines are among our last resort to end the raging COVID-19 pandemic. They, however, are prone to over 5000 mutations on the spike (S) protein uncovered by a Mutation Tracker based on over 200 000 genome isolates. It is imperative to understand how mutations will impact vaccines and antibodies in development. In this work, we first study the mechanism, frequency, and ratio of mutations on the S protein which is the common target of most COVID-19 vaccines and antibody therapies. Additionally, we build a library of 56 antibody structures and analyze their 2D and 3D characteristics. Moreover, we predict the mutation-induced binding free energy (BFE) changes for the complexes of S protein and antibodies or ACE2. By integrating genetics, biophysics, deep learning, and algebraic topology, we reveal that most of the 462 mutations on the receptor-binding domain (RBD) will weaken the binding of S protein and antibodies and disrupt the efficacy and reliability of antibody therapies and vaccines. A list of 31 antibody disrupting mutants is identified, while many other disruptive mutations are detailed as well. We also unveil that about 65% of the existing RBD mutations, including those variants recently found in the United Kingdom (UK) and South Africa, will strengthen the binding between the S protein and human angiotensin-converting enzyme 2 (ACE2), resulting in more infectious COVID-19 variants. We discover the disparity between the extreme values of RBD mutation-induced BFE strengthening and weakening of the bindings with antibodies and angiotensin-converting enzyme 2 (ACE2), suggesting that SARS-CoV-2 is at an advanced stage of evolution for human infection, while the human immune system is able to produce optimized antibodies. This discovery, unfortunately, implies the vulnerability of current vaccines and antibody drugs to new mutations. Our predictions were validated by comparison with more than 1400 deep mutations on the S protein RBD. Our results show the urgent need to develop new mutation-resistant vaccines and antibodies and to prepare for seasonal vaccinations.

Antibody therapeutics and vaccines are among our last resort to end the raging COVID-19 pandemic.     

In silico analysis of nonsynonymous single nucleotide polymorphisms of the human adiponectin receptor 2 (ADIPOR2) gene

Polymorphisms of the ADIPOR2 gene are frequently linked to a higher risk of developing diseases including obesity, type 2 diabetes and cardiovascular diseases. Though mutations of the ADIPOR2 gene are detrimental, there is a lack of comprehensive in silico analyses of the functional and structural impacts at the protein level. Considering the involvement of ADIPOR2 in glucose uptake and fatty acid oxidation, an in silico functional analysis was conducted to explore the possible association between genetic mutations and phenotypic variations. A genomic analysis of 82 nonsynonymous SNPs in ADIPOR2 was initiated using SIFT followed by the SNAP2, nsSNPAnalyzer, PolyPhen-2, SNPs&GO, FATHMM and PROVEAN servers. A total of 10 mutations (R126W, L160Q, L195P, F201S, L235R, L235P, L256R, Y328H, E334K and Q349H) were predicted to have deleterious effects on the ADIPOR2 protein and were therefore selected for further analysis. Theoretical models of the variants were generated by comparative modeling via MODELLER 9.16. A protein structural analysis of these amino acid variants was performed using SNPeffect, I-Mutant, ConSurf, Swiss-PDB Viewer and NetSurfP to explore their solvent accessibility, molecular dynamics and energy minimization calculations. In addition, FTSite was used to predict the ligand binding sites, while NetGlycate, NetPhos2.0, UbPerd and SUMOplot were used to predict post-translational modification sites. All of the variants showed increased free energy, though F201S exhibited the highest energy increase. The root mean square deviation values of the modeled mutants strongly indicated likely pathogenicity. Remarkably, three binding sites were detected on ADIPOR2, and two mutations at positions 328 and 201 were found in the first and second binding pockets, respectively. Interestingly, no mutations were found at the post-translational modification sites. These genetic variants can provide a better understanding of the wide range of disease susceptibility associated with ADIPOR2 and aid the development of new molecular diagnostic markers for these diseases. The findings may also facilitate the development of novel therapeutic elements for associated diseases.     

In Silico Discovery of Antimicrobial Peptides as an Alternative to Control SARS-CoV-2

A serious pandemic has been caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The interaction between spike surface viral protein (Sgp) and the angiotensin-converting enzyme 2 (ACE2) cellular receptor is essential to understand the SARS-CoV-2 infectivity and pathogenicity. Currently, no drugs are available to treat the infection caused by this coronavirus and the use of antimicrobial peptides (AMPs) may be a promising alternative therapeutic strategy to control SARS-CoV-2. In this study, we investigated the in silico interaction of AMPs with viral structural proteins and host cell receptors. We screened the antimicrobial peptide database (APD3) and selected 15 peptides based on their physicochemical and antiviral properties. The interactions of AMPs with Sgp and ACE2 were performed by docking analysis. The results revealed that two amphibian AMPs, caerin 1.6 and caerin 1.10, had the highest affinity for Sgp proteins while interaction with the ACE2 receptor was reduced. The effective AMPs interacted particularly with Arg995 located in the S2 subunits of Sgp, which is key subunit that plays an essential role in viral fusion and entry into the host cell through ACE2. Given these computational findings, new potentially effective AMPs with antiviral properties for SARS-CoV-2 were identified, but they need experimental validation for their therapeutic effectiveness.     

Calculation of hot spots for protein–protein interaction in p53/PMI-MDM2/MDMX complexes

The recently developed MM/GBSA_IE method is applied to computing hot and warm spots in p53/PMI-MDM2/MDMX protein–protein interaction systems. Comparison of the calculated hot (>2 kcal/mol) and warm spots (>1 kcal/mol) in P53 and PMI proteins interacting with MDM2 and MDMX shows a good quantitative agreement with the available experimental data. Further, our calculation predicted hot spots in MDM2 and MDMX proteins in their interactions with P53 and PMI and they help elucidate the interaction mechanism underlying this important PPI system. In agreement with the experimental result, the present calculation shows that PMI has more hot and warm spots and binds stronger to MDM2/MDMX. The analysis of these hot and warm spots helps elucidate the fundamental difference in binding between P53 and PMI to the MDM2/MDMX systems. Specifically, for p53/PMI-MDM2 systems, p53 and PMI use essentially the same residues (L54, I61, Y67, Q72, V93, H96, and I99) of MDM2 for binding. However, PMI enhanced interactions with residues L54, Y67, and Q72 of MDM2. For the p53/PMI-MDMX system, p53 and PMI use similar residues (M53, I60, Y66, Q71, V92, and Y99) of MDMX for binding. However, PMI exploited three extra residues (M61, K93, and L98) of MDMX for enhanced binding. In addition, PMI enhanced interaction with four residues (M53, Y66, Q71, and Y99) of MDMX. These results gave quantitative explanation on why the binding affinities of PMI-MDM2/MDMX interactions are stronger than that of p53-MDM2/MDMX although their binding modes are similar. © 2018 Wiley Periodicals, Inc.     

Computational simulation of interactions between SARS coronavirus spike mutants and host species-specific receptors

As a critical adaptive mechanism, amino acid replacements on the severe acute respiratory syndrome coronavirus (SARS-CoV) spike protein could alter the receptor-binding specificity of this envelope glycoprotein and in turn lead to the emergence or reemergence of this viral zoonosis. Based on the X-ray structures of SARS-CoV spike receptor-binding domain (RBD) in complex with its functional receptor (angiotensin-converting enzyme 2, ACE2), we perform computational simulations of interactions between three representative RBD mutants and four host species-specific receptors. The comparisons between computational predictions and experimental evidences validate our structural bioinformatics approaches. And the predictions further indicate that some viral prototypes might utilize the rat ACE2 while rats might serve as a vector or reservoir of SARS-CoV.     

Site Density Functional Theory and Structural Bioinformatics Analysis of the SARS-CoV Spike Protein and hACE2 Complex

The entry of the SARS-CoV-2, a causative agent of COVID-19, into human host cells is mediated by the SARS-CoV-2 spike (S) glycoprotein, which critically depends on the formation of complexes involving the spike protein receptor-binding domain (RBD) and the human cellular membrane receptor angiotensin-converting enzyme 2 (hACE2). Using classical site density functional theory (SDFT) and structural bioinformatics methods, we investigate binding and conformational properties of these complexes and study the overlooked role of water-mediated interactions. Analysis of the three-dimensional reference interaction site model (3DRISM) of SDFT indicates that water mediated interactions in the form of additional water bridges strongly increases the binding between SARS-CoV-2 spike protein and hACE2 compared to SARS-CoV-1-hACE2 complex. By analyzing structures of SARS-CoV-2 and SARS-CoV-1, we find that the homotrimer SARS-CoV-2 S receptor-binding domain (RBD) has expanded in size, indicating large conformational change relative to SARS-CoV-1 S protein. Protomer with the up-conformational form of RBD, which binds with hACE2, exhibits stronger intermolecular interactions at the RBD-ACE2 interface, with differential distributions and the inclusion of specific H-bonds in the CoV-2 complex. Further interface analysis has shown that interfacial water promotes and stabilizes the formation of CoV-2/hACE2 complex. This interaction causes a significant structural rigidification of the spike protein, favoring proteolytic processing of the S protein for the fusion of the viral and cellular membrane. Moreover, conformational dynamics simulations of RBD motions in SARS-CoV-2 and SARS-CoV-1 point to the role in modification of the RBD dynamics and their impact on infectivity.     

Glycyrrhizic Acid Inhibits SARS-CoV-2 Infection by Blocking Spike Protein-Mediated Cell Attachment

Glycyrrhizic acid (GA), also known as glycyrrhizin, is a triterpene glycoside isolated from plants of Glycyrrhiza species (licorice). GA possesses a wide range of pharmacological and antiviral activities against enveloped viruses including severe acute respiratory syndrome (SARS) virus. Since the S protein (S) mediates SARS coronavirus 2 (SARS-CoV-2) cell attachment and cell entry, we assayed the GA effect on SARS-CoV-2 infection using an S protein-pseudotyped lentivirus (Lenti-S). GA treatment dose-dependently blocked Lenti-S infection. We showed that incubation of Lenti-S virus, but not the host cells with GA prior to the infection, reduced Lenti-S infection, indicating that GA targeted the virus for infection. Surface plasmon resonance measurement showed that GA interacted with a recombinant S protein and blocked S protein binding to host cells. Autodocking analysis revealed that the S protein has several GA-binding pockets including one at the interaction interface to the receptor angiotensin-converting enzyme 2 (ACE2) and another at the inner side of the receptor-binding domain (RBD) which might impact the close-to-open conformation change of the S protein required for ACE2 interaction. In addition to identifying GA antiviral activity against SARS-CoV-2, the study linked GA antiviral activity to its effect on virus cell binding.     

Human Y‐chromosome SNP characterization by multiplex amplified product‐length polymorphism analysis

We designed an allele‐specific amplification protocol to optimize Y‐chromosome SNP typing, which is an unavoidable step for defining the phylogenetic status of paternal lineages. It allows the simultaneous highly specific definition of up to six mutations in a single reaction by amplification fragment length polymorphism (AFLP) without the need of specialized equipment, at a considerably lower cost than that based on single‐base primer extension (SNaPshot?) technology or PCR‐RFLP systems, requiring as little as 0.5 ng DNA and compatible with the small fragments characteristic of low‐quality DNA. By designation of two primers recognizing the derived and ancestral state for each SNP, which can be differentiated by size by the addition of a noncomplementary nucleotide tail, we could define major Y clades E, F, K, R, Q, and subhaplogroups R1, R1a, R1b, R1b1b, R1b1c, J1, J2, G1, G2, I1, Q1a3, and Q1a3a1 through amplification fragments that ranged between 60 and 158bp.     

Coronaviruses Use ACE2 Monomers as Entry-Receptors

The angiotensin-converting enzyme 2 (ACE2) has been identified as entry receptor on cells enabling binding and infection with the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) via trimeric spike (S) proteins protruding from the viral surface. It has been suggested that trimeric S proteins preferably bind to plasma membrane areas with high concentrations of possibly multimeric ACE2 receptors to achieve a higher binding and infection efficiency. Here we used direct stochastic optical reconstruction microscopy (dSTORM) in combination with different labeling approaches to visualize the distribution and quantify the expression of ACE2 on different cells. Our results reveal that endogenous ACE2 receptors are present as monomers in the plasma membrane with densities of only 1–2 receptors μm⁻². In addition, binding of trimeric S proteins does not induce the formation of ACE2 oligomers in the plasma membrane. Supported by infection studies using vesicular stomatitis virus (VSV) particles bearing S proteins our data demonstrate that a single S protein interaction per virus particle with a monomeric ACE2 receptor is sufficient for infection, which provides SARS-CoV-2 a high infectivity.     

Native Structure-Based Peptides as Potential Protein–Protein Interaction Inhibitors of SARS-CoV-2 Spike Protein and Human ACE2 Receptor

Severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2) is a positive-strand RNA virus that causes severe respiratory syndrome in humans, which is now referred to as coronavirus disease 2019 (COVID-19). Since December 2019, the new pathogen has rapidly spread globally, with over 65 million cases reported to the beginning of December 2020, including over 1.5 million deaths. Unfortunately, currently, there is no specific and effective treatment for COVID-19. As SARS-CoV-2 relies on its spike proteins (S) to bind to a host cell-surface receptor angiotensin-converting enzyme-2(ACE2), and this interaction is proved to be responsible for entering a virus into host cells, it makes an ideal target for antiviral drug development. In this work, we design three very short peptides based on the ACE2 sequence/structure fragments, which may effectively bind to the receptor-binding domain (RBD) of S protein and may, in turn, disrupt the important virus-host protein–protein interactions, blocking early steps of SARS-CoV-2 infection. Two of our peptides bind to virus protein with affinity in nanomolar range, and as very short peptides have great potential for drug development.     

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.     

Atomistic De-novo Inhibitor Generation-Guided Drug Repurposing for SARS-CoV-2 Spike Protein with Free-Energy Validation by Well-Tempered Metadynamics

Computational drug design is increasingly becoming important with new and unforeseen diseases like COVID-19. In this study, we present a new computational de novo drug design and repurposing method and applied it to find plausible drug candidates for the receptor binding domain (RBD) of SARS-CoV-2 (COVID-19). Our study comprises three steps: atom-by-atom generation of new molecules around a receptor, structural similarity mapping to existing approved and investigational drugs, and validation of their binding strengths to the viral spike proteins based on rigorous all-atom, explicit-water well-tempered metadynamics free energy calculations. By choosing the receptor binding domain of the viral spike protein, we showed that some of our new molecules and some of the repurposable drugs have stronger binding to RBD than hACE2. To validate our approach, we also calculated the free energy of hACE2 and RBD, and found it to be in an excellent agreement with experiments. These pool of drugs will allow strategic repurposing against COVID-19 for a particular prevailing conditions.     

Studying molecular interactions via capillary electrophoresis and microscale thermophoresis: A review

The process of choosing the most proper technique for studying the molecular interactions is based on critical factors such as instrumentation complexity, automation, experimental procedures, analysis time, consumables, and cost-value. This review has tracked the use of affinity capillary electrophoresis (ACE) and microscale thermophoresis (MST) techniques in the evaluation of molecular binding among different molecules during the 5 years 2016–2021. ACE has proved to be an attractive technique for biomolecular characterization with high resolution efficiency where small variations in several controlling factors can much improve such efficiency compared to other analytical techniques. Meanwhile, MST has proved its higher sensitivity for smaller amounts of complex non-purified biosamples without affecting its robustness while providing high through output. However, the main motivation to review both techniques in the proposed review was their capability to carry out all experiments without the need for immobilizing one interacting partner, besides a great flexibility in the use of buffering systems. The proposed review demonstrates the importance of both techniques in different areas of life sciences. Moreover, the recent advances in exploiting ACE and MST in other research interests have been discussed.