Many proteins have small molecule-binding pockets that are not easily detectable in the ligand-free structures. These cryptic sites require a conformational change to become apparent; a cryptic site can therefore be defined as a site that forms a pocket in a holo structure, but not in the apo structure. Because many proteins appear to lack druggable pockets, understanding and accurately identifying cryptic sites could expand the set of drug targets. Previously, cryptic sites were identified experimentally by fragment-based ligand discovery, and computationally by long molecular dynamics simulations and fragment docking. Here, we begin by constructing a set of structurally defined apo-holo pairs with cryptic sites. Next, we comprehensively characterize the cryptic sites in terms of their sequence, structure, and dynamics attributes. We find that cryptic sites tend to be as conserved in evolution as traditional binding pockets, but are less hydrophobic and more flexible. Relying on this characterization, we use machine learning to predict cryptic sites with relatively high accuracy (for our benchmark, the true positive and false positive rates are 73% and 29%, respectively). We then predict cryptic sites in the entire structurally characterized human proteome (11,201 structures, covering 23% of all residues in the proteome). CryptoSite increases the size of the potentially “druggable” human proteome from ~40% to ~78% of disease-associated proteins. Finally, to demonstrate the utility of our approach in practice, we experimentally validate a cryptic site in protein tyrosine phosphatase 1B using a covalent ligand and NMR spectroscopy. The CryptoSite web server is available at http://salilab.org/cryptosite.

Determining the interconverting conformations of dynamic proteins in atomic detail is a major challenge for structural biology. Conformational heterogeneity in the active site of the dynamic enzyme cyclophilin A (CypA) has been previously linked to its catalytic function, but the extent to which the different conformations of these residues are correlated is unclear. Here we compare the conformational ensembles of CypA by multitemperature synchrotron crystallography and fixed-target X-ray free-electron laser (XFEL) crystallography. The diffraction-before-destruction nature of XFEL experiments provides a radiation-damage-free view of the functionally important alternative conformations of CypA, confirming earlier synchrotron-based results. We monitored the temperature dependences of these alternative conformations with eight synchrotron datasets spanning 100-310 K. Multiconformer models show that many alternative conformations in CypA are populated only at 240 K and above, yet others remain populated or become populated at 180 K and below. These results point to a complex evolution of conformational heterogeneity between 180-–240 K that involves both thermal deactivation and solvent-driven arrest of protein motions in the crystal. The lack of a single shared conformational response to temperature within the dynamic active-site network provides evidence for a conformation shuffling model, in which exchange between rotamer states of a large aromatic ring in the middle of the network shifts the conformational ensemble for the other residues in the network. Together, our multitemperature analyses and XFEL data motivate a new generation of temperature- and time-resolved experiments to structurally characterize the dynamic underpinnings of protein function.DOI: http://dx.doi.org/10.7554/eLife.07574.001

The matrix 2 (M2) protein from influenza A virus is a proton channel that uses His37 as a selectivity filter. Here we report high-resolution (1.10 Å) cryogenic crystallographic structures of the transmembrane domain of M2 at low and high pH. These structures reveal that waters within the pore form hydrogen-bonded networks or “water wires” spanning 17 Å from the channel entrance to His37. Pore-lining carbonyl groups are well situated to stabilize hydronium via second-shell interactions involving bridging water molecules. In addition, room temperature crystallographic structures indicate that water becomes increasingly fluid with increasing temperature and decreasing pH, despite the higher electrostatic field. Complementary molecular dynamics simulations reveal a collective switch of hydrogen bond orientations that can contribute to the directionality of proton flux as His37 is dynamically protonated and deprotonated in the conduction cycle.

