A community effort in SARS-CoV-2 drug discovery

The COVID-19 pandemic continues to pose a substantial threat to human lives and is likely to do so for years to come. Despite the availability of vaccines, searching for efficient small-molecule drugs that are widely available, including in low- and middle-income countries, is an ongoing challenge. In this work, we report the results of an open science community effort, the "Billion molecules against COVID-19 challenge", to identify small-molecule inhibitors against SARS-CoV-2 or relevant human receptors. Participating teams used a wide variety of computational methods to screen a minimum of 1 billion virtual molecules against 6 protein targets. Overall, 31 teams participated, and they suggested a total of 639,024 molecules, which were subsequently ranked to find 'consensus compounds'. The organizing team coordinated with various contract research organizations (CROs) and collaborating institutions to synthesize and test 878 compounds for biological activity against proteases (Nsp5, Nsp3, TMPRSS2), nucleocapsid N, RdRP (only the Nsp12 domain), and (alpha) spike protein S. Overall, 27 compounds with weak inhibition/binding were experimentally identified by binding-, cleavage-, and/or viral suppression assays and are presented here. Open science approaches such as the one presented here contribute to the knowledge base of future drug discovery efforts in finding better SARS-CoV-2 treatments.

Identification and validation of fusidic acid and flufenamic acid as inhibitors of SARS-CoV-2 replication using DrugSolver CavitomiX

In this work, we present DrugSolver CavitomiX, a novel computational pipeline for drug repurposing and identifying ligands and inhibitors of target enzymes. The pipeline is based on cavity point clouds representing physico-chemical properties of the cavity induced solely by the protein. To test the pipeline's ability to identify inhibitors, we chose enzymes essential for SARS-CoV-2 replication as a test system. The active-site cavities of the viral enzymes main protease (M^pro) and papain-like protease (Pl^pro), as well as of the human transmembrane serine protease 2 (TMPRSS2), were selected as target cavities. Using active-site point-cloud comparisons, it was possible to identify two compounds-flufenamic acid and fusidic acid-which show strong inhibition of viral replication. The complexes from which fusidic acid and flufenamic acid were derived would not have been identified using classical sequence- and structure-based methods as they show very little structural (TM-score: 0.1 and 0.09, respectively) and very low sequence (~ 5%) identity to M^pro and TMPRSS2, respectively. Furthermore, a cavity-based off-target screening was performed using acetylcholinesterase (AChE) as an example. Using cavity comparisons, the human carboxylesterase was successfully identified, which is a described off-target for AChE inhibitors.

Optimizing variant-specific therapeutic SARS-CoV-2 decoys using deep-learning-guided molecular dynamics simulations

Treatment of COVID-19 with a soluble version of ACE2 that binds to SARS-CoV-2 virions before they enter host cells is a promising approach, however it needs to be optimized and adapted to emerging viral variants. The computational workflow presented here consists of molecular dynamics simulations for spike RBD-hACE2 binding affinity assessments of multiple spike RBD/hACE2 variants and a novel convolutional neural network architecture working on pairs of voxelized force-fields for efficient search-space reduction. We identified hACE2-Fc K31W and multi-mutation variants as high-affinity candidates, which we validated in vitro with virus neutralization assays. We evaluated binding affinities of these ACE2 variants with the RBDs of Omicron BA.3, Omicron BA.4/BA.5, and Omicron BA.2.75 in silico. In addition, candidates produced in Nicotiana benthamiana, an expression organism for potential large-scale production, showed a 4.6-fold reduction in half-maximal inhibitory concentration (IC50) compared with the same variant produced in CHO cells and an almost six-fold IC50 reduction compared with wild-type hACE2-Fc.

Structural bioinformatics analysis of SARS-CoV-2 variants reveals higher hACE2 receptor binding affinity for Omicron B.1.1.529 spike RBD compared to wild type reference

To date, more than 263 million people have been infected with SARS-CoV-2 during the COVID-19 pandemic. In many countries, the global spread occurred in multiple pandemic waves characterized by the emergence of new SARS-CoV-2 variants. Here we report a sequence and structural-bioinformatics analysis to estimate the effects of amino acid substitutions on the affinity of the SARS-CoV-2 spike receptor binding domain (RBD) to the human receptor hACE2. This is done through qualitative electrostatics and hydrophobicity analysis as well as molecular dynamics simulations used to develop a high-precision empirical scoring function (ESF) closely related to the linear interaction energy method and calibrated on a large set of experimental binding energies. For the latest variant of concern (VOC), B.1.1.529 Omicron, our Halo difference point cloud studies reveal the largest impact on the RBD binding interface compared to all other VOC. Moreover, according to our ESF model, Omicron achieves a much higher ACE2 binding affinity than the wild type and, in particular, the highest among all VOCs except Alpha and thus requires special attention and monitoring.

A multi-pronged approach targeting SARS-CoV-2 proteins using ultra-large virtual screening

