Effects of Xylanase A double mutation on substrate specificity and structural dynamics.

While protein activity is traditionally studied with a major focus on the active site, the activity of enzymes has been hypothesized to be linked to the flexibility of adjacent regions, warranting more exploration into how the dynamics in these regions affects catalytic turnover. One such enzyme is Xylanase A (XylA), which cleaves hemicellulose xylan polymers by hydrolysis at internal ß-1,4-xylosidic linkages. It contains a "thumb" region whose flexibility has been suggested to affect the activity. The double mutation D11F/R122D was previously found to affect activity and potentially bias the thumb region to a more open conformation. We find that the D11F/R122D double mutation shows substrate-dependent effects, increasing activity on the non-native substrate ONPX2 but decreasing activity on its native xylan substrate. To characterize how the double mutant causes these kinetics changes, nuclear magnetic resonance (NMR) and molecular dynamics (MD) simulations were used to probe structural and flexibility changes. NMR chemical shift perturbations revealed structural changes in the double mutant relative to the wild-type, specifically in the thumb and fingers regions. Increased slow-timescale dynamics in the fingers region was observed as intermediate-exchange line broadening. Lipari-Szabo order parameters show negligible changes in flexibility in the thumb region in the presence of the double mutation. To help understand if there is increased energetic accessibility to the open state upon mutation, alchemical free energy simulations were employed that indicated thumb opening is more favorable in the double mutant. These studies aid in further characterizing how flexibility in adjacent regions affects the function of XylA.

In Vitro and In Silico Explorations of the Protein Conformational Changes of Corynebacterium glutamicum MshA, a Model Retaining GT-B Glycosyltransferase.

MshA is a GT-B glycosyltransferase catalyzing the first step in the biosynthesis of mycothiol. While many GT-B enzymes undergo an open-to-closed transition, MshA is unique because its 97° rotation is beyond the usual range of 10-25°. Molecular dynamics (MD) simulations were carried out for MshA in both ligand bound and unbound states to investigate the effect of ligand binding on localized protein dynamics and its conformational free energy landscape. Simulations showed that both the unliganded "opened" and liganded "closed" forms of the enzyme sample a wide degree of dihedral angles and interdomain distances with relatively low overlapping populations. Calculation of the free energy surface using replica exchange MD for the apo "opened" and an artificial generated apo "closed" structure revealed overlaps in the geometries sampled, allowing calculation of a barrier of 2 kcal/mol for the open-to-closed transition in the absence of ligands. MD simulations of fully liganded MshA revealed a smaller sampling of the dihedral angles. The localized protein fluctuation changes suggest that UDP-GlcNAc binding activates the motions of loops in the 1-l-myo-inositol-1-phosphate (I1P)-binding site despite little change in the interactions with UDP-GlcNAc. Circular dichroism, intrinsic fluorescence spectroscopy, and mutagenesis studies were used to confirm the ligand-induced structural changes in MshA. The results support a proposed mechanism where UDP-GlcNAc binds with rigid interactions to the C-terminal domain of MshA and activates flexible loops in the N-terminal domain for binding and positioning of I1P. This model can be used for future structure-based drug development of inhibitors of the mycothiol biosynthetic pathway.

Shared and distinct mechanisms of UBA1 inactivation across different diseases.

Most cellular ubiquitin signaling is initiated by UBA1, which activates and transfers ubiquitin to tens of E2 enzymes. Clonally acquired UBA1 missense mutations cause an inflammatory-hematologic overlap disease called VEXAS (vacuoles, E1, X-linked, autoinflammatory, somatic) syndrome. Despite extensive clinical investigation into this lethal disease, little is known about the underlying molecular mechanisms. Here, by dissecting VEXAS-causing UBA1 mutations, we discovered that p.Met41 mutations alter cytoplasmic isoform expression, whereas other mutations reduce catalytic activity of nuclear and cytoplasmic isoforms by diverse mechanisms, including aberrant oxyester formation. Strikingly, non-p.Met41 mutations most prominently affect transthioesterification, revealing ubiquitin transfer to cytoplasmic E2 enzymes as a shared property of pathogenesis amongst different VEXAS syndrome genotypes. A similar E2 charging bottleneck exists in some lung cancer-associated UBA1 mutations, but not in spinal muscular atrophy-causing UBA1 mutations, which instead, render UBA1 thermolabile. Collectively, our results highlight the precision of conformational changes required for faithful ubiquitin transfer, define distinct and shared mechanisms of UBA1 inactivation in diverse diseases, and suggest that specific E1-E2 modules control different aspects of tissue differentiation and maintenance.

Shared and Distinct Mechanisms of UBA1 Inactivation Across Different Diseases.

