3V: cavity, channel and cleft volume calculator and extractor

The structures of Arabidopsis Deg5 and Deg8 reveal new insights into HtrA proteases

Plant Deg5 and Deg8 are two members of the HtrA proteases, a family of oligomeric serine endopeptidases that are involved in a variety of protein quality-control processes. These two HtrA proteases are located in the thylakoid lumen and participate in high-light stress responses by collaborating with other chloroplast proteins. Deg5 and Deg8 degrade photodamaged D1 protein of the photosystem II reaction centre, allowing its in situ replacement. Here, the crystal structures of Arabidopsis thaliana Deg5 (S266A) and Deg8 (S292A) are reported at 2.6 and 2.0 Å resolution, respectively. The Deg5 trimer contains two calcium ions in a central channel, suggesting a link between photodamage control and calcium ions in chloroplasts. Previous structures of HtrA proteases have indicated that their regulation usually requires C-terminal PDZ domain(s). Deg5 is unique in that it contains no PDZ domain and the trimeric structure of Deg5 (S266A) reveals a novel catalytic triad conformation. A similar triad conformation is observed in the hexameric structure of the single PDZ-domain-containing Deg8 (S292A). These findings suggest a novel activation mechanism for plant HtrA proteases and provide structural clues to their function in light-stress response.

Detergent release prolongs the lifetime of native-like membrane protein conformations in the gas-phase

Recent studies have suggested that detergents can protect the structure of membrane proteins during their transition from solution to the gas-phase. Here we provide mechanistic insights into this process by interrogating the structures of membrane protein-detergent assemblies in the gas-phase using ion mobility mass spectrometry. We show a clear correlation between the population of native-like protein conformations and the degree of detergent attachment to the protein in the gas-phase. Interrogation of these protein-detergent assemblies, by tandem mass spectrometry, enables us to define the mechanism by which detergents preserve native-like protein conformations in a solvent free environment. We show that the release of detergent is more central to the survival of these conformations than the physical presence of detergent bound to the protein. We propose that detergent release competes with structural collapse for the internal energy of the ion and permits the observation of transient native-like membrane protein conformations that are otherwise lost to structural rearrangement in the gas-phase.

Substrate placement influences reactivity in non-heme Fe(II) halogenases and hydroxylases

We employ error-corrected density functional theory methods to map out the dependence of reactivity on substrate position for SyrB2, a member of a family of non-heme iron halogenases and hydroxylases that are only reactive toward amino acid substrates delivered via prosthetic phosphopantetheine arms. For the initial hydrogen abstraction step, the inherent flexibility of the phosphopantetheine molecule weakens the position dependence for both the native substrate (threonine for SyrB2) and alternative substrates. Over a 5 Å window of substrate positions, the tethered hydrogen abstraction step proceeds with nearly identical activation energies and donor-acceptor distances in the transition state. The propensity of a particular substrate toward halogenation or hydroxylation is found to depend strongly on the substrate placement following hydrogen abstraction, with deeper substrate delivery into the active (for native substrates) site favoring halogenation and shallower substrate delivery favoring hydroxylation.

Nonspecific recognition is achieved in Pot1pC through the use of multiple binding modes

Pot1 is the protein responsible for binding to and protecting the 3' single-stranded DNA (ssDNA) overhang at most eukaryotic telomeres. Here, we present the crystal structure of one of the two oligonucleotide/oligosaccharide-binding folds (Pot1pC) that make up the ssDNA-binding domain in S. pombe Pot1. Comparison with the homologous human domain reveals unexpected structural divergence in the mode of ligand binding that explains the differing ligand requirements between species. Despite the presence of apparently base-specific hydrogen bonds, Pot1pC is able to bind a wide range of ssDNA sequences with thermodynamic equivalence. To address how Pot1pC binds ssDNA with little to no specificity, multiple structures of Pot1pC bound to noncognate ssDNA ligands were solved. These structures reveal that this promiscuity is implemented through new binding modes that thermodynamically compensate for base-substitutions through alternate stacking interactions and new H-bonding networks.

Structures of human acetylcholinesterase in complex with pharmacologically important ligands

Human acetylcholinesterase (AChE) is a significant target for therapeutic drugs. Here we present high resolution crystal structures of human AChE, alone and in complexes with drug ligands; donepezil, an Alzheimer's disease drug, binds differently to human AChE than it does to Torpedo AChE. These crystals of human AChE provide a more accurate platform for further drug development than previously available.

