Identification of tunnels in proteins, nucleic acids, inorganic materials and molecular ensembles

Modifying the amino acids in conformational motion pathway of the α-amylase of Geobacillus stearothermophilus improved its activity and stability

Amino acids along the conformational motion pathway of the enzyme molecule correlated to its flexibility and rigidity. To enhance the enzyme activity and thermal stability, the motion pathway of Geobacillus stearothermophilus α-amylase has been identified and molecularly modified by using the neural relational inference model and deep learning tool. The significant differences in substrate specificity, enzymatic kinetics, optimal temperature, and thermal stability were observed among the mutants with modified amino acids along the pathway. Mutants especially the P44E demonstrated enhanced hydrolytic activity and catalytic efficiency (kcat/KM) than the wild-type enzyme to 95.0% and 93.8% respectively, with the optimum temperature increased to 90°C. This mutation from proline to glutamic acid has increased the number and the radius of the bottleneck of the channels, which might facilitate transporting large starch substrates into the enzyme. The mutation could also optimize the hydrogen bonding network of the catalytic center, and diminish the spatial hindering to the substrate entry and exit from the catalytic center.

Lobaric acid prevents the adverse effects of tetramethrin on the estrous cycle of female albino Wistar rats

Tetramethrin (Tm) is a commonly used pesticide that has been reported to exert estrogen-antagonistic effects selectively on female rats. The present study was undertaken to assess the protective role of lobaric acid (La) on estrous cycle in Tm-treated female Wistar rats. Female rats were exposed to Tm (50 mg/kg b.w/day) only or in combination with La at low (50 mg/kg b.w/day) or high (100 mg/kg b.w/day) dose for 30 days. The results showed that Tm altered the estrous cycle of female rats by decreasing the levels of luteinizing hormone, follicular-stimulating hormone, progesterone, estrone, and estradiol while increasing testosterone level. The morphology of vaginal smears of Tm-treated female rats showed the presence of abnormal cells and/or structures at different phases of estrus cycle. Strikingly, in (Tm + La)-treated rats, all the observed adverse effects of Tm on the hormonal parameters, cell morphology, and the length of each phase of estrous cycle were significantly diminished in a dose-dependent manner. The docking results showed that La competes with Tm for Gonadotropin-Releasing Hormone (GnRH) receptor, thereby reducing the toxicity of Tm but did not cancel the response of GnRH receptor completely. In conclusion, our results designated that La could be used as a potential candidate in the management of insecticide-induced alterations of the reproductive cycle of rodents.

Potential of mean force and umbrella sampling simulation for the transport of 5-oxazolidinone in heterotetrameric sarcosine oxidase

The structure of heterotetrameric sarcosine oxidase (HSO) contains a highly complex system composed of a large cavity and tunnels, which are essential for the reaction and migration of the reactants, products, and intermediates. Previous geometrical analysis using the CAVER program has predicted that there are three possible tunnels, T1, T2, and T3, for the exit pathway of the iminium intermediate, 5-oxazolidinone (5-OXA), of the enzyme reaction. Previous molecular dynamics (MD) simulation of HSO has identified the regions containing the water channels from the density distribution of water. The simulation indicated that tunnel T3 is the most probable exit pathway of 5-OXA. In the present study, the potential of mean force (PMF) for the transport of 5-OXA through tunnels T1, T2, and T3 was calculated using umbrella sampling (US) MD simulations and the weighted histogram analysis method. The PMF profiles for the three tunnels support the notion that tunnel T3 is the exit pathway of 5-OXA, and that 5-OXA tends to stay at the middle of the tunnel. The maximum errors of the calculated PMF for the predicted exit pathway, tunnel T3, were estimated by repeating the US simulations using different sets of initial positions. The PMF profile was also calculated for the transport of glycine within T3. The PMF profiles from the US simulations were in good agreement with the previous predictions that 5-OXA escape through tunnel T3 and how glycine is released to the outside of HSO was discussed.

Estimation of the ligand-binding free energy of checkpoint kinase 1 via non-equilibrium MD simulations

