MOLE: A Voronoi Diagram-Based Explorer of Molecular Channels, Pores, and Tunnels

Crystal structure analysis of a fatty acid double-bond hydratase from Lactobacillus acidophilus

Bacteria have evolved mechanisms for the hydrogenation of unsaturated fatty acids. Hydroxy fatty acid formation may be the first step in such a process; however, knowledge of the structural and mechanistic aspects of this reaction is scarce. Recently, myosin cross-reactive antigen was shown to be a bacterial FAD-containing hydratase which acts on the 9Z and 12Z double bonds of C16 and C18 non-esterified fatty acids, with the formation of 10-hydroxy and 10,13-dihydroxy fatty acids. These fatty acid hydratases form a large protein family which is conserved across Gram-positive and Gram-negative bacteria with no sequence similarity to any known protein apart from the FAD-binding motif. In order to shed light on the substrate recognition and the mechanism of the hydratase reaction, the crystal structure of the hydratase from Lactobacillus acidophilus (LAH) was determined by single-wavelength anomalous dispersion. Crystal structures of apo LAH and of LAH with bound linoleic acid were refined at resolutions of 2.3 and 1.8 Å, respectively. LAH is a homodimer; each protomer consists of four intricately connected domains. Three of them form the FAD-binding and substrate-binding sites and reveal structural similarity to three domains of several flavin-dependent enzymes, including amine oxidoreductases. The additional fourth domain of LAH is located at the C-terminus and consists of three α-helices. It covers the entrance to the hydrophobic substrate channel leading from the protein surface to the active site. In the presence of linoleic acid, the fourth domain of one protomer undergoes conformational changes and opens the entrance to the substrate-binding channel of the other protomer of the LAH homodimer. The linoleic acid molecule is bound at the entrance to the substrate channel, suggesting movement of the lid domain triggered by substrate recognition.

Free energy landscape for the binding process of Huperzine A to acetylcholinesterase

Drug-target residence time (t = 1/k(off), where k(off) is the dissociation rate constant) has become an important index in discovering better- or best-in-class drugs. However, little effort has been dedicated to developing computational methods that can accurately predict this kinetic parameter or related parameters, k(off) and activation free energy of dissociation (ΔG(off)≠). In this paper, energy landscape theory that has been developed to understand protein folding and function is extended to develop a generally applicable computational framework that is able to construct a complete ligand-target binding free energy landscape. This enables both the binding affinity and the binding kinetics to be accurately estimated. We applied this method to simulate the binding event of the anti-Alzheimer's disease drug (-)-Huperzine A to its target acetylcholinesterase (AChE). The computational results are in excellent agreement with our concurrent experimental measurements. All of the predicted values of binding free energy and activation free energies of association and dissociation deviate from the experimental data only by less than 1 kcal/mol. The method also provides atomic resolution information for the (-)-Huperzine A binding pathway, which may be useful in designing more potent AChE inhibitors. We expect this methodology to be widely applicable to drug discovery and development.

Structure of a bacterial voltage-gated sodium channel pore reveals mechanisms of opening and closing

Sodium-gated ion channels open and close in response to the flow of ions. Here, McCusker et al. report the open structure of a sodium-gated ion channel pore from a bacterial homologue, and show, by comparison with the closed structure, that the movement of a C-terminal helix is sufficient to open the channel.

Voltage-gated sodium channels are vital membrane proteins essential for electrical signalling; in humans, they are key targets for the development of pharmaceutical drugs. Here we report the crystal structure of an open-channel conformation of NavMs, the bacterial channel pore from the marine bacterium Magnetococcus sp. (strain MC-1). It differs from the recently published crystal structure of a closed form of a related bacterial sodium channel (NavAb) by having its internal cavity accessible to the cytoplasmic surface as a result of a bend/rotation about a central residue in the carboxy-terminal transmembrane segment. This produces an open activation gate of sufficient diameter to allow hydrated sodium ions to pass through. Comparison of the open and closed structures provides new insight into the features of the functional states present in the activation cycles of sodium channels and the mechanism of channel opening and closing.

Alpha shape and Delaunay triangulation in studies of protein-related interactions

In recent years, more 3D protein structures have become available, which has made the analysis of large molecular structures much easier. There is a strong demand for geometric models for the study of protein-related interactions. Alpha shape and Delaunay triangulation are powerful tools to represent protein structures and have advantages in characterizing the surface curvature and atom contacts. This review presents state-of-the-art applications of alpha shape and Delaunay triangulation in the studies on protein-DNA, protein-protein, protein-ligand interactions and protein structure analysis.

Crystal structure of a voltage-gated K+ channel pore module in a closed state in lipid membranes

