Substrate Docking Research Articles

The crystal structure of Grindelia robusta 7,13-copalyl diphosphate synthase reveals active site features controlling catalytic specificity

Diterpenoid natural products serve critical functions in plant development and ecological adaptation and many diterpenoids have economic value as bioproducts. The family of class II diterpene synthases catalyzes the committed reactions in diterpenoid biosynthesis, converting a common geranylgeranyl diphosphate precursor into different bicyclic prenyl diphosphate scaffolds. Enzymatic rearrangement and modification of these precursors generates the diversity of bioactive diterpenoids. We report the crystal structure of Grindelia robusta 7,13-copalyl diphosphate synthase, GrTPS2, at 2.1 Å of resolution. GrTPS2 catalyzes the committed reaction in the biosynthesis of grindelic acid, which represents the signature metabolite in species of gumweed (Grindelia spp., Asteraceae). Grindelic acid has been explored as a potential source for drug leads and biofuel production. The GrTPS2 crystal structure adopts the conserved three-domain fold of class II diterpene synthases featuring a functional active site in the γβ-domain and a vestigial ɑ-domain. Substrate docking into the active site of the GrTPS2 apo protein structure predicted catalytic amino acids. Biochemical characterization of protein variants identified residues with impact on enzyme activity and catalytic specificity. Specifically, mutagenesis of Y457 provided mechanistic insight into the position-specific deprotonation of the intermediary carbocation to form the characteristic 7,13 double bond of 7,13-copalyl diphosphate.

Rigidifying Flexible Sites: A Promising Strategy to Improve Thermostability of Lysophospholipase From Pyrococcus abyssi.

High thermostability of the enzymes is one of the distinguishing characteristics that increase their industrial utility. In the current research work, rigidifying the flexible amino acid residues of a lysophospholipase (Pa-LPL) from Pyrococcus abyssi was used as a protein engineering approach to improve its thermostability. A truncated variant of Pa-LPL (t-LPL∆12) was constructed by trimming its 12 amino acid residues (50-61) through overlap extension PCR. The truncated enzyme worked optimally at 65°C and pH 6.5 with remarkable thermostability at 65°C-85°C. In comparison to wild-type Pa-LPL, 5.8 and 1.2-fold increase in half-life (t1/2) of t-LPL∆12 was observed at 65 (optimum temperature) and 95°C, respectively. The activity of t-LPL∆12 was stimulated by 1 mM Cu2+ followed by Ca2+, Ni2+, Co2+, and Mg2+. Both substrate docking and experimental results indicated that the truncated enzyme could hydrolyze a variety of p-nitrophenyl esters. Km, Vmax, and Kcat for enzymatic hydrolysis of p-nitrophenyl butyrate were calculated to be 1 ± 0.087 mM, 1456 ± 36.474 U/mg, and 1.397 × 1011 min-1, respectively. In short, broad substrate specificity and thermostability of t-LPL∆12 are some of the distinctive features that make it an ideal candidate for degumming of vegetable oils.

Structural Analysis of Mammalian Sialic Acid Esterase

Sialic acid esterase (SIAE) catalyzes the removal of O-acetyl groups from sialic acids found on cell surface glycoproteins to regulate cellular processes such as B cell receptor signalling and apoptosis. Loss-of-function mutations in SIAE are associated with several common autoimmune diseases including Crohn’s, ulcerative colitis, and arthritis. To gain a better understanding of the function and regulation of this protein, we determined crystal structures of SIAE from three mammalian homologs, including an acetate bound structure. The structures reveal that the catalytic domain adopts the fold of the SGNH hydrolase superfamily. The active site is composed of a catalytic dyad, as opposed to the previously reported catalytic triad. Attempts to determine a substrate-bound structure yielded only the hydrolyzed product acetate in the active site. Rigid docking of complete substrates followed by molecular dynamics simulations revealed that the active site does not form specific interactions with substrates, rather it appears to be broadly specific to accept sialoglycans with diverse modifications. Based on the acetate bound structure, a catalytic mechanism is proposed. Structural mapping of disease mutations reveals that most are located on the surface of the enzyme and would only cause minor disruptions to the protein fold, suggesting that these mutations likely affect binding to other factors. These results improve our understanding of SIAE biology and may aid in the development of therapies for autoimmune diseases and cancer.