Checkpoint kinase 1 (CHK1) is a serine/threonine-protein kinase that is involved in cell cycle regulation in eukaryotes. Inhibition of CHK1 is thus considered as a promising approach in cancer therapy. In this study, the fast pulling of ligand (FPL) process was applied to predict the relative binding affinities of CHK1 inhibitors using non-equilibrium molecular dynamics (MD) simulations. The work of external harmonic forces to pull the ligand out of the binding cavity strongly correlated with the experimental binding affinity of CHK1 inhibitors with the correlation coefficient of R = -0.88 and an overall root mean square error (RMSE) of 0.99 kcal/mol. The data indicate that the FPL method is highly accurate in predicting the relative binding free energies of CHK1 inhibitors with an affordable CPU time. A new set of molecules were designed based on the molecular modeling of interactions between the known inhibitor and CHK1 as inhibitory candidates. Molecular docking and FPL results exhibited that the binding affinities of developed ligands were similar to the known inhibitor in interaction with the catalytic site of CHK1, producing very potential CHK1 inhibitors of that the inhibitory activities should be further evaluated in vitro.

Enzyme discovery beyond homology: a unique hydroxynitrile lyase in the Bet v1 superfamily

Homology and similarity based approaches are most widely used for the identification of new enzymes for biocatalysis. However, they are not suitable to find truly novel scaffolds with a desired function and this averts options and diversity. Hydroxynitrile lyases (HNLs) are an example of non-homologous isofunctional enzymes for the synthesis of chiral cyanohydrins. Due to their convergent evolution, finding new representatives is challenging. Here we show the discovery of unique HNL enzymes from the fern Davallia tyermannii by coalescence of transcriptomics, proteomics and enzymatic screening. It is the first protein with a Bet v1-like protein fold exhibiting HNL activity, and has a new catalytic center, as shown by protein crystallography. Biochemical properties of D. tyermannii HNLs open perspectives for the development of a complementary class of biocatalysts for the stereoselective synthesis of cyanohydrins. This work shows that systematic integration of -omics data facilitates discovery of enzymes with unpredictable sequences and helps to extend our knowledge about enzyme diversity.

Ligand-based virtual screening interface between PyMOL and LiSiCA

Ligand-based virtual screening of large small-molecule databases is an important step in the early stages of drug development. It is based on the similarity principle and is used to reduce the chemical space of large databases to a manageable size where chosen ligands can be experimentally tested. Ligand-based virtual screening can also be used to identify bioactive molecules with different basic scaffolds compared to already known bioactive molecules, thus having the potential to increase the structural variability of compounds. Here, we present an interface between the popular molecular graphics system PyMOL and the ligand-based virtual screening software LiSiCA available at http://insilab.org/lisica-plugin and demonstrate how this interface can be used in the early stages of drug discovery process.Graphical Abstract

CAVER: Algorithms for Analyzing Dynamics of Tunnels in Macromolecules

The biological function of a macromolecule often requires that a small molecule or ion is transported through its structure. The transport pathway often leads through void spaces in the structure. The properties of transport pathways change significantly in time; therefore, the analysis of a trajectory from molecular dynamics rather than of a single static structure is needed for understanding the function of pathways. The identification and analysis of transport pathways are challenging because of the high complexity and diversity of macromolecular shapes, the thermal motion of their atoms, and the large amount of conformations needed to properly describe conformational space of protein structure. In this paper, we describe the principles of the CAVER 3.0 algorithms for the identification and analysis of properties of transport pathways both in static and dynamic structures. Moreover, we introduce the improved clustering solution for finding tunnels in macromolecules, which is included in the latest CAVER 3.02 version. Voronoi diagrams are used to identify potential pathways in each snapshot of a molecular dynamics trajectory and clustering is then used to find the correspondence between tunnels from different snapshots. Furthermore, the geometrical properties of pathways and their evolution in time are computed and visualized.

Analysis of water channels by molecular dynamics simulation of heterotetrameric sarcosine oxidase

A precise 100-ns molecular dynamics simulation in aquo was performed for the heterotetrameric sarcosine oxidase bound with a substrate analogue, dimethylglycine. The spatial region including the protein was divided into small rectangular cells. The average number of the water molecules locating within each cell was calculated based on the simulation trajectory. The clusters of the cells filled with water molecules were used to determine the water channels. The narrowness of the channels, the average hydropathy indices of the residues of the channels, and the number of migration events of water molecules through the channels were consistent with the selective transport hypothesis whereby tunnel T3 is the pathway for the exit of the iminium intermediate of the enzyme reaction.

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.

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

A single mutation in a tunnel to the active site changes the mechanism and kinetics of product release in haloalkane dehalogenase LinB

