Protein structural ensembles are revealed by redefining X-ray electron density noise.

To increase the power of X-ray crystallography to determine not only the structures but also the motions of biomolecules, we developed methods to address two classic crystallographic problems: putting electron density maps on the absolute scale of e(-)/Å(3) and calculating the noise at every point in the map. We find that noise varies with position and is often six to eight times lower than thresholds currently used in model building. Analyzing the rescaled electron density maps from 485 representative proteins revealed unmodeled conformations above the estimated noise for 45% of side chains and a previously hidden, low-occupancy inhibitor of HIV capsid protein. Comparing the electron density maps in the free and nucleotide-bound structures of three human protein kinases suggested that substrate binding perturbs distinct intrinsic allosteric networks that link the active site to surfaces that recognize regulatory proteins. These results illustrate general approaches to identify and analyze alternative conformations, low-occupancy small molecules, solvent distributions, communication pathways, and protein motions.

DOCK 6: Impact of new features and current docking performance.

This manuscript presents the latest algorithmic and methodological developments to the structure-based design program DOCK 6.7 focused on an updated internal energy function, new anchor selection control, enhanced minimization options, a footprint similarity scoring function, a symmetry-corrected root-mean-square deviation algorithm, a database filter, and docking forensic tools. An important strategy during development involved use of three orthogonal metrics for assessment and validation: pose reproduction over a large database of 1043 protein-ligand complexes (SB2012 test set), cross-docking to 24 drug-target protein families, and database enrichment using large active and decoy datasets (Directory of Useful Decoys [DUD]-E test set) for five important proteins including HIV protease and IGF-1R. Relative to earlier versions, a key outcome of the work is a significant increase in pose reproduction success in going from DOCK 4.0.2 (51.4%) → 5.4 (65.2%) → 6.7 (73.3%) as a result of significant decreases in failure arising from both sampling 24.1% → 13.6% → 9.1% and scoring 24.4% → 21.1% → 17.5%. Companion cross-docking and enrichment studies with the new version highlight other strengths and remaining areas for improvement, especially for systems containing metal ions. The source code for DOCK 6.7 is available for download and free for academic users at http://dock.compbio.ucsf.edu/.

Accessing protein conformational ensembles using room-temperature X-ray crystallography.

Modern protein crystal structures are based nearly exclusively on X-ray data collected at cryogenic temperatures (generally 100 K). The cooling process is thought to introduce little bias in the functional interpretation of structural results, because cryogenic temperatures minimally perturb the overall protein backbone fold. In contrast, here we show that flash cooling biases previously hidden structural ensembles in protein crystals. By analyzing available data for 30 different proteins using new computational tools for electron-density sampling, model refinement, and molecular packing analysis, we found that crystal cryocooling remodels the conformational distributions of more than 35% of side chains and eliminates packing defects necessary for functional motions. In the signaling switch protein, H-Ras, an allosteric network consistent with fluctuations detected in solution by NMR was uncovered in the room-temperature, but not the cryogenic, electron-density maps. These results expose a bias in structural databases toward smaller, overpacked, and unrealistically unique models. Monitoring room-temperature conformational ensembles by X-ray crystallography can reveal motions crucial for catalysis, ligand binding, and allosteric regulation.

DOCK 6: combining techniques to model RNA-small molecule complexes.

With an increasing interest in RNA therapeutics and for targeting RNA to treat disease, there is a need for the tools used in protein-based drug design, particularly DOCKing algorithms, to be extended or adapted for nucleic acids. Here, we have compiled a test set of RNA-ligand complexes to validate the ability of the DOCK suite of programs to successfully recreate experimentally determined binding poses. With the optimized parameters and a minimal scoring function, 70% of the test set with less than seven rotatable ligand bonds and 26% of the test set with less than 13 rotatable bonds can be successfully recreated within 2 A heavy-atom RMSD. When DOCKed conformations are rescored with the implicit solvent models AMBER generalized Born with solvent-accessible surface area (GB/SA) and Poisson-Boltzmann with solvent-accessible surface area (PB/SA) in combination with explicit water molecules and sodium counterions, the success rate increases to 80% with PB/SA for less than seven rotatable bonds and 58% with AMBER GB/SA and 47% with PB/SA for less than 13 rotatable bonds. These results indicate that DOCK can indeed be useful for structure-based drug design aimed at RNA. Our studies also suggest that RNA-directed ligands often differ from typical protein-ligand complexes in their electrostatic properties, but these differences can be accommodated through the choice of potential function. In addition, in the course of the study, we explore a variety of newly added DOCK functions, demonstrating the ease with which new functions can be added to address new scientific questions.

Fragment-based discovery of selective inhibitors of the Mycobacterium tuberculosis protein tyrosine phosphatase PtpA.

The development of low muM inhibitors of the Mycobacterium tuberculosis phosphatase PtpA is reported. The most potent of these inhibitors (K(i)=1.4+/-0.3 microM) was found to be selective when tested against a panel of human tyrosine and dual-specificity phosphatases (11-fold vs the highly homologous HCPtpA, and >70-fold vs all others tested). Modeling the inhibitor-PtpA complexes explained the structure-activity relationships observed in vitro and revealed further possibilities for compound development.