Docking interactions determine substrate specificity of members of a widespread family of protein phosphatases

How protein phosphatases achieve specificity for their substrates is a major outstanding question. PPM family serine/threonine phosphatases are widespread in bacteria and eukaryotes, where they dephosphorylate target proteins with a high degree of specificity. In bacteria, PPM phosphatases control diverse transcriptional responses by dephosphorylating anti-anti-sigma factors of the STAS domain family, exemplified by B. subtilis phosphatases SpoIIE, which controls cell-fate during endospore formation, and RsbU, which initiates the General Stress Response. Using a combination of forward genetics, biochemical reconstitution, and AlphaFold2 structure prediction, we identified a conserved, tripartite substrate docking interface comprised of three variable loops on the surface of the PPM phosphatase domains of SpoIIE and RsbU that recognize the three-dimensional structure of the substrate protein. Non-conserved amino acids in these loops facilitate the accommodation of the cognate substrate and prevent dephosphorylation of the non-cognate substrate. Together, single-amino acid substitutions in these three elements cause an over five-hundred fold change in specificity. Our data additionally suggest that substrate-docking interactions regulate phosphatase specificity through a conserved allosteric switch element that controls the catalytic efficiency of the phosphatase by positioning the metal cofactor and substrate. We hypothesize that this is a generalizable mechanistic model for PPM family phosphatase substrate specificity. Importantly, the substrate docking interface with the phosphatase is only partially overlapping with the much more extensive interface with the upstream kinase, suggesting the possibility that kinase and phosphatase specificity evolved independently.

FAM122A ensures cell cycle interphase progression and checkpoint control by inhibiting B55α/PP2A through helical motifs

The Ser/Thr protein phosphatase 2 A (PP2A) regulates the dephosphorylation of many phosphoproteins. Substrate recognition are mediated by B regulatory subunits. Here, we report the identification of a substrate conserved motif [RK]-V-x-x-[VI]-R in FAM122A, an inhibitor of B55α/PP2A. This motif is necessary for FAM122A binding to B55α, and computational structure prediction suggests the motif, which is helical, blocks substrate docking to the same site. In this model, FAM122A also spatially constrains substrate access by occluding the catalytic subunit. Consistently, FAM122A functions as a competitive inhibitor as it prevents substrate binding and dephosphorylation of CDK substrates by B55α/PP2A in cell lysates. FAM122A deficiency in human cell lines reduces the proliferation rate, cell cycle progression, and hinders G1/S and intra-S phase cell cycle checkpoints. FAM122A-KO in HEK293 cells attenuates CHK1 and CHK2 activation in response to replication stress. Overall, these data strongly suggest that FAM122A is a short helical motif (SHeM)-dependent, substrate-competitive inhibitor of B55α/PP2A that suppresses multiple functions of B55α in the DNA damage response and in timely progression through the cell cycle interphase.

The biochemical characterization of a TatD nuclease from Thermus thermophilus

Nucleases play pivotal roles in DNA repair and apoptosis. Moreover, they have various applications in biotechnology and industry. Among nucleases, TatD has been characterized as an exonuclease with various biological functions in different organisms. Here, we biochemically characterized the potential TatD nuclease from Thermus thermophilus. The tatD gene from T. thermophilus was cloned, then the recombinant TatD nuclease was expressed and purified. Our results revealed that the TthTatD nuclease could degrade both single-stranded and double-stranded DNA, and its activity is dependent on the divalent metal ions Mg2+ and Mn2+. Remarkably, the activity of TthTatD nuclease is highest at 37 °C and decreases with increasing temperature. TthTatD is not a thermostable enzyme, even though it is from a thermophilic bacterium. Based on the sequence similarity and molecular docking of the DNA substrate into the modeled TthTatD structure, several key conserved residues were identified and their roles were confirmed by analyzing the enzymatic activities of the site-directed mutants. The residues E86 and H149 play key roles in binding metal ions, residues R124/K126 and K211/R212 had a critical role in binding DNA substrate. Our results confirm the enzymatic properties of TthTatD and provide a primary basis for its possible application in biotechnology.

