scholarly journals Conserved Conformational Hierarchy across Functionally Divergent Glycosyltransferases of the GT-B Structural Superfamily as Determined from Microsecond Molecular Dynamics

2021 ◽  
Vol 22 (9) ◽  
pp. 4619
Author(s):  
Carlos A. Ramirez-Mondragon ◽  
Megin E. Nguyen ◽  
Jozafina Milicaj ◽  
Bakar A. Hassan ◽  
Frank J. Tucci ◽  
...  

It has long been understood that some proteins undergo conformational transitions en route to the Michaelis Complex to allow chemistry. Examination of crystal structures of glycosyltransferase enzymes in the GT-B structural class reveals that the presence of ligand in the active site triggers an open-to-closed conformation transition, necessary for their catalytic functions. Herein, we describe microsecond molecular dynamics simulations of two distantly related glycosyltransferases that are part of the GT-B structural superfamily, HepI and GtfA. Simulations were performed using the open and closed conformations of these unbound proteins, respectively, and we sought to identify the major dynamical modes and communication networks that interconnect the open and closed structures. We provide the first reported evidence within the scope of our simulation parameters that the interconversion between open and closed conformations is a hierarchical multistep process which can be a conserved feature of enzymes of the same structural superfamily. Each of these motions involves of a collection of smaller molecular reorientations distributed across both domains, highlighting the complexities of protein dynamic involved in the interconversion process. Additionally, dynamic cross-correlation analysis was employed to explore the potential effect of distal residues on the catalytic efficiency of HepI. Multiple distal nonionizable residues of the C-terminal domain exhibit motions anticorrelated to positively charged residues in the active site in the N-terminal domain involved in substrate binding. Mutations of these residues resulted in a reduction in negatively correlated motions and an altered enzymatic efficiency that is dominated by lower Km values with kcat effectively unchanged. The findings suggest that residues with opposing conformational motions involved in the opening and closing of the bidomain HepI protein can allosterically alter the population and conformation of the “closed” state, essential to the formation of the Michaelis complex. The stabilization effects of these mutations likely equally influence the energetics of both the ground state and the transition state of the catalytic reaction, leading to the unaltered kcat. Our study provides new insights into the role of conformational dynamics in glycosyltransferase’s function and new modality to modulate enzymatic efficiency.

2020 ◽  
Author(s):  
Carlos A. Ramirez-Mondragon ◽  
Megin E. Nguyen ◽  
Jozafina Milicaj ◽  
Frank J. Tucci ◽  
Ramaiah Muthyala ◽  
...  

AbstractIt has long been understood that some proteins to undergo conformational transitions enroute to the Michaelis Complex to allow chemistry. Examination of crystal structures of glycosyltransferase enzymes in the GT-B structural class reveals that the presence of ligand in the active site is necessary for the protein to crystalize in the closed conformation. Herein we describe microsecond molecular dynamics simulations of two evolutionarily unrelated glycosytransferases, HepI and GtfA. Simulations were performed using these proteins in the open and closed conformations, (respectively,) and we sought to identify the major dynamical modes and communication networks which allow conformational transition between the open and closed structures. We provide the first reported evidence (within the scope of our experimental parameters) that conformational hierarchy/directionality of the interconversion between open and closed conformations is a conserved feature of enzymes of the same structural superfamily. Additionally, residues previously identified to be important for substrate binding in HepI were shown to have strong negative correlations with non-ionizable residues distal to the active site. Mutagenesis of these residues produced mutants with altered enzymatic efficiency exhibiting lower Km values, while the kcat is effectively unchanged. The negatively correlated motions of these residues are important for substrate binding and forming the Michaelis complex, without impacting the activation barrier for catalysis. This suggests that in the bi-domain HepI, the enzyme dynamics did not impact the transition state stabilization or chemistry, but rather earlier steps along the reaction coordinate, leading to the reorganization of the active site electrostatic environment required for catalysis.


Author(s):  
Balaji Selvam ◽  
Ya-Chi Yu ◽  
Liqing Chen ◽  
Diwakar Shukla

<p>The SWEET family belongs to a class of transporters in plants that undergoes large conformational changes to facilitate transport of sugar molecules across the cell membrane. However, the structures of their functionally relevant conformational states in the transport cycle have not been reported. In this study, we have characterized the conformational dynamics and complete transport cycle of glucose in OsSWEET2b transporter using extensive molecular dynamics simulations. Using Markov state models, we estimated the free energy barrier associated with different states as well as 1 for the glucose the transport mechanism. SWEETs undergoes structural transition to outward-facing (OF), Occluded (OC) and inward-facing (IF) and strongly support alternate access transport mechanism. The glucose diffuses freely from outside to inside the cell without causing major conformational changes which means that the conformations of glucose unbound and bound snapshots are exactly same for OF, OC and IF states. We identified a network of hydrophobic core residues at the center of the transporter that restricts the glucose entry to the cytoplasmic side and act as an intracellular hydrophobic gate. The mechanistic predictions from molecular dynamics simulations are validated using site-directed mutagenesis experiments. Our simulation also revealed hourglass like intermediate states making the pore radius narrower at the center. This work provides new fundamental insights into how substrate-transporter interactions actively change the free energy landscape of the transport cycle to facilitate enhanced transport activity.</p>