Voltage-gated K(+) channels underlie the electrical excitability of cells. Each subunit of the functional tetramer consists of the tandem fusion of two modules, an N-terminal voltage-sensor and a C-terminal pore. To investigate how sensor coupling to the pore generates voltage-dependent channel opening, we solved the crystal structure and characterized the function of a voltage-gated K(+) channel pore in a lipid membrane. The structure of a functional channel in a membrane environment at 3.1 Å resolution establishes an unprecedented connection between channel structure and function. The structure is unique in delineating an ion-occupied ready to conduct selectivity filter, a confined aqueous cavity, and a closed activation gate, embodying a dynamic entity trapped in an unstable closed state.

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.

RNA polymerase II transcription: structure and mechanism

A minimal RNA polymerase II (pol II) transcription system comprises the polymerase and five general transcription factors (GTFs) TFIIB, -D, -E, -F, and -H. The addition of Mediator enables a response to regulatory factors. The GTFs are required for promoter recognition and the initiation of transcription. Following initiation, pol II alone is capable of RNA transcript elongation and of proofreading. Structural studies reviewed here reveal roles of GTFs in the initiation process and shed light on the transcription elongation mechanism. This article is part of a Special Issue entitled: RNA Polymerase II Transcript Elongation.

Partial least-squares functional mode analysis: application to the membrane proteins AQP1, Aqy1, and CLC-ec1

We introduce an approach based on the recently introduced functional mode analysis to identify collective modes of internal dynamics that maximally correlate to an external order parameter of functional interest. Input structural data can be either experimentally determined structure ensembles or simulated ensembles, such as molecular dynamics trajectories. Partial least-squares regression is shown to yield a robust solution to the multidimensional optimization problem, with a minimal and controllable risk of overfitting, as shown by extensive cross-validation. Several examples illustrate that the partial least-squares-based functional mode analysis successfully reveals the collective dynamics underlying the fluctuations in selected functional order parameters. Applications to T4 lysozyme, the Trp-cage, the aquaporin channels Aqy1 and hAQP1, and the CLC-ec1 chloride antiporter are presented in which the active site geometry, the hydrophobic solvent-accessible surface, channel gating dynamics, water permeability (p(f)), and a dihedral angle are defined as functional order parameters. The Aqy1 case reveals a gating mechanism that connects the inner channel gating residues with the protein surface, thereby providing an explanation of how the membrane may affect the channel. hAQP1 shows how the p(f) correlates with structural changes around the aromatic/arginine region of the pore. The CLC-ec1 application shows how local motions of the gating Glu(148) couple to a collective motion that affects ion affinity in the pore.

Biochemical and structural characterization of the trans-enoyl-CoA reductase from Treponema denticola

The production of fatty acids is an important cellular pathway for both cellular function and the development of engineered pathways for the synthesis of advanced biofuels. Despite the conserved reaction chemistry of various fatty acid synthase systems, the individual isozymes that catalyze these steps are quite diverse in their structural and biochemical features and are important for controlling differences at the cellular level. One of the key steps in the fatty acid elongation cycle is the enoyl-ACP (CoA) reductase function that drives the equilibrium forward toward chain extension. In this work, we report the structural and biochemical characterization of the trans-enoyl-CoA reductase from Treponema denticola (tdTer), which has been utilized for the engineering of synthetic biofuel pathways with an order of magnitude increase in product titers compared to those of pathways constructed with other enoyl-CoA reductase components. The crystal structure of tdTer was determined to 2.00 Å resolution and shows that the Ter enzymes are distinct from members of the FabI, FabK, and FabL families but are highly similar to members of the FabV family. Further biochemical studies show that tdTer uses an ordered bi-bi mechanism initiated by binding of the NADH redox cofactor, which is consistent with the behavior of other enoyl-ACP (CoA) reductases. Mutagenesis of the substrate binding loop, characterization of enzyme activity with respect to crotonyl-CoA, hexenoyl-CoA, and dodecenoyl-CoA substrates, and product inhibition by lauroyl-CoA suggest that this region is important for controlling chain length specificity, with the major portal playing a more important role for longer chain length substrates.

Detecting pore-lining regions in transmembrane protein sequences

BACKGROUND

Alpha-helical transmembrane channel and transporter proteins play vital roles in a diverse range of essential biological processes and are crucial in facilitating the passage of ions and molecules across the lipid bilayer. However, the experimental difficulties associated with obtaining high quality crystals has led to their significant under-representation in structural databases. Computational methods that can identify structural features from sequence alone are therefore of high importance.

RESULTS

We present a method capable of automatically identifying pore-lining regions in transmembrane proteins from sequence information alone, which can then be used to determine the pore stoichiometry. By labelling pore-lining residues in crystal structures using geometric criteria, we have trained a support vector machine classifier to predict the likelihood of a transmembrane helix being involved in pore formation. Results from testing this approach under stringent cross-validation indicate that prediction accuracy of 72% is possible, while a support vector regression model is able to predict the number of subunits participating in the pore with 62% accuracy.

CONCLUSION

To our knowledge, this is the first tool capable of identifying pore-lining regions in proteins and we present the results of applying it to a data set of sequences with available crystal structures. Our method provides a way to characterise pores in transmembrane proteins and may even provide a starting point for discovering novel routes of therapeutic intervention in a number of important diseases. This software is freely available as source code from: http://bioinf.cs.ucl.ac.uk/downloads/memsat-svm/.