Pocket-based drug design: exploring pocket space

The identification and application of druggable pockets of targets play a key role in in silico drug design, which is a fundamental step in structure-based drug design. Herein, some recent progresses and developments of the computational analysis of pockets have been covered. Also, the pockets at the protein-protein interfaces (PPI) have been considered to further explore the pocket space for drug discovery. We have presented two case studies targeting the kinetic pockets generated by normal mode analysis and molecular dynamics method, respectively, in which we focus upon incorporating the pocket flexibility into the two-dimensional virtual screening with both affinity and specificity. We applied the specificity and affinity (SPA) score to quantitatively estimate affinity and evaluate specificity using the intrinsic specificity ratio (ISR) as a quantitative criterion. In one of two cases, we also included some applications of pockets located at the dimer interfaces to emphasize the role of PPI in drug discovery. This review will attempt to summarize the current status of this pocket issue and will present some prospective avenues of further inquiry.

CAVER 3.0: a tool for the analysis of transport pathways in dynamic protein structures

Tunnels and channels facilitate the transport of small molecules, ions and water solvent in a large variety of proteins. Characteristics of individual transport pathways, including their geometry, physico-chemical properties and dynamics are instrumental for understanding of structure-function relationships of these proteins, for the design of new inhibitors and construction of improved biocatalysts. CAVER is a software tool widely used for the identification and characterization of transport pathways in static macromolecular structures. Herein we present a new version of CAVER enabling automatic analysis of tunnels and channels in large ensembles of protein conformations. CAVER 3.0 implements new algorithms for the calculation and clustering of pathways. A trajectory from a molecular dynamics simulation serves as the typical input, while detailed characteristics and summary statistics of the time evolution of individual pathways are provided in the outputs. To illustrate the capabilities of CAVER 3.0, the tool was applied for the analysis of molecular dynamics simulation of the microbial enzyme haloalkane dehalogenase DhaA. CAVER 3.0 safely identified and reliably estimated the importance of all previously published DhaA tunnels, including the tunnels closed in DhaA crystal structures. Obtained results clearly demonstrate that analysis of molecular dynamics simulation is essential for the estimation of pathway characteristics and elucidation of the structural basis of the tunnel gating. CAVER 3.0 paves the way for the study of important biochemical phenomena in the area of molecular transport, molecular recognition and enzymatic catalysis. The software is freely available as a multiplatform command-line application at http://www.caver.cz.

High resolution crystal structures of the Cerebratulus lacteus mini-Hb in the unligated and carbomonoxy states

The nerve tissue mini-hemoglobin from Cerebratulus lacteus (CerHb) displays an essential globin fold hosting a protein matrix tunnel held to allow traffic of small ligands to and from the heme. CerHb heme pocket hosts the distal TyrB10/GlnE7 pair, normally linked to low rates of O₂ dissociation and ultra-high O₂ affinity. However, CerHb affinity for O₂ is similar to that of mammalian myoglobins, due to a dynamic equilibrium between high and low affinity states driven by the ability of ThrE11 to orient the TyrB10 OH group relative to the heme ligand. We present here the high resolution crystal structures of CerHb in the unligated and carbomonoxy states. Although CO binds to the heme with an orientation different from the O₂ ligand, the overall binding schemes for CO and O₂ are essentially the same, both ligands being stabilized through a network of hydrogen bonds based on TyrB10, GlnE7, and ThrE11. No dramatic protein structural changes are needed to support binding of the ligands, which can freely reach the heme distal site through the apolar tunnel. A lack of main conformational changes between the heme-unligated and -ligated states grants stability to the folded mini-Hb and is a prerequisite for fast ligand diffusion to/from the heme.

MOLEonline 2.0: interactive web-based analysis of biomacromolecular channels

Biomolecular channels play important roles in many biological systems, e.g. enzymes, ribosomes and ion channels. This article introduces a web-based interactive MOLEonline 2.0 application for the analysis of access/egress paths to interior molecular voids. MOLEonline 2.0 enables platform-independent, easy-to-use and interactive analyses of (bio)macromolecular channels, tunnels and pores. Results are presented in a clear manner, making their interpretation easy. For each channel, MOLEonline displays a 3D graphical representation of the channel, its profile accompanied by a list of lining residues and also its basic physicochemical properties. The users can tune advanced parameters when performing a channel search to direct the search according to their needs. The MOLEonline 2.0 application is freely available via the Internet at http://ncbr.muni.cz/mole or http://mole.upol.cz.