Enhancing Paenibacillus sp. Cold-Active Acetyl Xylan Esterase Activity through Semi-Rational Protein Engineering

Interest in protein engineering for the enzymatic production of valuable products, such as pharmaceutical compounds and biofuels, is growing rapidly. The cold-active acetyl xylan esterase from Paenibacillus sp. (PbAcE) presents unusually broad substrate specificity. Here, we engineered a hydrophobic substrate-binding pocket to enable the accommodation of relatively large alcohol substrates, such as linalyl acetate and α-terpinyl acetate. To identify candidate residues for engineering, we performed covalent docking of substrates to the Ser185 active site using the HCovDock program. Functional hotspots were analyzed using HotSpot Wizard 3.1. Lys91, His93, and Tyr182 were selected for site-saturation mutagenesis (SSM). After generating the SSM mutant library, a qualitative colorimetric assay was conducted to identify positive mutants. Three, two, and five single mutants were selected for Lys91, His93, and Tyr182, respectively. The best single mutants were then sequentially combined to generate double and triple mutants. Single mutants exhibited a 10–30% increase in activity compared to that of wild-type PbAcE, while no significant synergistic improvements were observed in the double and triple mutants. The increase in activity against both linalyl acetate and α-terpinyl acetate was similar. Mutation did not affect the acetyl binding and catalysis. Further research on the acetyl binding pocket will provide insights into substrate specificity and aid in efficient biocatalyst development for industrial applications.

Increasing Multienzyme Cascade Efficiency and Stability of MOF via Partitioning Immobilization.

Enhancing the stability of multienzyme cascade reactions in metal-organic frameworks (MOFs) is a challenging task in the fields of biotechnology and chemistry. However, addressing this challenge could yield far-reaching benefits across the application range in the biomedical, food, and environmental sectors. In this study, multienzyme partitioning immobilization that sequentially immobilizes cascade enzymes with hierarchical MOFs is proposed to reduce substrate diffusion resistance. Conversion results of ginsenosides indicate that this strategy improves the cascade efficiency up to 1.26 times. The substrate diffusion model is used to investigate the dual-interenzyme mass transfer behavior of substrates in the restricted domain space and evaluate the substrate channeling effect under partitioning immobilization. Molecular docking and kinetic simulations reveal that the MOFs effectively limit the conformational changes of cascade enzymes at high temperatures and in organic solvents while maintaining a large pocket of active centers. This phenomenon increased efficient substrate docking to the enzyme molecules, further optimizing cascade efficiency. The results of the immobilization of GOX and horseradish peroxidase as model enzymes indicate that the partitioned MOF immobilization strategy could be used for universal adaptation of cascade enzymes.

Regiodivergent biosynthesis of bridged bicyclononanes

Medicinal compounds from plants include bicyclo[3.3.1]nonane derivatives, the majority of which are polycyclic polyprenylated acylphloroglucinols (PPAPs). Prototype molecules are hyperforin, the antidepressant constituent of St. John’s wort, and garcinol, a potential anticancer compound. Their complex structures have inspired innovative chemical syntheses, however, their biosynthesis in plants is still enigmatic. PPAPs are divided into two subclasses, named type A and B. Here we identify both types in Hypericum sampsonii plants and isolate two enzymes that regiodivergently convert a common precursor to pivotal type A and B products. Molecular modelling and substrate docking studies reveal inverted substrate binding modes in the two active site cavities. We identify amino acids that stabilize these alternative binding scenarios and use reciprocal mutagenesis to interconvert the enzymatic activities. Our studies elucidate the unique biochemistry that yields type A and B bicyclo[3.3.1]nonane cores in plants, thereby providing key building blocks for biotechnological efforts to sustainably produce these complex compounds for preclinical development.