The structure of monoacylglycerol lipase from Bacillus sp. H257 reveals unexpected conservation of the cap architecture between bacterial and human enzymes

Monoacylglycerol lipases (MGLs) catalyse the hydrolysis of monoacylglycerol into free fatty acid and glycerol. MGLs have been identified throughout all genera of life and have adopted different substrate specificities depending on their physiological role. In humans, MGL plays an integral part in lipid metabolism affecting energy homeostasis, signalling processes and cancer cell progression. In bacteria, MGLs degrade short-chain monoacylglycerols which are otherwise toxic to the organism. We report the crystal structures of MGL from the bacterium Bacillus sp. H257 (bMGL) in its free form at 1.2Å and in complex with phenylmethylsulfonyl fluoride at 1.8Å resolution. In both structures, bMGL adopts an α/β hydrolase fold with a cap in an open conformation. Access to the active site residues, which were unambiguously identified from the protein structure, is facilitated by two different channels. The larger channel constitutes the highly hydrophobic substrate binding pocket with enough room to accommodate monoacylglycerol. The other channel is rather small and resembles the proposed glycerol exit hole in human MGL. Molecular dynamics simulation of bMGL yielded open and closed states of the entrance channel and the glycerol exit hole. Despite differences in the number of residues, secondary structure elements, and low sequence identity in the cap region, this first structure of a bacterial MGL reveals striking structural conservation of the overall cap architecture in comparison with human MGL. Thus it provides insight into the structural conservation of the cap amongst MGLs throughout evolution and provides a framework for rationalising substrate specificities in each organism.

Allosterically gated enzyme dynamics in the cysteine synthase complex regulate cysteine biosynthesis in Arabidopsis thaliana

Plants and bacteria assimilate sulfur into cysteine. Cysteine biosynthesis involves a bienzyme complex, the cysteine synthase complex (CSC), which consists of serine-acetyl-transferase (SAT) and O-acetyl-serine-(thiol)-lyase (OAS-TL) enzymes. The activity of OAS-TL is reduced by formation of the CSC. Although this reduction is an inherent part of the self-regulation cycle of cysteine biosynthesis, there has until now been no explanation as to how OAS-TL loses activity in plants. Complexation of SAT and OAS-TL involves binding of the C-terminal tail of SAT in one of the active sites of the homodimeric OAS-TL. We here explore the flexibility of the unoccupied active site in Arabidopsis thaliana cytosolic and mitochondrial OAS-TLs. Our results reveal two gates in the OAS-TL active site that define its accessibility. The observed dynamics of the gates show allosteric closure of the unoccupied active site of OAS-TL in the CSC, which can hinder substrate binding, abolishing its turnover to cysteine.

Crystal structure of a voltage-gated sodium channel in two potentially inactivated states

In excitable cells, voltage-gated sodium (Na(V)) channels activate to initiate action potentials and then undergo fast and slow inactivation processes that terminate their ionic conductance. Inactivation is a hallmark of Na(V) channel function and is critical for control of membrane excitability, but the structural basis for this process has remained elusive. Here we report crystallographic snapshots of the wild-type Na(V)Ab channel from Arcobacter butzleri captured in two potentially inactivated states at 3.2 Å resolution. Compared to previous structures of Na(V)Ab channels with cysteine mutations in the pore-lining S6 helices (ref. 4), the S6 helices and the intracellular activation gate have undergone significant rearrangements: one pair of S6 helices has collapsed towards the central pore axis and the other S6 pair has moved outward to produce a striking dimer-of-dimers configuration. An increase in global structural asymmetry is observed throughout our wild-type Na(V)Ab models, reshaping the ion selectivity filter at the extracellular end of the pore, the central cavity and its residues that are analogous to the mammalian drug receptor site, and the lateral pore fenestrations. The voltage-sensing domains have also shifted around the perimeter of the pore module in wild-type Na(V)Ab, compared to the mutant channel, and local structural changes identify a conserved interaction network that connects distant molecular determinants involved in Na(V) channel gating and inactivation. These potential inactivated-state structures provide new insights into Na(V) channel gating and novel avenues to drug development and therapy for a range of debilitating Na(V) channelopathies.

Molecular determinants responsible for recognition of the single-stranded DNA regulatory sequence, χ, by RecBCD enzyme

