Open science discovery of potent noncovalent SARS-CoV-2 main protease inhibitors.

We report the results of the COVID Moonshot, a fully open-science, crowdsourced, and structure-enabled drug discovery campaign targeting the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main protease. We discovered a noncovalent, nonpeptidic inhibitor scaffold with lead-like properties that is differentiated from current main protease inhibitors. Our approach leveraged crowdsourcing, machine learning, exascale molecular simulations, and high-throughput structural biology and chemistry. We generated a detailed map of the structural plasticity of the SARS-CoV-2 main protease, extensive structure-activity relationships for multiple chemotypes, and a wealth of biochemical activity data. All compound designs (>18,000 designs), crystallographic data (>490 ligand-bound x-ray structures), assay data (>10,000 measurements), and synthesized molecules (>2400 compounds) for this campaign were shared rapidly and openly, creating a rich, open, and intellectual property-free knowledge base for future anticoronavirus drug discovery.

A Theoretical Mechanistic Study of the Asymmetric Desymmetrization of a Cyclic meso-Anhydride by a Bifunctional Quinine Sulfonamide Organocatalyst.

Cinchona alkaloids and their derivatives are widely used as organocatalysts in asymmetric synthesis. In particular, sulfonamide derivatives of cinchona alkaloids are highly enantioselective desymmetrization catalysts in the ring opening of a variety of cyclic anhydrides. To better understand the mechanism of catalysis, as well as to identify the basis for enantioselectivity by this catalyst, we have performed DFT calculations of this reaction with a cyclic meso anhydride. Herein, we report calculations for two reaction pathways, one concerted and one stepwise, for the production of each enantiomer of the desymmetrized product using the complete sulfonamide catalyst I. Our results are consistent with both the enantioselectivity of this transformation and the catalytic role of the quinuclidine moiety. We find that the stepwise pathway is the relevant pathway in the production of the major enantiomer. Our calculations highlight the role of differential distortion of the anhydride-methanol complex in the transition state as the factor leading to stereoselectivity.

The combined fragmentation and systematic molecular fragmentation methods.

Conspectus Chemistry, particularly organic chemistry, is mostly concerned with functional groups: amines, amides, alcohols, ketones, and so forth. This is because the reactivity of molecules can be categorized in terms of the reactions of these functional groups, and by the influence of other adjacent groups in the molecule. These simple truths ought to be reflected in the electronic structure and electronic energy of molecules, as reactivity is determined by electronic structure. However, sophisticated ab initio quantum calculations of the molecular electronic energy usually do not make these truths apparent. In recent years, several computational chemistry groups have discovered methods for estimating the electronic energy as a sum of the energies of small molecular fragments, or small sets of groups. By decomposing molecules into such fragments of adjacent functional groups, researchers can estimate the electronic energy to chemical accuracy; not just qualitative trends, but accurate enough to understand reactivity. In addition, this has the benefit of cutting down on both computational time and cost, as the necessary calculation time increases rapidly with an increasing number of electrons. Even with steady advances in computer technology, progress in the study of large molecules is slow. In this Account, we describe two related "fragmentation" methods for treating molecules, the combined fragmentation method (CFM) and systematic molecular fragmentation (SMF). In addition, we show how we can use the SMF approach to estimate the energy and properties of nonconducting crystals, by fragmenting the periodic crystal structure into relatively small pieces. A large part of this Account is devoted to simple overviews of how the methods work. We also discuss the application of these approaches to calculating reactivity and other useful properties, such as the NMR and vibrational spectra of molecules and crystals. These applications rely on the ability of these fragmentation methods to accurately estimate derivatives of the molecular and crystal energies. Finally, to provide some common applications of CFM and SMF, we present some specific examples of energy calculations for moderately large molecules. For computational chemists, this fragmentation approach represents an important practical advance. It reduces the computer time required to estimate the energies of molecules so dramatically, that accurate calculations of the energies and reactivity of very large organic and biological molecules become feasible.

Trouble with the Many-Body Expansion.

Longstanding conventional wisdom dictates that the widely used Many-Body Expansion (MBE) converges rapidly by the four-body term when applied to large chemical systems. We have found, however, that this is not true for calculations using many common, moderate-sized basis sets such as 6-311++G** and aug-cc-pVDZ. Energy calculations performed on water clusters using these basis sets showed a deceptively small error when the MBE was truncated at the three-body level, while inclusion of four- and five-body contributions drastically increased the error. Moreover, the error per monomer increases with system size, showing that the MBE is unsuitable to apply to large chemical systems when using these basis sets. Through a systematic study, we identified the cause of the poor MBE convergence to be a many-body basis set superposition effect exacerbated by diffuse functions. This was verified by analysis of MO coefficients and the behavior of the MBE with increasing monomer-monomer separation. We also found poor convergence of the MBE when applied to valence-bonded systems, which has implications for molecular fragmentation methods. The findings in this work suggest that calculations involving the MBE must be performed using the full-cluster basis set, using basis sets without diffuse functions, or using a basis set of at least aug-cc-pVTZ quality.