Most cellular ubiquitin signaling is initiated by UBA1, which activates and transfers ubiquitin to tens of E2 enzymes. Clonally acquired UBA1 missense mutations cause an inflammatory-hematologic overlap disease called VEXAS (vacuoles, E1, X-linked, autoinflammatory, somatic) syndrome. Despite extensive clinical investigation into this lethal disease, little is known about the underlying molecular mechanisms. Here, by dissecting VEXAS-causing UBA1 mutations, we discovered that p.Met41 mutations alter cytoplasmic isoform expression, whereas other mutations reduce catalytic activity of nuclear and cytoplasmic isoforms by diverse mechanisms, including aberrant oxyester formation. Strikingly, non-p.Met41 mutations most prominently affect transthioesterification, revealing ubiquitin transfer to cytoplasmic E2 enzymes as a shared property of pathogenesis amongst different VEXAS syndrome genotypes. A similar E2 charging bottleneck exists in some lung cancer-associated UBA1 mutations, but not in spinal muscular atrophy-causing UBA1 mutations, which instead, render UBA1 thermolabile. Collectively, our results highlight the precision of conformational changes required for faithful ubiquitin transfer, define distinct and shared mechanisms of UBA1 inactivation in diverse diseases, and suggest that specific E1-E2 modules control different aspects of tissue differentiation and maintenance.

Structural and photophysical characterization of the small ultra-red fluorescent protein.

The small Ultra-Red Fluorescent Protein (smURFP) represents a new class of fluorescent protein with exceptional photostability and brightness derived from allophycocyanin in a previous directed evolution. Here, we report the smURFP crystal structure to better understand properties and enable further engineering of improved variants. We compare this structure to the structures of allophycocyanin and smURFP mutants to identify the structural origins of the molecular brightness. We then use a structure-guided approach to develop monomeric smURFP variants that fluoresce with phycocyanobilin but not biliverdin. Furthermore, we measure smURFP photophysical properties necessary for advanced imaging modalities, such as those relevant for two-photon, fluorescence lifetime, and single-molecule imaging. We observe that smURFP has the largest two-photon cross-section measured for a fluorescent protein, and that it produces more photons than organic dyes. Altogether, this study expands our understanding of the smURFP, which will inform future engineering toward optimal FPs compatible with whole organism studies.

How Single Amino Acid Substitutions Can Disrupt a Protein Hetero-Dimer Interface: Computational and Experimental Studies of the LigAB Dioxygenase from Sphingobium sp. Strain SYK-6.

Protocatechuate 4,5-dioxygenase (LigAB) is a heterodimeric enzyme that catalyzes the dioxygenation of multiple lignin derived aromatic compounds. The active site of LigAB is at the heterodimeric interface, with specificity conferred by the alpha subunit and catalytic residues contributed by the beta subunit. Previous research has indicated that the phenylalanine at the 103 position of the alpha subunit (F103α) controls selectivity for the C5 position of the aromatic substrates, and mutations of this residue can enhance the rate of catalysis for substrates with larger functional groups at this position. While several of the mutations to this position (Valine, V; Threonine, T; Leucine, L; and Histidine, H) were catalytically active, other mutations (Alanine, A; and Serine, S) were found to have reduced dimer interface affinity, leading to challenges in copurifing the catalytically active enzyme complex under high salt conditions. In this study, we aimed to experimentally and computationally interrogate residues at the dimer interface to discern the importance of position 103α for maintaining the integrity of the heterodimer. Molecular dynamic simulations and electrophoretic mobility assays revealed a preference for nonpolar/aromatic amino acids in this position, suggesting that while substitutions to polar amino acids may produce a dioxygenase with a useful substrate utilization profile, those considerations may be off-set by potential destabilization of the catalytically active oligomer. Understanding the dimerization of LigAB provides insight into the multimeric proteins within the largely uncharacterized superfamily and characteristics to consider when engineering proteins that can degrade lignin efficiently. These results shed light on the challenges associated with engineering proteins for broader substrate specificity.

Kinetic Characterization and Computational Modeling of Escherichia coli Heptosyltransferase II: Exploring the Role of Protein Dynamics in Catalysis for GT-B Glycosyltransferase.

