A microscopic model for gas diffusion dynamics in a [NiFe]-hydrogenase

We describe and apply a microscopic model for the calculation of gas diffusion rates in a [NiFe]-hydrogenase. This enzyme has attracted much interest for use as a H(2) oxidising catalyst in biofuel cells, but a major problem is their inhibition by CO and O(2). In our model, the diffusive hopping of gas molecules in the protein interior is coarse grained using a master equation approach with transition rates estimated from equilibrium and non-equilibrium pulling simulations. Propagating the rate matrix in time, we find that the probability for a gas molecule to reach the enzyme active site follows a mono-exponential increase. Fits to a phenomenological rate law give an effective diffusion rate constant for CO that is in very good agreement with experimental measurements. We find that CO prefers to move along the canonical 'hydrophobic' main channel towards the active site, in contrast to O(2) and H(2), which were previously shown to explore larger fractions of the protein. Differences in the diffusion of the three gases are discussed in light of recent efforts to engineer a gas selectivity filter in the enzyme.     

Transient pockets as mediators of gas molecules routes inside proteins: The case study of dioxygen pathway in homogentisate 1,2-dioxygenase and its implication in Alkaptonuria development

Alkaptonuria (AKU) is an ultra-rare disease caused by mutations in homogentisate 1,2-dioxygenase (HGD) enzyme, characterized by the loss of enzymatic activity and the accumulation of its substrate, homogentisic acid (HGA) in different tissues, leading to ochronosis and organ degeneration. Although the pathological effects of HGD mutations are largely studied, less is known about the structure of the enzyme, in particular the pathways for dioxygen diffusion to the active site, required for the enzymatic reaction, are still uninvestigated. In the present project, the combination of two in silico techniques, Molecular Dynamics (MD) simulation and Implicit Ligand Sampling (ILS), was used to delineate gas diffusion routes in HGD enzyme. A route from the central opening of the hexameric structure of the enzyme to the back of the active site trough the protein moiety was identified as the path for dioxygen diffusion, also overlapping with a transient pocket, which then assumes an important role in dioxygen diffusion. Along the route the sequence location of the missense variant E401Q, responsible for AKU development, was also found, suggesting such mutation to be conducive of enzymatic activity loss by altering the flow dynamics of dioxygen. Our in silico approach allowed also to delineate the route of HGA substrate to the active site, until now only supposed.     

Carbonic anhydrase inhibitors: X-ray and molecular modeling study for the interaction of a fluorescent antitumor sulfonamide with isozyme II and IX

The X-ray crystal structure of the fluorescent antitumor sulfonamide carbonic anhydrase (CA, EC, 4.2.1.1) inhibitor (4-sulfamoylphenylethyl)thioureido fluorescein (1) in complex with the cytosolic isoform hCA II is reported, together with a modeling study of the adduct of 1 with the tumor-associated isoform hCA IX. Its binding to hCA II is similar to that of other benzesulfonamides, with the ionized sulfonamide coordinated to the Zn2+ ion within the enzyme active site, and also participating in a network of hydrogen bonds with residues Thr199 and Glu106. The scaffold of 1 did not establish polar interactions within the enzyme active site but made hydrophobic contacts (<4.5 A) with Gln92, Val121, Phe131, Val135, Leu198, Thr199, Thr200, and Pro202. The substituted 3-carboxy-amino-phenyl functionality was at van der Waals distance from Phe131, Gly132, and Val135. The bulky tricyclic fluorescein moiety was located at the rim of the active site, on the protein surface, and strongly interacted with the alpha-helix formed by residues Asp130-Val135. All these interactions were preserved in the hCA IX-1 adduct, but the carbonyl moiety of the fluorescein tail of 1 participates in a strong hydrogen bond with the guanidine moiety of Arg130, an amino acid characteristic of the hCA IX active site. This may account for the roughly 2 times higher affinity of 1 for hCA IX over hCA II and may explain why in vivo the compound specifically accumulates only in hypoxic tumors overexpressing CA IX and not in the normal tissues. The compound is in clinical studies as an imaging tool for acute hypoxic tumors.     

Structural analysis of charge discrimination in the binding of inhibitors to human carbonic anhydrases I and II

