Alternate structural conformations of Streptococcus pneumoniae hyaluronan lyase: insights into enzyme flexibility and underlying molecular mechanism of action

Streptococcus pneumoniae hyaluronan lyase is a surface enzyme of this Gram-positive bacterium. The enzyme degrades several biologically important, information-rich linear polymeric glycans: hyaluronan, unsulfated chondroitin, and some chondroitin sulfates. This degradation facilitates spreading of bacteria throughout the host tissues and presumably provides energy and a carbon source for pneumococcal cells. Its beta-elimination catalytic mechanism is an acid/base process termed proton acceptance and donation leading to cleavage of beta-1,4 linkages of the substrates. The degradation of hyaluronan occurs in two stages, initial endolytic cuts are followed by processive exolytic cleavage of one disaccharide at a time. In contrast, the degradation of chondroitins is purely endolytic. Structural studies together with flexibility analyses of two streptococcal enzymes, from S.pneumoniae and Streptococcus agalactiae, allowed for insights into this enzyme's molecular mechanism. Here, two new X-ray crystal structures of the pneumococcal enzyme in novel conformations are reported. These new conformations, complemented by molecular dynamics simulation results, directly confirm the predicted domain motions presumed to facilitate the processive degradative process. One of these new structures resembles the S.agalactiae enzyme conformation, and provides evidence of a uniform mechanistic/dynamic behavior of this protein across different bacteria.     

Structures of Streptococcus pneumoniae hyaluronate lyase in complex with chondroitin and chondroitin sulfate disaccharides. Insights into specificity and mechanism of action

Streptococcus pneumoniae hyaluronate lyase is a surface enzyme of this Gram-positive bacterium. The enzyme degrades hyaluronan and chondroitin/chondroitin sulfates by cleaving the beta1,4-glycosidic linkage between the glycan units of these polymeric substrates. This degradation helps spreading of this bacterial organism throughout the host tissues and facilitates the disease process caused by pneumococci. The mechanism of this degradative process is based on beta-elimination, is termed proton acceptance and donation, and involves selected residues of a well defined catalytic site of the enzyme. The degradation of hyaluronan alone is thought to proceed through a processive mode of action. The structures of complexes between the enzyme and chondroitin as well as chondroitin sulfate disaccharides allowed for the first detailed insights into these interactions and the mechanism of action on chondroitins. This degradation of chondroitin/chondroitin sulfates is nonprocessive and is selective for the chondroitin sulfates only with certain sulfation patterns. Chondroitin sulfation at the 4-position on the nonreducing site of the linkage to be cleaved or 2-sulfation prevent degradation due to steric clashes with the enzyme. Evolutionary studies suggest that hyaluronate lyases evolved from chondroitin lyases and still retained chondroitin/chondroitin sulfate degradation abilities while being specialized in the degradation of hyaluronan. The more efficient processive degradation mechanism has come to be preferred for the unsulfated substrate hyaluronan.     

Subtle functional collective motions in pancreatic-like ribonucleases: from ribonuclease A to angiogenin

The analysis of the dynamic behavior of enzymes is fundamental to structural biology. A direct relationship between protein flexibility and biological function has been shown for bovine pancreatic ribonuclease (RNase A) (Rasmussen et al., Nature 1992;357:423-424). More recently, crystallographic studies have shown that functional motions in RNase A involve the enzyme beta-sheet regions that move concertedly on substrate binding and release (Vitagliano et al., Proteins 2002;46:97-104). These motions have been shown to correspond to intrinsic dynamic properties of the native enzyme by molecular dynamics (MD) simulations. To unveil the occurrence of these collective motions in other members of pancreatic-like superfamily, we carried out MD simulations on human angiogenin (Ang). Essential dynamics (ED) analyses performed on the trajectories reveal that Ang exhibits collective motions similar to RNase A, despite the limited sequence identity (33%) of the two proteins. Furthermore, we show that these collective motions are also present in ensembles of experimentally determined structures of both Ang and RNase A. Finally, these subtle concerted beta-sheet motions were also observed for other two members of the pancreatic-like superfamily by comparing the ligand-bound and ligand-free structures of these enzymes. Taken together, these findings suggest that pancreatic-like ribonucleases share an evolutionary conserved dynamic behavior consisting of subtle beta-sheet motions, which are essential for substrate binding and release.     