The unparalleled global effort to combat the continuing severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic over the last year has resulted in promising prophylactic measures. However, a need still exists for cheap, effective therapeutics, and targeting multiple points in the viral life cycle could help tackle the current, as well as future, coronaviruses. Here, we leverage our recently developed, ultra-large-scale in silico screening platform, VirtualFlow, to search for inhibitors that target SARS-CoV-2. In this unprecedented structure-based virtual campaign, we screened roughly 1 billion molecules against each of 40 different target sites on 17 different potential viral and host targets. In addition to targeting the active sites of viral enzymes, we also targeted critical auxiliary sites such as functionally important protein-protein interactions.

A Multi-Pronged Approach Targeting SARS-CoV-2 Proteins Using Ultra-Large Virtual Screening

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), previously known as 2019 novel coronavirus (2019-nCoV), has spread rapidly across the globe, creating an unparalleled global health burden and spurring a deepening economic crisis. As of July 7th, 2020, almost seven months into the outbreak, there are no approved vaccines and few treatments available. Developing drugs that target multiple points in the viral life cycle could serve as a strategy to tackle the current as well as future coronavirus pandemics. Here we leverage the power of our recently developed in silico screening platform, VirtualFlow, to identify inhibitors that target SARS-CoV-2. VirtualFlow is able to efficiently harness the power of computing clusters and cloud-based computing platforms to carry out ultra-large scale virtual screens. In this unprecedented structure-based multi-target virtual screening campaign, we have used VirtualFlow to screen an average of approximately 1 billion molecules against each of 40 different target sites on 17 different potential viral and host targets in the cloud. In addition to targeting the active sites of viral enzymes, we also target critical auxiliary sites such as functionally important protein-protein interaction interfaces. This multi-target approach not only increases the likelihood of finding a potent inhibitor, but could also help identify a collection of anti-coronavirus drugs that would retain efficacy in the face of viral mutation. Drugs belonging to different regimen classes could be combined to develop possible combination therapies, and top hits that bind at highly conserved sites would be potential candidates for further development as coronavirus drugs. Here, we present the top 200 in silico hits for each target site. While in-house experimental validation of some of these compounds is currently underway, we want to make this array of potential inhibitor candidates available to researchers worldwide in consideration of the pressing need for fast-tracked drug development.

A nontoxic pain killer designed by modeling of pathological receptor conformations

Indiscriminate activation of opioid receptors provides pain relief but also severe central and intestinal side effects. We hypothesized that exploiting pathological (rather than physiological) conformation dynamics of opioid receptor-ligand interactions might yield ligands without adverse actions. By computer simulations at low pH, a hallmark of injured tissue, we designed an agonist that, because of its low acid dissociation constant, selectively activates peripheral μ-opioid receptors at the source of pain generation. Unlike the conventional opioid fentanyl, this agonist showed pH-sensitive binding, heterotrimeric guanine nucleotide-binding protein (G protein) subunit dissociation by fluorescence resonance energy transfer, and adenosine 3',5'-monophosphate inhibition in vitro. It produced injury-restricted analgesia in rats with different types of inflammatory pain without exhibiting respiratory depression, sedation, constipation, or addiction potential.

Effects of carbamazepine and two of its metabolites on the non-biting midge Chironomus riparius in a sediment full life cycle toxicity test

The antiepileptic drug carbamazepine (CBZ) and its main metabolites carbamazepine-10,11-epoxide (EP-CBZ) and 10,11-dihydro-10,11-dihydroxy-carbamazepine (DiOH-CBZ) were chosen as test substances to assess chronic toxicity on the non-biting midge Chironomus riparius. All the three substances were tested in a 40-day sediment full life cycle test (according to OECD 233) in which mortality, emergence, fertility, and clutch size were evaluated. In addition, these parameters were considered to calculate the population growth rate which represents an integrated measure to assess population relevant effects. With an LC50 of 0.20 mg/kg (time-weighted mean), the metabolite EP-CBZ was significantly more toxic than the parent substance CBZ (LC50: 1.1 mg/kg). Especially mortality, emergence, and fertility showed to be sensitive parameters under the exposure to CBZ and EP-CBZ. By using classical molecular dynamics (MD) simulations, the binding of CBZ to the ecdysone receptor was investigated as one possible mode of action (MoA) but appeared to be unlikely. The second metabolite DiOH-CBZ did not cause any effects within the tested concentration rage (0.17-1.2 mg/kg). Even though CBZ was less toxic compared to EP-CBZ, CBZ is found in the environment at much higher concentrations and therefore causes a higher potential risk for sediment dwelling organisms compared to its metabolites. Nevertheless, the current study illustrates the importance of including commonly found metabolites into the risk assessment of parent substances.

Peptide-polymer ligands for a tandem WW-domain, an adaptive multivalent protein-protein interaction: lessons on the thermodynamic fitness of flexible ligands

