Revelation of a catalytic calcium-binding site elucidates unusual metal dependence of a human apyrase

Intrinsic Apyrase-Like Activity of Cerium-Based Metal-Organic Frameworks (MOFs): Dephosphorylation of Adenosine Tri- and Diphosphate

Apyrase is an important family of extracellular enzymes that catalyse the hydrolysis of high-energy phosphate bonds (HEPBs) in ATP and ADP, thereby modulating many physiological processes and driving life activities. Herein, we report an unexpected discovery that cerium-based metal-organic frameworks (Ce-MOFs) of UiO-66(Ce) have intrinsic apyrase-like activity for ATP/ADP-related physiological processes. The abundant Ce^III /Ce^IV couple sites of Ce-MOFs endow them with the ability to selectively catalyse the hydrolysis of HEPBs of ATP and ADP under physiological conditions. Compared to natural enzymes, they could resist extreme pH and temperature, and present a broad range of working conditions. Based on this finding, a significant inhibitory effect on ADP-induced platelet aggregation was observed upon exposing the platelet-rich plasma (PRP) to the biomimetic UiO-66(Ce) films, prefiguring their wide application potentials in medicine and biotechnology.

Linking inhibitor motions to proteolytic stability of sunflower trypsin inhibitor-1

The remarkable capability of an enzyme isn't only determined by its active site but also controlled by the environment. To unravel the environment role in catalysis, the dynamic motions as well as the static mechanism need to be studied. In this work, QM/MM MD simulations were employed to study the proteolysis process of SFTI-1 and BiKF, which revealed that a combination of static non-bonded interactions and dynamic motions along the reaction coordinate can account for the different hydrolysis rates between them. A comparison among SFTI-1 and three analogs with similar non-bonded interactions further revealed a positive correlation between the mobility of inhibitors and the hydrolysis rates. Apart from the cyclic backbone and disulfide bond, intramolecular hydrogen bonds also increase the rigidity of the backbone of inhibitors, and therefore hinder inhibitor motions to resist proteolysis. These new detailed mechanistic insights suggest the need to consider inhibitor motions in the rational design of peptide inhibitors.

Besides the non-bonded interactions, inhibitor motions especially rotation of the scissile bond also influence proteolytic stability.

Molecular insight into chymotrypsin inhibitor 2 resisting proteolytic degradation

Chymotrypsin inhibitor 2 (CI2) is a special serine protease inhibitor which can resist hydrolysis for several days with a rapid equilibrium between the Michaelis complex and acyl-enzyme intermediate. The energies and conformational changes for subtilisin-catalyzed proteolysis of CI2 were examined in this paper for the first time by employing pseudo bond ab initio QM/MM MD simulations. In the acylation reaction, a low-barrier hydrogen bond between His64 and Asp32 in the transition state together with the lack of covalent backbone constraints makes the peptide bonds of CI2 break more easily than in other serine protease inhibitors. After acyl-enzyme formation, molecular dynamics simulations showed that the access of hydrolytic water to the active site requires partial dissociation of the leaving group. However, retention of the leaving group mainly by the intra- and inter-molecular H-bonding networks hinders the access of water and retards the deacylation reaction. Instead of the dissociation constant of inhibitors, we suggest employing the free energy at the acyl-enzyme state to predict the relative hydrolysis rates of CI2 mutants, which are testified by the experimental relative hydrolysis rates.

Two symmetric arginine residues play distinct roles in Thermus thermophilus Argonaute DNA guide strand-mediated DNA target cleavage

Bacterium Thermus thermophilus Argonaute (Ago; TtAgo) is a prokaryotic Ago (pAgo) that acts as the host defense against the uptake and propagation of foreign DNA by catalyzing the DNA cleavage reaction. The TtAgo active site consists of a plugged-in glutamate finger with two arginine residues (R545 and R486) located symmetrically around it. An interesting challenge is to understand how they can collaboratively facilitate enzymatic catalysis. In Kluyveromyces polysporus Ago, a eukaryotic Ago, the evolutionarily symmetrical residues are arginine and histidine, both of which function to stabilize the plugged-in catalytic tetrad conformation. Surprisingly, our simulation results indicated that, in TtAgo, only R545 is involved in the cleavage reaction by serving as a critical structural anchor to stabilize the catalytic tetrad Asp-Glu-Asp-Asp that is completed by the insertion of the glutamate finger, whereas R486 is not involved in target cleavage. The TtAgo-mediated target DNA cleavage occurs in a substrate-assisted mechanism, in which the pro-Rp (Rp, a tetrahedral phosphorus center with "R-type" chirality) oxygen of scissile phosphate acts as a general base to activate the nucleophilic water. Our unexpected theoretical findings on distinct roles played by R545 and R486 in TtAgo catalysis have been validated by single-point site-mutagenesis experiments, wherein the target cleavage is abolished for all mutants of R545. In sharp contrast, the cleavage activity is maintained for all mutants of R486. Our work provides mechanistic insights on the catalytic specificity of Ago proteins and could facilitate the design of new gene-editing tools in the long term.

