The periplasmic Cu+/Ag+ chaperone CusF features a novel cation-π interaction between a Cu+/Ag+ ion and Trp44 at the metal binding site. The nature and strength of the Cu+/Ag+-Trp44 interactions are investigated using computational methodologies. Quantum mechanical (QM) calculations show that Cu+ and Ag+ interactions with Trp44 are both of similar strength (~14 kcal/mol) and bond order. Quantum-mechanical/molecular mechanical (QM/MM) calculations show that Cu+ binds in a distorted tetrahedral coordination environment in the cation-π interaction-lacking Trp44Met mutant. Molecular dynamics (MD) simulations of CusF in the apo and Cu+ bound states emphasize the importance of the Cu+-Trp44 interactions in protecting Cu+ from water oxidation. The protein structure does not change over the time-scale of hundreds of nanoseconds in the metal bound state. The metal recognition site exhibits small motions in the apo state, but remains largely preorganized towards metal binding. Trp44 remains oriented to form the cation-π interaction in the apo state, and faces an energetic penalty to move away from the metal ion. Cu+ binding quenches the protein’s internal motions in regions linked to binding CusB, suggesting protein motions play an essential role in Cu+ transfer to CusB.

Accurately computing the free energy for biological processes like protein folding or protein-ligand association remains a challenging problem. Both describing the complex intermolecular forces involved and sampling the requisite configuration space make understanding these processes innately difficult. Herein, we address the sampling problem using a novel methodology we term “movable type”. Conceptually it can be understood by analogy with the evolution of printing and, hence, the name movable type. For example, a common approach to the study of protein-ligand complexation involves taking a database of intact drug-like molecules and exhaustively docking them into a binding pocket. This is reminiscent of early woodblock printing where each page had to be laboriously created prior to printing a book. However, printing evolved to an approach where a database of symbols (letters, numerals, etc.) was created and then assembled using a movable type system, which allowed for the creation of all possible combinations of symbols on a given page, thereby, revolutionizing the dissemination of knowledge. Our movable type (MT) method involves the identification of all atom pairs seen in protein-ligand complexes and then creating two databases: one with their associated pairwise distant dependent energies and another associated with the probability of how these pairs can combine in terms of bonds, angles, dihedrals and non-bonded interactions. Combining these two databases coupled with the principles of statistical mechanics allows us to accurately estimate binding free energies as well as the pose of a ligand in a receptor. This method, by its mathematical construction, samples all of configuration space of a selected region (the protein active site here) in one shot without resorting to brute force sampling schemes involving Monte Carlo, genetic algorithms or molecular dynamics simulations making the methodology extremely efficient. Importantly, this method explores the free energy surface eliminating the need to estimate the enthalpy and entropy components individually. Finally, low free energy structures can be obtained via a free energy minimization procedure yielding all low free energy poses on a given free energy surface. Besides revolutionizing the protein-ligand docking and scoring problem this approach can be utilized in a wide range of applications in computational biology which involve the computation of free energies for systems with extensive phase spaces including protein folding, protein-protein docking and protein design.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has killed more than 4 million humans globally, but there is no bona fide Food and Drug Administration–approved drug-like molecule to impede the COVID-19 pandemic. The sluggish pace of traditional therapeutic discovery is poorly suited to producing targeted treatments against rapidly evolving viruses. Here, we used an affinity-based screen of 4 billion DNA-encoded molecules en masse to identify a potent class of virus-specific inhibitors of the SARS-CoV-2 main protease (Mpro) without extensive and time-consuming medicinal chemistry. CDD-1714, the initial three-building-block screening hit (molecular weight [MW] = 542.5 g/mol), was a potent inhibitor (inhibition constant [Ki] = 20 nM). CDD-1713, a smaller two-building-block analog (MW = 353.3 g/mol) of CDD-1714, is a reversible covalent inhibitor of Mpro (Ki = 45 nM) that binds in the protease pocket, has specificity over human proteases, and shows in vitro efficacy in a SARS-CoV-2 infectivity model. Subsequently, key regions of CDD-1713 that were necessary for inhibitory activity were identified and a potent (Ki = 37 nM), smaller (MW = 323.4 g/mol), and metabolically more stable analog (CDD-1976) was generated. Thus, screening of DNA-encoded chemical libraries can accelerate the discovery of efficacious drug-like inhibitors of emerging viral disease targets.