Enzymatic one-step synthesis of natural 2-pyrones and new-to-nature derivatives from coenzyme A esters

The 2-pyrone moiety is present in a wide range of structurally diverse natural products with various biological activities. The plant biosynthetic routes towards these compounds mainly depend on the activity of either type III polyketide synthase-like 2-pyrone synthases or hydroxylating 2-oxoglutarate dependent dioxygenases. In the present study, the substrate specificity of these enzymes is investigated by a systematic screening using both natural and artificial substrates with the aims of efficiently forming (new) products and understanding the underlying catalytic mechanisms. In this framework, we focused on the in vitro functional characterization of a 2-pyrone synthase Gh2PS2 from Gerbera x hybrida and two dioxygenases AtF6’H1 and AtF6’H2 from Arabidopsis thaliana using a set of twenty aromatic and aliphatic CoA esters as substrates. UHPLC-ESI-HRMSn based analyses of reaction intermediates and products revealed a broad substrate specificity of the enzymes, enabling the facile "green" synthesis of this important class of natural products and derivatives in a one-step/one-pot reaction in aqueous environment without the need for halogenated or metal reagents and protective groups. Using protein modeling and substrate docking we identified amino acid residues that seem to be important for the observed product scope.

Crystal structure of activating sulfotransferase SgdX2 involved in biosynthesis of secondary metabolite sungeidine

Microorganisms synthesize a plethora of complex secondary metabolites, many of which are beneficial to human health, such as anticancer agents and antibiotics. Among these, the Sungeidines are a distinct class of secondary metabolites known for their bulky and intricate structures. They are produced by a specific biosynthetic gene cluster within the genome of the soil-dwelling actinomycete Micromonospora sp. MD118. A notable enzyme in the Sungeidine biosynthetic pathway is the activating sulfotransferase SgdX2. In this pathway, SgdX2 mediates a key sulfation step, after which the product undergoes spontaneous dehydration to yield a Sungeidine compound. To delineate the structural basis for SgdX2's substrate recognition and catalytic action, we have determined the crystal structure of SgdX2 in complex with its sulfate donor product, 3′-phosphoadenosine 5′-phosphate (PAP), at a resolution of 1.6 Å. Although SgdX2 presents a compact overall structure, its core elements are conserved among other activating sulfotransferases. Our structural analysis reveals a unique substrate-binding pocket that accommodates bulky, complex substrates, suggesting a specialized adaptation for Sungeidine synthesis. Moreover, we have constructed a substrate docking model that provides insights into the molecular interactions between SgdX2 and Sungeidine F, enhancing our understanding of the enzyme's specificity and catalytic mechanism. The model supports a general acid-base catalysis mechanism, akin to other sulfotransferases, and underscores the minor role of disordered regions in substrate recognition. This integrative study of crystallography and computational modeling advances our knowledge of microbial secondary metabolite biosynthesis and may facilitate the development of novel biotechnological applications.

CDK activity at the centrosome regulates the cell cycle

In human cells and yeast, an intact "hydrophobic patch" substrate docking site is needed for mitotic cyclin centrosomal localization. A hydrophobic patch mutant (HPM) of the fission yeast mitotic cyclin Cdc13 cannot enter mitosis, but whether this is due to defective centrosomal localization or defective cyclin-substrate docking more widely is unknown. Here, we show that artificially restoring Cdc13-HPM centrosomal localization promotes mitotic entry and increases CDK (cyclin-dependent kinase) substrate phosphorylation at the centrosome and in the cytoplasm. We also show that the S-phase B-cyclin hydrophobic patch is required for centrosomal localization but not for S phase. We propose that the hydrophobic patch is essential for mitosis due to its requirement for the local concentration of cyclin-CDK with CDK substrates and regulators at the centrosome. Our findings emphasize the central importance of the centrosome as a hub coordinating cell-cycle control and explain why the cyclin hydrophobic patch is essential for mitosis.