Glycosyltransferase (GT) enzymes promote the formation of glycosidic bonds between a sugar molecule and a diversity of substrates. Heptosyltransferase II (HepII) is a GT involved in the lipopolysaccharide (LPS) biosynthetic pathway that transfers the seven-carbon sugar (l-glycero-d-manno-heptose, Hep) onto a lipid-anchored glycopolymer (heptosylated Kdo2-Lipid A, Hep-Kdo2-Lipid A, or HLA). LPS plays a key role in Gram-negative bacterial sepsis, biofilm formation, and host colonization, and as such, LPS biosynthetic enzymes are targets for novel antimicrobial therapeutics. Three heptosyltransferases are involved in the inner-core LPS biosynthesis, with Escherichia coli HepII being the last to be quantitatively characterized in vivo. HepII shares modest sequence similarity with heptosyltransferase I (HepI) while maintaining a high degree of structural homology. Here, we report the first kinetic and biophysical characterization of HepII and demonstrate the properties of HepII that are shared with HepI, including sugar donor promiscuity and sugar acceptor-induced secondary structural changes, which results in significant thermal stabilization. HepII also has an increased catalytic efficiency and a significantly tighter binding affinity for both of its substrates compared to HepI. A structural model of the HepII ternary complex, refined by molecular dynamics simulations, was developed to probe the potentially important substrate-protein contacts. Ligand binding-induced changes in Trp fluorescence in HepII enabled the determination of substrate dissociation constants. Combined, these efforts meaningfully enhance our understanding of the heptosyltransferase family of enzymes and will aid in future efforts to design novel, potent, and specific inhibitors for this family of enzymes.

Discovery of first-in-class nanomolar inhibitors of heptosyltransferase I reveals a new aminoglycoside target and potential alternative mechanism of action.

A clinically relevant inhibitor for Heptosyltransferase I (HepI) has been sought after for many years because of its critical role in the biosynthesis of lipopolysaccharides on bacterial cell surfaces. While many labs have discovered or designed novel small molecule inhibitors, these compounds lacked the bioavailability and potency necessary for therapeutic use. Extensive characterization of the HepI protein has provided valuable insight into the dynamic motions necessary for catalysis that could be targeted for inhibition. Structural inspection of Kdo2-lipid A suggested aminoglycoside antibiotics as potential inhibitors for HepI. Multiple aminoglycosides have been experimentally validated to be first-in-class nanomolar inhibitors of HepI, with the best inhibitor demonstrating a Ki of 600 ± 90 nM. Detailed kinetic analyses were performed to determine the mechanism of inhibition while circular dichroism spectroscopy, intrinsic tryptophan fluorescence, docking, and molecular dynamics simulations were used to corroborate kinetic experimental findings. While aminoglycosides have long been described as potent antibiotics targeting bacterial ribosomes' protein synthesis leading to disruption of the stability of bacterial cell membranes, more recently researchers have shown that they only modestly impact protein production. Our research suggests an alternative and novel mechanism of action of aminoglycosides in the inhibition of HepI, which directly leads to modification of LPS production in vivo. This finding could change our understanding of how aminoglycoside antibiotics function, with interruption of LPS biosynthesis being an additional and important mechanism of aminoglycoside action. Further research to discern the microbiological impact of aminoglycosides on cells is warranted, as inhibition of the ribosome may not be the sole and primary mechanism of action. The inhibition of HepI by aminoglycosides may dramatically alter strategies to modify the structure of aminoglycosides to improve the efficacy in fighting bacterial infections.

Ligand-Induced Conformational and Dynamical Changes in a GT-B Glycosyltransferase: Molecular Dynamics Simulations of Heptosyltransferase I Complexes.

Understanding the dynamical motions and ligand recognition motifs of heptosyltransferase I (HepI) can be critical to discerning the behavior of other glycosyltransferase (GT) enzymes. Prior studies in our lab have demonstrated that GTs in the GT-B structural class, which are characterized by their connection of two Rossman-like domains by a linker region, have conserved structural fold and dynamical motions, despite low sequence homology, therefore making discoveries found in HepI transferable to other GT-B enzymes. Through molecular dynamics simulations and ligand binding free energy analysis of HepI in the apo and bound complexes (for all kinetically relevant combinations of the native substrates/products), we have determined the energetically favored enzymatic pathway for ligand binding and release. Our principal component, dynamic cross correlation, and network analyses of the simulations have revealed correlated motions involving residues within the N-terminal domain communicating with C-terminal domain residues via both proximal amino acid residues and also functional groups of the bound substrates. Analyses of the structural changes, energetics of substrate/product binding, and changes in pKa have elucidated a variety of inter and intradomain interactions that are critical for enzyme catalysis. These data corroborate our experimental observations of protein conformational changes observed in both presteady state kinetic and circular dichroism analyses of HepI. These simulations provided invaluable structural insights into the regions involved in HepI conformational rearrangement upon ligand binding. Understanding the specific interactions governing conformational changes is likely to enhance our efforts to develop novel dynamics disrupting inhibitors against GT-B structural enzymes in the future.

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

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.