The RecBCD enzyme is important for both restriction of foreign DNA and recombinational DNA repair. Switching enzyme function from the destructive antiviral state to the productive recombinational state is regulated by the recombination hotspot, χ (5'-GCTGGTGG-3'). Recognition of χ is unique in that it is recognized as a specific sequence within single-stranded DNA (ssDNA) during DNA translocation and unwinding by RecBCD. The molecular determinants of χ recognition and the subsequent alteration in function are unknown. Consequently, we mutated residues within the RecC subunit that comprise a channel where ssDNA is thought to be scanned for a χ sequence. These mutants were characterized in vivo with regard to χ recognition, UV-sensitivity, phage degradation, and recombination proficiency. Of 38 residues mutated, 11 were previously undescribed mutations that altered χ recognition. The mutants fell into two classes: five that failed to respond to χ, and six that suggested a relaxed specificity for χ recognition. The location of the first set of mutations defines a recognition structure responsible for sequence-specific binding of ssDNA. The second set defines a highly conserved structure, linked to the recognition structure, which we hypothesize regulates conversion of RecBCD from a molecular machine that destroys DNA to one that repairs it. These findings offer insight into the evolution of enzymes with alternate χ recognition specificities.

Structure of the Yersinia pestis FabV enoyl-ACP reductase and its interaction with two 2-pyridone inhibitors

The recently discovered FabV enoyl-ACP reductase, which catalyzes the last step of the bacterial fatty acid biosynthesis (FAS-II) pathway, is a promising but unexploited drug target against the reemerging pathogen Yersinia pestis. The structure of Y. pestis FabV in complex with its cofactor reveals that the enzyme features the common architecture of the short-chain dehydrogenase reductase superfamily, but contains additional structural elements that are mostly folded around the usually flexible substrate-binding loop, thereby stabilizing it in a very tight conformation that seals the active site. The structures of FabV in complex with NADH and two newly developed 2-pyridone inhibitors provide insights for the development of new lead compounds, and suggest a mechanism by which the substrate-binding loop opens to admit the inhibitor, a motion that could also be coupled to the interaction of FabV with the acyl-carrier protein substrate.

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.

Investigation of indazole unbinding pathways in CYP2E1 by molecular dynamics simulations

Human microsomal cytochrome P450 2E1 (CYP2E1) can oxidize not only low molecular weight xenobiotic compounds such as ethanol, but also many endogenous fatty acids. The crystal structure of CYP2E1 in complex with indazole reveals that the active site is deeply buried into the protein center. Thus, the unbinding pathways and associated unbinding mechanisms remain elusive. In this study, random acceleration molecular dynamics simulations combined with steered molecular dynamics and potential of mean force calculations were performed to identify the possible unbinding pathways in CYP2E1. The results show that channel 2c and 2a are most likely the unbinding channels of CYP2E1. The former channel is located between helices G and I and the B-C loop, and the latter resides between the region formed by the F-G loop, the B-C loop and the β1 sheet. Phe298 and Phe478 act as the gate keeper during indazole unbinding along channel 2c and 2a, respectively. Previous site-directed mutagenesis experiments also supported these findings.

Computational prediction of metabolism: sites, products, SAR, P450 enzyme dynamics, and mechanisms

Metabolism of xenobiotics remains a central challenge for the discovery and development of drugs, cosmetics, nutritional supplements, and agrochemicals. Metabolic transformations are frequently related to the incidence of toxic effects that may result from the emergence of reactive species, the systemic accumulation of metabolites, or by induction of metabolic pathways. Experimental investigation of the metabolism of small organic molecules is particularly resource demanding; hence, computational methods are of considerable interest to complement experimental approaches. This review provides a broad overview of structure- and ligand-based computational methods for the prediction of xenobiotic metabolism. Current computational approaches to address xenobiotic metabolism are discussed from three major perspectives: (i) prediction of sites of metabolism (SOMs), (ii) elucidation of potential metabolites and their chemical structures, and (iii) prediction of direct and indirect effects of xenobiotics on metabolizing enzymes, where the focus is on the cytochrome P450 (CYP) superfamily of enzymes, the cardinal xenobiotics metabolizing enzymes. For each of these domains, a variety of approaches and their applications are systematically reviewed, including expert systems, data mining approaches, quantitative structure–activity relationships (QSARs), and machine learning-based methods, pharmacophore-based algorithms, shape-focused techniques, molecular interaction fields (MIFs), reactivity-focused techniques, protein–ligand docking, molecular dynamics (MD) simulations, and combinations of methods. Predictive metabolism is a developing area, and there is still enormous potential for improvement. However, it is clear that the combination of rapidly increasing amounts of available ligand- and structure-related experimental data (in particular, quantitative data) with novel and diverse simulation and modeling approaches is accelerating the development of effective tools for prediction of in vivo metabolism, which is reflected by the diverse and comprehensive data sources and methods for metabolism prediction reviewed here. This review attempts to survey the range and scope of computational methods applied to metabolism prediction and also to compare and contrast their applicability and performance.

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.

Insights into Chi recognition from the structure of an AddAB-type helicase-nuclease complex

Homologous recombination DNA repair requires double-strand break resection by helicase–nuclease enzymes. The crystal structure of bacterial AddAB in complex with DNA substrates shows that it employs an inactive helicase site to recognize ‘Chi' recombination hotspot sequences that regulate resection.

