A molecular dynamics study of Beta-Glucosidase B upon small substrate binding

Paenibacillus polymyxa β-glucosidase B (BglB), belongs to a GH family 1, is a monomeric enzyme that acts as an exo-β-glucosidase hydrolysing cellobiose and cellodextrins of higher degree of polymerization using retaining mechanism. A molecular dynamics (MD) simulation was performed at 300 K under periodic boundary condition for 5 ns using the complexes structure obtained from previous docking study, namely BglB-Beta-d-glucose and BglB-Cellobiose. From the root-mean-square deviation analysis, both enzyme complexes were reported to deviate from the initial structure in the early part of the simulation but it was stable afterwards. The root-mean-square fluctuation analysis revealed that the most flexible regions comprised of the residues from 26 to 29, 43 to 53, 272 to 276, 306 to 325 and 364 to 367. The radius of gyration analysis had shown the structure of BglB without substrate became more compact towards the end of the simulation compare to other two complexes. The residues His122 and Trp410 were observed to form stable hydrogen bond with occupancy higher than 10%. In conclusion, the behaviour of BglB enzyme towards the substrate binding was successfully explored via MD simulation approaches.     

Conformational stability of PCID2 upon DSS1 binding with molecular dynamics simulation

DSS1 is a small acidic intrinsically disordered protein (IDP) that can fold upon binding with PCID2 TREX-2. The resulting complex plays a key role in mRNA export. However, the binding mechanism between DSS1 and PCID2 is unsolved. Here, three independent 500-ns molecular dynamics (MD) simulations were performed to study the DSS1–PCID2 binding mechanism by comparing apo-PCID2 and bound PCID2. The results show that the conformational variation of bound PCID2 is smaller than that of apo-PCID2, especially in the binding domain of two helices (helix IV and VIII). The probability of coil formation between helix III and helix IV of bound PCID2 increases, and a short anti-parallel β-sheet forms upon DSS1 binding. The decomposition of binding free energy into protein and residue pairs suggests that electrostatic and hydrophobic interactions play key roles in the recognition between DSS1 and PCID2. There is a hydrophobic core of seven residues in DSS1 favorable to the binding of PCID2. These analytical methods can be used to reveal the recognition mechanisms of other IDPs and their partners.

Open image in new window

Figure Differences of second structure of PCID2 in bound and unbound states. The interaction surface between the helix VIII of PCID2 and helix of DSS1

     

Conformational dynamics of Tetracenomycin aromatase/cyclase regulate polyketide binding and enzyme aggregation propensity

BackgroundThe N-terminal domain of Tetracenomycin aromatase/cyclase (TcmN), an enzyme derived from Streptomyces glaucescens, is involved in polyketide cyclization, aromatization, and folding. Polyketides are a diverse class of secondary metabolites produced by certain groups of bacteria, fungi, and plants with various pharmaceutical applications. Examples include antibiotics, such as tetracycline, and anticancer drugs, such as doxorubicin. Because TcmN is a promising enzyme for in vitro production of polyketides, it is important to identify conditions that enhance its thermal resistance and optimize its function.MethodsTcmN unfolding, stability, and dynamics were evaluated by fluorescence spectroscopy, circular dichroism, nuclear magnetic resonance ¹⁵N relaxation experiments, and microsecond molecular dynamics (MD) simulations.ResultsTcmN thermal resistance was enhanced at low protein and high salt concentrations, was pH-dependent, and denaturation was irreversible. Conformational dynamics on the μs-ms timescale were detected for residues in the substrate-binding cavity, and two predominant conformers representing opened and closed cavity states were observed in the MD simulations.ConclusionBased on the results, a mechanism was proposed in which the thermodynamics and kinetics of the TcmN conformational equilibrium modulate enzyme function by favoring ligand binding and avoiding aggregation.General significanceUnderstanding the principles underlying TcmN stability and dynamics may help in designing mutants with optimal properties for biotechnological applications.     

Protection of substrate against enzymatic action by binding to proteins. Dependence upon enzyme and binding protein