Despite the similarity in the active site pockets of carbonic anhydrase (CA) isozymes I and II, the binding affinities of benzenesulfonamide inhibitors are invariably higher with CA II as compared to CA I. To explore the structural basis of this molecular recognition phenomenon, we have designed and synthesized simple benzenesulfonamide inhibitors substituted at the para position with positively charged, negatively charged, and neutral functional groups, and we have determined the affinities and X-ray crystal structures of their enzyme complexes. The para-substituents are designed to bind in the midsection of the 15 A deep active site cleft, where interactions with enzyme residues and solvent molecules are possible. We find that a para-substituted positively charged amino group is more poorly tolerated in the active site of CA I compared with CA II. In contrast, a para-substituted negatively charged carboxylate substituent is tolerated equally well in the active sites of both CA isozymes. Notably, enzyme-inhibitor affinity increases upon neutralization of inhibitor charged groups by amidation or esterification. These results inform the design of short molecular linkers connecting the benzenesulfonamide group and a para-substituted tail group in "two-prong" CA inhibitors: an optimal linker segment will be electronically neutral, yet capable of engaging in at least some hydrogen bond interactions with protein residues and/or solvent. Microcalorimetric data reveal that inhibitor binding to CA I is enthalpically less favorable and entropically more favorable than inhibitor binding to CA II. This contrasting behavior may arise in part from differences in active site desolvation and the conformational entropy of inhibitor binding to each isozyme active site.     

喹啉加氧酶中氧气扩散的通道分析

通过生物信息学分析、量化计算优化、CAVER和MDpocket预测、随机加速分子动力学及伞状抽样动力学模拟等方法,对喹啉加氧酶(HOD)中的氧气扩散途径进行了计算预测.结果表明,氧气在HOD中的反应位点包埋在蛋白内部,而HOD中有数条可能的通道供氧气进出,其中长度最短的通道具有最高的优先度,不仅在随机加速动力学模拟中具有最高的氧气逸出概率,而且伞状抽样方法计算得到的自由能也最低.此通道的内端位于底物Re面的氧气结合位点,较好地解释了HOD的相关实验数据.     

Intricate role of water in proton transport through cytochrome c oxidase

Cytochrome c oxidase (CytcO), the final electron acceptor in the respiratory chain, catalyzes the reduction of O(2) to H(2)O while simultaneously pumping protons across the inner mitochondrial or bacterial membrane to maintain a transmembrane electrochemical gradient that drives, for example, ATP synthesis. In this work mutations that were predicted to alter proton translocation and enzyme activity in preliminary computational studies are characterized with extensive experimental and computational analysis. The mutations were introduced in the D pathway, one of two proton-uptake pathways, in CytcO from Rhodobacter sphaeroides . Serine residues 200 and 201, which are hydrogen-bonded to crystallographically resolved water molecules halfway up the D pathway, were replaced by more bulky hydrophobic residues (Ser200Ile, Ser200Val/Ser201Val, and Ser200Val/Ser201Tyr) to query the effects of changing the local structure on enzyme activity as well as proton uptake, release, and intermediate transitions. In addition, the effects of these mutations on internal proton transfer were investigated by blocking proton uptake at the pathway entrance (Asp132Asn replacement in addition to the above-mentioned mutations). Even though the overall activities of all mutant CytcO's were lowered, both the Ser200Ile and Ser200Val/Ser201Val variants maintained the ability to pump protons. The lowered activities were shown to be due to slowed oxidation kinetics during the P(R) → F and F → O transitions (P(R) is the "peroxy" intermediate formed at the catalytic site upon reaction of the four-electron-reduced CytcO with O(2), F is the oxoferryl intermediate, and O is the fully oxidized CytcO). Furthermore, the P(R) → F transition is shown to be essentially pH independent up to pH 12 (i.e., the apparent pK(a) of Glu286 is increased from 9.4 by at least 3 pK(a) units) in the Ser200Val/Ser201Val mutant. Explicit simulations of proton transport in the mutated enzymes revealed that the solvation dynamics can cause intriguing energetic consequences and hence provide mechanistic insights that would never be detected in static structures or simulations of the system with fixed protonation states (i.e., lacking explicit proton transport). The results are discussed in terms of the proton-pumping mechanism of CytcO.     

Enzymatic methyl transfer: role of an active site residue in generating active site compaction that correlates with catalytic efficiency

