Insights into the reaction mechanism of 3-O-sulfotransferase through QM/MM calculations

Polysaccharide sulfotransferases: the identification of putative sequences and respective functional characterisation

The vast structural diversity of sulfated polysaccharides demands an equally diverse array of enzymes known as polysaccharide sulfotransferases (PSTs). PSTs are present across all kingdoms of life, including algae, fungi and archaea, and their sulfation pathways are relatively unexplored. Sulfated polysaccharides possess anti-inflammatory, anticoagulant and anti-cancer properties and have great therapeutic potential. Current identification of PSTs using Pfam has been predominantly focused on the identification of glycosaminoglycan (GAG) sulfotransferases because of their pivotal roles in cell communication, extracellular matrix formation and coagulation. As a result, our knowledge of non-GAG PSTs structure and function remains limited. The major sulfotransferase families, Sulfotransfer_1 and Sulfotransfer_2, display broad homology and should enable the capture of a wide assortment of sulfotransferases but are limited in non-GAG PST sequence annotation. In addition, sequence annotation is further restricted by the paucity of biochemical analyses of PSTs. There are now high-throughput and robust assays for sulfotransferases such as colorimetric PAPS (3'-phosphoadenosine 5'-phosphosulfate) coupled assays, Europium-based fluorescent probes for ratiometric PAP (3'-phosphoadenosine-5'-phosphate) detection, and NMR methods for activity and product analysis. These techniques provide real-time and direct measurements to enhance the functional annotation and subsequent analysis of sulfated polysaccharides across the tree of life to improve putative PST identification and characterisation of function. Improved annotation and biochemical analysis of PST sequences will enhance the utility of PSTs across biomedical and biotechnological sectors.

From Steroid and Drug Metabolism to Glycobiology, Using Sulfotransferase Structures to Understand and Tailor Function

Sulfotransferases are ubiquitous enzymes that transfer a sulfo group from the universal cofactor donor 3'-phosphoadenosine 5'-phosphosulfate to a broad range of acceptor substrates. In humans, the cytosolic sulfotransferases are involved in the sulfation of endogenous compounds such as steroids, neurotransmitters, hormones, and bile acids as well as xenobiotics including drugs, toxins, and environmental chemicals. The Golgi associated membrane-bound sulfotransferases are involved in post-translational modification of macromolecules from glycosaminoglycans to proteins. The sulfation of small molecules can have profound biologic effects on the functionality of the acceptor, including activation, deactivation, or enhanced metabolism and elimination. Sulfation of macromolecules has been shown to regulate a number of physiologic and pathophysiological pathways by enhancing binding affinity to regulatory proteins or binding partners. Over the last 25 years, crystal structures of these enzymes have provided a wealth of information on the mechanisms of this process and the specificity of these enzymes. This review will focus on the general commonalities of the sulfotransferases, from enzyme structure to catalytic mechanism as well as providing examples into how structural information is being used to either design drugs that inhibit sulfotransferases or to modify the enzymes to improve drug synthesis. SIGNIFICANCE STATEMENT: This manuscript honors Dr. Masahiko Negishi's contribution to the understanding of sulfotransferase mechanism, specificity, and roles in biology by analyzing the crystal structures that have been solved over the last 25 years.

QM/MM study on enzymatic mechanism in sinigrin biosynthesis

As the major and abundant type of glucosinolates (GL) in plants, sinigrin has potential functions in promoting health and insect defense. The final step in the biosynthesis of sinigrin core structure is highly representative in GL compounds, which corresponds to the process from 3-methylthiopropyl ds-GL to 3-methylthiopropyl GL catalyzed by sulfotransferase (SOT). However, due to the lack of the crystallographic structure of SOT complexed with the 3-methylthiopropyl GL, little is known about this sulfonation process. Fortunately, the crystal structure of SOT 18 from Arabidopsis thaliana (AtSOT18) containing the substance (sinigrin) similar to 3-methylthiopropyl GL has been determined. To understand the enzymatic mechanism, we employed molecular dynamics (MD) simulation and quantum mechanics combined with molecular mechanics (QM/MM) methods to study the conversion from ds-sinigrin to sinigrin catalyzed by AtSOT18. The calculated results demonstrate that the reaction occurs through a concerted dissociative mechanism. Moreover, Lys93, Thr96, Thr97, Tyr130, His155, and two enzyme peptide chains (Pro92-Lys93 and Gln95-Thr96-Thr97) play a role in positioning the substrates and promoting the catalytic reaction by stabilizing the transition state geometry. Particularly, His155 acts as a catalytic base while Lys93 acts as a catalytic acid in the reaction process. The presently proposed concerted dissociative mechanism explains the role of AtSOT18 in sinigrin biosynthesis, and could be instructive for the study of GL biosynthesis catalyzed by other SOTs.

Heparan Sulfate Proteoglycans Biosynthesis and Post Synthesis Mechanisms Combine Few Enzymes and Few Core Proteins to Generate Extensive Structural and Functional Diversity

Glycosylation is a common and widespread post-translational modification that affects a large majority of proteins. Of these, a small minority, about 20, are specifically modified by the addition of heparan sulfate, a linear polysaccharide from the glycosaminoglycan family. The resulting molecules, heparan sulfate proteoglycans, nevertheless play a fundamental role in most biological functions by interacting with a myriad of proteins. This large functional repertoire stems from the ubiquitous presence of these molecules within the tissue and a tremendous structural variety of the heparan sulfate chains, generated through both biosynthesis and post synthesis mechanisms. The present review focusses on how proteoglycans are “gagosylated” and acquire structural complexity through the concerted action of Golgi-localized biosynthesis enzymes and extracellular modifying enzymes. It examines, in particular, the possibility that these enzymes form complexes of different modes of organization, leading to the synthesis of various oligosaccharide sequences.