Motional timescale predictions by molecular dynamics simulations: Case study using proline and hydroxyproline sidechain dynamics

We propose a new approach for force field optimizations which aims at reproducing dynamics characteristics using biomolecular MD simulations, in addition to improved prediction of motionally averaged structural properties available from experiment. As the source of experimental data for dynamics fittings, we use ¹³C NMR spin‐lattice relaxation times T₁ of backbone and sidechain carbons, which allow to determine correlation times of both overall molecular and intramolecular motions. For structural fittings, we use motionally averaged experimental values of NMR J couplings. The proline residue and its derivative 4‐hydroxyproline with relatively simple cyclic structure and sidechain dynamics were chosen for the assessment of the new approach in this work. Initially, grid search and simplexed MD simulations identified large number of parameter sets which fit equally well experimental J couplings. Using the Arrhenius‐type relationship between the force constant and the correlation time, the available MD data for a series of parameter sets were analyzed to predict the value of the force constant that best reproduces experimental timescale of the sidechain dynamics. Verification of the new force‐field (termed as AMBER99SB‐ILDNP) against NMR J couplings and correlation times showed consistent and significant improvements compared to the original force field in reproducing both structural and dynamics properties. The results suggest that matching experimental timescales of motions together with motionally averaged characteristics is the valid approach for force field parameter optimization. Such a comprehensive approach is not restricted to cyclic residues and can be extended to other amino acid residues, as well as to the backbone. Proteins 2014; 82:195–215. © 2013 Wiley Periodicals, Inc.     

Crystal structures of a family 8 polysaccharide lyase reveal open and highly occluded substrate-binding cleft conformations

Bacterial enzymatic degradation of glycosaminoglycans such as hyaluronan and chondroitin is facilitated by polysaccharide lyases. Family 8 polysaccharide lyase (PL8) enzymes contain at least two domains: one predominantly composed of α-helices, the α-domain, and another predominantly composed of β-sheets, the β-domain. Simulation flexibility analyses indicate that processive exolytic cleavage of hyaluronan, by PL8 hyaluronate lyases, is likely to involve an interdomain shift, resulting in the opening/closing of the substrate-binding cleft between the α- and β-domains, facilitating substrate translocation. Here, the Streptomyces coelicolor A3(2) PL8 enzyme was recombinantly expressed in and purified from Escherichia coli and biochemically characterized as a hyaluronate lyase. By using X-ray crystallography its structure was solved in complex with hyaluronan and chondroitin disaccharides. These findings show key catalytic interactions made by the different substrates, and on comparison with all other PL8 structures reveals that the substrate-binding cleft of the S. coelicolor enzyme is highly occluded. A third structure of the enzyme, harboring a mutation of the catalytic tyrosine, created via site-directed mutagenesis, interestingly revealed an interdomain shift that resulted in the opening of the substrate-binding cleft. These results add further support to the proposed processive mechanism of action of PL8 hyaluronate lyases and may indicate that the mechanism of action is likely to be universally used by PL8 hyaluronate lyases.     

Molecular Dynamics Study of the Structure,Flexibility and Dynamics of Thermostable L1 Lipase at High Temperatures

Molecular Dynamics (MD) simulations have been used to understand how protein structure, dynamics, and flexibility are affected by adaptation to high temperature for several years. We report here the results of the high temperature MD simulations of Bacillus stearothermophilus L1 (L1 lipase). We found that the N-terminal moiety of the enzyme showed a high flexibility and dynamics during high temperature simulations which preceded and followed by clear structural changes in two specific regions; the small domain and the main catalytic domain or core domain of the enzyme. These two domains interact with each other through a Zn²⁺-binding coordination with Asp-61 and Asp-238 from the core domain and His-81 and His-87 from the small domain. Interestingly, the His-81 and His-87 were among the highly fluctuated and mobile residues at high temperatures. The results appear to suggest that tight interactions of Zn²⁺-binding coordination with specified residues became weak at high temperature which suggests the contribution of this region to the thermostability of the enzyme. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