In bacterial cells, processing of double-stranded DNA breaks for repair by homologous recombination is dependent upon the recombination hotspot sequence Chi and is catalysed by either an AddAB- or RecBCD-type helicase–nuclease. Here, we report the crystal structure of AddAB bound to DNA. The structure allows identification of a putative Chi-recognition site in an inactivated helicase domain of the AddB subunit. By generating mutant protein complexes that do not respond to Chi, we show that residues responsible for Chi recognition are located in positions equivalent to the signature motifs of a conventional helicase. Comparison with the related RecBCD complex, which recognizes a different Chi sequence, provides further insight into the structural basis for sequence-specific ssDNA recognition. The structure suggests a simple mechanism for DNA break processing, explains how AddAB and RecBCD can accomplish the same overall reaction with different sets of functional modules and reveals details of the role of an Fe–S cluster in protein stability and DNA binding.

Molecular modeling of lipase binding to a substrate-water interface

Interactions of lipases with hydrophobic substrate-water interfaces are of great interest to design improved lipase variants and engineer reaction conditions. This chapter describes the necessary steps to carry out molecular dynamics simulations of Candida antarctica lipase B at tributyrin-water interface using the GROMACS simulation software. Special attention is drawn to the preparation of the protein and the substrate-water interface and to the analysis of the obtained trajectory.

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/.

MDpocket: open-source cavity detection and characterization on molecular dynamics trajectories

MOTIVATION

A variety of pocket detection algorithms are now freely or commercially available to the scientific community for the analysis of static protein structures. However, since proteins are dynamic entities, enhancing the capabilities of these programs for the straightforward detection and characterization of cavities taking into account protein conformational ensembles should be valuable for capturing the plasticity of pockets, and therefore allow gaining insight into structure-function relationships.

RESULTS

This article describes a new method, called MDpocket, providing a fast, free and open-source tool for tracking small molecule binding sites and gas migration pathways on molecular dynamics (MDs) trajectories or other conformational ensembles. MDpocket is based on the fpocket cavity detection algorithm and a valuable contribution to existing analysis tools. The capabilities of MDpocket are illustrated for three relevant cases: (i) the detection of transient subpockets using an ensemble of crystal structures of HSP90; (ii) the detection of known xenon binding sites and migration pathways in myoglobin; and (iii) the identification of suitable pockets for molecular docking in P38 Map kinase.

AVAILABILITY

MDpocket is free and open-source software and can be downloaded at http://fpocket.sourceforge.net.

CONTACT

pschmidtke@ub.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Structure and dynamics of the membrane-bound cytochrome P450 2C9

The microsomal, membrane-bound, human cytochrome P450 (CYP) 2C9 is a liver-specific monooxygenase essential for drug metabolism. CYPs require electron transfer from the membrane-bound CYP reductase (CPR) for catalysis. The structural details and functional relevance of the CYP-membrane interaction are not understood. From multiple coarse grained molecular simulations started with arbitrary configurations of protein-membrane complexes, we found two predominant orientations of CYP2C9 in the membrane, both consistent with experiments and conserved in atomic-resolution simulations. The dynamics of membrane-bound and soluble CYP2C9 revealed correlations between opening and closing of different tunnels from the enzyme's buried active site. The membrane facilitated the opening of a tunnel leading into it by stabilizing the open state of an internal aromatic gate. Other tunnels opened selectively in the simulations of product-bound CYP2C9. We propose that the membrane promotes binding of liposoluble substrates by stabilizing protein conformations with an open access tunnel and provide evidence for selective substrate access and product release routes in mammalian CYPs. The models derived here are suitable for extension to incorporate other CYPs for oligomerization studies or the CYP reductase for studies of the electron transfer mechanism, whereas the modeling procedure is generally applicable to study proteins anchored in the bilayer by a single transmembrane helix.

We describe the first atomic-detail models and simulations of a full-length, membrane-bound mammalian cytochrome P450. To date, all the structural studies of microsomal, drug-metabolizing cytochrome P450s have been performed using engineered, solubilized forms of the enzymes and it is not yet understood how the membrane influences their structure, dynamics, and ability to bind substrates. We focused on CYP2C9, the second most abundant cytochrome P450 in the human liver which oxidizes 20% of all marketed drugs. Here, we have derived models of CYP2C9-membrane complexes with a modeling procedure based on molecular dynamics simulations started with arbitrary configurations of the protein in the membrane and performed using both coarse grained and atomic-detail molecular representations. This procedure is expected to be generally applicable to proteins that are anchored in the membrane with a single transmembrane helix. The models and simulations provide evidence for selective substrate access and product release routes in membrane-bound CYPs. This observation may contribute to new strategies to manipulate the activity of CYPs and other enzymes with buried active sites. We anticipate that this study will bring about a paradigm shift towards studying microsomal CYPs as dynamic structures in their natural, lipid environment rather than in artificially solubilized forms.

Membrane position of ibuprofen agrees with suggested access path entrance to cytochrome P450 2C9 active site