Born-Oppenheimer Ab Initio QM/MM Molecular Dynamics Simulations of Enzyme Reactions

There are two key requirements for reliably simulating enzyme reactions: one is a reasonably accurate potential energy surface to describe the bond-forming/breaking process as well as to adequately model the heterogeneous enzyme environment; the other is to perform extensive sampling since an enzyme system consists of at least thousands of atoms and its energy landscape is very complex. One attractive approach to meet both daunting tasks is Born-Oppenheimer ab initio QM/MM molecular dynamics (aiQM/MM-MD) simulation with umbrella sampling. In this chapter, we describe our recently developed pseudobond Q-Chem-Amber interface, which employs a combined electrostatic-mechanical embedding scheme with periodic boundary condition and the particle mesh Ewald method for long-range electrostatics interactions. In our implementation, Q-Chem and the sander module of Amber are combined at the source code level without using system calls, and all necessary data communications between QM and MM calculations are achieved via computer memory. We demonstrate the applicability of this pseudobond Q-Chem-Amber interface by presenting two examples, one reaction in aqueous solution and one enzyme reaction. Finally, we describe our established aiQM/MM-MD enzyme simulation protocol, which has been successfully applied to study more than a dozen enzymes.

Thiol versus hydroxamate as zinc binding group in HDAC inhibition: An ab initio QM/MM molecular dynamics study

Zinc-dependent histone deacetylases (HDACs) play a critical role in transcriptional repression and gene silencing, and are among the most attractive targets for the development of new therapeutics against cancer and various other diseases. Two HDAC inhibitors have been approved by FDA as anti-cancer drugs: one is SAHA whose hydroxamate is directly bound to zinc, the other is FK228 whose active form may use thiol as the zinc binding group. In spite of extensive studies, it remains to be ambiguous regarding how thiol and hydroxamate are bound to the zinc active site of HDACs. In this work, our computational approaches center on Born-Oppenheimer ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics with umbrella sampling, which allow for modeling of the zinc active site with reasonable accuracy while properly including dynamics and effects of protein environment. Meanwhile, an improved short-long effective function (SLEF2) to describe non-bonded interactions between zinc and other atoms has been employed in initial MM equilibrations. Our ab initio QM/MM MD simulations have confirmed that hydroxamate is neutral when it is bound to HDAC8, and found that thiol is deprotonated when directly bound to zinc in the HDAC active site. By comparing thiol and hydroxamate, our results elucidated the differences in their binding environment in the HDAC active sites, and emphasized the importance of the linker design to achieve more specific binding toward class IIa HDACs.

Extracellular Streptomyces lividans vesicles: composition, biogenesis and antimicrobial activity

We selected Streptomyces lividans to elucidate firstly the biogenesis and antimicrobial activities of extracellular vesicles that a filamentous and highly differentiated Gram-positive bacterium produces. Vesicle types range in diameter from 110 to 230 nm and 20 to 60 nm, respectively; they assemble to clusters, and contain lipids and phospholipids allowing their in situ imaging by specific fluorescent dyes. The presence of the identified secondary metabolite undecylprodigiosin provokes red fluorescence of a portion of the heterogeneous vesicle populations facilitating in vivo monitoring. Protuberances containing vesicles generate at tips, and alongside of substrate hyphae, and enumerate during late vegetative growth to droplet-like exudates. Owing to in situ imaging in the presence and absence of a green fluorescent vancomycin derivative, we conclude that protuberances comprising vesicles arise at sites with enhanced levels of peptidoglycan subunits [pentapeptide of lipid II (C55)-linked disaccharides], and reduced levels of polymerized and cross-linked peptidoglycan within hyphae. These sites correlate with enhanced levels of anionic phospholipids and lipids. Vesicles provoke pronounced damages of Aspergillus proliferans, Verticillium dahliae and induced clumping and distortion of Escherichia coli. These harmful effects are likely attributable to the action of the identified vesicular compounds including different enzyme types, components of signal transduction cascades and undecylprodigiosin. Based on our pioneering findings, we highlight novel clues with environmental implications and application potential.

Mechanistic insights into a classic wonder drug--aspirin

Aspirin, one of the oldest and most common anti-inflammatory agents, has recently been shown to reduce cancer risks. The principal pharmacological effects of aspirin are known to arise from its covalent modification of cyclooxygenase-2 (COX-2) through acetylation of Ser530, but the detailed mechanism of its biochemical action and specificity remains to be elucidated. In this work, we have filled this gap by employing a state-of-the-art computational approach, Born-Oppenheimer molecular dynamics simulations with ab initio quantum mechanical/molecular mechanical potential and umbrella sampling. Our studies have characterized a substrate-assisted inhibition mechanism for aspirin acetylating COX: it proceeds in two successive stages with a metastable tetrahedral intermediate, in which the carboxyl group of aspirin serves as the general base. The computational results confirmed that aspirin would be 10-100 times more potent against COX-1 than against COX-2, and revealed that this inhibition specificity between the two COX isoforms can be attributed mainly to the difference in kinetics rate of the covalent inhibition reaction, not the aspirin-binding step. The structural origin of this differential inhibition of the COX enzymes by aspirin has also been elucidated.