Molecular Modeling of the Multiple-Substrate Activity of the Human Recombinant Intra-Melanosomal Domain of Tyrosinase and Its OCA1B-Related Mutant Variant P406L.

Tyrosinase serves as the key enzyme in melanin biosynthesis, catalyzing the initial steps of the pathway, the hydroxylation of the amino acid L-tyrosine into L-3,4-dihydroxyphenylalanine (L-DOPA), followed by the subsequent oxidation of L-DOPA into dopaquinone (DQ), and it facilitates the conversion of 5,6-dihydroxyindole-2-carboxylic acid (DHICA) into 5,6-indolequinone-2-carboxylic acid (IQCA) and 5,6-dihydroxy indole (DHI) into indolequinone (IQ). Despite its versatile substrate capabilities, the precise mechanism underlying tyrosinase's multi-substrate activity remains unclear. Previously, we expressed, purified, and characterized the recombinant intra-melanosomal domain of human tyrosinase (rTyr). Here, we demonstrate that rTyr mimics native human tyrosinase's catalytic activities in vitro and in silico. Molecular docking and molecular dynamics (MD) simulations, based on rTyr's homology model, reveal variable durability and binding preferences among tyrosinase substrates and products. Analysis of root mean square deviation (RMSD) highlights the significance of conserved residues (E203, K334, F347, and V377), which exhibit flexibility during the ligands' binding. Additionally, in silico analysis demonstrated that the OCA1B-related P406L mutation in tyrosinase substantially influences substrate binding, as evidenced by the decreased number of stable ligand conformations. This correlation underscores the mutation's impact on substrate docking, which aligns with the observed reduction in rTyr activity. Our study highlights how rTyr dynamically adjusts its structure to accommodate diverse substrates and suggests a way to modulate rTyr ligand plasticity.

Adapting an acyl CoA ligase from Metallosphaera sedula for lactam formation by structure-guided protein engineering

The CoA ligase from Metallosphaera sedula (MsACL) can be used for the chemoenzymatic synthesis of amides from carboxylic acids. In this CoA-independent conversion, the enzyme catalyzes the adenylation of a carboxylic acid with the help of ATP, followed by the uncatalyzed cleavage of acyl-AMP by a nucleophilic amine to yield an amide. With ω-amino acids as substrates this reaction may result in formation of lactams, but unfortunately the substrate preference of the wild-type enzyme is rather limited. To allow structure-based protein engineering and expand the substrate scope of the enzyme, crystal structures of MsACL were solved in the thioesterification conformational state with AMP, CoA and with the reaction intermediate acetyl-AMP bound in the active site. Using substrate docking and by comparing the crystals structures and sequence of MsACL to those of related CoA ligases, mutations were predicted which increase the affinity in the carboxylic acid binding pocket for ω-amino acids. The resulting mutations transformed a non-active enzyme into an active enzyme for ε-caprolactam synthesis, highlighting the potential of the thermophilic CoA ligase for this synthetic and biotechnologically relevant reaction.

Regiodivergent Radical Termination for Intermolecular Biocatalytic C-C Bond Formation.

Radical hydrofunctionalizations of electronically unbiased dienes are challenging to render regioselective, because the products are nearly identical in energy. Here, we report two engineered FMN-dependent "ene"-reductases (EREDs) that catalyze regiodivergent hydroalkylations of cyclic and linear dienes. While previous studies focused exclusively on the stereoselectivity of alkene hydroalkylation, this work highlights that EREDs can control the regioselectivity of hydrogen atom transfer, providing a method for selectively preparing constitutional isomers that would be challenging to prepare using traditional synthetic methods. Engineering the ERED from Gluconabacter sp. (GluER) furnished a variant that favors the γ,δ-unsaturated ketone, while an engineered variant from a commercial ERED panel favors the δ,ε-unsaturated ketone. The effect of beneficial mutations has been investigated using substrate docking studies and the mechanism probed by isotope labeling experiments. A variety of α-bromo ketones can be coupled with cyclic and linear dienes. These interesting building blocks can also be further modified to generate difficult-to-access heterocyclic compounds.