Many enzymes have buried active sites. The properties of the tunnels connecting the active site with bulk solvent affect ligand binding and unbinding and also the catalytic properties. Here, we investigate ligand passage in the haloalkane dehalogenase enzyme LinB and the effect of replacing leucine by a bulky tryptophan at a tunnel-lining position. Transient kinetic experiments show that the mutation significantly slows down the rate of product release. Moreover, the mechanism of bromide ion release is changed from a one-step process in the wild type enzyme to a two-step process in the mutant. The rate constant of bromide ion release corresponds to the overall steady-state turnover rate constant, suggesting that product release became the rate-limiting step of catalysis in the mutant. We explain the experimental findings by investigating the molecular details of the process computationally. Analysis of trajectories from molecular dynamics simulations with a tunnel detection software reveals differences in the tunnels available for ligand egress. Corresponding differences are seen in simulations of product egress using a specialized enhanced sampling technique. The differences in the free energy barriers for egress of a bromide ion obtained using potential of mean force calculations are in good agreement with the differences in rates obtained from the transient kinetic experiments. Interactions of the bromide ion with the introduced tryptophan are shown to affect the free energy barrier for its passage. The study demonstrates how the mechanism of an enzymatic catalytic cycle and reaction kinetics can be engineered by modification of protein tunnels.

Different functions of MdtB and MdtC subunits in the heterotrimeric efflux transporter MdtB(2)C complex of Escherichia coli

In contrast to homotrimeric transporters of the RND (resistance-nodulation-division) superfamily, which often conduct efflux transport of a wide range of substrates by the functionally rotating mechanism, the mechanism utilized by the heterotrimeric members of this family, which also perform multidrug efflux, is unclear. We examined one heterotrimeric transporter, the MdtB(2)C complex of Escherichia coli, by an extensive cysteine scanning mutagenesis of residues likely involved in ligand transport. Many such mutations in MdtC strongly decreased the level of cloxacillin transport, whereas mutations of corresponding residues in MdtB did not affect transport. Furthermore, many such residues in MdtC were covalently modified by fluorescein maleimide, which acted as a substrate and presumably produced labeling of the residues in the substrate path. In contrast, few residues in MdtB were labeled. Together with the previous data showing that the inactivation of proton translocation channel in MdtC has an only modest effect on transport yet in MdtB totally inactivated the activity, these results suggest that the two subunits, MdtB and MdtC, play very different roles, MdtC likely involved in substrate binding and transport and MdtB presumably inducing the conformational change needed for transport through proton translocation. Three-dimensional models of MdtB and MdtC, based on sequence homology with the AcrB transporter, also support this interpretation.

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.

Structural and kinetic isotope effect studies of nicotinamidase (Pnc1) from Saccharomyces cerevisiae

Nicotinamidases catalyze the hydrolysis of nicotinamide to nicotinic acid and ammonia. Nicotinamidases are absent in mammals but function in NAD(+) salvage in many bacteria, yeast, plants, protozoa, and metazoans. We have performed structural and kinetic investigations of the nicotinamidase from Saccharomyces cerevisiae (Pnc1). Steady-state product inhibitor analysis revealed an irreversible reaction in which ammonia is the first product released, followed by nicotinic acid. A series of nicotinamide analogues acting as inhibitors or substrates were examined, revealing that the nicotinamide carbonyl oxygen and ring nitrogen are critical for binding and reactivity. X-ray structural analysis revealed a covalent adduct between nicotinaldehyde and Cys167 of Pnc1 and coordination of the nicotinamide ring nitrogen to the active-site zinc ion. Using this structure as a guide, the function of several residues was probed via mutagenesis and primary (15)N and (13)C kinetic isotope effects (KIEs) on V/K for amide bond hydrolysis. The KIE values of almost all variants were increased, indicating that C-N bond cleavage is at least partially rate limiting; however, a decreased KIE for D51N was indicative of a stronger commitment to catalysis. In addition, KIE values using slower alternate substrates indicated that C-N bond cleavage is at least partially rate limiting with nicotinamide to highly rate limiting with thionicotinamide. A detailed mechanism involving nucleophilic attack of Cys167, followed by elimination of ammonia and then hydrolysis to liberate nicotinic acid, is discussed. These results will aid in the design of mechanism-based inhibitors to target pathogens that rely on nicotinamidase activity.

Structural differences between soluble and membrane bound cytochrome P450s

