Brownian dynamics simulations of glycolytic enzyme subsets with F-actin.

Previous Brownian dynamics (BD) simulations identified specific basic residues on fructose-1,6-bisphophate aldolase (aldolase) (I. V. Ouporov et al., Biophysical Journal, 1999, Vol. 76, pp. 17-27) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) (I. V. Ouporov et al., Journal of Molecular Recognition, 2001, Vol. 14, pp. 29-41) involved in binding F-actin, and suggested that the quaternary structure of the enzymes may be important. Herein, BD simulations of F-actin binding by enzyme dimers or peptides matching particular sequences of the enzyme and the intact enzyme triose phosphate isomerase (TIM) are compared. BD confirms the experimental observation that TIM has little affinity for F-actin. For aldolase, the critical residues identified by BD are found in surface grooves, formed by subunits A/D and B/C, where they face like residues of the neighboring subunit enhancing their electrostatic potentials. BD simulations between F-actin and aldolase A/D dimers give results similar to the native tetramer. Aldolase A/B dimers form complexes involving residues that are buried in the native structure and are energetically weaker; these results support the importance of quaternary structure for aldolase. GAPDH, however, placed the critical residues on the corners of the tetramer so there is no enhancement of the electrostatic potential between the subunits. Simulations using GAPDH dimers composed of either S/H or G/H subunits show reduced binding energetics compared to the tetramer, but for both dimers, the sets of residues involved in binding are similar to those found for the native tetramer. BD simulations using either aldolase or GAPDH peptides that bind F-actin experimentally show complex formation. The GAPDH peptide bound to the same F-actin domain as did the intact tetramer; however, unlike the tetramer, the aldolase peptide lacked specificity for binding a single F-actin domain.

Characterization of antigenic properties of horseradish peroxidase with monoclonal antibodies.

Monoclonal antibodies (mcAbs) specific to alkaline isoenzymes of horseradish peroxidase were used to characterize the antigenic properties of horseradish peroxidase. The results of a competitive binding assay indicated that monoclonal antibodies can be divided into three groups directed against distinct parts of the protein. The interaction of monoclonal antibodies with native and modified horseradish peroxidase showed also three different patterns of reactivity. Antibodies from groups I and II are directed against epitopes which are conformational and formed by tertiary structure elements. Epitopes recognized by these antibodies are sensitive to heme removal or partial denaturation of peroxidase. Antibodies from group III bind specifically with epitopes consisting of primary or secondary structure elements. The antigenic determinants recognized by antibodies from group III PO(1) and 36F(9) were shown to be linear (continuous) and formed by amino acid residues 261-267 and 271-277, respectively, as determined by the peptide scanning method (PEPSCAN). The location of revealed linear antigenic determinants in the molecular structure of peroxidase is analyzed.

Brownian dynamics simulations of aldolase binding glyceraldehyde 3-phosphate dehydrogenase and the possibility of substrate channeling.

Brownian dynamics (BD) simulations test for channeling of the substrate, glyceraldehyde 3-phosphate (GAP), as it passes between the enzymes fructose-1,6-bisphosphate aldolase (aldolase) and glyceraldehyde 3-phosphate dehydrogenase (GAPDH). First, BD simulations determined the favorable complexes between aldolase and GAPDH; two adjacent subunits of GAPDH form salt bridges with two subunits of aldolase. These intermolecular contacts provide a strong electrostatic interaction between the enzymes. Second, BD simulates GAP moving out of the active site of the A or D aldolase subunit and entering any of the four active sites of GAPDH. The efficiency of transfer is determined as the relative number of BD trajectories that reached any active site of GAPDH. The distribution functions of the transfer time were calculated based on the duration of successful trajectories. BD simulations of the GAP binding from solution to aldolase/GAPDH complex were compared to the channeling simulations. The efficiency of transfer of GAP within an aldolase/GAPDH complex was 2 to 3% compared to 1.3% when GAP was binding to GAPDH from solution. There is a preference for GAP channeling between aldolase and GAPDH when compared to binding from solution. However, this preference is not large enough to be considered as a theoretical proof of channeling between these proteins.