Dirigent isoflavene-forming PsPTS2: 3D structure, stereochemical, and kinetic characterization comparison with pterocarpan-forming PsPTS1 homolog in pea

Pea phytoalexins (-)-maackiain and (+)-pisatin have opposite C6a/C11a configurations, but biosynthetically how this occurs is unknown. Pea dirigent-protein (DP) PsPTS2 generates 7,2'-dihydroxy-4',5'-methylenedioxyisoflav-3-ene (DMDIF), and stereoselectivity toward four possible 7,2'-dihydroxy-4',5'-methylenedioxyisoflavan-4-ol (DMDI) stereoisomers was investigated. Stereoisomer configurations were determined using NMR spectroscopy, electronic circular dichroism, and molecular orbital analyses. PsPTS2 efficiently converted cis-(3R,4R)-DMDI into DMDIF 20-fold faster than the trans-(3R,4S)-isomer. The 4R-configured substrate's near β-axial OH orientation significantly enhanced its leaving group abilities in generating A-ring mono-quinone methide (QM), whereas 4S-isomer's α-equatorial-OH was a poorer leaving group. Docking simulations indicated that the 4R-configured β-axial OH was closest to Asp51, whereas 4S-isomer's α-equatorial OH was further away. Neither cis-(3S,4S)- nor trans-(3S,4R)-DMDIs were substrates, even with the former having C3/C4 stereochemistry as in (+)-pisatin. PsPTS2 used cis-(3R,4R)-7,2'-dihydroxy-4'-methoxyisoflavan-4-ol [cis-(3R,4R)-DMI] and C3/C4 stereoisomers to give 2',7-dihydroxy-4'-methoxyisoflav-3-ene (DMIF). DP homologs may exist in licorice (Glycyrrhiza pallidiflora) and tree legume Bolusanthus speciosus, as DMIF occurs in both species. PsPTS1 utilized cis-(3R,4R)-DMDI to give (-)-maackiain 2200-fold more efficiently than with cis-(3R,4R)-DMI to give (-)-medicarpin. PsPTS1 also slowly converted trans-(3S,4R)-DMDI into (+)-maackiain, reflecting the better 4R configured OH leaving group. PsPTS2 and PsPTS1 provisionally provide the means to enable differing C6a and C11a configurations in (+)-pisatin and (-)-maackiain, via identical DP-engendered mono-QM bound intermediate generation, which PsPTS2 either re-aromatizes to give DMDIF or PsPTS1 intramolecularly cyclizes to afford (-)-maackiain. Substrate docking simulations using PsPTS2 and PsPTS1 indicate cis-(3R,4R)-DMDI binds in the anti-configuration in PsPTS2 to afford DMDIF, and the syn-configuration in PsPTS1 to give maackiain.

Characterization of the Baeyer-Villiger monooxygenase in the pathway of the bacterial pyrrolizidine alkaloids, legonmycins.

The Baeyer-Villiger monooxygenase (BVMO), LgnC, plays a crucial role in the biosynthesis of bacterial pyrrolizidine alkaloids, legonmycins. It processes bicyclic indolizidine substrates generated from the coordinative action of two non-ribosomal peptide synthetases (LgnB and LgnD) and the standalone type II thioesterase-like enzyme (LgnA). It has been demonstrated that the enzyme selectively inserts molecular oxygen into the carbon-carbon bond adjacent to the carbonyl group in legonindolizidines to form bicyclic 1,3-oxazepine carbamate intermediates. After ring opening and contraction, the most advanced products, prelegonmycins, are formed. However, factors controlling the final hydroxylation step and how the enzyme handles the substrates have remained elusive. In this study, we show that the final hydroxylation at the activated carbon of the electron-rich pyrrole system is attributed to either spontaneous oxidation or the action of an endogenous redox reagent. Substrate docking on the structural model of LgnC combined with site-directed mutagenesis allows the identification of several key amino acids that are essential for substrate/intermediate binding and a mechanism of LgnC-catalysed transformation is proposed.