Synthesis and testing of a focused phenothiazine library for binding to HIV-1 TAR RNA.

We have synthesized a series of phenothiazine derivatives, which were used to test the structure-activity relationship of binding to HIV-1 TAR RNA. Variations from our initial compound, 2-acetylphenothiazine, focused on two moieties: ring substitutions and n-alkyl substitutions. Binding characteristics were ascertained via NMR, principally by saturation transfer difference spectra of the ligand and imino proton resonance shifts of the RNA. Both ring and alkyl substitutions manifested NMR changes upon binding. In general, the active site, while somewhat flexible, has regions that can be capitalized for increased binding through van der Waals interactions and others that can be optimized for solubility in subsequent stages of development. However, binding can be nontrivially enhanced several-fold through optimization of van der Waals and hydrophilic sites of the scaffold.

Cyclic AMP regulation of protein lysine acetylation in Mycobacterium tuberculosis.

Protein lysine acetylation networks can regulate central processes such as carbon metabolism and gene expression in bacteria. In Escherichia coli, cyclic AMP (cAMP) regulates protein lysine acetyltransferase (PAT) activity at the transcriptional level, but in Mycobacterium tuberculosis, fusion of a cyclic nucleotide-binding domain to a Gcn5-like PAT domain enables direct cAMP control of protein acetylation. Here we describe the allosteric activation mechanism of M. tuberculosis PAT. The crystal structures of the autoinhibited and cAMP-activated PAT reveal that cAMP binds to a cryptic site in the regulatory domain that is over 32 Å from the catalytic site. An extensive conformational rearrangement relieves this autoinhibition by means of a substrate-mimicking lid that covers the protein-substrate binding surface. A steric double latch couples the domains by harnessing a classic, cAMP-mediated conformational switch. The structures suggest general features that enable the evolution of long-range communication between linked domains.

Backdoor opening mechanism in acetylcholinesterase based on X-ray crystallography and molecular dynamics simulations.

The transient opening of a backdoor in the active-site wall of acetylcholinesterase, one of nature's most rapid enzymes, has been suggested to contribute to the efficient traffic of substrates and products. A crystal structure of Torpedo californica acetylcholinesterase in complex with the peripheral-site inhibitor aflatoxin is now presented, in which a tyrosine at the bottom of the active-site gorge rotates to create a 3.4-Å wide exit channel. Molecular dynamics simulations show that the opening can be further enlarged by movement of Trp84. The crystallographic and molecular dynamics simulation data thus point to the interface between Tyr442 and Trp84 as the key element of a backdoor, whose opening permits rapid clearance of catalysis products from the active site. Furthermore, the crystal structure presented provides a novel template for rational design of inhibitors and reactivators, including anti-Alzheimer drugs and antidotes against organophosphate poisoning.

Automated electron-density sampling reveals widespread conformational polymorphism in proteins.

Although proteins populate large structural ensembles, X-ray diffraction data are traditionally interpreted using a single model. To search for evidence of alternate conformers, we developed a program, Ringer, which systematically samples electron density around the dihedral angles of protein side chains. In a diverse set of 402 structures, Ringer identified weak, nonrandom electron-density features that suggest of the presence of hidden, lowly populated conformations for >18% of uniquely modeled residues. Although these peaks occur at electron-density levels traditionally regarded as noise, statistically significant (P < 10(-5)) enrichment of peaks at successive rotameric chi angles validates the assignment of these features as unmodeled conformations. Weak electron density corresponding to alternate rotamers also was detected in an accurate electron density map free of model bias. Ringer analysis of the high-resolution structures of free and peptide-bound calmodulin identified shifts in ensembles and connected the alternate conformations to ligand recognition. These results show that the signal in high-resolution electron density maps extends below the traditional 1 sigma cutoff, and crystalline proteins are more polymorphic than current crystallographic models. Ringer provides an objective, systematic method to identify previously undiscovered alternate conformations that can mediate protein folding and function.

Development and validation of a modular, extensible docking program: DOCK 5.

We report on the development and validation of a new version of DOCK. The algorithm has been rewritten in a modular format, which allows for easy implementation of new scoring functions, sampling methods and analysis tools. We validated the sampling algorithm with a test set of 114 protein-ligand complexes. Using an optimized parameter set, we are able to reproduce the crystal ligand pose to within 2 A of the crystal structure for 79% of the test cases using our rigid ligand docking algorithm with an average run time of 1 min per complex and for 72% of the test cases using our flexible ligand docking algorithm with an average run time of 5 min per complex. Finally, we perform an analysis of the docking failures in the test set and determine that the sampling algorithm is generally sufficient for the binding pose prediction problem for up to 7 rotatable bonds; i.e. 99% of the rigid ligand docking cases and 95% of the flexible ligand docking cases are sampled successfully. We point out that success rates could be improved through more advanced modeling of the receptor prior to docking and through improvement of the force field parameters, particularly for structures containing metal-based cofactors.