How does it really move? Recent progress in the investigation of protein nanosecond dynamics by NMR and simulation

Nuclear magnetic resonance (NMR) spin relaxation experiments currently probe molecular motions on timescales from picoseconds to nanoseconds. The detailed interpretation of these motions in atomic detail benefits from complementarity with the results from molecular dynamics (MD) simulations. In this mini-review, we describe the recent developments in experimental techniques to study the backbone dynamics from ¹⁵N relaxation and side-chain dynamics from ¹³C relaxation, discuss the different analysis approaches from model-free to dynamics detectors, and highlight the many ways that NMR relaxation experiments and MD simulations can be used together to improve the interpretation and gain insights into protein dynamics.     

Asymmetric structure and domain binding interfaces of human tyrosyl‐tRNA synthetase studied by molecular dynamics simulations

Human tyrosyl‐tRNA synthetase (HsTyrRS) is composed of two structural modules: N‐terminal catalytic core and an EMAP II‐like C‐terminal domain. The structures of these modules are known, but no crystal structure of the full‐length HsTyrRS is currently available. An all‐atom model of the full‐length HsTyrRS was developed in this work. The structure, dynamics, and domain binding interfaces of HsTyrRS were investigated by extensive molecular dynamics (MD) simulations. Our data suggest that HsTyrRS in solution consists of a number of compact asymmetric conformations, which differ significantly by their rigidity, internal mobility, orientation of C‐terminal modules, and the strength of interdomain binding. Interfaces of domain binding obtained in MD simulations are in perfect agreement with our previous coarse‐grained hierarchical rotations technique simulations. Formation of the hydrogen bonds between R93 residue of the ELR cytokine motif and the residues A340 and E479 in the C‐module was observed. This observation supports the idea that the lack of cytokine activity in the full‐length HsTyrRS is explained by interactions between N‐modules and C‐modules, which block the ELR motif. Copyright © 2013 John Wiley & Sons, Ltd.     

Role of conformational fluctuations in the enzymatic reaction of HIV-1 protease

The emergence of compensatory drug-resistant mutations in HIV-1 protease challenges the common view of the reaction mechanism of this enzyme. Here, we address this issue by performing classical and ab initio molecular dynamics simulations (MD) on a complex between the enzyme and a peptide substrate. The classical MD calculation reveals large-scale protein motions involving the flaps and the cantilever. These motions modulate the conformational properties of the substrate at the cleavage site. The ab initio calculations show in turn that substrate motion modulates the activation free energy barrier of the enzymatic reaction dramatically. Thus, the catalytic power of the enzyme does not arise from the presence of a pre-organized active site but from the protein mechanical fluctuations. The implications of this finding for the emergence of drug-resistance are discussed.     

Insight Derived from Molecular Dynamics Simulation into Substrate-Induced Changes in Protein Motions of Proteinase K

Abstract

Because of the significant industrial, agricultural and biotechnological importance of serine protease proteinase K, it has been extensively investigated using experimental approaches such as X-ray crystallography, site-directed mutagenesis and kinetic measurement. However, detailed aspects of enzymatic mechanism such as substrate binding, release and relevant regulation remain unstudied. Molecular dynamics (MD) simulations of the proteinase K alone and in complex with the peptide substrate AAPA were performed to investigate the effect of substrate binding on the dynamics/molecular motions of proteinase K. The results indicate that during simulations the substrate-complexed proteinase K adopt a more compact and stable conformation than the substrate-free form. Further essential dynamics (ED) analysis reveals that the major internal motions are confined within a subspace of very small dimension. Upon substrate binding, the overall flexibility of the protease is reduced; and the noticeable displacements are observed not only in substrate-binding regions but also in regions opposite the substrate-binding groove/pockets. The dynamic pockets caused by the large concerted motions are proposed to be linked to the substrate recognition, binding, orientation and product release; and the significant displacements in regions opposite the binding groove/pockets are considered to play a role in modulating the dynamics of enzyme-substrate interaction. Our simulation results complement the biochemical and structural studies, highlighting the dynamic mechanism of the functional properties of proteinase K.     