Three polymers, poly(N-(2-hydroxypropyl)methacrylamide) (pHPMA), hyperbranched polyglycerol (hPG), and dextran were investigated as carriers for multivalent ligands targeting the adaptive tandem WW-domain of formin-binding protein (FBP21). Polymer carriers were conjugated with 3–9 copies of the proline-rich decapeptide GPPPRGPPPR-NH₂ (P1). Binding of the obtained peptide–polymer conjugates to the tandem WW-domain was investigated employing isothermal titration calorimetry (ITC) to determine the binding affinity, the enthalpic and entropic contributions to free binding energy, and the stoichiometry of binding for all peptide–polymer conjugates. Binding affinities of all multivalent ligands were in the µM range, strongly amplified compared to the monovalent ligand P1 with a K_D > 1 mM. In addition, concise differences were observed, pHPMA and hPG carriers showed moderate affinity and bound 2.3–2.8 peptides per protein binding site resulting in the formation of aggregates. Dextran-based conjugates displayed affinities down to 1.2 µM, forming complexes with low stoichiometry, and no precipitation. Experimental results were compared with parameters obtained from molecular dynamics simulations in order to understand the observed differences between the three carrier materials. In summary, the more rigid and condensed peptide–polymer conjugates based on the dextran scaffold seem to be superior to induce multivalent binding and to increase affinity, while the more flexible and dendritic polymers, pHPMA and hPG are suitable to induce crosslinking upon binding.

Markov model-based polymer assembly from force field-parameterized building blocks

A conventional by hand construction and parameterization of a polymer model for the purpose of molecular simulations can quickly become very work-intensive and time-consuming. Using the example of polyglycerol, I present a polymer decomposition strategy yielding a set of five monomeric residues that are convenient for an instantaneous assembly and subsequent force field simulation of a polyglycerol polymer model. Force field parameters have been developed in accordance with the classical Amber force field. Partial charges of each unit were fitted to the electrostatic potential using quantum-chemical methods and slightly modified in order to guarantee a neutral total polymer charge. In contrast to similarly constructed models of amino acid and nucleotide sequences, the glycerol building blocks may yield an arbitrary degree of bifurcations depending on the underlying probabilistic model. The iterative development of the overall structure as well as the relation of linear to branching units is controlled by a simple Markov model which is presented with few algorithmic details. The resulting polymer is highly suitable for classical explicit water molecular dynamics simulations on the atomistic level after a structural relaxation step. Moreover, the decomposition strategy presented here can easily be adopted to many other (co)polymers.

[Medicine from the computer]

Small molecules can have a significant effect on human metabolic processes. Computational drug design aims at constructing specialized small molecules that selectively and efficiently address specific proteins. The basic ideas of computational molecular design are presented and it will be shown how a virtual protein can be computer designed. This virtual protein can be used to predict the binding affinity of given small molecules without having to synthesize them in a laboratory. Modern computational drug design goes far beyond the lock and key principle. Possible future developments are discussed and a current successful example of computational drug design in the field of painkiller medication is demonstrated.

[Medicine from the computer]

Small molecules can have a significant effect on human metabolic processes. Computational drug design aims at constructing specialized small molecules that selectively and efficiently address specific proteins. The basic ideas of computational molecular design are presented and it will be shown how a virtual protein can be computer designed. This virtual protein can be used to predict the binding affinity of given small molecules without having to synthesize them in a laboratory. Modern computational drug design goes far beyond the lock and key principle. Possible future developments are discussed and a current successful example of computational drug design in the field of painkiller medication is demonstrated.

Photochemical trans-/cis-isomerization and quantitation of zearalenone in edible oils

The emphasis of the present work was to investigate the photochemical conversion of trans- to cis-zearalenone in edible oils under real-life conditions. For quantitation purposes a cis-zearalenone standard was synthesized and characterized for its identity and purity (≥95%) by (1)H NMR, X-ray crystallography, HPLC fluorescence and mass spectrometric detection. In a sample survey of 12 edible oils (9 corn oils, 3 hempseed oils) from local supermarkets all corn oils contained trans-zearalenone (median 194 μg/kg), but no cis-zearalenone was detected. For alteration studies trans-zearalenone contaminated corn oils were exposed to sunlight over 4 and 30 weeks, revealing an obvious shift toward cis-zearalenone up to a cis/trans ratio of 9:1 by storage in colorless glass bottles. Irradiation experiments of trans-zearalenone in different organic solvents confirmed the preferred formation of cis-zearalenone possibly caused by entropic effects rather than by enthalpic entities as investigated by quantum chemical and classical force field simulations.

Evaluation of the apical leakage of isobutyl cyanoacrylate when used as a root canal sealant

The sealing properties of isobutyl cyanoacrylate were evaluated in extracted single-rooted human teeth. In the experimental group, instrumented root canals were obturated with isobutyl cyanoacrylate and gutta-percha by using the lateral condensation technique. In the control group instrumented root canals were obturated with gutta-percha only. The penetration of radioactive 131I was evaluated by an autoradiographic technique. No leakage was observed in the experimental group. This study shows that isobutyl cyanoacrylate is a good apical sealant.