COVID Moonshot: Open Science Discovery of SARS-CoV-2 Main Protease Inhibitors by Combining Crowdsourcing, High-Throughput Experiments, Computational Simulations, and Machine Learning

The COVID-19 pandemic is a stark reminder that a barren global antiviral pipeline has grave humanitarian consequences. Future pandemics could be prevented by accessible, easily deployable broad-spectrum oral antivirals and open knowledge bases that derisk and accelerate novel antiviral discovery and development. Here, we report the results of the COVID Moonshot, a fully open-science structure-enabled drug discovery campaign targeting the SARS-CoV-2 main protease. We discovered a novel chemical scaffold that is differentiated from current clinical candidates in terms of toxicity, resistance, and pharmacokinetics liabilities, and developed it into noncovalent orally-bioavailable nanomolar inhibitors with clinical potential. Our approach leveraged crowdsourcing, high-throughput structural biology, machine learning, and exascale molecular simulations. In the process, we generated a detailed map of the structural plasticity of the main protease, extensive structure-activity relationships for multiple chemotypes, and a wealth of biochemical activity data. In a first for a structure-based drug discovery campaign, all compound designs (>18,000 designs), crystallographic data (>500 ligand-bound X-ray structures), assay data (>10,000 measurements), and synthesized molecules (>2,400 compounds) for this campaign were shared rapidly and openly, creating a rich open and IP-free knowledgebase for future anti-coronavirus drug discovery.

GPathFinder: Identification of Ligand-Binding Pathways by a Multi-Objective Genetic Algorithm.

Protein-ligand docking is a widely used method to generate solutions for the binding of a small molecule with its target in a short amount of time. However, these methods provide identification of physically sound protein-ligand complexes without a complete view of the binding process dynamics, which has been recognized to be a major discriminant in binding affinity and ligand selectivity. In this paper, a novel piece of open-source software to approach this problem is presented, called GPathFinder. It is built as an extension of the modular GaudiMM platform and is able to simulate ligand diffusion pathways at atomistic level. The method has been benchmarked on a set of 20 systems whose ligand-binding routes were studied by other computational tools or suggested from experimental "snapshots". In all of this set, GPathFinder identifies those channels that were already reported in the literature. Interestingly, the low-energy pathways in some cases indicate novel possible binding routes. To show the usefulness of GPathFinder, the analysis of three case systems is reported. We believe that GPathFinder is a software solution with a good balance between accuracy and computational cost, and represents a step forward in extending protein-ligand docking capacities, with implications in several fields such as drug or enzyme design.

Differential response to Cadmium exposure by expression of a two and a three-domain metallothionein isoform in the land winkle Pomatias elegans: Valuating the marine heritage of a land snail.

Through evolution, marine snails have adapted several times independently to terrestrial life. A prime example for such transitions is the adaptation to terrestrial conditions in members of the gastropod clade of Littorinoidea (Caenogastropoda). Some species of this lineage like the periwinkle (Littorina littorea), live in intertidal habitats, where they are intermittently exposed to semi-terrestrial conditions. Pomatias elegans is a close relative of Littorina littorea that has successfully colonized terrestrial habitats. Evolutionary transitions from marine to terrestrial conditions have often been fostered in marine ancestors by acquisition of physiological pre-adaptations to terrestrial life. Such pre-adaptations are based, among others, on the optimization of a wide repertoire of stress resistance mechanisms, such as the expression of metal inactivating metallothioneins (MTs). The objective of our study was to explore the Cd handling strategy in the terrestrial snail Pomatias elegans in comparison to that observed previously in Littorina littorea. After Cd exposure, the metal is accumulated mainly in the midgut gland of Pomatias elegans, in a similar way as in its marine relative. Upon Cd exposure, Pomatias elegans expresses Cd-specific MTs, as also described from Littorina littorea. In contrast to the latter species, however, the detoxification of Cd in Pomatias elegans is mediated by two different MT isoforms, one two-domain and one three-domain MT. Although the MT proteins of both species are homologous and clearly originate from one common ancestor, the three-domain MT isoform of Pomatias elegans has evolved independently from the three-domain MT of its marine counterpart, probably by addition of a third domain to the pre-existing two-domain MT. Obviously, the occurrence of homologous MT structures in both species is a hereditary character, whereas the differentiation into two distinct MT isoforms with different upregulation capacities in Pomatias elegans is an adaptive feature that probably emerged upon transition to life on land.

GARLEEK: Adding an extra flavor to ONIOM.

The ONIOM method, developed in the group of Keiji Morokuma, is one of the most successful examples of quantum mechanics/molecular mechanics (QM/MM) treatments, and of multilayer methods in general. Its implementation in the Gaussian program package is in particular widely used. This implementation has access to the wide variety of QM methods available in Gaussian, but is limited to only three specific force fields. The current article presents the GARLEEK interface, which expands the availability of molecular mechanics methods to the wide variety of force fields available in MM packages. The focus is in the simple installation and use. Two examples of the performance of the interface with selected systems are provided. GARLEEK is MIT-licensed and freely available at https://github.com/insilichem/garleek. © 2018 Wiley Periodicals, Inc.