Common functionally important motions of the nucleotide‐binding domain of Hsp70

The 70 kDa heat shock proteins (Hsp70) are a family of molecular chaperones involved in protein folding, aggregate prevention, and protein disaggregation. They consist of the substrate‐binding domain (SBD) that binds client substrates, and the nucleotide‐binding domain (NBD), whose cycles of nucleotide hydrolysis and exchange underpin the activity of the chaperone. To characterize the structure–function relationships that link the binding state of the NBD to its conformational behavior, we analyzed the dynamics of the NBD of the Hsp70 chaperone from Bos taurus (PDB 3C7N:B) by all‐atom canonical molecular dynamics simulations. It was found that essential motions within the NBD fall into three major classes: the mutual class, reflecting tendencies common to all binding states, and the ADP‐ and ATP‐unique classes, which reflect conformational trends that are unique to either the ADP‐ or ATP‐bound states, respectively. “Mutual” class motions generally describe “in‐plane” and/or “out‐of‐plane” (scissor‐like) rotation of the subdomains within the NBD. This result is consistent with experimental nuclear magnetic resonance data on the NBD. The “unique” class motions target specific regions on the NBD, usually surface loops or sites involved in nucleotide binding and are, therefore, expected to be involved in allostery and signal transmission. For all classes, and especially for those of the “unique” type, regions of enhanced mobility can be identified; these are termed “hot spots,” and their locations generally parallel those found by NMR spectroscopy. The presence of magnesium and potassium cations in the nucleotide‐binding pocket was also found to influence the dynamics of the NBD significantly. Proteins 2015; 83:282–299. © 2014 Wiley Periodicals, Inc.     

Computational study of conformational changes in human 3-hydroxy-3-methylglutaryl coenzyme reductase induced by substrate binding

Abstract

The enzyme 3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGR) is mainly involved in the regulation of cholesterol biosynthesis. HMGR catalyses the reduction of 3-hydroxy-3-methylglutaryl coenzyme A (HMG-CoA) to mevalonate at the expense of two NADPH molecules in a two-step reversible reaction. In the present study, we constructed a model of human HMGR (hHMGR) to explore the conformational changes of HMGR in complex with HMG-CoA and NADPH. In addition, we analysed the complete sequence of the Flap domain using molecular dynamics (MD) simulations and principal component analysis (PCA). The simulations revealed that the Flap domain plays an important role in catalytic site activation and substrate binding. The apo form of hHMGR remained in an open state, while a substrate-induced closure of the Flap domain was observed for holo hHMGR. Our study also demonstrated that the phosphorylation of Ser872 induces significant conformational changes in the Flap domain that lead to a complete closure of the active site, suggesting three principal conformations for the first stage of hHMGR catalysis. Our results were consistent with previous proposed models for the catalytic mechanism of hHMGR.

Communicated by Ramaswamy H. Sarma     

Molecular dynamics of leucine and dopamine transporter proteins in a model cell membrane lipid bilayer

