Pathways and Mechanisms for Product Release in the Engineered Haloalkane Dehalogenases Explored Using Classical and Random Acceleration Molecular Dynamics Simulations

Eight mutants of the DhaA haloalkane dehalogenase carrying mutations at the residues lining two tunnels, previously observed by protein X-ray crystallography, were constructed and biochemically characterized. The mutants showed distinct catalytic efficiencies with the halogenated substrate 1,2,3-trichloropropane. Release pathways for the two dehalogenation products, 2,3-dichloropropane-1-ol and the chloride ion, and exchange pathways for water molecules, were studied using classical and random acceleration molecular dynamics simulations. Five different pathways, denoted p1, p2a, p2b, p2c, and p3, were identified. The individual pathways showed differing selectivity for the products: the chloride ion releases solely through p1, whereas the alcohol releases through all five pathways. Water molecules play a crucial role for release of both products by breakage of their hydrogen-bonding interactions with the active-site residues and shielding the charged chloride ion during its passage through a hydrophobic tunnel. Exchange of the chloride ions, the alcohol product, and the waters between the buried active site and the bulk solvent can be realized by three different mechanisms: (i) passage through a permanent tunnel, (ii) passage through a transient tunnel, and (iii) migration through a protein matrix. We demonstrate that the accessibility of the pathways and the mechanisms of ligand exchange were modified by mutations. Insertion of bulky aromatic residues in the tunnel corresponding to pathway p1 leads to reduced accessibility to the ligands and a change in mechanism of opening from permanent to transient. We propose that engineering the accessibility of tunnels and the mechanisms of ligand exchange is a powerful strategy for modification of the functional properties of enzymes with buried active sites.     

The structural dynamics of myoglobin

Conformational fluctuations in proteins were initially invoked to explain the observation that diffusion of small ligands through the matrix is a global phenomenon. Small globular proteins contain internal cavities that play a role not only in matrix dynamics but also in controlling function, tracing a pathway for the diffusion of the ligand to and from the active site. This is the main point addressed in this Review, which presents pertinent information obtained on myoglobin (Mb). Mb, a simple globular heme protein which binds reversibly oxygen and other ligands. The bond between the heme Fe(II) and gaseous ligands can be photodissociated by a laser pulse, generating a non-equilibrium population of protein structures that relaxes on a picosecond to millisecond time range. This process is associated with migration of the ligand to internal cavities of the protein, which are known to bind xenon. Some of the results obtained by laser photolysis, molecular dynamics simulations, and X-ray diffraction of intermediate states of wild-type and mutant myoglobins are summarized. The extended relaxation of the globin moiety directly observed by Laue crystallography reflects re-equilibration among conformational substates known to play an essential role in controlling protein function.     

Hydration energy landscape of the active site cavity in cytochrome P450cam

Hydration of protein cavities influences protein stability, dynamics, and function. Protein active sites usually contain water molecules that, upon ligand binding, are either displaced into bulk solvent or retained to mediate protein–ligand interactions. The contribution of water molecules to ligand binding must be accounted for to compute accurate values of binding affinities. This requires estimation of the extent of hydration of the binding site. However, it is often difficult to identify the water molecules involved in the binding process when ligands bind on the surface of a protein. Cytochrome P450cam is, therefore, an ideal model system because its substrate binds in a buried active site, displacing partially disordered solvent, and the protein is well characterized experimentally. We calculated the free energy differences for having five to eight water molecules in the active site cavity of the unliganded enzyme from molecular dynamics simulations by thermodynamic integration employing a three-stage perturbation scheme. The computed free energy differences between the hydration states are small (within 12 kJ mol⁻¹) but distinct. Consistent with the crystallographic determination and studies employing hydrostatic pressure, we calculated that, although ten water molecules could in principle occupy the volume of the active site, occupation by five to six water molecules is thermodynamically most favorable. Proteins 32:381–396, 1998. © 1998 Wiley-Liss, Inc.     

The particle concept: placing discrete water molecules during protein-ligand docking predictions.

