Polypeptide matrices as active catalytic supports for stereoselective electron-transfer reactions

The oxidation of L -adrenaline (epinephrine) in the presence of [Fe(tetpy)(OH)₂]⁺ ions bound to poly(L -glutamate) or to poly(D-glutamate) has been studied at pH 7 (tetpy = 2,2′:6′,2″:6″,2?-tetrapyridyl). Electron transfer from the substrate to the central metal ion, which is rate-determining, proceeds stereoselectively only when extensive and possibly specific interactions between adrenaline and the peptidic residues of the ordered polymer in the close environment of the active sites occur. This ensures different steric constraints for the two diastereomeric precursor complexes, which are thought to affect the separation and orientation of the redox centers differently, leading to the observed phenomena. Some data on the catalytic oxidation of L -dopa(3,4-dihydroxyphenylalanine) are also presented, showing stereoselective effects similar to those observed with L -adrenaline, despite the diverse distance of the chiral center from the reacting OH groups. A mechanistic interpretation of the results is discussed in the light of a few general considerations concerning the structural features of the catalytic systems. Possible explanations for the finding that stereoselectivity occurs at the expense of the efficiency of catalysis are also considered.     

Energetics of formation of electron-transfer complexes between asymmetric species

The energetics of formation of diastereomeric electron-transfer complexes between L-dopa or L-adrenaline and iron(III) complex ions bound to poly(L-glutamate) or poly(D-glutamate) were measured by microcalorimetry at 25°C and pH 7. When the association of substrates by Fe-polypeptide systems was virtualy complete, ΔH_ll = 1.3 ± 0.1 and ΔH_dl = 0.9 ± 0.1 kcal/mol with both substrates. A diastereomeric discrimination energy of 400 ± 100 cal/mol is thus observed, which compares satisfactorily with that directely obtained by differential microcalorimetric measurements. These results are fully consistent with previous findings indicating that thermodynamic effects are of minor importance in the overall stereoselectivity of the electron transfer reactions investigated.     

Chiral discrimination in the energetics of formation of diastereomeric adducts involving polypeptides

Chiral discrimination in the energetics of formation of diastereomeric pairs between iron (III)-complex ions bound to poly(L -glutamate) or poly (D -glutamate) and L -catecholamines was measured by differential microcalorimetry at 25°C. When the association of substrates by Fe–polypeptide systems was virtually complete, diastereomeric discrimination energies of a few hundred calories per mol were observed, while diastereomeric discrimination entropies were found to be negligibly small. These results are consistent with the finding that the overall stereoselectivity in the electron transfer reaction between the optically active catecholamines and enantiomeric iron(III) materials—corresponding to Δ(ΔG^≠) values around 900 cal/mol—is largely controlled by transition state effects in the kinetics.     

Stereospecific oxidation of ascorbic acid by chiral association complexes between polypeptides and iron(III) chelate ions

2,2,′,2″,2?-Tetrapyridineiron(III) complex ions anchored to poly(L -glutamate) (FeL) or poly(D -glutamate) (FeD) were used as catalysts for the H₂O₂ oxidation of L (+)-ascorbic acid at pH 7 and varying complex:polymer-residue molar ratios [C]/[P]. Evidence is produced that the reaction is a composite process reflecting contributions from parallel pathways, one of which corresponds to a catalytic route and is [H₂O₂]-independent and the other to an uncatalyzed electron-transfer process between the ascorbate anion and hydrogen peroxide. Stereospecific effects in the catalysis are observed on increasing the complex:polymer ratio, which corresponds to an increase of the amount of α-helical fraction in the polypeptide supports (x_a). Thus, at [C]/[P] = 0.01 (x_a < 0.05), k_FeD/k_FeL = 1.0; but at [C]/[P] = 0.20 (x_a ≈ 0.70), k_FeD/k_FeL = 4.0 ± 0.5, k_FeD and k_FeL being the second-order rate constants of the electron-transfer reaction between the FeD or FeL isomer of the asymmetric catalyst and the L -ascorbate anion. The activation energies were found to increase markedly on going from the former to the latter complex:polymer ratio but, at the same time, to exhibit equal values with both enantiomeric catalysts. Stereoselectivity therefore appears to be an entropy-controlled process, arising from the conformational rigidity of the precursor complex, which very likely sees the substrate molecules bound to the chiral residues of the ordered polymer surrounding the active sites. The implications of the stereochemical features of the substrate–catalyst adduct on the mechanism of electron transfer are also discussed. Evidence is presented that the asymmetric [Fe(tetpy)(OH)₂]⁺–polyelectrolyte systems play the additional role of environmental controller of the uncatalyzed oxidation of the L -ascorbate anion.     