The dopamine transporter (DAT) operates via facilitated diffusion, harnessing an inward Na⁺ gradient to drive dopamine from the extracellular synaptic cleft to the neuron interior. The DAT is relevant to central nervous system disorders such as Parkinson disease and attention‐deficit hyperactivity disorder and is the primary site of action for the abused psychostimulants cocaine and amphetamines. Crystallization of a DAT homolog, the bacterial leucine transporter LeuT, provided the first reliable 3‐D DAT template. Here, the LeuT crystal structure and the DAT molecular model have been combined with their respective substrates, leucine and dopamine, in lipid bilayer molecular dynamics simulations toward tracking substrate movement along the protein's substrate/ion permeation pathway. Specifically, movement of residue pairs that comprise the “external gate” was followed as a function of substrate presence. The transmembrane (TM) 1 arginine‐TM 10 aspartate strut formed less readily in DAT compared with LeuT, with or without substrate present. For LeuT but not DAT, the addition of substrate enhanced the chances of forming the TM 1‐10 bridge. Also, movement of the fourth extracellular loop EL‐4 in the presence of substrate was more pronounced for DAT, the EL‐4 unwinding to a degree. The overall similarity between the LeuT and DAT molecular dynamics simulations indicated that LeuT was a legitimate model to guide DAT structure‐function predictions. There were, nevertheless, differences significant enough to allow for DAT‐unique insights, which may include how cocaine, methylphenidate (Ritalin, NIDA Drug Supply, Rockville, MD), and other DAT blockers are not recognized as substrates even though they can access the primary substrate binding pocket. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Reactive Oxygen Species and Hyaluronidase 2 Regulate Airway Epithelial Hyaluronan Fragmentation

Hyaluronidase 2 (Hyal2) is a hyaluronan (HA)-degrading enzyme found intracellularly or/and anchored to the plasma membrane through glycosylphosphatidylinositol (GPI). Normal human bronchial epithelial cells (NHBE) grown at the air-liquid interphase (ALI), treated with PI-specific phospholipase C (PI-PLC), exhibited increased Hyal activity in secretions and decreased protein and activity on the apical membrane, confirming that GPI-anchored Hyal2 is expressed in NHBE cells and it remains active in its soluble form. We have reported that HA degradation was mediated by reactive oxygen species (ROS) in human airways. Here we show that ROS increase Hyal2 expression and activity in NHBE cells and that the p38MAPK signaling pathway is involved in this effect. Hyal2 induction was confirmed by using small interfering RNA (siRNA) expressing lentivirus. These in vitro findings correlated in vivo with smokers, where increased Hyal2 immunoreactivity in the epithelium was associated with augmented levels of HA and the appearance of low molecular mass HA species in bronchial secretions. In summary, this work provides evidence that ROS induce Hyal2, suggesting that Hyal2 is likely responsible for the sustained HA fragmentation in the airway lumen observed in inflammatory conditions associated with oxidative stress.     

Substrate Binding Specifically Modulates Domain Arrangements in Adenylate Kinase

The enzyme adenylate kinase (ADK) features two substrate binding domains that undergo large-scale motions during catalysis. In the apo state, the enzyme preferentially adopts a globally open state with accessible binding sites. Binding of two substrate molecules (AMP + ATP or ADP + ADP) results in a closed domain conformation, allowing efficient phosphoryl-transfer catalysis. We employed molecular dynamics simulations to systematically investigate how the individual domain motions are modulated by the binding of substrates. Two-dimensional free-energy landscapes were calculated along the opening of the two flexible lid domains for apo and holo ADK as well as for all single natural substrates bound to one of the two binding sites of ADK. The simulations reveal a strong dependence of the conformational ensembles on type and binding position of the bound substrates and a nonsymmetric behavior of the lid domains. Altogether, the ensembles suggest that, upon initial substrate binding to the corresponding lid site, the opposing lid is maintained open and accessible for subsequent substrate binding. In contrast, ATP binding to the AMP-lid induces global domain closing, preventing further substrate binding to the ATP-lid site. This might constitute a mechanism by which the enzyme avoids the formation of a stable but enzymatically unproductive state.     

Terminal sialic acids on CD44 N‐glycans can block hyaluronan binding by forming competing intramolecular contacts with arginine sidechains