The superfamily of cytochrome P450s forms a large class of heme monooxygenases with more than 13,000 enzymes represented in organisms from all biological kingdoms. Despite impressive variability in sizes, sequences, location, and function, all cytochrome P450s from various organisms have very similar tertiary structures within the same fold. Here we show that systematic comparison of all available X-ray structures of cytochrome P450s reveals the presence of two distinct structural classes of cytochrome P450s. For all membrane bound enzymes, except the CYP51 family, the beta-domain and the A-propionate heme side chain are shifted towards the proximal side of the heme plane, which may result in an increase of the volume of the substrate binding pocket and an opening of a potential channel for the substrate access and/or product escape directly into the membrane. This structural feature is also observed in several soluble cytochrome P450s, such as CYP108, CYP151, and CYP158A2, which catalyze transformations of bulky substrates. Alternatively, both beta-domains and the A-propionate side chains in the soluble isozymes extend towards the distal site of the heme. This difference between the structures of soluble and membrane bound cytochrome P450s can be rationalized through the presence of several amino acid inserts in the latter class which are involved in direct interactions with the membrane, namely the F'- and G'-helices. Molecular dynamics using the most abundant human cytochrome P450, CYP3A4, incorporated into a model POPC bilayer reveals the facile conservation of a substrate access channel, directed into the membrane between the B-C loop and the beta domain, and the closure of the peripheral substrate access channel directed through the B-C loop. This is in contrast to the case when the same simulation is run in buffer, where no such channel closing occurs. Taken together, these results reveal a key structural difference between membrane bound and soluble cytochrome P450s with important functional implications induced by the lipid bilayer.

Interactions of cations with RNA loop-loop complexes

RNA loop-loop interactions are essential in many biological processes, including initiation of RNA folding into complex tertiary shapes, promotion of dimerization, and viral replication. In this article, we examine interactions of metal ions with five RNA loop-loop complexes of unique biological significance using explicit-solvent molecular-dynamics simulations. These simulations revealed the presence of solvent-accessible tunnels through the major groove of loop-loop interactions that attract and retain cations. Ion dynamics inside these loop-loop complexes were distinctly different from the dynamics of the counterion cloud surrounding RNA and depend on the number of basepairs between loops, purine sequence symmetry, and presence of unpaired nucleotides. The cationic uptake by kissing loops depends on the number of basepairs between loops. It is interesting that loop-loop complexes with similar functionality showed similarities in cation dynamics despite differences in sequence and loop size.

Why boys will be boys: two pathways of fetal testicular androgen biosynthesis are needed for male sexual differentiation

Human sexual determination is initiated by a cascade of genes that lead to the development of the fetal gonad. Whereas development of the female external genitalia does not require fetal ovarian hormones, male genital development requires the action of testicular testosterone and its more potent derivative dihydrotestosterone (DHT). The "classic" biosynthetic pathway from cholesterol to testosterone in the testis and the subsequent conversion of testosterone to DHT in genital skin is well established. Recently, an alternative pathway leading to DHT has been described in marsupials, but its potential importance to human development is unclear. AKR1C2 is an enzyme that participates in the alternative but not the classic pathway. Using a candidate gene approach, we identified AKR1C2 mutations with sex-limited recessive inheritance in four 46,XY individuals with disordered sexual development (DSD). Analysis of the inheritance of microsatellite markers excluded other candidate loci. Affected individuals had moderate to severe undervirilization at birth; when recreated by site-directed mutagenesis and expressed in bacteria, the mutant AKR1C2 had diminished but not absent catalytic activities. The 46,XY DSD individuals also carry a mutation causing aberrant splicing in AKR1C4, which encodes an enzyme with similar activity. This suggests a mode of inheritance where the severity of the developmental defect depends on the number of mutations in the two genes. An unrelated 46,XY DSD patient carried AKR1C2 mutations on both alleles, confirming the essential role of AKR1C2 and corroborating the hypothesis that both the classic and alternative pathways of testicular androgen biosynthesis are needed for normal human male sexual differentiation.

Role of the omega loop in specificity determination in subsite 2 of the D-alanine:D-alanine (D-lactate) ligase from Leuconostoc mesenteroides: a molecular docking study