2020 ◽  
Vol 16 (4) ◽  
pp. 451-459 ◽  
Author(s):  
Fortunatus C. Ezebuo ◽  
Ikemefuna C. Uzochukwu

Background: Sulfotransferase family comprises key enzymes involved in drug metabolism. Oxamniquine is a pro-drug converted into its active form by schistosomal sulfotransferase. The conformational dynamics of side-chain amino acid residues at the binding site of schistosomal sulfotransferase towards activation of oxamniquine has not received attention. Objective: The study investigated the conformational dynamics of binding site residues in free and oxamniquine bound schistosomal sulfotransferase systems and their contribution to the mechanism of oxamniquine activation by schistosomal sulfotransferase using molecular dynamics simulations and binding energy calculations. Methods: Schistosomal sulfotransferase was obtained from Protein Data Bank and both the free and oxamniquine bound forms were subjected to molecular dynamics simulations using GROMACS-4.5.5 after modeling it’s missing amino acid residues with SWISS-MODEL. Amino acid residues at its binding site for oxamniquine was determined and used for Principal Component Analysis and calculations of side-chain dihedrals. In addition, binding energy of the oxamniquine bound system was calculated using g_MMPBSA. Results: The results showed that binding site amino acid residues in free and oxamniquine bound sulfotransferase sampled different conformational space involving several rotameric states. Importantly, Phe45, Ile145 and Leu241 generated newly induced conformations, whereas Phe41 exhibited shift in equilibrium of its conformational distribution. In addition, the result showed binding energy of -130.091 ± 8.800 KJ/mol and Phe45 contributed -9.8576 KJ/mol. Conclusion: The results showed that schistosomal sulfotransferase binds oxamniquine by relying on hybrid mechanism of induced fit and conformational selection models. The findings offer new insight into sulfotransferase engineering and design of new drugs that target sulfotransferase.


2021 ◽  
Vol 22 (9) ◽  
pp. 4769
Author(s):  
Pablo Maturana ◽  
María S. Orellana ◽  
Sixto M. Herrera ◽  
Ignacio Martínez ◽  
Maximiliano Figueroa ◽  
...  

Agmatine is the product of the decarboxylation of L-arginine by the enzyme arginine decarboxylase. This amine has been attributed to neurotransmitter functions, anticonvulsant, anti-neurotoxic, and antidepressant in mammals and is a potential therapeutic agent for diseases such as Alzheimer’s, Parkinson’s, and cancer. Agmatinase enzyme hydrolyze agmatine into urea and putrescine, which belong to one of the pathways producing polyamines, essential for cell proliferation. Agmatinase from Escherichia coli (EcAGM) has been widely studied and kinetically characterized, described as highly specific for agmatine. In this study, we analyze the amino acids involved in the high specificity of EcAGM, performing a series of mutations in two loops critical to the active-site entrance. Two structures in different space groups were solved by X-ray crystallography, one at low resolution (3.2 Å), including a guanidine group; and other at high resolution (1.8 Å) which presents urea and agmatine in the active site. These structures made it possible to understand the interface interactions between subunits that allow the hexameric state and postulate a catalytic mechanism according to the Mn2+ and urea/guanidine binding site. Molecular dynamics simulations evaluated the conformational dynamics of EcAGM and residues participating in non-binding interactions. Simulations showed the high dynamics of loops of the active site entrance and evidenced the relevance of Trp68, located in the adjacent subunit, to stabilize the amino group of agmatine by cation-pi interaction. These results allow to have a structural view of the best-kinetic characterized agmatinase in literature up to now.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Nur Suhanawati Ashaari ◽  
Mohd Hairul Ab. Rahim ◽  
Suriana Sabri ◽  
Kok Song Lai ◽  
Adelene Ai-Lian Song ◽  
...  

AbstractLinalool and nerolidol are terpene alcohols that occur naturally in many aromatic plants and are commonly used in food and cosmetic industries as flavors and fragrances. In plants, linalool and nerolidol are biosynthesized as a result of respective linalool synthase and nerolidol synthase, or a single linalool/nerolidol synthase. In our previous work, we have isolated a linalool/nerolidol synthase (designated as PamTps1) from a local herbal plant, Plectranthus amboinicus, and successfully demonstrated the production of linalool and nerolidol in an Escherichia coli system. In this work, the biochemical properties of PamTps1 were analyzed, and its 3D homology model with the docking positions of its substrates, geranyl pyrophosphate (C10) and farnesyl pyrophosphate (C15) in the active site were constructed. PamTps1 exhibited the highest enzymatic activity at an optimal pH and temperature of 6.5 and 30 °C, respectively, and in the presence of 20 mM magnesium as a cofactor. The Michaelis–Menten constant (Km) and catalytic efficiency (kcat/Km) values of 16.72 ± 1.32 µM and 9.57 × 10–3 µM−1 s−1, respectively, showed that PamTps1 had a higher binding affinity and specificity for GPP instead of FPP as expected for a monoterpene synthase. The PamTps1 exhibits feature of a class I terpene synthase fold that made up of α-helices architecture with N-terminal domain and catalytic C-terminal domain. Nine aromatic residues (W268, Y272, Y299, F371, Y378, Y379, F447, Y517 and Y523) outlined the hydrophobic walls of the active site cavity, whilst residues from the RRx8W motif, RxR motif, H-α1 and J-K loops formed the active site lid that shielded the highly reactive carbocationic intermediates from the solvents. The dual substrates use by PamTps1 was hypothesized to be possible due to the architecture and residues lining the catalytic site that can accommodate larger substrate (FPP) as demonstrated by the protein modelling and docking analysis. This model serves as a first glimpse into the structural insights of the PamTps1 catalytic active site as a multi-substrate linalool/nerolidol synthase.