Water is known to play a significant role in the formation of protein-ligand complexes. In this paper, we focus on the influence of water molecules on the structure of protein-ligand complexes. We present an algorithmic approach, called the particle concept, for integrating the placement of single water molecules in the docking algorithm of FLEXX. FLEXX is an incremental construction approach to ligand docking consisting of three phases: the selection of base fragments, the placement of the base fragments, and the incremental reconstruction of the ligand inside the active site of a protein. The goal of the extension is to find water molecules at favorable places in the protein-ligand interface which may guide the placement of the ligand. In a preprocessing phase, favorable positions of water molecules inside the active site are calculated and stored in a list of possible water positions. During the incremental construction phase, water molecules are placed at the precomputed positions if they can form additional hydrogen bonds to the ligand. Steric constraints resulting from the water molecules as well as the geometry of the hydrogen bonds are used to optimize the ligand orientation in the active site during the reconstruction process. We have tested the particle concept on a series of 200 protein-ligand complexes. Although the average improvement of the prediction results is minor, we were able to predict water molecules between the protein and the ligand correctly in several cases. For instance in the case of HIV-1 protease, where a single water molecule between the protein and the ligand is known to be of importance in complex formation, significant improvements can be achieved.     

How catalase recognizes H2O2 in a sea of water

Monofunctional heme‐catalases have been studied for many decades but there is still an incomplete understanding of why such a large tetrameric protein with deeply buried active sites is required to accomplish such a simple reaction as H₂O₂ dismutation. Catalase accomplishes this reaction at a high rate although water at 55 M is expected to compete with H₂O₂ for the enzyme's active site. Using molecular dynamics simulations we addressed the question as to how catalase selects H₂O₂ in water. Selection is accomplished through different mechanisms: higher residence time of H₂O₂ in the vicinity of certain prevalent amino acid residues at the protein surface and substrate channel, coordinated motion of the main passage amino acids that is increased in the presence of H₂O₂, a gate valve mechanism consisting of the motion of two contiguous phenylalanine residues that drive water molecules out of the final section of the substrate channel, a hydrophobic barrier before the active site that was crossed more easily by H₂O₂ which kept most of its hydrogen bonds while passing, and finally an increased residence time for H₂O₂ at the active site. These mechanisms, based on the physicochemical differences between H₂O₂ and water, provide an explanation as to why such a large tetrameric protein with deeply buried active sites is required to accomplish efficient H₂O₂ dismutation. Proteins 2014; 82:45–56. © 2013 Wiley Periodicals, Inc.     

Role of a conserved glutamine residue in tuning the catalytic activity of Escherichia coli cytochrome c nitrite reductase

The pentaheme cytochrome c nitrite reductase (NrfA) of Escherichia coli is responsible for nitrite reduction during anaerobic respiration when nitrate is scarce. The NrfA active site consists of a hexacoordinate high-spin heme with a lysine ligand on the proximal side and water/hydroxide or substrate on the distal side. There are four further highly conserved active site residues including a glutamine (Q263) positioned 8 A from the heme iron for which the side chain, unusually, coordinates a conserved, essential calcium ion. Mutation of this glutamine to the more usual calcium ligand, glutamate, results in an increase in the K m for nitrite by around 10-fold, while V max is unaltered. Protein film voltammetry showed that lower potentials were required to detect activity from NrfA Q263E when compared with native enzyme, consistent with the introduction of a negative charge into the vicinity of the active site heme. EPR and MCD spectroscopic studies revealed the high spin state of the active site to be preserved, indicating that a water/hydroxide molecule is still coordinated to the heme in the resting state of the enzyme. Comparison of the X-ray crystal structures of the as-prepared, oxidized native and mutant enzymes showed an increased bond distance between the active site heme Fe(III) iron and the distal ligand in the latter as well as changes to the structure and mobility of the active site water molecule network. These results suggest that an important function of the unusual Q263-calcium ion pair is to increase substrate affinity through its role in supporting a network of hydrogen bonded water molecules stabilizing the active site heme distal ligand.     

Bromodomains识别并结合配体的机理研究

关于酵母重组蛋白内的Bromodomain识别乙酰化赖氨酸的研究近年来受到广泛的关注,但是其识别配体并与之相紧密结合的机理有待进一步的研究。本文采用2015年开发的结合位点拓扑学方法(FCTM)和分子动力学模拟的方式对Bromodomains识别并结合配体的机理进行了充分研究,其中分子动力学模拟时间达24 ns。通过FCTM方法发现结合位点的几何结构具有高度的凹性,且其alphaspace达到了131。分子动力学模拟的结果显示:在模拟的过程中结合位点表面的脯氨酸(Pro66)始终对配体保持着强的分子间相互作用,同时pocket内的水分子分布对配体的氢键网络也一直存在影响。以上结果表明Bromodomains识别并结合配体有两个重要因素:蛋白结构域自身的几何结构和配体受到来自于结合位点表面的氨基酸分子相互作用和pocket内水分子的氢键网络作用。     