The synthesis of D-ala-D-lactate in Leuconostoc mesenteroides is catalyzed by D-alanine:D-alanine (D-lactate) ligase (ADP). The ability to assemble this depsipeptide as well as D-ala-D-ala provides a mechanism for the organism's intrinsic resistance to vancomycin. Mutation of Phe261 to Tyr261 in the Ω-loop of this ligase showed a complete loss of the ability to make D-ala-D-lactate (Park and Walsh, J. Biol. Chem. 272 (1997) 9210-9214). Phe261 is a key specificity determinant in the α-helical cap of the Ω-loop when folded into the closed conformation. A molecular docking study of the closed ligase using AutoDock 4.2 defines additional specificity constraints promoted by the Ω-loop capping the catalytic center. Attaining productive orientations of D-lactate with favorable ligation chemistry requires the flexibilities of Phe261 and Arg301 in the docking protocol. These are in addition to the optimization of van der Waals contacts with Lys260, Met326, and Ser327. The location of Phe261 and Lys260 in the α-helical cap of the Ω-loop over subsite 2 is an essential part of the folding process ensuring depsipeptide formation in the hydrophobic environment of the catalytic center. The importance of the F261Y mutation suggests that the hydroxyl of Tyr261 plays an instrumental role in determining non-productive docking orientations of D-lactate. Two of these are presented: (A) D-lactate-OH as an H-bond donor to the Tyr261-OH; (B) D-lactate as an H-bond donor to the phosphoryl of the intermediate D-alanyl phosphate, and the D-lactate-COO- as an H-bond acceptor for the Tyr261-OH. Neither orientation, A or B, show the bifurcated H-bonding with Arg301 recently proposed for the activation of the nucleophilic D-lactate for D-ala-D-lactate formation. Insights into the role of the Ω-loop and its K(F/Y) signature provide additional background for inhibitor design targeted to subsite 2 of the D-alanine:D-alanine (D-X) ligases.

Channeling of electrons within SLAC, the small laccase from Streptomyces coelicolor

The reduction kinetics of the fluorescently labeled small laccase (SLAC) from Streptomyces coelicolor was studied by stopped flow kinetic measurements. The tryptophan fluorescence and the emission from a covalently attached label were used to selectively follow the progress of the reduction of the trinuclear copper center (TNC) and the type-1 (T1) Cu site in the enzyme as a function of time. A numerical analysis of the kinetic traces provided new insight into the midpoint potential difference between the T1 and the TNC site as the TNC becomes stepwise charged with electrons. The change in fluorescence of the TNC as the reduction of the TNC proceeds provided evidence that the type-3 dinuclear part of the TNC becomes charged prior to the reduction of the type-2 (T2) center of the TNC. The rate of reduction of the enzyme by dithionite (DT) appeared proportional to the square root of the DT concentration with a rate constant of k(red) = 0.28 +/- 0.02 microM(-1/2) s(-1).

Rational engineering of L-asparaginase reveals importance of dual activity for cancer cell toxicity

Using proteins in a therapeutic context often requires engineering to modify functionality and enhance efficacy. We have previously reported that the therapeutic antileukemic protein macromolecule Escherichia coli L-asparaginase is degraded by leukemic lysosomal cysteine proteases. In the present study, we successfully engineered L-asparaginase to resist proteolytic cleavage and at the same time improve activity. We employed a novel combination of mutant sampling using a genetic algorithm in tandem with flexibility studies using molecular dynamics to investigate the impact of lid-loop and mutations on drug activity. Applying these methods, we successfully predicted the more active L-asparaginase mutants N24T and N24A. For the latter, a unique hydrogen bond network contributes to higher activity. Furthermore, interface mutations controlling secondary glutaminase activity demonstrated the importance of this enzymatic activity for drug cytotoxicity. All selected mutants were expressed, purified, and tested for activity and for their ability to form the active tetrameric form. By introducing the N24A and N24A R195S mutations to the drug L-asparaginase, we are a step closer to individualized drug design.

Molecular basis of reduced glucosylceramidase activity in the most common Gaucher disease mutant, N370S

Gaucher disease is caused by the defective activity of the lysosomal hydrolase, glucosylceramidase. Although the x-ray structure of wild type glucosylceramidase has been resolved, little is known about the structural features of any of the >200 mutations. Various treatments for Gaucher disease are available, including enzyme replacement and chaperone therapies. The latter involves binding of competitive inhibitors at the active site to enable correct folding and transport of the mutant enzyme to the lysosome. We now use molecular dynamics, a set of structural analysis tools, and several statistical methods to determine the flexible behavior of the N370S Gaucher mutant at various pH values, with and without binding the chaperone, N-butyl-deoxynojirimycin. We focus on the effect of the chaperone on the whole protein, on the active site, and on three important structural loops, and we demonstrate how the chaperone modifies the behavior of N370S in such a way that it becomes more active at lysosomal pH. Our results suggest a mechanism whereby the binding of N-butyl-deoxynojirimycin helps target correctly folded glucosylceramidase to the lysosome, contributes to binding with saposin C, and explains the initiation of the substrate-enzyme complex. Such analysis provides a new framework for determination of the structure of other Gaucher disease mutants and suggests new approaches for rational drug design.