In fetal or adult rats estradiol is carried in the plasma by alpha-fetoprotein or albumin. The protection of the carriers toward enzymatic oxidation by 17 beta-hydroxysteroid dehydrogenase from rat liver has been studied. Concentrations of carrier protein and estradiol were adjusted to give free estradiol concentrations varying from Km/10 to Km/100 and the ratio of the catalytic velocity to that observed for the same concentration of free estradiol in the absence of carrier protein were recorded. With alpha-fetoprotein the ratio fell as expected as the carrier concentration was increased, but with serum albumin the ratio was close to unity when the concentration of carrier protein was increased from 70 to 700 microM. Thus, a fetoprotein but not albumin protected the steroid against catalytic oxidation. Similar experiments were carried out replacing the mammalian enzyme by the dehydrogenase from Pseudomonas testosteroni, and in this case, both carrier proteins protected the substrate. The lack of protection by albumin against the dehydrogenation by the mammalian enzyme is discussed.     

Conformational changes of colicin Ia channel-forming domain upon membrane binding: a solid-state NMR study

Channel-forming colicins are bactericidal proteins that spontaneously insert into hydrophobic lipid bilayers. We have used magic-angle spinning solid-state nuclear magnetic resonance spectroscopy to examine the conformational differences between the water-soluble and the membrane-bound states of colicin Ia channel domain, and to study the effect of bound colicin on lipid bilayer structure and dynamics. We detected (13)C and (15)N isotropic chemical shift differences between the two forms of the protein, which indicate structural changes of the protein due to membrane binding. The Val C(alpha) signal, unambiguously assigned by double-quantum experiments, gave a 0.6 ppm downfield shift in the isotropic position and a 4 ppm reduction in the anisotropic chemical shift span after membrane binding. These suggest that the alpha-helices in the membrane-bound colicin adopt more ideal helical torsion angles as they spread onto the membrane. Colicin binding significantly reduced the lipid chain order, as manifested by (2)H quadrupolar couplings. These results are consistent with the model that colicin Ia channel domain forms an extended helical array at the membrane-water interface upon membrane binding.     

Conformational changes of colicin Ia channel-forming domain upon membrane binding: a solid-state NMR study

Channel-forming colicins are bactericidal proteins that spontaneously insert into hydrophobic lipid bilayers. We have used magic-angle spinning solid-state nuclear magnetic resonance spectroscopy to examine the conformational differences between the water-soluble and the membrane-bound states of colicin Ia channel domain, and to study the effect of bound colicin on lipid bilayer structure and dynamics. We detected ¹³C and ¹⁵N isotropic chemical shift differences between the two forms of the protein, which indicate structural changes of the protein due to membrane binding. The Val Cα signal, unambiguously assigned by double-quantum experiments, gave a 0.6 ppm downfield shift in the isotropic position and a 4 ppm reduction in the anisotropic chemical shift span after membrane binding. These suggest that the α-helices in the membrane-bound colicin adopt more ideal helical torsion angles as they spread onto the membrane. Colicin binding significantly reduced the lipid chain order, as manifested by ²H quadrupolar couplings. These results are consistent with the model that colicin Ia channel domain forms an extended helical array at the membrane-water interface upon membrane binding.     

Conformational change upon ligand binding and dynamics of the PDZ domain from leukemia-associated Rho guanine nucleotide exchange factor

Leukemia-associated Rho guanine nucleotide exchange factor (LARG) is a RhoA-specific guanine nucleotide exchange factor (GEF) that can activate RhoA. The PDZ (PSD-95/Disc-large/ZO-1 homology) domain of LARG interacts with membrane receptors, which can relay extracellular signals to RhoA signal transduction pathways. Until now there is no structural and dynamic information about these interactions. Here we report the NMR structures of the LARG PDZ in the apo form and in complex with the plexin-B1 C-terminal octapeptide. Unobservable resonances of the residues in betaB/betaC and betaE/alphaB loops in apo state were observed in the complex state. A distinct region of the binding groove in the LARG PDZ was found to undergo conformational change compared with other PDZs. Analysis of the (15)N relaxation data using reduced spectral density mapping shows that the apo LARG PDZ (especially its ligand-binding groove) is flexible and exhibits internal motions on both picosecond to nanosecond and microsecond to millisecond timescales. Mutagenesis and thermodynamic studies indicate that the conformation of the betaB/betaC and betaE/alphaB loops affects the PDZ-peptide interaction. It is suggested that the conformational flexibility could facilitate the change of structures upon ligand binding.     