Interactions of glyceraldehyde-3-phosphate dehydrogenase with G- and F-actin predicted by Brownian dynamics.

Brownian dynamics (BD) was used to simulate the binding of glyceraldehyde-3-phosphate dehydrogenase (GAPDH) to G- and F-actin. High-resolution three-dimensional models (X-ray and homology built) of the proteins were used in the simulations. The electrostatic potential about each protein was predicted by solving the linearized Poisson-Boltzmann equation for use in BD simulations. The BD simulations resulted in complexes of GAPDH with G- or F-actin involving positively charged surface patches on GAPDH (Lyses 24, 69, 110 and 114) and negatively charged residues of the N- and C-termini (Asps 1, 25 and 363 and Glus 2, 4, 224 and 364) of actin. The actin residues all belong to subdomain 1. Although the positively charged surface patches of GAPDH are not close enough to each other to enhance their electrostatic potential, occasionally two subunits of the GAPDH tetramer may simultaneously interact with two neighboring monomers of F-actin. These results are different from those of fructose-1,6-bisphosphate aldolase, where quaternary structure directly influenced binding by two subunits combining their electrostatic potentials (see previous study, Ouporov et al., 1999, Biophys. J. 76: 17-27). Instead, GAPDH uses its quaternary structure to span the distance between two different actin subunits so that it can interact with two different actin subunits simultaneously.

Computer simulations of glycolytic enzyme interactions with F-actin.

Muscle actin and fructose-1,6-bisphosphate aldolase (aldolase) were chemically crosslinked to produce an 80 kDa product representing one subunit of aldolase linked to one subunit of actin. Hydroxylamine digestion of the crosslinked product resulted in two 40.5 kDa fragments, one that was aldolase linked to the 12 N-terminal residues of actin. Brownian dynamics simulations of muscle aldolase and GAPDH with F-actin (muscle, yeast, and various mutants) estimated the association free energy. Mutations of residues 1-4 of muscle actin to Ala individually or two in combination of the first four residues reduced the estimated binding free energy. Simulations showed that muscle aldolase binds with the same affinity to the yeast actin as to the double mutated muscle actin; these mutations make the N-terminal of muscle actin identical to yeast, supporting the conclusion that the actin N-terminus participates in binding. Because the depth of free energy wells for yeast and the double mutants is less than for native rabbit actin, the simulations support experimental findings that muscle aldolase and GAPDH have a higher affinity for muscle actin than for yeast actin. Furthermore, Brownian dynamics revealed that the lower affinity of yeast actin for aldolase and GAPDH compared to muscle actin, was directly related to the acidic residues at the N-terminus of actin.

Brownian dynamics simulations of interactions between aldolase and G- or F-actin.

Compartmentation of proteins in cells is important to proper cell function. Interactions of F-actin and glycolytic enzymes is one mechanism by which glycolytic enzymes can compartment. Brownian dynamics (BD) simulations of the binding of the muscle form of the glycolytic enzyme fructose-1,6-bisphosphate aldolase (aldolase) to F- or G-actin provide first-encounter snapshots of these interactions. Using x-ray structures of aldolase, G-actin, and three-dimensional models of F-actin, the electrostatic potential about each protein was predicted by solving the linearized Poisson-Boltzmann equation for use in BD simulations. The BD simulations provided solution complexes of aldolase with F- or G-actin. All complexes demonstrate the close contacts between oppositely charged regions of the protein surfaces. Positively charged surface regions of aldolase (residues Lys 13, 27, 288, 293, and 341 and Arg 257) are attracted to the negatively charged amino terminus (Asp 1 and Glu 2 and 4) and other patches (Asp 24, 25, and 363 and Glu 361, 364, 99, and 100) of actin subunits. According to BD results, the most important factor for aldolase binding to actin is the quaternary structure of aldolase and actin. Two pairs of adjacent aldolase subunits greatly add to the positive electrostatic potential of each other creating a region of attraction for the negatively charged subdomain 1 of the actin subunit that is exposed to solvent in the quaternary F-actin structure.

Evidence for indole-3-acetic acid binding site in plant peroxidases. Structural similarity between peroxidases and auxin-binding proteins.