Cytochrome P450 2C9 (CYP2C9) is a membrane-anchored human microsomal protein involved in the drug metabolism in liver. CYP2C9 consists of an N-terminal transmembrane anchor and a catalytic cytoplasmic domain. While the structure of the catalytic domain is well-known from X-ray experiments, the complete structure and its incorporation into the membrane remains unsolved. We constructed an atomistic model of complete CYP2C9 in a dioleoylphosphatidylcholine membrane and evolved it by molecular dynamics simulations in explicit water on a 100+ ns time-scale. The model agrees well with known experimental data about membrane positioning of cytochromes P450. The entry to the substrate access channel is proposed to be facing the membrane interior while the exit of the product egress channel is situated above the interface pointing toward the water phase. The positions of openings of the substrate access and product egress channels correspond to free energy minima of CYP2C9 substrate ibuprofen and its metabolite in the membrane, respectively.

The crystal structure of a voltage-gated sodium channel

Voltage-gated sodium (Na(V)) channels initiate electrical signalling in excitable cells and are the molecular targets for drugs and disease mutations, but the structural basis for their voltage-dependent activation, ion selectivity and drug block is unknown. Here we report the crystal structure of a voltage-gated Na(+) channel from Arcobacter butzleri (NavAb) captured in a closed-pore conformation with four activated voltage sensors at 2.7 Å resolution. The arginine gating charges make multiple hydrophilic interactions within the voltage sensor, including unanticipated hydrogen bonds to the protein backbone. Comparisons to previous open-pore potassium channel structures indicate that the voltage-sensor domains and the S4-S5 linkers dilate the central pore by pivoting together around a hinge at the base of the pore module. The NavAb selectivity filter is short, ∼4.6 Å wide, and water filled, with four acidic side chains surrounding the narrowest part of the ion conduction pathway. This unique structure presents a high-field-strength anionic coordination site, which confers Na(+) selectivity through partial dehydration via direct interaction with glutamate side chains. Fenestrations in the sides of the pore module are unexpectedly penetrated by fatty acyl chains that extend into the central cavity, and these portals are large enough for the entry of small, hydrophobic pore-blocking drugs. This structure provides the template for understanding electrical signalling in excitable cells and the actions of drugs used for pain, epilepsy and cardiac arrhythmia at the atomic level.

Encoding molecular motions in voxel maps

This paper builds on the combination of robotic path planning algorithms and molecular modeling methods for computing large-amplitude molecular motions, and introduces voxel maps as a computational tool to encode and to represent such motions. We investigate several applications and show results that illustrate the interest of such representation.

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.

Structural characterization of the complex between alpha-naphthoflavone and human cytochrome P450 1B1

The atomic structure of human P450 1B1 was determined by x-ray crystallography to 2.7 Å resolution with α-naphthoflavone (ANF) bound in the active site cavity. Although the amino acid sequences of human P450s 1B1 and 1A2 have diverged significantly, both enzymes exhibit narrow active site cavities, which underlie similarities in their substrate profiles. Helix I residues adopt a relatively flat conformation in both enzymes, and a characteristic distortion of helix F places Phe(231) in 1B1 and Phe(226) in 1A2 in similar positions for π-π stacking with ANF. ANF binds in a distinctly different orientation in P450 1B1 from that observed for 1A2. This reflects, in part, divergent conformations of the helix B'-C loop that are stabilized by different hydrogen-bonding interactions in the two enzymes. Additionally, differences between the two enzymes for other amino acids that line the edges of the cavity contribute to distinct orientations of ANF in the two active sites. Thus, the narrow cavity is conserved in both P450 subfamily 1A and P450 subfamily 1B with sequence divergence around the edges of the cavity that modify substrate and inhibitor binding. The conservation of these P450 1B1 active site amino acid residues across vertebrate species suggests that these structural features are conserved.

Plant Lessons: Exploring ABCB Functionality Through Structural Modeling

In contrast to mammalian ABCB1 proteins, narrow substrate specificity has been extensively documented for plant orthologs shown to catalyze the transport of the plant hormone, auxin. Using the crystal structures of the multidrug exporters Sav1866 and MmABCB1 as templates, we have developed structural models of plant ABCB proteins with a common architecture. Comparisons of these structures identified kingdom-specific candidate substrate-binding regions within the translocation chamber formed by the transmembrane domains of ABCBs from the model plant Arabidopsis. These results suggest an early evolutionary divergence of plant and mammalian ABCBs. Validation of these models becomes a priority for efforts to elucidate ABCB function and manipulate this class of transporters to enhance plant productivity and quality.

Dynamics and structural changes induced by ATP binding in SAV1866, a bacterial ABC exporter