Elastic network normal mode dynamics reveal the GPCR activation mechanism

G‐protein‐coupled receptors (GPCR) are a family of membrane‐embedded metabotropic receptors which translate extracellular ligand binding into an intracellular response. Here, we calculate the motion of several GPCR family members such as the M2 and M3 muscarinic acetylcholine receptors, the A_2A adenosine receptor, the β₂‐adrenergic receptor, and the CXCR4 chemokine receptor using elastic network normal modes. The normal modes reveal a dilation and a contraction of the GPCR vestibule associated with ligand passage, and activation, respectively. Contraction of the vestibule on the extracellular side is correlated with cavity formation of the G‐protein binding pocket on the intracellular side, which initiates intracellular signaling. Interestingly, the normal modes of rhodopsin do not correlate well with the motion of other GPCR family members. Electrostatic potential calculation of the GPCRs reveal a negatively charged field around the ligand binding site acting as a siphon to draw‐in positively charged ligands on the membrane surface. Altogether, these results expose the GPCR activation mechanism and show how conformational changes on the cell surface side of the receptor are allosterically translated into structural changes on the inside. Proteins 2014; 82:579–586. © 2013 Wiley Periodicals, Inc.     

Structural dynamics of myoglobin

Conformational fluctuations have been invoked to explain the observation that the diffusion of small ligands through a protein is a global phenomenon, as suggested (for example) by the oxygen induced fluorescence quenching of buried tryptophans. In enzymes processing large substrates, a channel to the catalytic site is often seen in the crystal structure; on the other hand in small globular proteins, it is not known if the cavities identified in the interior space are important in controlling their function by defining specific pathways in the diffusion to the active site. This point is addressed in this paper, which reports some relevant results obtained on myoglobin, the hydrogen atom of molecular biology. Protein conformational relaxations have been extensively investigated with myoglobin because the photosensivity of the adduct with CO, O2 and NO allows us to follow events related to the migration of the ligand through the matrix. Results obtained by laser photolysis, molecular dynamics simulations, X-ray diffraction of intermediate states of wt type and mutant myoglobins are briefly summarized. Crystallographic data on the photochemical intermediate of a new triple mutant of sperm whale myoglobin (Mb-YQR) show, for the first time, the photolyzed CO* sitting in one of the Xe-binding cavities, removed from the heme group. These results support the viewpoint that pre-existing 'packing defects' in the protein interior play a major role in controlling the dynamics of ligand binding, including oxygen, and thereby acquire a survival value.     

Strong and weak hydrogen bonds in protein-ligand complexes of kinases: a comparative study

Strong and weak hydrogen bonds between protein and ligand are analyzed in a group of 233 X-ray crystal structures of the kinase family. These kinases are from both eukaryotic and prokaryotic organisms. The dataset comprises of 44 sub-families, out of which 35 are of human origin and the rest belong to other organisms. Interaction analysis was carried out in the active sites, defined here as a sphere of 10 A radius around the ligand. A majority of the interactions are observed between the main chain of the protein and the ligand atoms. As a donor, the ligand frequently interacts with amino acid residues like Leu, Glu and His. As an acceptor, the ligand interacts often with Gly, and Leu. Strong hydrogen bonds N-H...O, O-H...O, N-H...N and weak bonds C-H...O, C-H...N are common between the protein and ligand. The hydrogen bond donor capacity of Gly in N-H...O and C-H...O interactions is noteworthy. Similarly, the acceptor capacity of main chain Glu is ubiquitous in several kinase sub-families. Hydrogen bonds between protein and ligand form characteristic hydrogen bond patterns (supramolecular synthons). These synthon patterns are unique to each sub-family. The synthon locations are conserved across sub-families due to a higher percentage of conserved sequences in the active sites. The nature of active site water molecules was studied through a novel classification scheme, based on the extent of exposure of water molecules. Water which is least exposed usually participates in hydrogen bond formation with the ligand. These findings will help structural biologists, crystallographers and medicinal chemists to design better kinase inhibitors.     

Molecular dynamics simulation of Axillaridine–A: A potent natural cholinesterase inhibitor