Specific sugar residues and their linkages form the basis of molecular recognition for interactions of glycoproteins with other biomolecules. Seemingly small changes, like the addition of a single monosaccharide in the covalently attached glycan component of glycoproteins, can greatly affect these interactions. For instance, the sialic acid capping of glycans affects protein‐ligand binding involved in cell–cell and cell–matrix interactions. CD44 is a single‐pass transmembrane glycoprotein whose binding with its carbohydrate ligand hyaluronan (HA), an extracellular matrix component, mediates processes such as leukocyte homing, cell adhesion, and tumor metastasis. This binding is highly regulated by glycosylation of the N‐terminal extracellular hyaluronan‐binding domain (HABD); specifically, sialic acid capped N‐glycans of HABD inhibit ligand binding. However, the molecular mechanism behind this sialic acid mediated regulation has remained unknown. Two of the five N‐glycosyation sites of HABD have been previously identified as having the greatest inhibitory effect on HA binding, but only if the glycans contain terminal sialic acid residues. These two sites, Asn25 and Asn120, were chosen for in silico glycosylation in this study. Here, from extensive standard molecular dynamics simulations and biased simulations, we propose a molecular mechanism for this behavior based on spontaneously‐formed charge‐paired hydrogen bonding interactions between the negatively‐charged sialic acid residues and positively‐charged Arg sidechains known to be critically important for binding to HA, which itself is negatively charged. Such intramolecular hydrogen bonds would preclude associations critical to hyaluronan binding. This observation suggests how CD44 and related glycoprotein binding is regulated by sialylation as cellular environments fluctuate. Proteins 2014; 82:3079–3089. © 2014 Wiley Periodicals, Inc.     

Thermal stability of single‐domain antibodies estimated by molecular dynamics simulations

Single‐domain antibodies (sdAbs) function like regular antibodies, however, consist of only one domain. Because of their low molecular weight, sdAbs have advantages with respect to production and delivery to their targets and for applications such as antibody drugs and biosensors. Thus, sdAbs with high thermal stability are required. In this work, we chose seven sdAbs, which have a wide range of melting temperature (T_m) values and known structures. We applied molecular dynamics (MD) simulations to estimate their relative stability and compared them with the experimental data. High‐temperature MD simulations at 400 K and 500 K were executed with simulations at 300 K as a control. The fraction of native atomic contacts, Q, measured for the 400 K simulations showed a fairly good correlation with the T_m values. Interestingly, when the residues were classified by their hydrophobicity and size, the Q values of hydrophilic residues exhibited an even better correlation, suggesting that stabilization is correlated with favorable interactions of hydrophilic residues. Measuring the Q value on a per‐residue level enabled us to identify residues that contribute significantly to the instability and thus demonstrating how our analysis can be used in a mutant case study.     

Molecular dynamics simulations reveal that apo‐HisJ can sample a closed conformation

The Escherichia coli histidine binding protein HisJ is a type II periplasmic binding protein (PBP) that preferentially binds histidine and interacts with its cytoplasmic membrane ABC transporter, HisQMP₂, to initiate histidine transport. HisJ is a bilobal protein where the larger Domain 1 is connected to the smaller Domain 2 via two linking strands. Type II PBPs are thought to undergo “Venus flytrap” movements where the protein is able to reversibly transition from an open to a closed conformation. To explore the accessibility of the closed conformation to the apo state of the protein, we performed a set of all‐atom molecular dynamics simulations of HisJ starting from four different conformations: apo‐open, apo‐closed, apo‐semiopen, and holo‐closed. The simulations reveal that the closed conformation is less dynamic than the open one. HisJ experienced closing motions and explored semiopen conformations that reverted to closed forms resembling those found in the holo‐closed state. Essential dynamics analysis of the simulations identified domain closing/opening and twisting as main motions. The formation of specific inter‐hinge strand and interdomain polar interactions contributed to the adoption of the closed apo‐conformations although they are up to 2.5‐fold less prevalent compared with the holo‐closed simulations. The overall sampling of the closed form by apo‐HisJ provides a rationale for the binding of unliganded PBPs with their cytoplasmic membrane ABC transporters. Proteins 2014; 82:386–398. © 2013 Wiley Periodicals, Inc.