Multidrug transporters of the ATP-binding cassette family export a wide variety of compounds across membranes in both prokaryotes and eukaryotes, using ATP hydrolysis as energy source. Several of these membrane proteins are of clinical importance. Although biochemical and structural studies have provided insights into the mechanism underlying substrate transport, many key questions subsist regarding the molecular and structural nature of this mechanism. In particular, the detailed conformational changes occurring during the catalytic cycle are still elusive. We explored the conformational changes occurring upon ATP/Mg(2+) binding using molecular dynamics simulations starting from the nucleotide-bound structure of SAV1866 embedded in an explicit lipid bilayer. The removal of nucleotide revealed a major rearrangement in the outer membrane leaflet portion of the transmembrane domain (TMD) resulting in the closure of the central cavity at the extracellular side. This closure is similar to that observed in the crystal nucleotide-free structures. The interface of the nucleotide-binding domain dimer (NDB) is significantly more hydrated in the nucleotide-free trajectory though it is not disrupted. This finding suggests that the TMD closure could occur as a first step preceding the dissociation of the dimer. The transmission pathway of the signal triggered by the removal of ATP/Mg(2+) mainly involves the conserved Q-loop and X-loop as well as TM6.

Molecular modeling on inhibitor complexes and active-site dynamics of cytochrome P450 C17, a target for prostate cancer therapy

A molecular model for the P450 enzyme cytochrome P450 C17 (CYP17) is presented based on sequence alignments of multiple template structures and homology modeling. This enzyme plays a central role in the biosynthesis of testosterone and is emerging as a major target in prostate cancer, with the recently developed inhibitor abiraterone currently in advanced clinical trials. The model is described in detail, together with its validation, by providing structural explanations to available site-directed mutagenesis data. The CYP17 molecule in this model is in the form of a triangular prism, with an edge of approximately 55 A and a thickness of approximately 37 A. It is predominantly helical, comprising 13 alpha helices interspersed by six 3(10) helices and 11 beta-sheets. Multinanosecond molecular dynamics simulations in explicit solvent have been carried out, and principal components analysis has been used to reveal the details of dynamics around the active site. Coarse-grained methods have also been used to verify low-frequency motions, which have been correlated with active-site gating. The work also describes the results of docking synthetic inhibitors, including the drug abiraterone and the natural substrate pregnenolone, in the CYP17 active site together with molecular dynamics simulations on the complexes.

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

As larger macromolecular structures become available, there is a growing need to understand their ‘internal’ volumes—such as deep clefts, channels and cavities—as these often play critical roles in their function. The 3V web server can automatically extract and comprehensively analyze all the internal volumes from input RNA and protein structures. It rapidly finds internal volumes by taking the difference between two rolling-probe solvent-excluded surfaces, one with as large as possible a probe radius and the other with a solvent radius (typically 1.5 Å for water). The outputs are volumetric representations, both as images and downloadable files, which can be used for further analysis.

The 3V server and source code are available from http://3vee.molmovdb.org.

Molecular mechanisms underlying the early stage of protein translocation through the Sec translocon

The Sec translocon, a protein-conducting channel, consists of a heterotrimeric complex (SecYEG in bacteria and Sec61 alpha beta gamma in eukaryotes) that provides a pathway for secretary proteins to cross membranes, or for membrane proteins to integrate into the membrane. The Sec translocon alone is a passive channel, and association with channel partners, including the ribosome or SecA ATPase in bacteria, is needed for protein translocation. Three recently published crystal structures of SecY are considered to represent the closed (resting state), pre-open (transitional state determined with the bound Fab fragment mimicking SecA interaction), and SecA-bound forms. To elucidate mechanisms of transition between closed and pre-open forms, we performed all-atom molecular dynamics simulations for the pre-open form of Thermus thermophilus SecYE and the closed form of Methanococcus janaschii SecYE beta in explicit solvent and membranes. We found that the closed form of SecY is stable, while the pre-open form without the Fab fragment undergoes large conformational changes toward the closed form. The pre-open form of SecY with Fab remains unchanged, suggesting that the cytosolic interaction mimicking SecA binding stabilizes the pre-open form of SecY. Importantly, a lipid molecule at the lateral gate region appears to be required to maintain the pre-open form in the membrane. We propose that the conformational transition from closed to pre-open states of SecY upon association with SecA facilitates intercalation of phospholipids at the lateral gate, inducing initial entry of the positively charged signal peptide into the channel.

Crystal structure of the bifunctional proline utilization A flavoenzyme from Bradyrhizobium japonicum

The bifunctional proline catabolic flavoenzyme, proline utilization A (PutA), catalyzes the oxidation of proline to glutamate via the sequential activities of FAD-dependent proline dehydrogenase (PRODH) and NAD(+)-dependent Delta(1)-pyrroline-5-carboxylate dehydrogenase (P5CDH) domains. Although structures for some of the domains of PutA are known, a structure for the full-length protein has not previously been solved. Here we report the 2.1 A resolution crystal structure of PutA from Bradyrhizobium japonicum, along with data from small-angle x-ray scattering, analytical ultracentrifugation, and steady-state and rapid-reaction kinetics. PutA forms a ring-shaped tetramer in solution having a diameter of 150 A. Within each protomer, the PRODH and P5CDH active sites face each other at a distance of 41 A and are connected by a large, irregularly shaped cavity. Kinetics measurements show that glutamate production occurs without a lag phase, suggesting that the intermediate, Delta(1)-pyrroline-5-carboxylate, is preferably transferred to the P5CDH domain rather than released into the bulk medium. The structural and kinetic data imply that the cavity serves both as a microscopic vessel for the hydrolysis of Delta(1)-pyrroline-5-carboxylate to glutamate semialdehyde and a protected conduit for the transport of glutamate semialdehyde to the P5CDH active site.