Molecular Dynamics (MD) simulations were carried out for human acetylcholinesterase (hAChE) and its complex with Axillaridine–A, in order to dynamically explore the active site of the protein and the behaviour of the ligand at the peripheral binding site. Simulation of the enzyme alone showed that the active site of AChE is located at the bottom of a deep and narrow cavity whose surface is lined with rings of aromatic residues while Tyr72 is almost perpendicular to the Trp286, which is responsible for stable π -π interactions. The complexation of AChE with Axillaridine-A, results in the reduction of gorge size due to interaction between the ligand and the active site residues. The gorge size was determined by the distance between the center of mass of Glu81 and Trp286. As far as the geometry of the active site is concerned, the presence of ligand in the active site alters its specific conformation, as revealed by stable hydrogen bondings established between amino acids. With the increasing interaction between ligand and the active amino acids, size of the active site of the complex decreases with respect to time. Axillaridine-A, forms stable π -π interactions with the aromatic ring of Tyr124 that results in inhibition of catalytic activity of the enzyme. This π -π interaction keeps the substrate stable at the edge of the catalytic gorge by inhibiting its catalytic activity. The MD results clearly provide an explanation for the binding pattern of bulky steroidal alkaloids at the active site of AChE.     

Characterization of the channel constriction allowing the access of the substrate to the active site of yeast oxidosqualene cyclase

In oxidosqualene cyclases (OSCs), an enzyme which has been extensively studied as a target for hypocholesterolemic or antifungal drugs, a lipophilic channel connects the surface of the protein with the active site cavity. Active site and channel are separated by a narrow constriction operating as a mobile gate for the substrate passage. In Saccharomyces cerevisiae OSC, two aminoacidic residues of the channel/constriction apparatus, Ala525 and Glu526, were previously showed as critical for maintaining the enzyme functionality. In this work sixteen novel mutants, each bearing a substitution at or around the channel constrictions, were tested for their enzymatic activity. Modelling studies showed that the most functionality-lowering substitutions deeply alter the H-bond network involving the channel/constriction apparatus. A rotation of Tyr239 is proposed as part of the mechanism permitting the access of the substrate to the active site. The inhibition of OSC by squalene was used as a tool for understanding whether the residues under study are involved in a pre-catalytic selection and docking of the substrate oxidosqualene.     

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.     

Effects of ligand binding on the mechanical stability of protein GB1 studied by steered molecular dynamics simulation

Regulation of the mechanical properties of proteins plays an important role in many biological processes, and sheds light on the design of biomaterials comprised of protein. At present, strategies to regulate protein mechanical stability focus mainly on direct modulation of the force-bearing region of the protein. Interestingly, the mechanical stability of GB1 can be significantly enhanced by the binding of Fc fragments of human IgG antibody, where the binding site is distant from the force-bearing region of the protein. The mechanism of this long-range allosteric control of protein mechanics is still elusive. In this work, the impact of ligand binding on the mechanical stability of GB1 was investigated using steered molecular dynamics simulation, and a mechanism underlying the enhanced protein mechanical stability is proposed. We found that the external force causes deformation of both force-bearing region and ligand binding site. In other words, there is a long-range coupling between these two regions. The binding of ligand restricts the distortion of the binding site and reduces the deformation of the force-bearing region through a long-range allosteric communication, which thus improves the overall mechanical stability of the protein. The simulation results are very consistent with previous experimental observations. Our studies thus provide atomic-level insights into the mechanical unfolding process of GB1, and explain the impact of ligand binding on the mechanical properties of the protein through long-range allosteric regulation, which should facilitate effective modulation of protein mechanical properties.     

A core process in receptor function, general anesthesia, sleep, and aging

The diffusion of ligands and proteins was proposed to be guided by chreodes in water organized by protein-surface side chains with varying hydropathic states. These chreodes are proposed to be the target of volatile general anesthetic agents. The similarity between this effect and sleep deprivation leads to a proposal of an external agent responsible for sleep. This agent is elemental nitrogen. An extension of this effect is the concept that elemental nitrogen is a core factor in aging.     

A salt-bridge switch in the molecular recognition between RS receptor and RGD ligand from the ABEEM σπ molecular dynamics simulations

Abstract

How the receptor and ligand recognise each other is a challenging subject in explaining the mechanism of recognition at the molecular level. As a starting point, here, a synthesised RS receptor and its RGD ligand were investigated as a proper model to simulate their recognition process in terms of ABEEMσπ/MM polarisable force field. It is found that a switch of forming up a salt bridge in the ligand triggers the recognition of the receptor and ligand. Through the salt-bridge switch that undergoes several cycles from on-state with parallel hydrogen bonds to off-state with bifurcated hydrogen bonds, the active site of ligand can flex easily to interact with the active site of the receptor. In addition, the water molecules form a decisive bridge connecting the active sites of the bound system. The salt-bridge switch and water-mediated movement are cooperative as the important factors for the receptor-ligand recognition. In addition, the properties, such as binding free energy, conformational flexibility and solvent accessible surface area have been calculated to provide adequate evidence for the whole recognition process. According to the simulation, a detailed mechanism was derived involving diffusion, a switch triggered cooperative water-mediated movement, and conformational folding, for the flexible recognition.     