Conformational dynamics of DnaB helicase upon DNA and nucleotide binding: analysis by intrinsic tryptophan fluorescence quenching

DnaB helicase of E. coli unwinds duplex DNA in the replication fork using the energy of ATP hydrolysis. We have analyzed structural and conformational changes in the DnaB protein in various nucleotides and DNA bound intermediate states by fluorescence quenching analysis of intrinsic fluorescence of native tryptophan (Trp) residues in DnaB. Fluorescence quenching analysis indicated that Trp48 in domain alpha is in a hydrophobic environment and resistant to fluorescence quenchers such as potassium iodide (KI). In domain beta, Trp294 was found to be in a partially hydrophobic environment, whereas Trp456 in domain gamma appeared to be in the least hydrophobic environment. Binding of oligonucleotides to DnaB helicase resulted in a significant attenuation of the fluorescence quenching profile, indicating a change in conformation. ATPgammaS or ATP binding appeared to lead to a conformation in which Trp residues had a higher degree of solvent exposure and fluorescence quenching. However, the most dramatic increase of Trp fluorescence quenching was observed with ADP binding with a possible conformational relaxation. Site-specific Trp --> Cys mutants of DnaB helicase demonstrated that conformational change upon ADP binding could be attributed exclusively to a conformational transition in the alpha domain leading to an increase in the solvent exposure of Trp48. However, formation of DnaB.ATPgammaS.DNA ternary complex led to a conformation with a fluorescence quenching profile similar to that observed with DnaB alone. The DnaB.ADP.DNA ternary complex produced a quenching curve similar to that of DnaB.ADP complex pointing to a change in conformation due to ATP hydrolysis. There are at least four identifiable structural/conformational states of DnaB helicase that are likely important in the helicase activity. The noncatalytic alpha domain in the N-terminus appeared to undergo the most significant conformational changes during nucleotide binding and hydrolysis. This is the first reported elucidation of the putative role of domain alpha, which is essential for DNA helicase action. We have correlated these results with partial structural models of alpha, beta, and gamma domains     

Differentiation of enzyme and substrate binding in the prothrombinase complex

The Ca2+ dependence of factor Xa binding to phospholipid vesicles was measured in the presence and absence of factor Va. The increase in polarization of a fluorescently labeled derivative of factor Xa, [5-(dimethylamino)-1-naphthalenesulfonyl] glutamylglycylarginyl factor Xa (Dns-EGR-Xa), was used as a probe to measure the interaction of factor Xa with phospholipid. The Ca2+ concentration required for half-maximal binding of Dns-EGR-Xa to phospholipid vesicles was 3.5 X 10(-4) M in the presence of factor Va and 9.5 X 10(-4) M in the absence of factor Va. At a Ca2+ concentration of 5 X 10(-4) M, the binding of Dns-EGR-Xa to phospholipid-bound factor Va was near maximal, whereas there was no detectable interaction of Dns-EGR-Xa with phospholipid alone at this Ca2+ concentration as detected by fluorescence polarization. These results were qualitatively confirmed by high-performance liquid chromatography. The rate of hydrolysis of the factor Xa synthetic substrate, benzoylisoleucylglutamylglycylarginine p-nitroanilide, by factor Xa in the presence of factor Va and phospholipid decreased in a Ca2+-dependent manner. These data were analyzed as fraction of factor Xa bound to the phospholipid. A Ca2+ concentration of 2.7 X 10(-4) M resulted in half-maximal binding by this technique. The relationship observed between rates of prothrombin activation and Ca2+ concentration could be predicted quantitatively from calculations of local enzyme and substrate concentrations.     