Evidence for a bicarbonate "escort" site in Haemophilus influenzae beta-carbonic anhydrase

The Haemophilus influenzae beta-carbonic anhydrase (HICA) allosteric site variants V47A and G41A were overexpressed and purified to homogeneity. These variants have k(cat)/K(m) values similar to that of the wild-type enzyme and exhibit a similar dramatic decrease in catalytic activity at pH <8.0. However, both HICA-G41A and -V47A were serendipitously found to bind sulfate ion or bicarbonate ion near pairs of Glu50 and Arg64 residues located on the dimerization interface. In the case of HICA-V47A, bicarbonate ions simultaneously bind to both the dimerization interface and the allosteric sites. For HICA-G41A, two of 12 chains in the asymmetric unit bind bicarbonate ion exclusively at the dimerization interface, while the remaining 10 chains bind bicarbonate ion exclusively at the allosteric site. We propose that the new anion binding site along the dimerization interface of HICA is an "escort" site that represents an intermediate along the ingress and egress route of bicarbonate ion to and from the allosteric binding site, respectively. The structural evidence for sulfate binding at the escort site suggests that the mechanism of sulfate activation of HICA is the result of sulfate ion competing for bicarbonate at the escort site, preventing passage of bicarbonate from the bulk solution to its allosteric site.

Ligand docking and binding site analysis with PyMOL and Autodock/Vina

Docking of small molecule compounds into the binding site of a receptor and estimating the binding affinity of the complex is an important part of the structure-based drug design process. For a thorough understanding of the structural principles that determine the strength of a protein/ligand complex both, an accurate and fast docking protocol and the ability to visualize binding geometries and interactions are mandatory. Here we present an interface between the popular molecular graphics system PyMOL and the molecular docking suites Autodock and Vina and demonstrate how the combination of docking and visualization can aid structure-based drug design efforts.

Mechanism of recognition of compounds of diverse structures by the multidrug efflux pump AcrB of Escherichia coli

The AcrB trimeric multidrug efflux transporter of Escherichia coli pumps out a very wide spectrum of compounds. Although minocycline and doxorubicin have been cocrystallized within the large binding pocket in the periplasmic domain of the binding protomer, nothing is known about the binding of many other ligands to this protein. We used computer docking to evaluate the interaction of about 30 compounds with the binding protomer and found that many of them are predicted to bind to a narrow groove at one end of the pocket whereas some others prefer to bind to a wide cave at the other end. Competition assays using nitrocefin efflux and covalent labeling of Phe615Cys mutant AcrB with fluorescein-5-maleimide showed that presumed groove-binders competed against each other, but cave-binders did not compete against groove-binders, although the number of compounds tested was limited. These results give us at least a hypothesis to be tested by more biochemical and genetic experiments in the future.

Finding molecular dioxygen tunnels in homoprotocatechuate 2,3-dioxygenase: implications for different reactivity of identical subunits

Extradiol dioxygenases facilitate microbial aerobic degradation of catechol and its derivatives by activating molecular dioxygen and incorporating both oxygen atoms into their substrates. Experimental and theoretical studies have focused on the mechanism of the reaction at the active site. However, whether the catalytic rate is limited by O(2) access to the active site has not yet been explored. Here, we choose a recently solved X-ray structure of homoprotocatechuate 2,3-dioxygenase as a typical example to determine potential pathways for O(2) migration from the solvent into the enzyme center. On the basis of the trajectories of two 10-ns molecular dynamics simulations, implicit ligand sampling was used to calculate the 3D free energy map for O(2) inside the protein. The energetically optimal routes for O(2) diffusion were identified for each subunit of the homotetrameric protein structure. The O(2) tunnels formed because of thermal fluctuations were also characterized by connecting elongated cavities inside the protein. By superimposing the favorable O(2) tunnels on to the free energy map, both energetically and geometrically preferred O(2) pathways were determined, as also were the amino acids that may be critical for O(2) passage along these paths. Our results demonstrate that identical subunits possess quite distinct O(2) tunnels. The order of O(2) affinity of these tunnels is generally consistent with the order of the catalytic rate of each subunit. As a consequence, the probability of finding the reaction product is highest in the subunit containing the highest O(2) affinity pathway.