Human catechol-O-methyltransferase (COMT) catalyzes a methyl transfer from S-adenosylmethionine (AdoMet) to dopamine. Site-specific mutants at three positions (Tyr68, Trp38, and Val108) have been characterized with regard to product distribution, catalytic efficiency, and secondary kinetic isotope effects. The series of mutations at Tyr68 within wild-type protein and the common polymorphic variant (Val108Met) yields a linear correlation between the catalytic efficiency and the size of the secondary kinetic isotope effect. We conclude that active site compaction in COMT is modulated by a proximal side chain residing behind the sulfur-bearing methyl group of AdoMet. These findings are discussed in the context of the active site compression that has been postulated to accompany enzyme-supported hydrogen tunneling.     

Focusing mutations into the P. fluorescens esterase binding site increases enantioselectivity more effectively than distant mutations

Rational design of enzymes with improved properties, such as enantioselectivity, usually focuses mutations within the substrate binding site. On the other hand, directed evolution of enzymes usually targets the entire protein and discovers beneficial mutations far from the substrate binding site. In this paper, we propose an explanation for this discrepancy and show that a combined approach--random mutagenesis within the substrate binding site--is better. To increase the enantioselectivity (E) of a Pseudomonas fluorescens esterase (PFE) toward methyl 3-bromo-2-methylpropionate, we focused mutagenesis into the substrate binding site at Trp28, Val121, Phe198, and Val225. Five of the catalytically active mutants (13%) showed better enantioselectivity than wild-type PFE. The increases in enantioselectivity were higher (up to 5-fold, reaching E = 61) than with mutants identified by random mutagenesis of the entire enzyme.     

Diffusion-influenced ligand binding to buried sites in macromolecules and transmembrane channels

We consider diffusion-influenced binding to a buried binding site that is connected to the surface by a narrow tunnel. Under the single assumption of an equilibrium distribution of ligands over the tunnel cross section, we reduce the calculation of the time-dependent rate coefficient to the solution of a one-dimensional diffusion equation with appropriate boundary conditions. We obtain a simple analytical expression for the steady-state rate that depends on the potential of mean force in the tunnel and the diffusion-controlled rate of binding to the tunnel entrance. Potential applications of our theory include substrate binding to a buried active site of an enzyme and permeant ion binding to an internal site in a transmembrane channel.     

Numerical Simulation of Nonequilibrium Species Diffusion in a Low-Power Nitrogen–Hydrogen Arcjet Thruster

Species distributions in a low-power arcjet thruster are investigated using a two-dimensional thermal and chemical nonequilibrium numerical model that incorporates the self-consistent effective binary diffusion coefficient approximation treatment of diffusion. Plasma flows in arcjet thruster with different input mole ratios of nitrogen to hydrogen are modelled. It is found that species separation due to nonequilibrium chemical kinetic processes occurs mainly in the regions where the dissociation and ionization of nitrogen and hydrogen species take place. The enrichment of nitrogen molecules at the fringes of the arc and hydrogen molecules near the anode wall of the thruster occurs mainly because the recombination processes of these two gases occur in different temperature ranges. In the expansion portion of the thruster nozzle, the gas residence times are of the same order as some chemical kinetic processes. Comparison between the nitrogen and hydrogen species profiles at the constrictor and thruster exit shows that the recombination of hydrogen ions and atoms are dominant kinetic processes near the thruster centreline, while the chemical reactions of nitrogen species are almost frozen in the high speed flow. The effects of temperature and pressure gradients on the species diffusion inside the arcjet thruster are also presented, with thermal diffusion found to have a much larger influence than pressure diffusion.     

Theoretical study of the mechanisms of substrate recognition by catalase

A variety of theoretical methods including classical molecular interaction potentials, classical molecular dynamics, and activated molecular dynamics have been used to analyze the substrate recognition mechanisms of peroxisomal catalase from Saccharomyces cerevisiae. Special attention is paid to the existence of channels connecting the heme group with the exterior of the protein. On the basis of these calculations a rationale is given for the unique catalytic properties of this enzyme, as well as for the change in enzyme efficiency related to key mutations. According to our calculations the water is expected to be a competitive inhibitor of the enzyme, blocking the access of hydrogen peroxide to the active site. The main channel is the preferred route for substrate access to the enzyme and shows a cooperative binding to hydrogen peroxide. However, the overall affinity of the main channel for H(2)O(2) is only slightly larger than that for H(2)O. Alternative channels connecting the heme group with the monomer interface and the NADP(H) binding site are detected. These secondary channels might be important for product release.     