Binding Mode Prediction and Identification of New Lead Compounds from Natural Products as 3-OST Enzyme Inhibitors

Background and Introduction:

The availability of antiviral medicines for the treatment of viral diseases is limited, hence the discovery of novel bioactive molecules is required. The present investigation has been carried out to develop novel 3-O-sulfotransferase enzyme inhibitors to treat viral diseases.

Method:

Virtual screening study (QSAR, docking and pharmacophore analysis) and binding mode analysis have been performed on a dataset collected from the literature (synthetic and natural compounds).

Results:

The docking studies showed that Glu184, His186, Lys215 and Lys368 residues established the most important hydrogen bonding with several hit compounds. The QSAR results explained that the presence of electronegative atoms/groups in the aromatic or heteroaromatic rings confer increased activity. Furthermore, the flexibility and the aromatic rings with less polar groups have better activity than the compounds connected to purine rings. Finally, the structurebased pharmacophore studies illustrated that the ligand has many polar interaction sites, and the projected acceptor and donor groups in the molecules make a significant contribution to the pharmacophore model building.

Conclusion:

These studies identified two compounds, Phomoidride B and Barceloneic acid A, as potential 3-OST inhibitors.

Revisiting the catalytic mechanism of Mo-Cu carbon monoxide dehydrogenase using QM/MM and DFT calculations

Previous density functional theory (DFT) studies have shown that the release of the produced carbon dioxide (CO2) from an active-site cluster is a thermodynamically or kinetically difficult step in the enzymatic carbon monoxide (CO) oxidation catalyzed by Mo-Cu carbon monoxide dehydrogenase (Mo-Cu CODH). To better understand the effect of the protein environment on this difficult CO2 release step as well as other reaction steps, we applied hybrid quantum mechanics and molecular mechanics (QM/MM) calculations to the Mo-Cu CODH enzyme. The results show that in the first step, the equatorial Mo[double bond, length as m-dash]O group in the active-site cluster attacks the nearby CO molecule bound to the Cu site. Afterward, a stable thiocarbonate intermediate is formed in which the CO2 molecule is embedded and the copper-S(μ-sulfido) bond is broken. A free CO2 molecule, i.e., the final product, is then released from the active-site cluster, not directly from the thiocarbonate intermediate but via a previously formed intermediate that also contains CO2 but retains the Cu-S(μ-sulfido) bond. In contrast to the previous DFT results, the calculated barrier for this process was low in our QM/MM calculations. An additional QM/MM analysis of the barrier height showed that the effect of the protein environment on this barrier lowering is not very large. We found that the reason for the low barrier obtained by QM/MM is that the barrier for CO2 release is already not high at the DFT level. These results allow us to conclude that the CO oxidation reaction passes through the formation of a thiocarbonate intermediate, and that the subsequent CO2 release is kinetically not difficult. Nevertheless, the protein environment has an important role to play in making the latter process thermodynamically favored. No low-barrier pathway for the product release could be obtained for the reaction of n-butylisocyanide, which is consistent with the experimental fact that n-butylisocyanide inhibits Mo-Cu CODH.

Mechanistic Insights on Human Phosphoglucomutase Revealed by Transition Path Sampling and Molecular Dynamics Calculations

Human α-phosphoglucomutase 1 (α-PGM) catalyzes the isomerization of glucose-1-phosphate into glucose-6-phosphate (G6P) through two sequential phosphoryl transfer steps with a glucose-1,6-bisphosphate (G16P) intermediate. Given that the release of G6P in the gluconeogenesis raises the glucose output levels, α-PGM represents a tempting pharmacological target for type 2 diabetes. Here, we provide the first theoretical study of the catalytic mechanism of human α-PGM. We performed transition-path sampling simulations to unveil the atomic details of the two catalytic chemical steps, which could be key for developing transition state (TS) analogue molecules with inhibitory properties. Our calculations revealed that both steps proceed through a concerted S_N 2-like mechanism, with a loose metaphosphate-like TS. Even though experimental data suggests that the two steps are identical, we observed noticeable differences: 1) the transition state ensemble has a well-defined TS region and a late TS for the second step, and 2) larger coordinated protein motions are required to reach the TS of the second step. We have identified key residues (Arg23, Ser117, His118, Lys389), and the Mg²⁺ ion that contribute in different ways to the reaction coordinate. Accelerated molecular dynamics simulations suggest that the G16P intermediate may reorient without leaving the enzymatic binding pocket, through significant conformational rearrangements of the G16P and of specific loop regions of the human α-PGM.

What regulates the catalytic activities in AGE catalysis? An answer from quantum mechanics/molecular mechanics simulations

The AGE superfamily (AGEs) is made up of kinds of isomerase which are very important both physiologically and industrially. One of the most intriguing aspects of AGEs has to do with the mechanism that regulates their activities in single conserved active pocket. In order to clarify the relationship among single conserved active pocket and two activities in AGEs, results for the epimerization activity catalyzed by RaCE and the isomerization activity catalyzed by SeYihS were obtained by using QM/MM umbrella sampling simulations and 2D-FES calculations. Our results show that both of them have similar enzyme-substrate combination mode for inner pyranose ring in single conserved active pocket even though they have different substrate specificity. This means that the pathways of ring opening catalyzed by them are similar. However, one non-conserved residue (Leu183 in RaCE, Met175 in SeYihS) in the active site, which has different steric hindrance, causes a small but effective change in the direction of ring opening in stage 1. And then this change will induce a fundamentally different catalytic activity for RaCE and SeYihS in stage 2. Our results give a novel viewpoint about the regulatory mechanism between CE and YihS in AGEs, and may be helpful for further experiments of rational enzyme design based on the (α/α)₆-barrel basic scaffold.