The structure of neuroglobin at high Xe and Kr pressure reveals partial conservation of globin internal cavities

Neuroglobin (Ngb) is a hexacoordinate globin expressed in the brain of vertebrates. Ferrous Ngb binds dioxygen with high affinity and the O(2) adduct is able to scavenge NO. Convincing in vitro and in vivo data indicate that Ngb is involved in neuroprotection during hypoxia and ischemia. The 3D structure of Ngb reveals the presence of a wide internal cavity connecting its heme active site with the bulk. To explore the role of this "tunnel" in the control of ligand binding, we determined the structure of metNgb and NgbCO equilibrated with Xe or Kr. We show four docking sites for Xe (only two for Kr); two of the four Xe sites are within the large cavity. They are only partially conserved in globins, since the two proximal Xe sites identified in myoglobin (Xe1 and Xe2) are absent in Ngb, as well as in cytoglobin. The Xe docking sites in Ngb map a pathway within the protein matrix, leading to the heme, which becomes more accessible in the ligand-bound species. This may be of significance in connection with the redox chemistry that may be the primary function of this hexacoordinate globin.

A substrate channel in the nitrogenase MoFe protein

Nitrogenase catalyzes the six electron/six proton reduction of N2 to two ammonia molecules at a complex organometallocluster called "FeMo cofactor." This cofactor is buried within the alpha-subunit of the MoFe protein, with no obvious access for substrates. Examination of high-resolution X-ray crystal structures of MoFe proteins from several organisms has revealed the existence of a water-filled channel that extends from the solvent-exposed surface to a specific face of FeMo cofactor. This channel could provide a pathway for substrate and product access to the active site. In the present work, we examine this possibility by substituting four different amino acids that line the channel with other residues and analyze the impact of these substitutions on substrate reduction kinetic parameters. Each of the MoFe protein variants was purified and kinetic parameters were established for the reduction of the substrates N2, acetylene, azide, and propyne. For each MoFe protein, V (max) values for the different substrates were found to be nearly unchanged when compared with the values for the wild-type MoFe protein, indicating that electron delivery to the active site is not compromised by the various substitutions. In contrast, the K(m) values for these substrates were found to increase significantly (up to 22-fold) in some of the MoFe protein variants compared with the wild-type MoFe protein values. Given that each of the amino acids that were substituted is remote from the active site, these results are consistent with the water-filled channel functioning as a substrate channel in the MoFe protein.

Computational tools for designing and engineering biocatalysts

Current computational tools to assist experimentalists for the design and engineering of proteins with desired catalytic properties are reviewed. The applications of these tools for de novo design of protein active sites, optimization of substrate access and product exit pathways, redesign of protein-protein interfaces, identification of neutral/advantageous/deleterious mutations in the libraries from directed evolution and stabilization of protein structures are described. Remarkable progress is seen in de novo design of enzymes catalyzing a chemical reaction for which a natural biocatalyst does not exist. Yet, constructed biocatalysts do not match natural enzymes in their efficiency, suggesting that more research is needed to capture all the important features of natural biocatalysts in theoretical designs.

Finding and characterizing tunnels in macromolecules with application to ion channels and pores

We describe a new algorithm, CHUNNEL, to automatically find, characterize, and display tunnels or pores in proteins. The correctness and accuracy of the algorithm is verified on a constructed set of proteins and used to analyze large sets of real proteins. The verification set contains proteins with artificially created pores of known path and width profile. The previous benchmark algorithm, HOLE, is compared with the new algorithm. Results show that the major advantage of the new algorithm is that it can successfully find and characterize tunnels with no a priori guidance or clues about the location of the tunnel mouth, and it will successfully find multiple tunnels if present. CHUNNEL can also be used in conjunction with HOLE, with the former used toprime HOLE and the latter to track and characterize the pores. Analysis was conducted on families of membrane protein structures culled from the Protein Data Bank as well as on a set of transmembrane proteins with predicted membrane-aqueous phase interfaces, yielding the first completely automated examination of tunnels through membrane proteins, including tunnels that exit in the membrane bilayer.