Application of computer methods allowed us to demonstrate that plant peroxidases and auxin-binding proteins contain structurally similar fragments. The mapping of the fragments was done using a model structure of horseradish peroxidase. Five of six structurally similar fragments belong to the distal domain and form a subdomain in plant peroxidases that includes the distal heme-coordinating sequence, LHFHDC (amino acid residues 39-44 in horseradish peroxidase). The existence of a substrate-binding site for indole-3-acetic acid in the distal subdomain comprising helices A (whole), B (middle), C (beginning), and D (whole) and the loop between helices D and D' is discussed.

Effect of pH on tobacco anionic peroxidase stability and its interaction with hydrogen peroxide.

The effect of extremely acidic pH on the stability of tobacco peroxidase and lignin peroxidase holoenzymes has been studied. Stabilization of tobacco peroxidase holoenzyme in the presence of calcium cations at pH < 2 and stabilization of lignin peroxidase at pH > 2 in the presence of veratryl alcohol have been shown. The dependence of the reaction rate constant for hydrogen peroxide interaction with tobacco peroxidase on pH suggests that the reaction rate is under control of a group with pK of 2.5. A tobacco peroxidase model structure has been created by means of homology modeling on the basis of the tobacco peroxidase sequence and the coordinates of peanut peroxidase crystal structure. The model structure demonstrates the presence of the negatively charged Glu-141 at the entrance to the active site and its electrostatic repulsion from heme propionates and triad of Asp-76, -79, and -80 residues. The results on tobacco holoperoxidase stabilization at pH 1.8 in the presence of calcium cations and drop in reaction rate constant for the enzyme interaction with hydrogen peroxide are explained by a hypothetical formation of ionic bonds between Glu-141 and the triad of aspartic acid residues via calcium cation lowering the accessibility of the active site and stabilizing the holoenzyme.

Epitope mapping of horseradish peroxidase (isoenzyme C).

Peptide scanning (PEPSCAN) was used to determine linear antigenic determinants of horseradish peroxidase isoenzyme C (HRPC). For this purpose, we synthesized 303 overlapping hexapeptide fragments (with a step of one amino acid residue) of the protein primary structure and studied their interactions with anti-HRPC polyclonal antisera by ELISA. Experiments with various titers of antisera allowed us to determine linear antigenic determinants of HRPC; several such determinants were spatially located in regular elements of the secondary structure (alpha-helices) found both inside and outside the protein globule. A fraction of epitopes were located in loops and folds of the HRPC peptide chain with irregular shapes. These epitopes contained several functionally important residues: Arg 38, which is part of the active site of the enzyme, as well as Phe 142 and Phe 143, which form a channel allowing aromatic substrates to reach the active site. Amino acid residues that form calcium-binding sites or occur in the vicinity of disulfide bonds are not involved in these epitopes.

Formation and properties of dimeric recombinant horseradish peroxidase in a system of reversed micelles.

Wild-type recombinant horseradish peroxidase purified and refolded from Escherichia coli inclusion bodies has been studied in the system of bis(2-ethylhexyl)sulphosuccinate sodium salt (Aerosol OT)-reversed micelles in octane. In contrast with native horseradish peroxidase the wild-type recombinant enzyme forms dimeric structures as judged by sedimentation analysis. Peroxidase substrates affect the equilibrium between monomeric and dimeric enzyme forms. The dependence of the catalytic activity of recombinant peroxidase on the degree of hydration of the surfactant exhibits two maxima with pyrogallol, o-phenylene- diamine, guaiacol and o-dianisidine, with different ratios of activities for the first and second maxima. The differences in activities of monomeric and dimeric forms of the recombinant horseradish peroxidase provide evidence for active-site screening in dimeric forms. This has been used to model a dimeric structure of recombinant horseradish peroxidase with the screened entrance to the active site. In the model structure obtained, three of eight glycosylation sites were screened. This might explain the absence of dimeric structures in native enzyme peroxidase. The system of reversed micelles provides, for the first time, evidence for the formation of dimeric structures by recombinant plant peroxidase with an altered substrate specificity compared with the native enzyme. Thus one can assume that haem-containing peroxidases in general are able to form dimeric structures.