Background: Hypertrophic cardiomyopathy (HCM) is a complex disease partly explained by the effects of individual gene variants on sarcomeric protein biomechanics. At the cellular level, HCM mutations most commonly enhance force production, leading to higher energy demands. Despite significant advances in elucidating sarcomeric structure-function relationships, there is still much to be learned about the mechanisms that link altered cardiac energetics to HCM phenotypes. In this work, we test the hypothesis that changes in cardiac energetics represent a common pathophysiologic pathway in HCM. Methods: We performed a comprehensive multi-omics profile of the molecular (transcripts, metabolites, and complex lipids), ultrastructural, and functional components of HCM energetics using myocardial samples from 27 HCM patients and 13 normal controls (donor hearts). Results: Integrated omics analysis revealed alterations in a wide array of biochemical pathways with major dysregulation in fatty acid metabolism, reduction of acylcarnitines, and accumulation of free fatty acids. HCM hearts showed evidence of global energetic decompensation manifested by a decrease in high energy phosphate metabolites [ATP, ADP, and phosphocreatine (PCr)] and a reduction in mitochondrial genes involved in creatine kinase and ATP synthesis. Accompanying these metabolic derangements, electron microscopy showed an increased fraction of severely damaged mitochondria with reduced cristae density, coinciding with reduced citrate synthase (CS) activity and mitochondrial oxidative respiration. These mitochondrial abnormalities were associated with elevated reactive oxygen species (ROS) and reduced antioxidant defenses. However, despite significant mitochondrial injury, HCM hearts failed to upregulate mitophagic clearance. Conclusions: Overall, our findings suggest that perturbed metabolic signaling and mitochondrial dysfunction are common pathogenic mechanisms in patients with HCM. These results highlight potential new drug targets for attenuation of the clinical disease through improving metabolic function and reducing mitochondrial injury.

The M2 proton channel of influenza A is a drug target that is essential for the reproduction of the flu virus. It is also a model system for the study of selective, unidirectional proton transport across a membrane. Ordered water molecules arranged in "wires" inside the channel pore have been proposed to play a role in both the conduction of protons to the four gating His37 residues and the stabilization of multiple positive charges within the channel. To visualize the solvent in the pore of the channel at room temperature while minimizing the effects of radiation damage, data were collected to a resolution of 1.4 Å using an X-ray free-electron laser (XFEL) at three different pH conditions: pH 5.5, pH 6.5, and pH 8.0. Data were collected on the Inward open state, which is an intermediate that accumulates at high protonation of the His37 tetrad. At pH 5.5, a continuous hydrogen-bonded network of water molecules spans the vertical length of the channel, consistent with a Grotthuss mechanism model for proton transport to the His37 tetrad. This ordered solvent at pH 5.5 could act to stabilize the positive charges that build up on the gating His37 tetrad during the proton conduction cycle. The number of ordered pore waters decreases at pH 6.5 and 8.0, where the Inward open state is less stable. These studies provide a graphical view of the response of water to a change in charge within a restricted channel environment.W ater molecules in transmembrane protein pores participate in the transport of protons across the membrane bilayer. This process has been extensively studied experimentally and through computational simulations, particularly in small channels such as gramicidin A and influenza A M2. The movement of ions through channels is coupled to the diffusion of water through the pore, but protons are transported at a rate that is faster than the diffusion of H 3 O + (1, 2). Instead of diffusing through channels, protons move concertedly across networks of hydrogen-bonded waters through what is known as the Grotthuss mechanism (3-5). This mechanism of proton transport was initially discovered by the behavior of water in solution, and it has also been proposed to occur within membrane proteins containing water-filled pores (6-9). The matrix 2 (M2) protein of influenza A is one of the smallest proton-selective channels found in nature. This makes it an ideal system for studying the involvement of water in the selective transport of protons across the membrane. The M2 protein of influenza A is a tetramer made up of four 97-residue-long monomers. M2 is multifunctional, with different functions lying in different regions of the sequence. Residues 1-22 make up a conserved N-terminal domain that assists the incorporation of M2 into the virion (10) and is absent in influenza B viruses. The transmembrane domain of M2 (residues 22-46) is necessary for tetramerization (11) and forms a proton-selective channel (12-15) that is the target of the adamantane class of drugs, amantadine and rimantadine (11,(16)(17)(18). The transme...