Software tools for identification, visualization and analysis of protein tunnels and channels

Protein structures contain highly complex systems of voids, making up specific features such as surface clefts or grooves, pockets, protrusions, cavities, pores or channels, and tunnels. Many of them are essential for the migration of solvents, ions and small molecules through proteins, and their binding to the functional sites. Analysis of these structural features is very important for understanding of structure-function relationships, for the design of potential inhibitors or proteins with improved functional properties. Here we critically review existing software tools specialized in rapid identification, visualization, analysis and design of protein tunnels and channels. The strengths and weaknesses of individual tools are reported together with examples of their applications for the analysis and engineering of various biological systems. This review can assist users with selecting a proper software tool for study of their biological problem as well as highlighting possible avenues for further development of existing tools. Development of novel descriptors representing not only geometry, but also electrostatics, hydrophobicity or dynamics, is needed for reliable identification of biologically relevant tunnels and channels.

Structural anatomy of telomere OB proteins

Telomere DNA-binding proteins protect the ends of chromosomes in eukaryotes. A subset of these proteins are constructed with one or more OB folds and bind with G+T-rich single-stranded DNA found at the extreme termini. The resulting DNA-OB protein complex interacts with other telomere components to coordinate critical telomere functions of DNA protection and DNA synthesis. While the first crystal and NMR structures readily explained protection of telomere ends, the picture of how single-stranded DNA becomes available to serve as primer and template for synthesis of new telomere DNA is only recently coming into focus. New structures of telomere OB fold proteins alongside insights from genetic and biochemical experiments have made significant contributions towards understanding how protein-binding OB proteins collaborate with DNA-binding OB proteins to recruit telomerase and DNA polymerase for telomere homeostasis. This review surveys telomere OB protein structures alongside highly comparable structures derived from replication protein A (RPA) components, with the goal of providing a molecular context for understanding telomere OB protein evolution and mechanism of action in protection and synthesis of telomere DNA.

Mitochondrial localization and structure-based phosphate activation mechanism of Glutaminase C with implications for cancer metabolism

Glutamine is an essential nutrient for cancer cell proliferation, especially in the context of citric acid cycle anaplerosis. In this manuscript we present results that collectively demonstrate that, of the three major mammalian glutaminases identified to date, the lesser studied splice variant of the gene gls, known as Glutaminase C (GAC), is important for tumor metabolism. We show that, although levels of both the kidney-type isoforms are elevated in tumor vs. normal tissues, GAC is distinctly mitochondrial. GAC is also most responsive to the activator inorganic phosphate, the content of which is supposedly higher in mitochondria subject to hypoxia. Analysis of X-ray crystal structures of GAC in different bound states suggests a mechanism that introduces the tetramerization-induced lifting of a "gating loop" as essential for the phosphate-dependent activation process. Surprisingly, phosphate binds inside the catalytic pocket rather than at the oligomerization interface. Phosphate also mediates substrate entry by competing with glutamate. A greater tendency to oligomerize differentiates GAC from its alternatively spliced isoform and the cycling of phosphate in and out of the active site distinguishes it from the liver-type isozyme, which is known to be less dependent on this ion.

Voronoi-based extraction and visualization of molecular paths

Visual analysis is widely used to study the behavior of molecules. Of particular interest are the analysis of molecular interactions and the investigation of binding sites. For large molecules, however, it is difficult to detect possible binding sites and paths leading to these sites by pure visual inspection. In this paper, we present new methods for the computation and visualization of potential molecular paths. Using a novel filtering method, we extract the significant paths from the Voronoi diagram of spheres. For the interactive visualization of molecules and their paths, we present several methods using deferred shading and other state-of-the-art techniques. To allow for a fast overview of reachable regions of the molecule, we illuminate the molecular surface using a large number of light sources placed on the extracted paths. We also provide a method to compute the extension surface of selected paths and visualize it using the skin surface. Furthermore, we use the extension surface to clip the molecule to allow easy visual tracking of even deeply buried paths. The methods are applied to several proteins to demonstrate their usefulness.