Molecular docking analysis of a dermatan sulfate tetra-saccharide to human alpha-L-iduronidase.

Human alpha-L-iduronidase (IDUA) is a 653 amino acid protein involved in the sequential degradation of glycos-amino-glycans (GAG), heparan sulfate (HS), and dermatan sulfate (DS). Some variants in the IDUA gene produce a deficient enzyme that causes un-degraded DS and HS to accumulate in multiple tissues, leading to an organ dysfunction known as muco-poly-saccharidosis type I (MPS I). Molecular and catalytic activity assays of new or rare variants of IDUA do not predict the phenotype that a patient will develop. Therefore, it is of interest to describe the molecular docking analysis, to locate binding regions of DS to IDUA to better understand the effect of a variant on MPS I development. The results presented herein demonstrate the presence of a polar/acidic catalytic site and a basic region in the putative binding site of DS to IDUA. Further, synthetic substrate docking with the enzyme could help in the predictions of the MPS I phenotype.

Development of a Novel CYP3A4 Classifier Model via Site of Metabolism (SOM)-based Molecular Docking, Multivariate Analysis and Molecular Dynamics of Known Substrates and Inhibitors

CYP3A4 is a major hepatic enzyme essential for metabolizing diverse chemical entities. The development of CYP3A4 substrate and inhibitor classifier models is a valuable research strategy to prevent toxicokinetics. Currently, no molecular docking model is available to classify the substrates and inhibitors for a wide range of chemical entities. This study generated a CYP3A4 substrate-inhibitor classifier model using multivariate analysis of selected variables from site of metabolism (SOM)-based docking results and molecular descriptors. CYP3A4 SOM-based molecular docking experiment was performed on the substrates and inhibitors from the DrugBank database. The relevant descriptors were selected using the area under the receiving operating curve (AUROC) and boxplot analysis. The selected variables were used to generate a CYP3A4 substrate-inhibitor classifier model using multivariate analysis. Two complexes were selected for molecular dynamic simulations. A total of 326 substrates and 154 inhibitors were successfully docked. The SOM-based docking model predicted 78.5% SOM correctly of the substrates. The CDOCKER energy, improper dihedral energy, potential energy, initial RMS (root mean square) gradient, and RMS gradient demonstrated acceptable AUROC, sensitivity, specificity, and Youden’s index results. The selected variables were subjected to multivariate analysis, and the PLS-DA analysis classification model showed promising performance in predicting chemical entities’ behavior as substrates or inhibitors. This model classified the external samples correctly with over 70.0% prediction accuracy. Molecular dynamic simulations showed that the selected complexes obtained from the docking model produced stable complexes. The classifier model is applicable as a research tool for the early stages of drug discovery.

Role of Conformational Dynamics of Sulfotransferases SULT1A1 and SULT1A3 in Substrate Specificity

Sulfotransferases (SULTs) are phase II metabolizing enzymes catalyzing the sulfoconjugation from the co-factor 3′-Phosphoadenosine 5′-Phosphosulfate (PAPS) to a wide variety of endogenous compounds, drugs and natural products. Although SULT1A1 and SULT1A3 share 93% identity, SULT1A1, the most abundant SULT isoform in humans, exhibits a broad substrate range with specificity for small phenolic compounds, while SULT1A3 displays a high affinity toward monoamine neurotransmitters like dopamine. To elucidate the factors determining the substrate specificity of the SULT1 isoenzymes, we studied the dynamic behavior and structural specificities of SULT1A1 and SULT1A3 by using molecular dynamics (MD) simulations and ensemble docking of common and specific substrates of the two isoforms. Our results demonstrated that while SULT1A1 exhibits a relatively rigid structure by showing lower conformational flexibility except for the lip (loop L1), the loop L2 and the cap (L3) of SULT1A3 are extremely flexible. We identified protein residues strongly involved in the recognition of different substrates for the two isoforms. Our analyses indicated that being more specific and highly flexible, the structure of SULT1A3 has particularities in the binding site, which are crucial for its substrate selectivity.