Structural basis for cofactor-independent dioxygenation of N-heteroaromatic compounds at the alpha/beta-hydrolase fold

Enzymatic catalysis of oxygenation reactions in the absence of metal or organic cofactors is a considerable biochemical challenge. The CO-forming 1-H-3-hydroxy-4-oxoquinaldine 2,4-dioxygenase (HOD) from Arthrobacter nitroguajacolicus Rü61a and 1-H-3-hydroxy-4-oxoquinoline 2,4-dioxygenase (QDO) from Pseudomonas putida 33/1 are homologous cofactor-independent dioxygenases involved in the breakdown of N-heteroaromatic compounds. To date, they are the only dioxygenases suggested to belong to the alpha/beta-hydrolase fold superfamily. Members of this family typically catalyze hydrolytic processes rather than oxygenation reactions. We present here the crystal structures of both HOD and QDO in their native state as well as the structure of HOD in complex with its natural 1-H-3-hydroxy-4-oxoquinaldine substrate, its N-acetylanthranilate reaction product, and chloride as dioxygen mimic. HOD and QDO are structurally very similar. They possess a classical alpha/beta-hydrolase fold core domain additionally equipped with a cap domain. Organic substrates bind in a preorganized active site with an orientation ideally suited for selective deprotonation of their hydroxyl group by a His/Asp charge-relay system affording the generation of electron-donating species. The "oxyanion hole" of the alpha/beta-hydrolase fold, typically employed to stabilize the tetrahedral intermediate in ester hydrolysis reactions, is utilized here to host and control oxygen chemistry, which is proposed to involve a peroxide anion intermediate. Product release by proton back transfer from the catalytic histidine is driven by minimization of intramolecular charge repulsion. Structural and kinetic data suggest a nonnucleophilic general-base mechanism. Our analysis provides a framework to explain cofactor-independent dioxygenation within a protein architecture generally employed to catalyze hydrolytic reactions.

Visualizing the mapped ion pathway through the Na,K-ATPase pump

The Na(+),K(+)-ATPase pump achieves thermodynamically uphill exchange of cytoplasmic Na(+) ions for extracellular K(+) ions by using ATP-mediated phosphorylation, followed by autodephosphorylation, to power conformational changes that allow ion access to the pump's binding sites from only one side of the membrane at a time. Formally, the pump behaves like an ion channel with two tightly coupled gates that are constrained to open and close alternately. The marine agent palytoxin disrupts this coupling, allowing both gates to sometimes be open, so temporarily transforming a pump into an ion channel. We made a cysteine scan of Na(+),K(+)-ATPase transmembrane (TM) segments TM1 to TM6, and used recordings of Na(+) current flow through palytoxin-bound pump-channels to monitor accessibility of introduced cysteine residues via their reaction with hydrophilic methanethiosulfonate (MTS) reagents. To visualize the open-channel pathway, the reactive positions were mapped onto a homology model of Na(+),K(+)-ATPase based on the structure of the related sarcoplasmicand endoplasmic-reticulum (SERCA) Ca(2+)-ATPase in a BeF(3)(-)-trapped state,(1,2) in which the extra-cytoplasmic gate is wide open (although the cytoplasmic access pathway is firmly shut). The results revealed a single unbroken chain of reactive positions that traverses the pump from the extracellular surface to the cytoplasm, comprises residues from TM1, TM2, TM4 and TM6, and passes through the equivalent of cation binding site II in SERCA, but not through site I. Cavity search analysis of the homology model validated its use for mapping the data by yielding a calculated extra-cytoplasmic pathway surrounded by MTS-reactive residues. As predicted by previous experimental results, that calculated extra-cytoplasmic pathway abruptly broadens above residue T806, at the outermost end of TM6 that forms the floor of the extracellular-facing vestibule. These findings provide a structural basis for further understanding cation translocation by the Na(+),K(+)-ATPase and by other P-type pumps like the Ca(2+)- and H(+),K(+)-ATPases.

Conformational flexibility in immunoglobulin E-Fc 3-4 revealed in multiple crystal forms

The structure of immunoglobulin E (IgE)-Fc(3-4) has been solved in three new crystal forms, providing 13 snapshots of the Fc conformation and revealing a diverse range of open-closed motions among subunit chains and dimers. A more detailed analysis of the open-to-closed motion of IgE-Fc(3-4) was possible with so many structures, and the new structures allow a more thorough examination of the flexibility of IgE-Fc and its implications for receptor binding. The existence of a hydrophobic pocket at the elbow region of the Fc appears to be conformation dependent and suggests a means of regulating the IgE-Fc conformation (and potentially receptor binding) with small molecules.