Parallel DNA double helices. II. Conformational analysis of regular helices having a second order axis of symmetry

Conformational analysis of double helices of DNA with parallel arranged sugar-phosphate chains connected by twofold symmetry has been performed. Homopolymers poly(dA).poly(dA), poly(dC).poly(dC), poly(dG).poly(dG) and poly(dT).poly(dT) were studied. For each of the homopolymers all variants of H-bond pairing were checked. The maps of closing of sugar-phosphate backbone were previously computed. By the optimization of potential energy the dihedral angles and helix parameters of relatively stable conformations of parallel stranded polynucleotides were calculated. The dependence of conformational energy on the nucleic base character and the base pair type were studied. Two main conformational regions for favourable "parallel" helix of polynucleotides were found. The former of these two regions coincide with the region of typical conformational parameters of B-DNA. On an average the conformational energy of "parallel" DNA is close to the energy of canonic "antiparallel" B-DNA.     

Parallel double helices of DNA. Conformational analysis of regular helices with the second order symmetry axis

We have performed a conformational analysis of DNA double helices with parallel directed backbone strands connected with the second order symmetry axis being at the same time the helix axis. The calculations were made for homopolymers poly(dA).poly(dA), poly(dC).poly(dC), poly(dG) poly(dG), and poly(dT).poly(dT). All possible variants of hydrogen bonding of base pairs of the same name were studied for each polymer. The maps of backbone chain geometrical existence were constructed. Conformational and helical parameters corresponding to local minima of conformational energy of "parallel" DNA helices, calculated at atom-atom approximation, were determined. The dependence of conformational energy on the base pair and on the hydrogen bond type was analysed. Two major conformational advantageous for "parallel" DNA's do not depend much on the hydrogen-bonded base pair type were indicated. One of them coincided with the conformational region typical for "antiparallel" DNA, in particular for the B-form DNA. Conformational energy of "parallel" DNA depends on the base pair type and for the most part is similar to the conformational energy of "antiparallel" B-DNA.     

Volume and hydration changes of DNA-ligand interactions

We report the volumetric and other thermodynamic properties of ethidium bromide (EB), propidium iodide (PI) and daunomycin (DAU) intercalating with poly(dA).poly(dT), poly[d(A-T)].poly[d(A-T)], and poly[d(G-C)].poly[d(G-C)], respectively, as well as minor groove binder Hoechst 33258 binding with poly[d(A-T)].poly[d(A-T)]. The data were obtained using fluorescence titration and hydrostatic pressure measurements. Our thermodynamic data are combined with enthalpies from literature reports to analyze the thermodynamic characteristics of the different interactions. The differences are interpreted based on three processes related to hydration: I. burial of non-polar hydrophobic solvent accessible surface, II. burial of polar surface and formation of solute-solute H-bonds, and III. disruption of "structural" hydration. Sequence dependent conformational changes may also be important when comparing ligand binding to different DNA sequences. We conclude that a combination of different thermodynamic parameters, especially volume change, is essential in order to understand the role of hydration in the energetics of DNA-ligand interactions.     

Equilibrium analyses of the active-site asymmetry in enterococcal NADH oxidase: role of the cysteine-sulfenic acid redox center