Exploring the binding sites of the haloalkane dehalogenase DhlA from Xanthobacter autotrophicus GJ10

The catalytic site of haloalkane dehalogenase DhlA is buried more than 10 A from the protein surface. While potential access channels to this site have been reported, the precise mechanism of substrate import and product export is still unconfirmed. We used computational methods to examine surface pockets and their putative roles in ligand access to and from the catalytic site. Computational solvent mapping moves small organic molecule as probes over the protein surface in order to identify energetically favorable sites, that is, regions that tend to bind a variety of molecules. The mapping of three DhlA structures identifies seven such regions, some of which have been previously suggested to be involved in the binding and the import/export of substrates or products. These sites are the active site, the putative entrance of the channel leading to the active site, two pockets that bind Br- ions, a pocket in the slot region, and two additional sites between the main domain and the cap of DhlA. We also performed mapping and free energy analysis of the DhlA structures using the substrate, 1,2-dichloroethane, and halide ions as probes. The findings were compared to crystallographic data and to results obtained by CAVER, a program developed for finding routes from protein clefts and cavities to the surface. Solvent mapping precisely reproduced all three Br- binding sites identified by protein crystallography and the openings to four channels found by CAVER. The analyses suggest that (i) the active site has the highest affinity for the substrate molecule, (ii) the substrate initially binds at the entrance of the main tunnel, (iii) the site Br2, close to the entrance, is likely to serve as an intermediate binding site in product export, (iv) the site Br3, induced in the structure at high concentrations of Br-, could be part of an auxiliary route for product release, and (v) three of the identified sites are likely to be entrances of water-access channels leading to the active site. For comparison, we also mapped haloalkane dehalogenases DhaA and LinB, both of which contain significantly larger and more solvent accessible binding sites than DhlA. The mapping of DhaA and LinB places the majority of probes in the active site, but most of the other six regions consistently identified in DhlA were not observed, suggesting that the more open active site eliminates the need for intermediate binding sites for the collision complex seen in DhlA.     

Allosteric regulation and catalysis emerge via a common route

Allosteric regulation of protein function is a mechanism by which an event in one place of a protein structure causes an effect at another site, much like the behavior of a telecommunications network in which a collection of transmitters, receivers and transceivers communicate with each other across long distances. For example, ligand binding or an amino acid mutation at an allosteric site can alter enzymatic activity or binding affinity in a distal region such as the active site or a second binding site. The mechanism of this site-to-site communication is of great interest, especially since allosteric effects must be considered in drug design and protein engineering. In this review, conformational mobility as the common route between allosteric regulation and catalysis is discussed. We summarize recent experimental data and the resulting insights into allostery within proteins, and we discuss the nature of future studies and the new applications that may result from increased understanding of this regulatory mechanism.     

Molecular mimicry of substrate oxygen atoms by water molecules in the beta-amylase active site

Soybean beta-amylase (EC 3.2.1.2) has been crystallized both free and complexed with a variety of ligands. Four water molecules in the free-enzyme catalytic cleft form a multihydrogen-bond network with eight strategic residues involved in enzyme-ligand hydrogen bonds. We show here that the positions of these four water molecules are coincident with the positions of four potential oxygen atoms of the ligands within the complex. Some of these waters are displaced from the active site when the ligands bind to the enzyme. How many are displaced depends on the shape of the ligand. This means that when one of the four positions is not occupied by a ligand oxygen atom, the corresponding water remains. We studied the functional/structural role of these four waters and conclude that their presence means that the conformation of the eight side chains is fixed in all situations (free or complexed enzyme) and preserved from unwanted or forbidden conformational changes that could hamper the catalytic mechanism. The water structure at the active pocket of beta-amylase is therefore essential for providing the ligand recognition process with plasticity. It does not affect the protein active-site geometry and preserves the overall hydrogen-bonding network, irrespective of which ligand is bound to the enzyme. We also investigated whether other enzymes showed a similar role for water. Finally, we discuss the potential use of these results for predicting whether water molecules can mimic ligand atoms in the active center.