Hydrogen bonding pathways in human dihydroorotate dehydrogenase

Dihydroorotate dehydrogenase (DHOD) catalyzes the only redox reaction in the pathway for pyrimidine biosynthesis. In this reaction, a proton is transferred from a carbon atom of the substrate to a serine residue, and a hydride is transferred from another carbon atom of the substrate to a cofactor. The deprotonation of the substrate is postulated to involve a proton relay mechanism along a hydrogen bonding pathway in the active site. In this paper, molecular dynamics simulations are used to identify and characterize potential hydrogen bonding pathways that could facilitate the redox reaction catalyzed by human DHOD. The observed pathways involve hydrogen bonding of the active base serine to a water molecule, which is hydrogen bonded to the substrate carboxylate group or a threonine residue. The threonine residue is positioned to enable proton transfer to another water molecule leading to the bulk solvent. The impact of mutating the active base serine to cysteine is also investigated. This mutation is found to increase the average donor-acceptor distances for proton and hydride transfer and to disrupt the hydrogen bonding pathways observed for the wild-type enzyme. These effects could lead to a significant decrease in enzyme activity, as observed experimentally for the analogous mutant in Escherichia coli DHOD.     

Computational studies for reduced graphene oxide in hydrogen-rich environment

We employ molecular dynamic simulations to study the reduction process of graphene oxide (GO) in a chemically active environment enriched with hydrogen. We examine the concentration and pressure of hydrogen gas as a function of temperature in which abstraction of oxygen is possible with minimum damage to C-sp(2) bonds, hence preserving the integrity of the graphene sheet. Through these studies we find chemical pathways that demonstrate beneficiary mechanisms for the quality of graphene including formation of water as well as suppression of carbonyl pair holes in favor of hydroxyl and epoxide formation facilitated by hydrogen gas in the environment.     

Diffusion of gas molecules in the polystyrene matrix

We have studied the diffusion of gas molecules inside an amorphous polystyrene matrix. The diffusion constant of several gases at T = 450 K in polystyrene was computed. Particular attention was given to CO₂ for temperatures between 300 and 800 K. The temperature dependence of the diffusion constant and the relationship between the diffusion constant and the diameter of the gas molecule were analysed. We further examined the motion of the gas molecules on the short time scale not readily accessible to experimental observation. Here we used the cage overlap function which gives information on the typical cage sizes and distribution times. On the short time scale the gas molecules show a hopping behaviour. The distribution of the time period between hopping events, the distance between the cages and the size of the cages in the polystyrene matrix in the presence of guest molecules were calculated. Through the analyzing, we get a clearer picture from the behaviour of gas molecules on short time scale in the polymer matrix.     

HIV蛋白酶抑制剂——利迪链菌素的分子对接研究

用分子对接的方法, 对利迪链菌素的抗HIV蛋白酶活性进行了研究. 为了更准确地反映利迪链菌素分子与酶蛋白结合的情况, 充分考虑受体活性部位的柔性, 采用了FlexX(初步对接)和Flexidock(精确对接)分两步将配体与受体进行对接. 在初步对接中, 设计了不同的受体活性部位来考察是否有结合水分子参与抑制剂与酶的结合. 对一种作用方式已知的非肽类HIV蛋白酶抑制剂Aha006进行的对接研究显示, 分子模拟的结果与实际情况吻合得较好, 证明了本文所采用的方法的可靠性. 利迪链菌素与蛋白酶活性部位的对接结果显示, 配体分子与受体之间的结合没有结合水分子的参与, 两者通过5对氢键作用结合成为稳定的复合物. 利迪链菌素占据结合腔, 覆盖了蛋白酶的活性三联体Asp25-Thr26-Gly27, 从而起到抑制其生物活性的作用.     

Gas pressure effects on the rates of catalytic H(2) oxidation by hydrogenases

During catalysis by hydrogenases, entities no larger than H(2) or H(+) reach and leave a deeply buried active site, by as yet unidentified pathways. Novel experiments, conducted mainly with the membrane-bound [NiFe]-hydrogenase from Ralstonia eutropha, explore why small excess gas pressures (H(2) or He) attenuate the rate of H(2) oxidation.     

Contribution of multicentered short hydrogen bond arrays to potency of active site-directed serine protease inhibitors