Recent studies [Mallett, T. C., and Claiborne, A. (1998) Biochemistry 37, 8790-8802] of the O2 reactivity of C42S NADH oxidase (O2 --> H2O2) revealed an asymmetric mechanism in which the two FADH2.NAD+ per reduced dimer display kinetic inequivalence. In this report we provide evidence indicating that the fully active, recombinant wild-type oxidase (O2 --> 2H2O) displays thermodynamic inequivalence between the two active sites per dimer. Using NADPH to generate the free reduced wild-type enzyme (EH2'/EH4), we have shown that NAD+ titrations lead to differential behavior as only one FADH2 per dimer binds NAD+ tightly to give the charge-transfer complex. The second FADH2, in contrast, transfers its electrons to the single Cys42-sulfenic acid (Cys42-SOH) redox center, which remains oxidized during the reductive titration. Titrations of the reduced NADH oxidase with oxidized 3-acetylpyridine and 3-aminopyridine adenine dinucleotides further support the conclusion that the two FADH2 per dimer in wild-type enzyme can be described as distinct "charge-transfer" and "electron-transfer" sites, with the latter site giving rise to either intramolecular (Cys42-SOH) or bimolecular (pyridine nucleotide) reduction. The reduced C42S mutant is not capable of intramolecular electron transfer on binding pyridine nucleotides, thus confirming that the Cys42-SOH center is in fact the source of the redox asymmetry observed with wild-type oxidase. These observations on the role of Cys42-SOH in the expression of thermodynamic inequivalence as observed in wild-type NADH oxidase complement the previously described kinetic inequivalence of the C42S mutant; taken together, these results provide the overlapping framework for an alternating sites cooperativity model of oxidase action.     

Stereoselective catalysis by Fe(III) complex ions supported on asymmetric polymers: Oxidation of ascorbic acid in aqueous solution

Quaterpyridyneiron (III) complex ions anchored to partially ordered poly (L-glutamate) or poly (D-glutamate) were used as (enantiomeric) catalysts for the H₂O₂-oxidation of L(+) ascorbic acid at pH 7. When the α-helical fraction of polypeptide matrices was low, the configuration dissymmetry of the active sites was unable to impart any stereoselective effect in the catalysis, i.e. k = 3.66 x 10³ M^?1?sec^?1 (25.9°C) with both catalysts. On the contrary, by increasing the amount of α-helix in the polymeric supports the stereoselectivity increases, the second-order rate constants k_FeD being definitely higher than k_FeL.Implications of the role played by the conformational dissymmetry of the active sites in the stereospecificity of the process are briefly discussed.     

Stereoselectivity of the epoxidation of 7,8-dihydrobenzo[a]pyrene by prostaglandin H synthase and cytochrome P-450 determined by the identification of polyguanylic acid adducts

The stereoselectivity of the oxidation of 7,8-dihydrobenzo[a]pyrene (H2BP) to 9,10-epoxy-7,8,9,10-tetrahydrobenzo[a]pyrene (H4BP-epoxide) by prostaglandin H (PGH) synthase and cytochrome P-450 has been studied using microsomal preparations from ram seminal vesicles and rat liver. Incubations were performed in the presence of polyguanylic acid and the adducts formed between H4BP-epoxide and guanosine were isolated following the recovery and hydrolysis of the poly(G). When (+/-)-H4BP-epoxide was reacted with poly(G), four diastereomeric adducts were formed by the cis and trans addition of the exocyclic amino group of guanine to the benzylic carbon of the epoxide enantiomers. Each diastereomer was identified by a combination of ultraviolet, nuclear magnetic resonance, circular dichroism, and mass spectroscopy. Under comparable conditions, ram seminal vesicle microsomes in the presence of arachidonic acid triggered the binding of H2BP to poly(G) to a greater extent than rat liver microsomes from untreated and phenobarbital- and methylcholanthrene pretreated animals in the presence of NADPH. Quantitation of the (-)-cis- and (+)-cis-guanosine adducts revealed the degree of stereoselectivity of epoxidation. The ratio of (-)/(+) adducts was 54:46 for PGH synthase and 89:11 (control), 62:38 (phenobarbital), and 69:31 (methylcholanthrene) for cytochrome P-450-catalyzed reactions. PGH synthase catalyzed the epoxidation of H2BP with little or no stereoselectivity in contrast to cytochrome P-450. The utility of the poly(G) binding technique for the elucidation of the stereoselective generation of chiral electrophiles is discussed along with the mechanistic implications of the results.     

On the use of chiral compounds for probing the DNA handedness: Z to B conversion in poly(dGm5dC) upon binding of Fe(phen)3(2+) and Ru(phen)3(2+)

In order to examine whether chiral metal complexes can be used to discriminate between right- and left-handed DNA conformational states we have studied the enantioselective interactions of Fe(phen)3(2+) and Ru(phen)3(2+) (phen = 1,10-phenanthroline) with poly(dGm5dC) under B- and Z-form conditions. With the inversion-labile Fe(phen)3(2+), enantioselectivity leads to shifts in the diastereomeric binding equilibria. This effect, known as the "Pfeiffer effect" (1-4), is monitored as slowly emerging circular dichroism of the solution, corresponding to a net excess of the favoured enantiomer. With Ru(phen)3(2+), which is stable to intramolecular inversion, the difference in DNA-binding strengths of the enantiomers results in an excess of the less favoured enantiomer in the bulk solution. This excess is detected in the dialysate of the DNA/metal complex solution. With both complexes we find that the delta-enantiomer is favoured when the polynucleotide adopts the B-form, as previously shown, but also when it initially adopts the Z-form conformational state. This observation, together with evidence from UV-circular dichroism and binding data, indicates that the binding of these metal complexes induces a Z- to B-form transition in Z-form poly(dGm5dC). Consequently, neither of the studied chiral DNA-binders can easily be used to discriminate the DNA handedness.     

Synchronous (synphaseous) additive models for mitochondrial ATP synthetase complex and for other energy transducing systems of cells (a hypothesis)

There exist some energy transducing enzymes containing immobile molecules of substrates which are not exchangeable during protein functioning. According to the proposed models the immobile substrates are localized at the "idle" (or "partial" or "imitational") catalytic sites, which differ from normal ("working") active sites of enzymes. Only some steps of a complete reaction sequence which take place at the "working" sites are carried out at the "idle" sites. On the other hand, cyclic conversion of the immobile substrate at an "idle" catalytic site may include some steps which are absent in the "working" site cycle. The occurrence of identical steps on the "idle" and "working" catalytic sites allows to synchronize their action through conformational interconversions of tightly packed and structurally related "idle" and "working" subunits of the enzyme. The presence of covaletly bound substrates or substrates localized in closed cavities of the "idle" sites allows to synchronize the action of many monomers containing such sites due to the absence of the rate-limiting step of simultaneous saturation of many catalytic sites by substrate molecules from solution, and due to the lack of substrate inhibition on the "idle" sites. The functions of the "idle" sites are miscellaneous e.g. in ion-transporting systems these sites are directly involved in ion translocation. In the actomyosin complex the "idle" sites imitate conformational alterations of "working" sites, thus allowing synchronous functioning of the polymeric structure. Variations in the number of the "idle" sites operating simultaneously with one "working" site allow to regulate some parameters of enzymatic processes, e.g. the stoichiometry (number of transported ions per ATP hydrolysed (or synthesised) or electron-transported, or hv-absorbed ones) for ion transported systems or the ratio (velocity of contraction to developed efforts) for the actomyosin complex.     

Four-stranded DNA helices: conformational analysis of regular poly(dT).poly(dA).poly(dA).poly(dT) helices with various types of base binding

The paper presents results obtained in conformational analysis of homopolymeric four-stranded poly(dT).poly(dA).poly(dA).poly(dT) DNA helices in which the pairs of strands with identical bases are parallel and have a two-fold symmetry axis. All possible models of base binding to yield a symmetric complex have been considered. The dihedral angles of sugar-phosphate backbones and helix parameters, which are consistent with the minima of conformational energy for four-stranded DNAs, have been determined using the results of optimization of conformational energy calculated at atom-atom approximation. Potential energy is shown to depend on the structure of base complexes and on the mutual orientation of unlike strands. Possible biological functions of four-stranded helices are discussed.     

Crystal structures of MAP kinase p38 complexed to the docking sites on its nuclear substrate MEF2A and activator MKK3b

The structures of the MAP kinase p38 in complex with docking site peptides containing a phi(A)-X-phi(B) motif, derived from substrate MEF2A and activating enzyme MKK3b, have been solved. The peptides bind to the same site in the C-terminal domain of the kinase, which is both outside the active site and distinct from the "CD" domain previously implicated in docking site interactions. Mutational analysis on the interaction of p38 with the docking sites supports the crystallographic models and has uncovered two novel residues on the docking groove that are critical for binding. The two peptides induce similar large conformational changes local to the peptide binding groove. The peptides also induce unexpected and different conformational changes in the active site, as well as structural disorder in the phosphorylation lip.     

Reconciling the "old" and "new" views of protein allostery: a molecular simulation study of chemotaxis Y protein (CheY)

A combination of thirty-two 10-ns-scale molecular dynamics simulations were used to explore the coupling between conformational transition and phosphorylation in the bacteria chemotaxis Y protein (CheY), as a simple but representative example of protein allostery. Results from these simulations support an activation mechanism in which the beta4-alpha4 loop, at least partially, gates the isomerization of Tyr106. The roles of phosphorylation and the conserved Thr87 are deemed indirect in that they stabilize the active configuration of the beta4-alpha4 loop. The indirect role of the activation event (phosphorylation) and/or conserved residues in stabilizing, rather than causing, specific conformational transition is likely a feature in many signaling systems. The current analysis of CheY also helps to make clear that neither the "old" (induced fit) nor the "new" (population shift) views for protein allostery are complete, because they emphasize the kinetic (mechanistic) and thermodynamic aspects of allosteric transitions, respectively. In this regard, an issue that warrants further analysis concerns the interplay of concerted collective motion and sequential local structural changes in modulating cooperativity between distant sites in biomolecules.     

Induced fit or conformational selection? The role of the semi-closed state in the maltose binding protein

A full characterization of the thermodynamic forces underlying ligand-associated conformational changes in proteins is essential for understanding and manipulating diverse biological processes, including transport, signaling, and enzymatic activity. Recent experiments on the maltose binding protein (MBP) have provided valuable data about the different conformational states implicated in the ligand recognition process; however, a complete picture of the accessible pathways and the associated changes in free energy remains elusive. Here we describe results from advanced accelerated molecular dynamics (aMD) simulations, coupled with adaptively biased force (ABF) and thermodynamic integration (TI) free energy methods. The combination of approaches allows us to track the ligand recognition process on the microsecond time scale and provides a detailed characterization of the protein's dynamic and the relative energy of stable states. We find that an induced-fit (IF) mechanism is most likely and that a mechanism involving both a conformational selection (CS) step and an IF step is also possible. The complete recognition process is best viewed as a "Pac Man" type action where the ligand is initially localized to one domain and naturally occurring hinge-bending vibrations in the protein are able to assist the recognition process by increasing the chances of a favorable encounter with side chains on the other domain, leading to a population shift. This interpretation is consistent with experiments and provides new insight into the complex recognition mechanism. The methods employed here are able to describe IF and CS effects and provide formally rigorous means of computing free energy changes. As such, they are superior to conventional MD and flexible docking alone and hold great promise for future development and applications to drug discovery.     

Stereochemical aspects of cholinesterase catalysis

Proceeding from the sterochemical regularities of the nucleophilic substitution reaction at the carbonyl group and the assumption that the spatial structure of the active center of cholinesterases is complementary to the molecule of the ester substrates for these enzymes, some general features of the stereoselectivity phenomena in the reactions of acetylcholinesterase (EC 3.1.1.7) and butyrylcholinesterase (EC 3.1.1.8) with organophosphorus inhibitors are discussed. For these enzymes the models of the active center are proposed in terms of different binding sites and the catalytic center. On the basis of this model, the stereochemical pecularities and the physicochemical background of the stereoselectivity effects on enzyme inhibition, reactivation, and “aging” reactions can be understood. Knowledge of the absolute configuration of several chiral organophosphorus inhibitors also makes it possible to determine the absolute spatial arrangement of the hydrophobic binding sites on the active surface of cholinesterases.     

Stereoselectivity of lipases. II. Stereoselective hydrolysis of triglycerides by gastric and pancreatic lipases

In the present study, porcine pancreatic lipase, rabbit gastric lipase, and human gastric lipase stereospecificity toward chemically alike, but sterically nonequivalent ester groups within one single triglyceride molecule was investigated. Lipolysis reactions were carried out on synthetic trioctanoin or triolein, which are homogenous, prochiral triglycerides, chosen as models for physiological lipase substrates. Diglyceride mixtures resulting from lipolysis were derivatized with optically active R-(+)-1-phenylethylisocyanate, to give diastereomeric carbamate mixtures, which were further separated by high performance liquid chromatography. Resolution of diastereomeric carbamates gave enantiomeric excess values, which reflect the lipases stereobias and clearly demonstrate the existence of a stereopreference by both gastric lipases for the sn-3 position. The stereoselectivity of human and rabbit gastric lipases, expressed as the enantiomeric excess percentage, was 54% and 70% for trioctanoin and 74% and 47% for triolein, respectively. The corresponding values with porcine pancreatic lipase were 3% in the case of trioctanoin and 8% in that of triolein. It is worth noting that rabbit gastric lipase, unlike human gastric lipase, became more stereoselective for the triglyceride with shorter acyl chains (trioctanoin). This is one of the most striking catalytic differences observed between these two gastric lipases.     

Thermodynamic analysis shows conformational coupling and dynamics confer substrate specificity in fructose-1,6-bisphosphate aldolase

Conformational flexibility is emerging as a central theme in enzyme catalysis. Thus, identifying and characterizing enzyme dynamics are critical for understanding catalytic mechanisms. Herein, coupling analysis, which uses thermodynamic analysis to assess cooperativity and coupling between distal regions on an enzyme, is used to interrogate substrate specificity among fructose-1,6-(bis)phosphate aldolase (aldolase) isozymes. Aldolase exists as three isozymes, A, B, and C, distinguished by their unique substrate preferences despite the fact that the structures of the active sites of the three isozymes are nearly identical. While conformational flexibility has been observed in aldolase A, its function in the catalytic reaction of aldolase has not been demonstrated. To explore the role of conformational dynamics in substrate specificity, those residues associated with isozyme specificity (ISRs) were swapped and the resulting chimeras were subjected to steady-state kinetics. Thermodynamic analyses suggest cooperativity between a terminal surface patch (TSP) and a distal surface patch (DSP) of ISRs that are separated by >8.9 A. Notably, the coupling energy (DeltaGI) is anticorrelated with respect to the two substrates, fructose 1,6-bisphosphate and fructose 1-phosphate. The difference in coupling energy with respect to these two substrates accounts for approximately 70% of the energy difference for the ratio of kcat/Km for the two substrates between aldolase A and aldolase B. These nonadditive mutational effects between the TSP and DSP provide functional evidence that coupling interactions arising from conformational flexibility during catalysis are a major determinant of substrate specificity.     

Trapping conformational intermediate states in the reaction center protein from photosynthetic bacteria

In protein, conformational changes are often crucial for function but not easy to observe. Two functionally relevant conformational intermediate states of photosynthetic reaction center protein (RCs) are trapped and characterized at low temperature. RCs frozen in the dark do not allow electron transfer from the reduced primary quinone, Q(A)(-), to the secondary quinone, Q(B). In contrast, RCs frozen under illumination in the product (P(+)Q(A)Q(B)(-)) state, with the oxidized electron donor, P(+), and reduced Q(B)(-), return to the ground state at cryogenic temperature in a conformation that allows a high yield of Q(B) reduction. Thus, RCs frozen under illumination are found to be trapped above the ground state in a conformation that allows product formation. When the temperature is raised above 120 K, the protein relaxes to an inactive conformation which is different from the RCs frozen in the dark. The activation energy for this change is 87 +/- 8 meV, and the active and inactive states differ in energy by only 16 +/- 3 meV. Thus, there are several conformational substates along the reaction coordinate with different transition temperatures. The ground state spectra of the RCs in active and inactive conformations report differences in the intraprotein electrostatic field, demonstrating that the dipole or charge distribution has changed. In addition, the electrochromic shift associated with the Q(A)(-) to Q(B) electron transfer at low temperature was characterized. The electron-transfer rate from Q(B)(-) to P(+) was measured at cryogenic temperature and is similar to the rate at room temperature, as expected for an exothermic, electron tunneling reaction in RCs.