Conformational dynamics of RNA-peptide binding: a molecular dynamics simulation study

Molecular dynamics simulations of the binding of the heterochiral tripeptide KkN to the transactivation responsive (TAR) RNA of HIV-1 is presented, using an all-atom force field with explicit water. To obtain starting structures for the TAR-KkN complex, semirigid docking calculations were performed that employ an NMR structure of free TAR RNA. The molecular dynamics simulations show that the starting structures in which KkN binds to the major groove of TAR (as it is the case for the Tat-TAR complex of HIV-1) are unstable. On the other hand, the minor-groove starting structures are found to lead to several binding modes, which are stabilized by a complex interplay of stacking, hydrogen bonding, and electrostatic interactions. Although the ligand does not occupy the binding position of Tat protein, it is shown to hinder the interhelical motion of free TAR RNA. The latter is presumably necessary to achieve the conformational change of TAR RNA to bind Tat protein. Considering the time evolution of the trajectories, the binding process is found to be ligand-induced and cooperative. That is, the conformational rearrangement only occurs in the presence of the ligand and the concerted motion of the ligand and a large part of the RNA binding site is necessary to achieve the final low-energy binding state.     

Changes in dynamics upon oligomerization regulate substrate binding and allostery in amino acid kinase family members

Oligomerization is a functional requirement for many proteins. The interfacial interactions and the overall packing geometry of the individual monomers are viewed as important determinants of the thermodynamic stability and allosteric regulation of oligomers. The present study focuses on the role of the interfacial interactions and overall contact topology in the dynamic features acquired in the oligomeric state. To this aim, the collective dynamics of enzymes belonging to the amino acid kinase family both in dimeric and hexameric forms are examined by means of an elastic network model, and the softest collective motions (i.e., lowest frequency or global modes of motions) favored by the overall architecture are analyzed. Notably, the lowest-frequency modes accessible to the individual subunits in the absence of multimerization are conserved to a large extent in the oligomer, suggesting that the oligomer takes advantage of the intrinsic dynamics of the individual monomers. At the same time, oligomerization stiffens the interfacial regions of the monomers and confers new cooperative modes that exploit the rigid-body translational and rotational degrees of freedom of the intact monomers. The present study sheds light on the mechanism of cooperative inhibition of hexameric N-acetyl-L-glutamate kinase by arginine and on the allosteric regulation of UMP kinases. It also highlights the significance of the particular quaternary design in selectively determining the oligomer dynamics congruent with required ligand-binding and allosteric activities.     

Conformational dynamics play important roles upon the function of N-acetylglutamate kinase

N-acetylglutamate kinase (NAGK) catalyzes the phosphorylation of N-acetylglutamate. In many bacteria, NAGK catalysis is the rate controlling step in the L-arginine biosynthesis pathway from glutamate to L-arginine and is allosterically inhibited by L-arginine. Many data show that conformational dynamics of NAGKs are essential for their function. The demonstration of the conformational mechanism provides a potential way to improve the yield of arginine. Due to the lack of NAGK catalysis step in arginine synthesis route of mammals, the elucidation of the dynamic mechanism can also provide a way to design a new antivirus drug. This paper reviews how the dynamics affect the activity of NAGKs and are controlled by the effectors. X-ray crystallography and modeling data have shown that in NAGKs, the structural elements required for inhibitor and substrate binding, catalysis and product release, are highly mobile. It is possible to eliminate the inhibition of the arginine and/or block the synthesis of arginine by disturbing the flexibility of the NAGKs. Amino acid kinase family is thought to share some common dynamic features; the flexible structural elements of NAGKs have been identified, but the details of the dynamics and the signal transfer pathways are yet to be elucidated.

     

Conformational changes in the multidrug transporter EmrE associated with substrate binding

