Chemical probes that form a covalent bond with a protein target often show enhanced selectivity, potency, and utility for biological studies. Despite these advantages, protein-reactive compounds are usually avoided in high-throughput screening campaigns. Here we describe a general method (DOCKovalent) for screening large virtual libraries of electrophilic small molecules. We apply this method prospectively to discover reversible covalent fragments that target distinct protein nucleophiles, including the catalytic serine of AmpC β-lactamase and noncatalytic cysteines in RSK2, MSK1, and JAK3 kinases. We identify submicromolar to low-nanomolar hits with high ligand efficiency, cellular activity and selectivity, including the first reported reversible covalent inhibitors of JAK3. Crystal structures of inhibitor complexes with AmpC and RSK2 confirm the docking predictions and guide further optimization. As covalent virtual screening may have broad utility for the rapid discovery of chemical probes, we have made the method freely available through an automated web server (http://covalent.docking.org).

Virtual and high-throughput screens (HTS) should have complementary strengths and weaknesses, but studies that prospectively and comprehensively compare them are rare. We undertook a parallel docking and HTS screen of 197861 compounds against cruzain, a thiol protease target for Chagas disease, looking for reversible, competitive inhibitors. On workup, 99% of the hits were eliminated as false positives, yielding 146 well-behaved, competitive ligands. These fell into five chemotypes: two were prioritized by scoring among the top 0.1% of the docking-ranked library, two were prioritized by behavior in the HTS and by clustering, and one chemotype was prioritized by both approaches. Determination of an inhibitor/cruzain crystal structure and comparison of the high-scoring docking hits to experiment illuminated the origins of docking false-negatives and false-positives. Prioritizing molecules that are both predicted by docking and are HTS-active yields well-behaved molecules, relatively unobscured by the false-positives to which both techniques are individually prone.

Adhesive type 1 pili from uropathogenic Escherichia coli are filamentous protein complexes that are attached to the assembly platform FimD in the outer membrane. During pilus assembly, FimD binds complexes between the chaperone FimC and type 1 pilus subunits in the periplasm and mediates subunit translocation to the cell surface. Here we report nuclear magnetic resonance and X-ray protein structures of the N-terminal substrate recognition domain of FimD (FimD N ) before and after binding of a chaperone-subunit complex. FimD N consists of a flexible N-terminal segment of 24 residues, a structured core with a novel fold, and a C-terminal hinge segment. In the ternary complex, residues 1-24 of FimD N specifically interact with both FimC and the subunit, acting as a sensor for loaded FimC molecules. Together with in vivo complementation studies, we show how this mechanism enables recognition and discrimination of different chaperonesubunit complexes by bacterial pilus assembly platforms.

Conformational change in protein–ligand complexes is widely modeled, but the protein accommodation expected on binding a congeneric series of ligands has received less attention. Given their use in medicinal chemistry, there are surprisingly few substantial series of congeneric ligand complexes in the Protein Data Bank (PDB). Here we determine the structures of eight alkyl benzenes, in single-methylene increases from benzene to n-hexylbenzene, bound to an enclosed cavity in T4 lysozyme. The volume of the apo cavity suffices to accommodate benzene but, even with toluene, larger cavity conformations become observable in the electron density, and over the series two other major conformations are observed. These involve discrete changes in main-chain conformation, expanding the site; few continuous changes in the site are observed. In most structures, two discrete protein conformations are observed simultaneously, and energetic considerations suggest that these conformations are low in energy relative to the ground state. An analysis of 121 lysozyme cavity structures in the PDB finds that these three conformations dominate the previously determined structures, largely modeled in a single conformation. An investigation of the few congeneric series in the PDB suggests that discrete changes are common adaptations to a series of growing ligands. The discrete, but relatively few, conformational states observed here, and their energetic accessibility, may have implications for anticipating protein conformational change in ligand design.

Adhesive type 1 pili from uropathogenic Escherichia coli strains are heat and denaturant resistant, filamentous protein complexes. Individual pilus subunits associate through "donor strand complementation," whereby the incomplete immunoglobulin-like fold of each subunit is completed by the N-terminal extension of a neighboring subunit. We show that antiparallel donor strand insertion generally causes nonequilibrium behavior in protein folding and extreme activation energy barriers for dissociation of subunit-subunit complexes. We identify the most kinetically stable, noncovalent protein complex known to date. The complex between the pilus subunit FimG and the donor strand peptide of the subunit FimF shows an extrapolated dissociation half-life of 3 x 10(9) years. The 15 residue peptide forms ideal intermolecular beta sheet H-bonds with FimG over 10 residues, and its hydrophobic side chains strongly interact with the hydrophobic core of FimG. The results show that kinetic stability and nonequilibrium behavior in protein folding confers infinite stability against dissociation in extracellular protein complexes.