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.

Dipole interaction model predicted pi-pi* circular dichroism of cyclo(L-Pro)3 using structures created by semi-empirical, ab initio, and molecular mechanics methods.

Cyclo(l-Pro)3 (CP3) is a synthetic peptide created to model cis and torsionally strained peptide bonds that also exhibits a strong distinctive UV circular dichroic (CD) spectrum. Circular dichroic spectra were computed for the amide pi-pi* transition using the dipole interaction model for various conformations of the peptide. Conformations of CP3 were created initially from crystal data, and followed by energy minimizations via molecular mechanics using the cvff force field; the effects of additional geometric optimizations by semi-empirical and ab initio quantum mechanics were investigated. The CD spectra for each conformation were calculated using a variety of different parameters, and each result was compared with the published experimental spectrum [Deber, C.M., Scatturin, A., Vaidya, V.M. & Blout, E.R. (1970) Small cyclic proline peptides: UV absorption and CD. In: Peptides: Chemistry and Biochemistry, Proceedings of the First American Peptide Symposium (Weinstein, B., ed.), Marcel Dekker, New York pp. 163-173]. Herein, two distinct conformations, a C3 symmetric and an asymmetric form, gave CD predictions that separately did not resemble the experimental spectrum. Energy differences were predicted at various theoretical levels, including MP2 and density functional theory. When the predicted CD spectra for each conformation were multiplied by Boltzmann weighting factors created using heats of formation determined by the AM1 optimizations, the weighted composite CD spectrum created did resemble experiment for the pi-pi* transition indicating that both conformations may exist simultaneously in solution.

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.

Force-field dependence of Boltzmann weighting factors on predicted pi-pi* circular dichroic spectra of cyclo (Gly-Pro-Gly)2.

Semi-empirical energy calculations were performed for published conformations of cyclo (Gly-Pro-Gly)2 using different force fields (DISCOVER cvff and cff91, AMBER, and CHARMM). The resulting potential energies were then used to create Boltzmann weighting factors for an ensemble of cyclo(Gly-Pro-Gly)2 structures. The dipole interaction model was used to predict pi-pi* circular dichroic spectra (CD) for the individual structures of cyclo(Gly-Pro-Gly)2. The Boltzmann weighting factors were applied to the individual spectra so that a composite spectrum was constructed to represent a CD arising from a collection of different structures in solution. Weighting factors determined from different force fields were compared. Boltzmann-weighted spectra better resembled the experimental CD than any calculated spectrum using only a single conformation of cyclo (Gly-Pro-Gly)2. The structures most heavily weighted contained at least one type I beta-turn.

Effects of charged amino acid mutations on the bimolecular kinetics of reduction of yeast iso-1-ferricytochrome c by bovine ferrocytochrome b5.

The reduction of wild-type yeast iso-1-ferricytochrome c (ycytc) and several mutants by trypsin-solubilized bovine liver ferrocytochrome b5 (cytb5) has been studied under conditions in which the electron-transfer reaction is bimolecular. The effect of electrostatic charge modifications and steric changes on the kinetics has been determined by experimental and theoretical observations of the electron-transfer rates of ycytc mutants K79A, K'72A, K79A/K'72A, and R38A (K' is used to signify trimethyllysine (Tml)). A structurally robust Brownian dynamics (BD) method simulating diffusional docking and electron transfer was employed to predict the mutation effect on the rate constants. A realistic model of the electron-transfer event embodied in an intrinsic unimolecular rate constant is used which varies exponentially with donor-acceptor distance. The BD method quantitatively predicts rate constants over a considerable range of ionic strengths. Semiquantitative agreement is obtained in predicting the perturbing influence of the mutations on the rate constants. Both the experimentally observed rate constants and those predicted by BD descend in the following order: native ycytc > K79A > K'72A > K79A/K'72A. Variant R38A was studied at a different ionic strength than this series of mutations, and the theory agreed with experiment in predicting a smaller rate constant for the mutant. In all cases the predicted effect of mutation was in the correct direction, but not as large as that observed. The BD simulations predict that the two proteins dock through essentially a single domain, with a distance of closest approach of the two heme groups in rigid body docking typically around 12 A. Two predominant classes of complexes were calculated, the most frequent involving the quartet of cytb5/ycytc interactions, Glu48-Arg13, Glu56-Lys87, Asp60-Lys86, and heme-Tml72, having an average electrostatic energy of -13.0 kcal/mol. The second most important complexes were of the type previously postulated (Salemme, 1976; Mauk et al., 1986; Rodgers et al., 1988) with interactions Glu44-Lys27, Glu48-Arg13, Asp60-Tml72, and heme-Lys79 and having an energy of -6.4 kcal/mol. The ionic strength dependence of the bimolecular reaction rate was well reproduced using a discontinuous dielectric model, but poorly so for a uniform dielectric model.

Effects of proline ring conformation on theoretical pi-pi* absorption and CD spectra of helical poly(L-proline) forms I and II.

Absorption and CD spectra of the pi-pi* transition near 200 nm are calculated for helical (Pro)10 forms I and II with a variable proline ring conformation characterized by torsion angle chi 2 in the range -60 degrees to 60 degrees. The spectra for poly(Pro) I are not sufficiently sensitive to chi 2 to suggest a preferred ring conformation. The spectra for poly(Pro) II are more sensitive to chi 2, and suggest preferred ring conformations near either or both of the chi 2 regions -50 +/- 10 degrees and 50 +/- 10 degrees.

Bond-optimized ring closure for proline: comparison of conformations and semiempirical energies with small molecule X-ray structures.

A method is described for generating proline ring structures by successive addition of atoms, wherein ring closure is achieved by optimizing the fit to known ring bond-angles and one closing bond-length ("bond-optimized ring closure"). Two ring torsion angles are fixed independently within broad, allowed ranges, and the remaining torsion angles are determined uniquely in most cases. The independent torsion angles are chosen as phi and chi 2, and ring closure is achieved without prohibitive strain through most of the ranges -130 degrees less than phi less than -20 degrees and -60 degrees less than chi 2 less than 60 degrees. Comparisons of predicted ring structures to 191 X-ray diffraction structures from the literature, starting with the known values of phi and chi 2, yielded root-mean-square deviations of 4.8 degrees in chi 1, 4.7 degrees in chi 3, 8.3 degrees in chi 4, and 0.3-2% in the ring bond angles and the N-C delta distance. Semiempirical energies were calculated for the optimized structures using three sets of energy parameters from the literature. The energy surfaces show broad minima coinciding with the torsion angle regions in which the highest concentrations of observed structures are found. Two of the sets of energy parameters produce double minima corresponding to the "up" and "down" puckered conformations.