EmrE is a bacterial multidrug transporter of the small multidrug resistance family, which extrudes large hydrophobic cations such as tetraphenylphosphonium (TPP(+)) out of the cell by a proton antiport mechanism. Binding measurements were performed on purified EmrE solubilized in dodecylmaltoside to determine the stoichiometry of TPP(+) binding; the data showed that one TPP(+) molecule bound per EmrE dimer. Reconstitution of purified EmrE at low lipid:protein ratios in either the presence or the absence of TPP(+) produced well ordered two-dimensional crystals. Electron cryo-microscopy was used to collect images of frozen hydrated EmrE crystals and projection maps were determined by image processing to 7A resolution. An average native EmrE projection structure was calculated from the c222 and p222(1) crystals, which was subsequently subtracted from the average of two independent p2 projection maps of EmrE with TPP(+) bound. The interpretation of the difference density image most consistent with biochemical data suggested that TPP(+) bound at the monomer-monomer interface in the centre of the EmrE dimer, and resulted in the movement of at least one transmembrane alpha-helix.     

Conformational rearrangement of beta-lactoglobulin upon interaction with an anionic membrane

The recent availability of the SHV-1 beta-lactamase crystal structure provides a framework for the understanding of the functional role of amino acid residues in this enzyme. To that end, we have constructed by site-directed mutagenesis 18 variants of the SHV beta-lactamase: an extended spectrum group: Gly238Ser, Gly238Ser-Glu240Lys, Asp104Lys-Gly238Ser, Asp104Lys-Thr235Ser-Gly238Ser, Asp179Asn, Arg164His, and Arg164Ser; an inhibitor resistant group: Arg244Ser, Met69Ile, Met69Leu, and Ser130Gly; mutants that are synergistic with those that confer resistance to oxyimino-cephalosporins: Asp104Glu, Asp104Lys, Glu240Lys, and Glu240Gln; and structurally conserved mutants: Thr235Ser, Thr235Ala and Glu166Ala. Among the extended spectrum group the combination of high-level ampicillin and cephalosporin resistance was demonstrated in the Escherichia coli DH10B strains possessing the Gly238Ser mutation: Gly238Ser, Gly238Ser-Glu240Lys, Asp104Lys-Gly238Ser, and Asp104Lys-Thr235Ser-Gly238Ser. Of the inhibitor resistant group, the Ser130Gly mutant was the most resistant to ampicillin/clavulanate. Using a polyclonal anti-SHV antibody, we assayed steady state protein expression levels of the SHV beta-lactamase variants. Mutants with the Gly238Ser substitution were among the most highly expressed. The Gly238Ser substitution resulted in an improved relative k(cat)/K(m) value for cephaloridine and oxyimino-cephalosporins compared to SHV-1 and Met69Ile. In our comparative survey, the Gly238Ser and extended spectrum beta-lactamase variants containing this substitution exhibited the greatest substrate versatility against penicillins and cephalosporins and greatest protein expression. This defines a unique role of Gly238Ser in broad-spectrum beta-lactam resistance in this family of class A beta-lactamases.     

Conformational change of erythroid alpha-spectrin at the tetramerization site upon binding beta-spectrin

We previously determined the solution structures of the first 156 residues of human erythroid alpha-spectrin (SpalphaI-1-156, or simply Spalpha). Spalpha consists of the tetramerization site of alpha-spectrin and associates with a model beta-spectrin protein (Spbeta) with an affinity similar to that of native alpha- and beta-spectrin. Upon alphabeta-complex formation, our previous results indicate that there is an increase in helicity in the complex, suggesting conformational change in either Spalpha or Spbeta or in both. We have now used isothermal titration calorimetry, circular dichroism, static and dynamic light scattering, and solution NMR methods to investigate properties of the complex as well as the conformation of Spalpha in the complex. The results reveal a highly asymmetric complex, with a Perrin shape parameter of 1.23, which could correspond to a prolate ellipsoid with a major axis of about five and a minor axis of about one. We identified 12 residues, five prior to and seven following the partial domain helix in Spalpha that moved freely relative to the structural domain in the absence of Spbeta but when in the complex moved with a mobility similar to that of the structural domain. Thus, it appears that the association with Spbeta induced an unstructured-to-helical conformational transition in these residues to produce a rigid and asymmetric complex. Our findings may provide insight toward understanding different association affinities of alphabeta-spectrin at the tetramerization site for erythroid and non-erythroid spectrin and a possible mechanism to understand some of the clinical mutations, such as L49F of alpha-spectrin, which occur outside the functional partial domain region.     