2013 ◽  
Vol 12 (08) ◽  
pp. 1341002 ◽  
Author(s):  
XIN ZHANG ◽  
MING LEI

The deamination process of isoxanthopterin catalyzed by isoxanthopterin deaminase was determined using the combined QM(PM3)/MM molecular dynamics simulations. In this paper, the updated PM3 parameters were employed for zinc ions and the initial model was built up based on the crystal structure. Proton transfer and following steps have been investigated in two paths: Asp336 and His285 serve as the proton shuttle, respectively. Our simulations showed that His285 is more effective than Aap336 in proton transfer for deamination of isoxanthopterin. As hydrogen bonds between the substrate and surrounding residues play a key role in nucleophilic attack, we suggested mutating Thr195 to glutamic acid, which could enhance the hydrogen bonds and help isoxanthopterin get close to the active site. The simulations which change the substrate to pterin 6-carboxylate also performed for comparison. Our results provide reference for understanding of the mechanism of deaminase and for enhancing the deamination rate of isoxanthopterin deaminase.


2018 ◽  
Vol 115 (52) ◽  
pp. E12192-E12200 ◽  
Author(s):  
Haoran Yu ◽  
Paul A. Dalby

The directed evolution of enzymes for improved activity or substrate specificity commonly leads to a trade-off in stability. We have identified an activity–stability trade-off and a loss in unfolding cooperativity for a variant (3M) of Escherichia coli transketolase (TK) engineered to accept aromatic substrates. Molecular dynamics simulations of 3M revealed increased flexibility in several interconnected active-site regions that also form part of the dimer interface. Mutating the newly flexible active-site residues to regain stability risked losing the new activity. We hypothesized that stabilizing mutations could be targeted to residues outside of the active site, whose dynamics were correlated with the newly flexible active-site residues. We previously stabilized WT TK by targeting mutations to highly flexible regions. These regions were much less flexible in 3M and would not have been selected a priori as targets using the same strategy based on flexibility alone. However, their dynamics were highly correlated with the newly flexible active-site regions of 3M. Introducing the previous mutations into 3M reestablished the WT level of stability and unfolding cooperativity, giving a 10.8-fold improved half-life at 55 °C, and increased midpoint and aggregation onset temperatures by 3 °C and 4.3 °C, respectively. Even the activity toward aromatic aldehydes increased up to threefold. Molecular dynamics simulations confirmed that the mutations rigidified the active-site via the correlated network. This work provides insights into the impact of rigidifying mutations within highly correlated dynamic networks that could also be useful for developing improved computational protein engineering strategies.


2019 ◽  
Author(s):  
Christian Curado-Carballada ◽  
Ferran Feixas ◽  
Sílvia Osuna

<p><b> </b><i>Aspergillus niger </i>Monoamine Oxidase (MAO-N) is a homodimeric enzyme responsible for the oxidation of amines into the corresponding imine. Laboratory evolved variants of MAO-N in combination with a non-selective chemical reductant represents a powerful strategy for the deracemisation of chiral amine mixtures and, thus, is of interest for obtaining chiral amine building blocks. MAO-N presents a rich conformational dynamics with a flexible ß-hairpin region that can adopt closed, partially closed and open states. Despite the ß-hairpin conformational dynamics is altered along the laboratory evolutionary pathway of MAO-N, the connection between the ß-hairpin conformational dynamics and active site catalysis still remains unclear. In this work, we use accelerated molecular dynamics to elucidate the potential interplay between the ß-hairpin conformational dynamics and catalytic activity in MAO-N wild type and its evolved D5 variant. Our study reveals a delicate communication between both MAO-N subunits that impacts the active site architecture, and thus its catalytic efficiency. In both MAO-N WT and the laboratory evolved D5 variant, the ß-hairpin conformation in one of the monomers affects the productive binding of the substrate in the active site of the other subunit. However, both MAO-N WT and D5 variants show a quite different behaviour due to the distal mutations introduced experimentally with Directed Evolution. </p>


Sign in / Sign up

Export Citation Format

Share Document