Effect of strong electric field on the conformational integrity of insulin

A series of molecular dynamics (MD) simulations up to 1 μs for bovine insulin monomer in different external electric fields were carried out to study the effect of external electric field on conformational integrity of insulin. Our results show that the secondary structure of insulin is kept intact under the external electric field strength below 0.15 V/nm, but disruption of secondary structure is observed at 0.25 V/nm or higher electric field strength. Although the starting time of secondary structure disruption of insulin is not clearly correlated with the strength of the external electric field ranging between 0.15 and 0.60 V/nm, long time MD simulations demonstrate that the cumulative effect of exposure time under the electric field is a major cause for the damage of insulin's secondary structure. In addition, the strength of the external electric field has a significant impact on the lifetime of hydrogen bonds when it is higher than 0.60 V/nm. The fast evolution of some hydrogen bonds of bovine insulin in the presence of the 1.0 V/nm electric field shows that different microwaves could either speed up protein folding or destroy the secondary structure of globular proteins deponding on the intensity of the external electric field.

Catalytic mechanism of histone acetyltransferase p300: from the proton transfer to acetylation reaction

The transcriptional coactivator and histone acetyltransferase (HAT) p300 acetylates the four core histones and other transcription factors to regulate a plethora of fundamental biological processes including cell growth, development, oncogenesis and apoptosis. Recent structural and biochemical studies on the p300 HAT domain revealed a Theorell-Chance, or "hit-and-run", catalytic mechanism. Nonetheless, the chemical mechanism of the entire reaction process including the proton transfer (PT) scheme and consequent acetylation reaction route remains unclear. In this study, a combined computational strategy consisting of molecular modeling, molecular dynamic (MD) simulation, and quantum mechanics/molecular mechanics (QM/MM) simulation was applied to elucidate these important issues. An initial p300/H3/Ac-CoA complex structure was modeled and optimized using a 100 ns MD simulation. Residues that play important roles in substrate binding and the acetylation reaction were comprehensively investigated. For the first time, these studies reveal a plausible PT scheme consisting of Y1394, D1507, and a conserved crystallographic water molecule, with all components of the scheme being stable during the MD simulation and the energy barrier low for PT to occur. The two-dimensional potential energy surface for the nucleophilic attack process was also calculated. The comparison of potential energies for two possible elimination half-reaction mechanisms revealed that Y1467 reprotonates the coenzyme-A leaving group to form product. This study provides new insights into the detailed catalytic mechanism of p300 and has important implications for the discovery of novel small molecule regulators for p300.

A QM/MM MD study of the pH-dependent ring-opening catalysis and lid motif flexibility in glucosamine 6-phosphate deaminase

QM/MM MD and MM MD simulations reveal pH-dependent proton-shuttle ring-opening mechanisms of GlcN6P and dynamical behavior of the lid motif inSmuNagB.

A full picture of enzymatic catalysis by hydroxynitrile lyases from Hevea brasiliensis: protonation dependent reaction steps and residue-gated movement of the substrate and the product

Hydroxynitrile lyases (HNLs) defend plants from herbivores and microbial attack by releasing cyanide from hydroxynitriles.

DNA cytosine methylation: structural and thermodynamic characterization of the epigenetic marking mechanism

DNA cytosine methyltransferases regulate the expression of the genome through the precise epigenetic marking of certain cytosines with a methyl group, and aberrant methylation is a hallmark of human diseases including cancer. Targeting these enzymes for drug design is currently a high priority. We have utilized ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations to investigate extensively the reaction mechanism of the representative DNA methyltransferase HhaI (M.HhaI) from prokaryotes, whose overall mechanism is shared with the mammalian enzymes. We obtain for the first time full free energy profiles for the complete reaction, together with reaction dynamics in atomistic detail. Our results show an energetically preferred mechanism in which nucleophilic attack of cytosine C5 on the S-adenosyl-L-methionine (AdoMet) methyl group is concerted with formation of the Michael adduct between a conserved Cys in the active site with cytosine C6. Spontaneous and reversible proton transfer between a conserved Glu in the active site and cytosine N3 at the transition state was observed in our simulations, revealing the chemical participation of this Glu residue in the catalytic mechanism. Subsequently, the β-elimination of the C5 proton utilizes as base an OH(-) derived from a conserved crystal water that is part of a proton wire water channel, and this syn β-elimination reaction is the rate-limiting step. Design of novel cytosine methylation inhibitors would be advanced by our structural and thermodynamic characterization of the reaction mechanism.