Conformational changes and complement activation induced upon antigen binding to antibodies

The interaction between antibodies specific for the loop region of lysozyme and this monovalent antigenic determinant was compared to their interaction with either the dimeric derivative of the loop (bis-loop) or the polyvalent antigen loop-A--L. This comparison was based on complement fixation capacity, and on spectroscopic changes in the circular polarized fluorescence (CPL) of the antibodies. The binding of the loop to its antibodies causes spectroscopic changes, assigned to conformational changes in the Fc region, which are not accompanied by complement activation. However, the binding of the bisloop led to distinctly different CPL changes and concomitantly to complement fixation comparable to that induced upon binding of the loop-A--L. Complement fixation and the CPL changes can be induced in the intact antibodies only; reduction of the interchain disulfide bridges, or cleavage by papain or pepsin digestion abolishes both activities.     

Conformational change of pseudouridine 55 synthase upon its association with RNA substrate

Pseudouridine 55 synthase (Ψ55S) catalyzes isomerization of uridine (U) to pseudouridine (Ψ) at position 55 in transfer RNA. The crystal structures of Thermotoga maritima Ψ55S, and its complex with RNA, have been determined at 2.9 and 3.0 Å resolutions, respectively. Structural comparisons with other families of pseudouridine synthases (ΨS) indicate that Ψ55S may acquire its ability to recognize a stem–loop RNA substrate by two insertions of polypeptides into the ΨS core. The structure of apo-Ψ55S reveals that these two insertions interact with each other. However, association with RNA substrate induces substantial conformational change in one of the insertions, resulting in disruption of interaction between insertions and association of both insertions with the RNA substrate. Specific interactions between two insertions, as well as between the insertions and the RNA substrate, account for the molecular basis of the conformational change.     

Stalk dynamics and lipid flow upon membrane hemifusion

A hydrodynamic theory describing the stalk dynamics and lipid flows upon BLM hemifusion was developed. The value of intermonolayer viscosity, etar, for membranes formed from azolectin mixed with lysophosphatidylcholine (7.10(-4) mg/ml) in n-decane, etar approximately 10(-9) g/s, was determined from comparison of the theoretical calculations and literature data. For membranes formed in squalene, the values etar approximately 10(-7) g/s for phosphatidylethanolamine and etar approximately 2-10(-7) g/s for azolectin were obtained. The calculated values are close to the published results of independent experiments which shows that the developed theory describes well the stalk growth and lipid flow.     

Structural changes of creatine kinase upon substrate binding.

Small-angle x-ray scattering was used to investigate structural changes upon binding of individual substrates or a transition state analog complex (TSAC; Mg-ADP, creatine, and KNO3) to creatine kinase (CK) isoenzymes (dimeric muscle-type (M)-CK and octameric mitochondrial (Mi)-CK) and monomeric arginine kinase (AK). Considerable changes in the shape and the size of the molecules occurred upon binding of Mg-nucleotide or TSAC. The radius of gyration of Mi-CK was reduced from 55.6 A (free enzyme) to 48.9 A (enzyme plus Mg-ATP) and to 48.2 A (enzyme plus TSAC). M-CK showed similar changes from 28.0 A (free enzyme) to 25.6 A (enzyme plus Mg-ATP) and to 25.5 A (enzyme plus TSAC). Creatine alone did not lead to significant changes in the radii of gyration, nor did free ATP or ADP. AK also showed a change of the radius of gyration from 21.5 A (free enzyme) to 19.7 A (enzyme plus Mg-ATP), whereas with arginine alone only a minor change could be observed. The primary change in structure as seen with monomeric AK seems to be a Mg-nucleotide-induced domain movement relative to each other, whereas the effect of substrate may be of local order only. In CK, however, additional movements have to be involved.