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.
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.
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 ⁻/Å ³ 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.
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 χ 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 σ 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.
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.
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.
The development of low μM inhibitors of the Mycobacterium tuberculosis phosphatase PtpA is reported.
The development of low μM inhibitors of the Mycobacterium tuberculosis phosphatase PtpA is ...reported. The most potent of these inhibitors (Ki=1.4±0.3μM) 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.
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.
PDB Code(s): 2XI4
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.
PDF