Dimerization of recombinant horseradish peroxidase in a reversed micellar system.

Recombinant horseradish peroxidase reactivated from E. coli inclusion bodies was studied in a reversed micellar system of AOT in octane. The ability of the recombinant enzyme, in contrast to native horseradish peroxidase, to form a dimeric structure was found. The existence of the dimer was proved by results of sedimentation analysis. Dimer/monomer ratio in the enzyme-containing micelles and dimer catalytic activity were found to depend on the substrate used (pyrogallol, guaiacol, o-dianisidine, o-phenylenediamine). Computer modelling was used to describe possible structures of the dimeric recombinant horseradish peroxidase.

Helical stacking in DNA three-way junctions containing two unpaired pyrimidines: proton NMR studies.

The proton NMR spectra of DNA three-way junction complexes (TWJ) having unpaired pyrimidines, 5'-TT- and 5'-TC- on one strand at the junction site were assigned from 2D NOESY spectra acquired in H2O and D2O solvents and homonuclear 3D NOESY-TOCSY and 3D NOESY-NOESY in D2O solvent. TWJ are the simplest branched structures found in biologically active nucleic acids. Unpaired nucleotides are common features of such structures and have been shown to stabilize junction formation. The NMR data confirm that the component oligonucleotides assemble to form conformationally homogeneous TWJ complexes having three double-helical, B-form arms. Two of the helical arms stack upon each other. The unpaired pyrimidine bases lie in the minor groove of one of the helices and are partly exposed to solvent. The coaxial stacking arrangement deduced is different from that determined by Rosen and Patel (Rosen, M.A., and D.J. Patel. 1993. Biochemistry. 32:6576-6587) for a DNA three-way junction having two unpaired cytosines, but identical to that suggested by Welch et al. (Welch, J. B., D. R. Duckett, D. M. J. Lilley. 1993. Nucleic Acids Res. 21:4548-4555) on the basis of gel electrophoretic studies of DNA three-way junctions containing unpaired adenosines and thymidines.

Refinement of the solution structure of a branched DNA three-way junction.

We have refined the structure of the DNA Three-Way Junction complex, TWJ-TC, described in the companion paper by quantitative analysis of two 2D NOESY spectra (mixing times 60 and 200 ms) obtained in D2O solution. NOESY crosspeak intensities extracted from these spectra were used in two kinds of refinement procedure: 1) distance-restrained energy minimization (EM) and molecular dynamics (MD) and 2) full relaxation matrix back calculation refinement. The global geometry of the refined model is very similar to that of a published, preliminary model (Leontis, 1993). Two of the helical arms of the junction are stacked. These are Helix 1, defined by basepairs S1-G1/S3-C12 through S1-C5/S3-G8 and Helix 2, which comprises basepairs S1-C6/S2-G5 through S1-G10/S2-G1. The third helical arm (Helix 3), comprised of basepairs S2-C6/S3-G5 through S2-C10/S3-G1 extends almost perpendicularly from the axis defined by Helices 1 and 2. The bases S1-C5 and S1-C6 of Strand 1 are continuously stacked across the junction region. The conformation of this strand is close to that of B-form DNA along its entire length, including the S1-C5 to S1-C6 dinucleotide step at the junction. The two unpaired bases S3-T6 and S3-C7 lie outside of the junction along the minor groove of Helix 1 and largely exposed to solvent. Analysis of the refined structure reveals that the glycosidic bond of S3-T6 exists in the syn conformation, allowing the methyl group of this residue to contact the hydrophobic surface of the minor groove of Helix 1, at S3-G11. The helical parameters of the three helical arms of the structure exhibit only weak deviations from typical values for right-handed B-form DNA. Unusual dihedral angles are only observed for the sugarphosphate backbone joining the "hinge" residues, S2-G5 and S2-C6, and S3-G5 through S3-G8. The glycosidic bond of S3-G8also lies within the syn range, allowing favorable Watson-Crick base-pairing interactions with Si -C5. The stability of this structure was checked in 39 ps molecular dynamic simulation at 330 K in water. The structure of TWJ-TC retained the geometrical features mentioned above at the end of the simulation period. The final R(1/6)-factor of the refined structure is 5%.