Identifying continuous pores in protein structures with PROPORES by computational repositioning of gating residues

Proteins containing concavities such as pockets, cavities, and tunnels or pores perform important functions in ligand-induced signal transduction, enzymatic catalysis, and in facilitating the permeation of small molecules through membranes. Computational algorithms for identifying such shapes are therefore of great use for studying the mechanisms of these reactions. We developed the novel toolkit PROPORES for pore identification and applied our program to the systems aquaporin, tryptophan synthase, leucine transporter, and acetylcholinesterase. As a novel feature, the program checks whether access to occluded ligand binding pockets or blocked channels can be achieved by systematically rotating side chains of the gating residues. In this way, we obtain a more flexible view of the putative structural adaptability of protein structures. For the four systems mentioned, the new method was able to identify connections between pores that are separated in the X-ray structures or to connect internal pores with the protein surrounding. The software is available from http://gepard.bioinformatik.uni-saarland.de/software/propores/.

The crystal structure of necrosis- and ethylene-inducing protein 2 from the causal agent of cacao's Witches' Broom disease reveals key elements for its activity

The necrosis- and ethylene-inducing peptide 1 (NEP1)-like proteins (NLPs) are proteins secreted from bacteria, fungi and oomycetes, triggering immune responses and cell death in dicotyledonous plants. Genomic-scale studies of Moniliophthora perniciosa, the fungus that causes the Witches' Broom disease in cacao, which is a serious economic concern for South and Central American crops, have identified five members of this family (termed MpNEP1-5). Here, we show by RNA-seq that MpNEP2 is virtually the only NLP expressed during the fungus infection. The quantitative real-time polymerase chain reaction results revealed that MpNEP2 has an expression pattern that positively correlates with the necrotic symptoms, with MpNEP2 reaching its highest level of expression at the advanced necrotic stage. To improve our understanding of MpNEP2's molecular mechanism of action, we determined the crystallographic structure of MpNEP2 at 1.8 Å resolution, unveiling some key structural features. The implications of a cation coordination found in the crystal structure were explored, and we show that MpNEP2, in contrast to another previously described member of the NLP family, NLP(Pya) from Pythium aphanidermatum, does not depend on an ion to accomplish its necrosis- and electrolyte leakage-promoting activities. Results of site-directed mutagenesis experiments confirmed the importance of a negatively charged cavity and an unforeseen hydrophobic β-hairpin loop for MpNEP2 activity, thus offering a platform for compound design with implications for disease control. Electron paramagnetic resonance and fluorescence assays with MpNEP2 performed in the presence of lipid vesicles of different compositions showed no sign of interaction between the protein and the lipids, implying that MpNEP2 likely requires other anchoring elements from the membrane to promote cytolysis or send death signals.

Structural insight of dopamine β-hydroxylase, a drug target for complex traits, and functional significance of exonic single nucleotide polymorphisms

Background

Human dopamine β-hydroxylase (DBH) is an important therapeutic target for complex traits. Several single nucleotide polymorphisms (SNPs) have also been identified in DBH with potential adverse physiological effect. However, difficulty in obtaining diffractable crystals and lack of a suitable template for modeling the protein has ensured that neither crystallographic three-dimensional structure nor computational model for the enzyme is available to aid rational drug design, prediction of functional significance of SNPs or analytical protein engineering.

Principal Findings

Adequate biochemical information regarding human DBH, structural coordinates for peptidylglycine alpha-hydroxylating monooxygenase and computational data from a partial model of rat DBH were used along with logical manual intervention in a novel way to build an in silico model of human DBH. The model provides structural insight into the active site, metal coordination, subunit interface, substrate recognition and inhibitor binding. It reveals that DOMON domain potentially promotes tetramerization, while substrate dopamine and a potential therapeutic inhibitor nepicastat are stabilized in the active site through multiple hydrogen bonding. Functional significance of several exonic SNPs could be described from a structural analysis of the model. The model confirms that SNP resulting in Ala318Ser or Leu317Pro mutation may not influence enzyme activity, while Gly482Arg might actually do so being in the proximity of the active site. Arg549Cys may cause abnormal oligomerization through non-native disulfide bond formation. Other SNPs like Glu181, Glu250, Lys239 and Asp290 could potentially inhibit tetramerization thus affecting function.

Conclusions