We describe and compare the pH dependencies of the potencies and of the bound structures of two inhibitor isosteres that form multicentered short hydrogen bond arrays at the active sites of trypsin, thrombin, and urokinase type plasminogen activator (urokinase or uPA) over certain ranges of pH. Depending on the pH, short hydrogen bond arrays at the active site are mediated by two waters, one in the oxyanion hole (H(2)O(oxy)) and one on the other (S2) side of the inhibitor (H(2)O(S2)), by one water (H(2)O(oxy)), or by no water. The dramatic variation in the length of the active site hydrogen bonds as a function of pH, of inhibitor, and of enzyme, along with the involvement or absence of ordered water, produces a large structural manifold of active site hydrogen bond motifs. Diverse examples of multicentered and two-centered short hydrogen bond arrays, both at and away from the active site, recently discovered in several protein crystal systems, suggest that short hydrogen bonds in proteins may be more common than has been recognized. The short hydrogen bond arrays resemble one another with respect to ionic nature, highly polar environment, multitude of associated ordinary hydrogen bonds, and disparate pK(a) values of participating groups. Comparison of structures and K(i) values of trypsin complexes at pH values where the multicentered short hydrogen bond arrays mediating inhibitor binding are present or absent indicate that these arrays have a minor effect on inhibitor potency. These features suggest little covalent nature within the short hydrogen bonds, despite their extraordinary shortness (as short as 2.0 A).     

Phoenix mars mission—The thermal evolved gas analyzer

The Phoenix spacecraft that was launched to Mars in August 2007 landed safely on the Martian northern arctic region on May 25, 2008. It carried six experiments to study the history of water on the planet and search for organic molecules in the icy subsurface Martian soil. The spacecraft is a lander with an arm and scoop designed to dig a trench though the top soil to reach an expected ice layer near the surface. One of the instruments on board is the thermal evolved gas analyzer (TEGA), which consists of two components, a set of eight very small ovens that will heat samples of the ice soil mixtures from the trench to release imbedded gases and mineral decomposition products, and a mass spectrometer that serves as the analysis tool for the evolved gases, and also for measurements of the composition and isotopic ratios of the gases that comprise the atmosphere of Mars. The mass spectrometer is a miniature magnetic sector instrument controlled by microprocessor-driven power supplies. One feature is the gas enrichment cell that will increase the partial pressures of the noble gases in an atmosphere sample by removing all the active gases, carbon dioxide, and nitrogen, to improve the accuracy of their isotopic ratio measurements.     

Mechanistic insight into the inactivation of carboxypeptidase A by alpha-benzyl-2-oxo-1,3-oxazolidine-4-acetic acid, a novel type of irreversible inhibitor for carboxypeptidase A with no stereospecificity

On the basis of the active site topology and enzymic catalytic mechanism of carboxypeptidase A (CPA), a prototypical zinc-containing proteolytic enzyme, alpha-benzyl-2-oxo-1,3-oxazolidine-4-acetic acid (1), was designed as a novel type of mechanism-based inactivator of the enzyme. All four possible stereoisomers of the inhibitor were synthesized in an enantiomerically pure form starting with optically active aspartic acid, and their CPA inhibitory activities were evaluated to find that surprisingly all of the four stereoisomers inhibit CPA in a time dependent manner. The inhibited enzyme did not regain its enzymic activity upon dialysis. The inactivations were prevented by 2-benzylsuccinic acid, a competitive inhibitor that is known to bind the active site of the enzyme. These kinetic results strongly support that the inactivators attach covalently to the enzyme at the active site. The analysis of ESI mass spectral data of the inactivated CPA ascertained the conclusion from the kinetic results. The values of second-order inhibitory rate constants (k(obs)/[I](o)) fall in the range of 1.7-3.6 M(-1) min(-1). The lack of stereospecificity shown in the inactivation led us to propose that the ring cleavage occurs by the nucleophilic attack at the 2-position rather than at the 5-position and the ring opening takes place in an addition-elimination mechanism. The tetrahedral transition state that would be generated in this pathway is thought to be stabilized by the active site zinc ion, which was supported by the PM3 semiemprical calculations. In addition, alpha-benzyl-2-oxo-1,3-oxazolidine-5-acetic acid (18), a structural isomer of 1 was also found to inactivate CPA in an irreversible manner, reinforcing the nucleophilic addition-elimination mechanism. The present study demonstrates that the transition state for the inactivation pathway plays a critical role in determining stereochemistry of the inactivation.