The first three-dimensional model of full-length human DBH protein was obtained in a novel manner with a set of experimental data as guideline for consistency of in silico prediction. Preliminary physicochemical tests validated the model. The model confirms, rationalizes and provides structural basis for several biochemical data and claims testable hypotheses regarding function. It provides a reasonable template for drug design as well.

A few low-frequency normal modes predominantly contribute to conformational responses of hen egg white lysozyme in the tetragonal crystal to variations of molecular packing controlled by environmental humidity

The structures of proteins in crystals are fixed by molecular interactions with neighboring molecules, except in non-contacting flexible regions. Thus, it is difficult to imagine what conformational changes occur in solution. However, if molecular interactions can be changed by manipulating molecular packing in crystals, it may be possible to visualize conformational responses of proteins at atomic resolution by diffraction experiments. For this purpose, it is suitable to control the molecular packing in protein crystals by changing the volume of solvent channels through variation of the environmental relative humidity. Here, we studied conformational responses of hen egg white lysozyme (HEWL) in the tetragonal crystal by X-ray diffraction experiments using a humidity-control apparatus, which provided air flow of 20-98%rh at 298 K. First, we monitored the lattice parameters and crystalline order during dehydration and rehydration of HEWL crystal between 61 and 94%rh at 300 K. Then two crystal structures at a resolution of 2.1 Å using diffraction data obtained at 84.2 and 71.9%rh were determined to discuss the conformational responses of HEWL against the external perturbation induced by changes in molecular packing. The structure at 71.9%rh displayed a closure movement that was likely induced by the molecular contacts formed during dehydration and could be approximated by ten low-frequency normal modes for the crystal structure obtained at 84.2%rh. In addition, we observed reorganization of hydration structures at the molecular interfaces between symmetry neighbors. These findings suggest that humidity-controlled X-ray crystallography is an effective tool to investigate the responses of inherent intramolecular motions of proteins to external perturbations.

dxTuber: detecting protein cavities, tunnels and clefts based on protein and solvent dynamics

Empty space in a protein structure can provide valuable insight into protein properties such as internal hydration, structure stabilization, substrate translocation, storage compartments or binding sites. This information can be visualized by means of cavity analysis. Numerous tools are available depicting cavities directly or identifying lining residues. So far, all available techniques base on a single conformation neglecting any form of protein and cavity dynamics. Here we report a novel, grid-based cavity detection method that uses protein and solvent residence probabilities derived from molecular dynamics simulations to identify (I) internal cavities, (II) tunnels or (III) clefts on the protein surface. Driven by a graphical user interface, output can be exported in PDB format where cavities are described as individually selectable groups of adjacent voxels representing regions of high solvent residence probability. Cavities can be analyzed in terms of solvent density, cavity volume and cross-sectional area along a principal axis. To assess dxTuber performance we performed test runs on a set of six example proteins representing the three main classes of protein cavities and compared our findings to results obtained with SURFNET, CAVER and PyMol.

Identification and characterization of two families of F420 H2-dependent reductases from Mycobacteria that catalyse aflatoxin degradation

Aflatoxins are polyaromatic mycotoxins that contaminate a range of food crops as a result of fungal growth and contribute to serious health problems in the developing world because of their toxicity and mutagenicity. Although relatively resistant to biotic degradation, aflatoxins can be metabolized by certain species of Actinomycetales. However, the enzymatic basis for their breakdown has not been reported until now. We have identified nine Mycobacterium smegmatis enzymes that utilize the deazaflavin cofactor F₄₂₀H₂ to catalyse the reduction of the α,β-unsaturated ester moiety of aflatoxins, activating the molecules for spontaneous hydrolysis and detoxification. These enzymes belong to two previously uncharacterized F₄₂₀H₂ dependent reductase (FDR-A and -B) families that are distantly related to the flavin mononucleotide (FMN) dependent pyridoxamine 5′-phosphate oxidases (PNPOxs). We have solved crystal structures of an enzyme from each FDR family and show that they, like the PNPOxs, adopt a split barrel protein fold, although the FDRs also possess an extended and highly charged F₄₂₀H₂ binding groove. A general role for these enzymes in xenobiotic metabolism is discussed, including the observation that the nitro-reductase Rv3547 from Mycobacterium tuberculosis that is responsible for the activation of bicyclic nitroimidazole prodrugs belongs to the FDR-A family.