Powering DNA repair through substrate electrostatic interactions

The reaction catalyzed by the DNA repair enzyme uracil DNA glycosylase (UDG) proceeds through an unprecedented stepwise mechanism involving a positively charged oxacarbenium ion sugar and uracil anion leaving group. Here we use a novel approach to evaluate the catalytic contribution of electrostatic interactions between four essential phosphodiester groups of the DNA substrate and the cationic transition state. Our strategy was to substitute each of these phosphate groups with an uncharged (R)- or (S)-methylphosphonate linkage (MeP). We then compared the damaging effects of these methylphosphonate substitutions on catalysis with their damaging effects on binding of a cationic 1-azadeoxyribose (1-aza-dR(+)) oxacarbenium ion analogue to the UDG-uracil anion binary complex. A plot of log k(cat)/K(m) for the series of MeP-substituted substrates against log K(D) for binding of the 1-aza-dR(+) inhibitors gives a linear correlation of unit slope, confirming that the electronic features of the transition state resemble that of the 1-aza-dR(+), and that the anionic backbone of DNA is used in transition state stabilization. We estimate that all of the combined phosphodiester interactions with the substrate contribute 6-8 kcal/mol toward lowering the activation barrier, a stabilization that is significant compared to the 16 kcal/mol catalytic power of UDG. However, unlike groups of the enzyme that selectively stabilize the charged transition state by an estimated 7 kcal/mol, these phosphodiester groups also interact strongly in the ground state. To our knowledge, these results provide the first experimental evidence for electrostatic stabilization of a charged enzymatic transition state and intermediate using the anionic backbone of DNA.     

The pH dependence of the thermodynamics of the interaction of 3'-cytidine monophosphate with ribonuclease A.

The apparent free energy (deltaGapp) and enthalpy changes (deltaHB) associated with the interaction of 3'-cytosine monophosphate (3'-CMP) and ribonuclease A (RNase) are reported for the pH range 4--9, T = 25 degrees, mu = 0.05. The pH dependence of deltaGapp and deltaHB has been interpreted in terms of coupled ionization of histidine residues 12, 48, and 119, assuming that only the dianionic form of the inhibitor is bound. The results of this analysis are consistent with the calorimetric and potentiometric titration results for the free enzyme and its 3'-CMP complex reported in the previous paper (M. Flogel and R. L. Biltonen ((1975), Biochemistry, preceding paper in this issue). This analysis allows the calculation of the thermodynamic quantities associated with hypothetical but clearly defined reactions (e.g., the reaction of the dianionic inhibitor with the completely protonated enzyme). It is concluded that the primary thermodynamic driving forces for the reaction are van der Waals interactions between the riboside moiety and the protein fabric and electrostatic interaction between the negatively charged phosphate group of the inhibitor and the positively charged histidine residues at the binding locus. It is also suggested that the binding reaction is weakly coupled (approximately to 0.5 kcal/mol) with a conformational change of the protein associated with protonation of residue 48. These results are consistent with the model originally proposed by G. G. Hammes ((1968), Adv. Protein Chem. 23, 1) and lend additional quantitative detail to the nature of the reaction.     

Structural determinants of trypsin affinity and specificity for cationic inhibitors

The binding free energies of four inhibitors to bovine beta-trypsin are calculated. The inhibitors use either ornithine, lysine, or arginine to bind to the S1 specificity site. The electrostatic contribution to binding free energy is calculated by solving the finite difference Poisson-Boltzmann equation, the contribution of nonpolar interactions is calculated using a free energy-surface area relationship and the loss of conformational entropy is estimated both for trypsin and ligand side chains. Binding free energy values are of a reasonable magnitude and the relative affinity of the four inhibitors for trypsin is correctly predicted. Electrostatic interactions are found to oppose binding in all cases. However, in the case of ornithine- and lysine-based inhibitors, the salt bridge formed between their charged group and the partially buried carboxylate of Asp189 is found to stabilize the complex. Our analysis reveals how the molecular architecture of the trypsin binding site results in highly specific recognition of substrates and inhibitors. Specifically, partially burying Asp189 in the inhibitor-free enzyme decreases the penalty for desolvation of this group upon complexation. Water molecules trapped in the binding interface further stabilize the buried ion pair, resulting in a favorable electrostatic contribution of the ion pair formed with ornithine and lysine side chains. Moreover, all side chains that form the trypsin specificity site are partially buried, and hence, relatively immobile in the inhibitor-free state, thus reducing the entropic cost of complexation. The implications of the results for the general problem of recognition and binding are considered. A novel finding in this regard is that like charged molecules can have electrostatic contributions to binding that are more favorable than oppositely charged molecules due to enhanced interactions with the solvent in the highly charged complex that is formed.     

Shift in nucleotide conformational equilibrium contributes to increased rate of catalysis of GpAp versus GpA in barnase

The microbial ribonuclease barnase exhibits low catalytic activity toward GpN dinucleotides, where G is guanosine, p is phosphate and N represents any nucleoside. When a phosphate is added to the 3'-end, as in GpNp, substrate affinity is enhanced by one order of magnitude, and the catalytic rate by two. In order to gain insight into this phenomenon, we analyzed the nucleotide conformations and protein-nucleotide interactions of 4 ns molecular dynamics (MD) trajectories of complexes of barnase with guanylyl(3'-5') adenosine (GpA) and guanylyl(3'-5') adenosine 3'-monophosphate (GpAp), respectively, in the presence of solvent and counter ions. We found that, in a majority of the bound GpA conformations, the guanine base was firmly bound to the recognition site. The phosphate and adenosine moieties pointed into the solvent, and interactions with key catalytic residues were absent. In contrast, the bound GpAp adopted conformations in which all of the nucleotide portions remained tightly bound to the enzyme and interactions with key catalytic residues were maintained. These observations indicate that, for GpA, a significant proportion of the bound nucleotide adopts non-productive conformations and that adding the terminal phosphate as in GpAp shifts the equilibrium of the bound conformations towards structures capable of undergoing catalysis. Incorporating this property into the kinetic equations yields an increase in both the apparent rate constant (kcat) and the apparent dissociation constant (K(M)) for GpAp versus GpA. The increase in K(M), caused by the presence of additional non-productive binding modes for GpA, should however be counterbalanced by the propensity of free GpA to adopt folded conformations in solution, which are unable to bind the enzyme and by the tighter binding of GpAp (Giraldo J, Wodak SJ, Van Belle D. Conformational analysis of GpA and GpAp in aqueous solution by molecular dynamics and statistical methods. J Mol Biol 1998; 283:863-882). Addition of the terminal phosphate is shown to significantly influence the collective motion of the enzyme in a manner that fosters interactions with key catalytic residues, representing a further likely contribution to the catalytic rate enhancement.     

Electrostatic interactions in catalytic centers of F1-ATPase

The first part of this paper is a brief review of works concerned with the mechanisms of functioning of F0F1-ATP synthases. F0F1-ATP syntheses operate as rotating molecular machines that provide the synthesis of ATP from ADP and inorganic phosphate (Pi) in mitochondria, chloroplasts, and bacteria at the expense of the energy of electrochemical gradient of hydrogen ions generated across energy-transducing mitochondrial, chloroplast or, bacterial membranes. A distinguishing feature of these enzymes is that they operate as rotary molecular motors. In the second part of the work, we calculated the contribution of electrostatic interactions between charged groups of a substrate (MgATP), reaction products (MgADP and Pi), and charged amino acid residues of the F1-ATPase molecule to energy changes associated with the binding of ATP and its chemical transformations in the catalytic centers located at the interface of the alpha- and beta-subunits of the enzyme (oligomer complex alpha 3 beta 3 gamma of bovine mitochondrial ATPase). The catalytic cycle of ATP hydrolysis considered in the work includes conformational changes of alpha- and beta-subunits caused by unidirectional rotations of the central gamma-subunit. The results of our calculations are consistent with the idea that the energetically favorable process of ATP binding to the "open" catalytic center of F1-ATPase initiates the rotation of the gamma-subunit followed by ATP hydrolysis in another ("closed") catalytic center of the enzyme.     

The roles of ATP4- and Mg2+ in control steps of phosphoglycerate kinase

1H-NMR measurements were made of solutions of yeast phosphoglycerate kinase containing the nucleotide substrate, ATP, and Mg2+ in varying concentrations in order to investigate the affect that the metal ion has on the mode of ATP binding to the enzyme. From the change in the chemical shifts of the 'basic-patch' histidine resonances (His62, His167 and His170) and the nucleotide C8H, C2H and C1'H resonances it is apparent that there are at least two ATP-binding sites on the enzyme. Downfield shifts observed for the above histidine resonances at low nucleotide/enzyme molar ratios indicates that the primary binding site involves electrostatic interactions between the nucleotide triphosphate chain and the basic-patch region of the N-terminal domain. The secondary binding site is shown to involve predominantly hydrophobic interactions between the adenosine moiety and the protein. Evidence from previous two-dimensional NMR experiments [Fairbrother et al. (1990) Eur. J. Biochem. 190, 161-169] suggests that the secondary site is equivalent to the crystallographically observed catalytic site. The affinity of the catalytic site is increased relative to the primary electrostatic site with increasing Mg2+ concentration. The possible importance of these observations in the regulation of this enzyme in vivo are discussed.     

Electrostatic attraction by surface charge does not contribute to the catalytic efficiency of acetylcholinesterase.

Acetylcholinesterases (AChEs) are characterized by a high net negative charge and by an uneven surface charge distribution, giving rise to a negative electrostatic potential extending over most of the molecular surface. To evaluate the contribution of these electrostatic properties to the catalytic efficiency, 20 single- and multiple-site mutants of human AChE were generated by replacing up to seven acidic residues, vicinal to the rim of the active-center gorge (Glu84, Glu285, Glu292, Asp349, Glu358, Glu389 and Asp390), by neutral amino acids. Progressive simulated replacement of these charged residues results in a gradual decrease of the negative electrostatic potential which is essentially eliminated by neutralizing six or seven charges. In marked contrast to the shrinking of the electrostatic potential, the corresponding mutations had no significant effect on the apparent bimolecular rate constants of hydrolysis for charged and non-charged substrates, or on the Ki value for a charged active center inhibitor. Moreover, the kcat values for all 20 mutants are essentially identical to that of the wild type enzyme, and the apparent bimolecular rate constants show a moderate dependence on the ionic strength, which is invariant for all the enzymes examined. These findings suggest that the surface electrostatic properties of AChE do not contribute to the catalytic rate, that this rate is probably not diffusion-controlled and that long-range electrostatic interactions play no role in stabilization of the transition states of the catalytic process.     

The flattened face of type II beta phosphatidylinositol phosphate kinase binds acidic phospholipid membranes

Type II beta phosphatidylinositol phosphate kinase is a representative phosphatidylinositol phosphate kinase that is active against membrane-bound substrates. The structure of the enzyme contains a flattened basic face that spans the crystallographic dimer interface and is adjacent to the active site. Analytical ultracentrifugation shows that phosphatidylinositol phosphate kinase is a dimer in solution. Modeling suggested that the flattened face binds to acidic phospholipids by electrostatic interactions. The enzyme binds to acidic vesicles containing phosphatidylserine, phosphatidic acid, or phosphoinositides mixed with phosphatidylcholine, but not to neutral phosphatidylcholine vesicles. Binding to acidic vesicles is abolished in the presence of 1.0 M NaCl, consistent with an essential electrostatic contribution to the free energy of binding. The +14 charge on the flattened face of the dimer was reduced to +2 in the triple mutant Lys72Glu/Lys76Glu/Lys78Glu. The mutation has no effect on dimerization, but reduces the apparent KA for 25% phosphatidylserine/75% phosphatidylcholine mixed vesicles by 16-fold. The reduction in the level of binding can be ascribed to a loss of electrostatic interactions based on the finite difference solution to the Poisson-Boltzmann equation. The mutant reduces catalytic activity toward phosphatidylinositol 5-phosphate by approximately 50-fold. The wild-type enzyme binds half-maximally to phosphatidylinositol 4,5-bisphosphate-containing vesicles at a mole fraction of 0.3% in a phosphatidylcholine background, as compared to a 22% mole fraction in phosphatidylserine. The binding to phosphatidylinositol 4,5-bisphosphate-containing membranes is less sensitive to salt and to the triple mutation than binding to phosphatidylserine-containing membranes, suggesting that at least part of phosphatidylinositol 4,5-bisphosphate's interaction with the enzyme is independent of the flattened face. It is concluded that the flattened face of type II beta phosphatidylinositol phosphate kinase binds to membranes through nonspecific interactions, and that this interaction is essential for efficient catalysis.     

The roles of ADP2- and Mg2+ in control steps of phosphoglycerate kinase

1H-NMR measurements were made of solutions of yeast phosphoglycerate kinase containing the nucleotide, ADP, and Mg2+ in varying concentrations in order to investigate the affect that the metal ion has on the mode of ADP binding to the enzyme. A preliminary study of adenosine binding to phosphoglycerate kinase was made in order to be sure of the nature of the adenine site. From the change in chemical shifts of the 'basic patch' histidine resonances (His62, 167 and 170), the nucleotide C8-H, C2-H and C1'-H resonances and resonances 40 and 41 (assigned to Thr373 and Thr375 in the hydrophobic, i.e. catalytic, site), it is apparent that there are at least two ADP binding sites on the enzyme: one at the hydrophobic (catalytic) site and one at the electrostatic site. A comparison of the results for ADP and ATP reveals differences due to the differential binding of the phosphate groups. The presence of Mg2+ results in further differences being observed. The data suggest that the primary binding site of ADP, in the absence of Mg2+, involves electrostatic interactions between the diphosphate chain of the substrate and the 'basic patch' region of the N-terminal domain. In the presence of greater than or equal to 1:1 ratio of Mg2+/ADP, however, the primary binding site involves predominantly hydrophobic interactions between the adenosine moiety and the catalytic site, with secondary binding occurring at the electrostatic site. Addition of Mg2+, therefore, tends to reduce the affinity of the electrostatic site (presumably by competing for ADP). It is suggested that alpha-helix XII, including residues 372, 373 and 375, moves differentially on binding ADP, Mg ADP, ATP or Mg . ATP, consistent with Mg2+ assisting the transfer of the gamma-phosphate of ATP to 3-phosphoglycerate during catalysis.     

Uracil DNA glycosylase: revisiting substrate-assisted catalysis by DNA phosphate anions

Uracil DNA glycosylase (UNG) is a powerful DNA repair enzyme that has been shown to stabilize a glycosyl cation reaction intermediate and a related tight binding inhibitor using electrostatic interactions with the +1 and -1, but not the +2, phosphodiester group of the single-stranded DNA substrate Ap (2+)Ap (1+)Up (1-)ApA. These experimental results differed considerably from computational findings using duplex DNA, where the +2 phosphate was found to stabilize the transition state by approximately 5 kcal/mol, suggesting that UNG uses different catalytic strategies with single-stranded and double-stranded DNA substrates. In addition, the computational studies indicated that the conserved and positively charged His148 (which hydrogen bonds to the +2 phosphate) destabilized the glycosyl cation intermediate by 6-8 kcal/mol through anticatalytic electrostatic interactions. To evaluate these interesting proposals, we measured the kinetic effects of neutral methylphosphonate (MeP) stereoisomers at the +1 and +2 positions of a 12-mer dsDNA substrate and also the catalytic contribution and ionization state of His148. For MeP substitutions at the +1 position, single-turnover kinetic studies showed that the activation barrier was increased by 9.8 and 3.1 kcal/mol, corresponding to a stereoselectivity of nearly 40000-fold for the respective MeP isomers. Identical to the findings with ssDNA, MeP substitutions at the +2 position resulted in only small changes in the activation barrier (+/-0.3 kcal/mol), with little stereoselectivity ( approximately 4-fold). However, the H148A mutation destabilizes both the ground state and transition states by 2.4 and 4.3 kcal/mol, respectively. Thus, His148 is catalytic because it stabilizes the transition state to a greater extent (1.9 kcal/mol) than the ground state. Heteronuclear NMR studies established that His148 was neutral in the free enzyme at neutral pH, and in conformational exchange in a specific DNA complex containing uracil. We conclude that the +1 and +2 phosphate esters play identical catalytic roles in the reactions of single-stranded and double-stranded DNA substrates, and that His148 serves a catalytic role by positioning the substrate and catalytic water, or by an environmental effect.     

Structurally and catalytically important residues in the phosphate binding loop of adenylate kinase of Escherichia coli

Amino acids in the phosphate binding loop of adenylate kinase of Escherichia coli were mutated by site-directed mutagenesis. The mutant proteins with a Pro-9----Gly (P9G) and with a Lys-13----Gln (K13Q) exchange were overexpressed and purified. They were characterized by steady-state kinetics, fluorescence binding, and structural studies, together with the phosphate binding loop mutants P9L and G10V prepared earlier [Reinstein, J., Brune, M., & Wittinghofer, A. (1988) Biochemistry 27, 4712-4720]. The results obtained show that all these mutations change the structure of the protein as evidenced by NMR spectroscopy and temperature-stability studies. All the mutant proteins have increased dissociation constants for substrates and inhibitors, but their catalytic activity, except for K13Q, is not reduced. The results obtained with K13Q suggest that this lysine residue, which is conserved in all guanine and many adenine nucleotide proteins, might have an important role in catalysis.     

Binding interactions between the core central domain of 16S rRNA and the ribosomal protein S15 determined by molecular dynamics simulations

The goal of the current study is to utilize molecular dynamic (MD) simulations to investigate the dynamic behavior of 16S rRNA in the presence and absence of S15 and to identify the binding interactions between these two molecules. The simulations show that: (i) 16S rRNA remains in a highly folded structure when it is bound to S15; (ii) in the absence of S15, 16S rRNA significantly alters its conformation and transiently forms conformations that are similar to the bound structure that make it available for binding with S15; (iii) the unbound rRNA spends the majority of its time in extended conformations. The formation of the extended conformations is a result of the molecule reaching a lower electrostatic energy and the formation of the highly folded, crystal-like conformation is a result of achieving a lower solvation energy. In addition, our MD simulations show that 16S rRNA and S15 bind across the major groove of helix 22 (H22) via electrostatic interactions. The negatively charged phosphate groups of G658, U740, G741 and G742 bind to the positively charged S15 residues Lys7, Arg34 and Arg37. The current study provides a dynamic view of the binding of 16S rRNA with S15.     

Automated computational framework for the analysis of electrostatic similarities of proteins

Charge plays an important role in protein-protein interactions. In the case of excessively charged proteins, their electrostatic potentials contribute to the processes of recognition and binding with other proteins or ligands. We present an automated computational framework for determining the contribution of each charged amino acid to the electrostatic properties of proteins, at atomic resolution level. This framework involves computational alanine scans, calculation of Poisson-Boltzmann electrostatic potentials, calculation of electrostatic similarity distances (ESDs), hierarchical clustering analysis of ESDs, calculation of solvation free energies of association, and visualization of the spatial distributions of electrostatic potentials. The framework is useful to classify families of mutants with similar electrostatic properties and to compare them with the parent proteins in the complex. The alanine scan mutants introduce perturbations in the local electrostatic properties of the proteins and aim in delineating the contribution of each mutated amino acid in the spatial distribution of electrostatic potential, and in biological function when electrostatics is a dominant contributing factor in protein-protein interactions. The framework can be used to design new proteins with tailored electrostatic properties, such as immune system regulators, inhibitors, and vaccines, and in guiding experimental studies. We present an example for the interaction of the immune system protein C3d (the d-fragment of complement protein C3) with its receptor CR2, and we discuss our data in view of a binding site controversy.     

Nucleotide binding to the isolated beta subunit of the chloroplast ATP synthase.

The beta subunit isolated from the chloroplast ATP synthase F1 (CF1) has a single dissociable nucleotide binding site, consistent with the proposed function of this subunit in nucleotide binding and catalysis. The beta subunit bound the nucleotide analogs trinitrophenyl-ATP (TNP-ATP) or trinitrophenyl-ADP (TNP-ADP) with nearly equal affinities (Kd = 1-2 microM) but did not bind trinitrophenyl-AMP. Both ATP and ADP effectively competed with TNP-ATP for binding. Other nucleoside triphosphates were also able to compete with TNP-ATP for binding to beta; their order of effectiveness (ATP greater than GTP, ITP greater than CTP) mimicked the normal substrate specificity of CF1. The single nucleotide binding site on the isolated beta subunit very closely resembles the low affinity catalytic site (site 3) of CF1 (Bruist, M.F., and Hammes, G. G. (1981) Biochemistry 20, 6298-6305), suggesting that tight nucleotide binding by other sites on the enzyme involves other CF1 subunits in addition to the beta subunit. The results are inconsistent with an earlier report (Frasch, W.D., Green, J., Caguial, J., and Mejia, A. (1989) J. Biol. Chem. 264, 5064-5069), which suggested more than one nucleotide binding site per beta subunit. Binding of nucleotides to the isolated beta subunit was eliminated by a brief heat treatment (40 degrees C for 10 min) of the protein. A small change in the circular dichroism spectrum of beta accompanied the heat treatment indicating that a localized (rather than global) change in the folding of beta, involving at least part of the nucleotide binding domain, had occurred. Also accompanying the loss of nucleotide binding was a loss of the reconstitutive capacity of the beta subunit. ATP protected against the effects of the heat treatment.     

Highly organized but pliant active site of DNA polymerase beta: compensatory mechanisms in mutant enzymes revealed by dynamics simulations and energy analyses

To link conformational transitions noted for DNA polymerases with kinetic results describing catalytic efficiency and fidelity, we investigate the role of key DNA polymerase beta residues on subdomain motion through simulations of five single-residue mutants: Arg-283-Ala, Tyr-271-Ala, Asp-276-Val, Arg-258-Lys, and Arg-258-Ala. Since a movement toward a closed state was only observed for R258A, we suggest that Arg(258) is crucial in modulating motion preceding chemistry. Analyses of protein/DNA interactions in the mutant active site indicate distinctive hydrogen bonding and van der Waals patterns arising from compensatory structural adjustments. By comparing closed mutant complexes with the wild-type enzyme, we interpret experimentally derived nucleotide binding affinities in molecular terms: R283A (decreased), Y271A (increased), D276V (increased), and R258A (decreased). Thus, compensatory interactions (e.g., in Y271A with adjacent residues Phe(272), Asn(279), and Arg(283)) increase the overall binding affinity for the incoming nucleotide although direct interactions may decrease. Together with energetic analyses, we predict that R258G might increase the rate of nucleotide insertion and maintain enzyme fidelity as R258A; D276L might increase the nucleotide binding affinity more than D276V; and R283A/K280A might decrease the nucleotide binding affinity and increase misinsertion more than R283A. The combined observations regarding key roles of specific residues (e.g., Arg(258)) and compensatory interactions echo the dual nature of polymerase active site, namely versatility (to accommodate various basepairs) and specificity (for preserving fidelity) and underscore an organized but pliant active site essential to enzyme function.     

Pseudomonas aeruginosa aspartate transcarbamoylase. Characterization of its catalytic and regulatory properties

Aspartate transcarbamoylase from Pseudomonadaceae is a class A enzyme consisting of six copies of a 36-kDa catalytic chain and six copies of a 45-kDa polypeptide of unknown function. The 45-kDa polypeptide is homologous to dihydroorotase but lacks catalytic activity. Pseudomonas aeruginosa aspartate transcarbamoylase was overexpressed in Escherichia coli. The homogeneous His-tagged protein isolated in high yield, 30 mg/liter of culture, by affinity chromatography and crystallized. Attempts to dissociate the catalytic and pseudo-dihydroorotase (pDHO) subunits or to express catalytic subunits only were unsuccessful suggesting that the pDHO subunits are required for the proper folding and assembly of the complex. As reported previously, the enzyme was inhibited by micromolar concentrations of all nucleotide triphosphates. In the absence of effectors, the aspartate saturation curves were hyperbolic but became strongly sigmoidal in the presence of low concentrations of nucleotide triphosphates. The inhibition was unusual in that only free ATP, not MgATP, inhibits the enzyme. Moreover, kinetic and binding studies with a fluorescent ATP analog suggested that ATP induces a conformational change that interferes with the binding of carbamoyl phosphate but has little effect once carbamoyl phosphate is bound. The peculiar allosteric properties suggest that the enzyme may be a potential target for novel chemotherapeutic agents designed to combat Pseudomonas infection.     

Inhibition of yeast mitochondrial F1-ATPase, F0F1-ATPase and submitochondrial particles by rhodamines and ethidium bromide

ATP hydrolysis by F1-ATPase is strongly inhibited by cationic rhodamines; neutral rhodamines are very poor inhibitors. Rhodamine 6G is a noncompetitive inhibitor of purified F0F1-ATPase and submitochondrial particles, however, an uncompetitive inhibitor of F1-ATPase (KI approximately equal to 2.4 microM for all three enzyme forms). Ethidium bromide is a noncompetitive inhibitor of F0F1-ATPase, submitochondrial particles and also F1-ATPase (KI approximately equal to 270 microM). Neither of the inhibitors affects the negative cooperativity (nH approximately equal to 0.7). The non-identical binding sites for rhodamine 6G and ethidium bromide are located on the F1-moiety and are topologically distinct from the catalytic site. Binding of the inhibitors prevents the conformational changes essential for energy transduction. It is concluded that the inhibitor binding sites are involved in proton translocation. In F1-ATPase, binding of MgATP at a catalytic site causes conformational changes, which allosterically induce the correct structure of the rhodamine 6G binding site. In F0F1-ATPase, this conformation of the F1-moiety exists a priori, due to allosteric interactions with F0-subunits. The binding site for ethidium bromide on F1-ATPase does not require substrate binding at the catalytic site and is not affected by F0F1-subunit interactions.     

Calculation of HyHel10-lysozyme binding free energy changes: Effect of ten point mutations

The change in free energy of binding of hen egg white lysozyme (HEL) to the antibody HyHel-10 arising from ten point mutations in HEL (D101K, D101G, K96M, K97D, K97G, K97G, R21E, R21K, W62Y, and W63Y) was calculated using a combination of the finite difference Poisson-Boltzmann method for the electrostatic contribution, a solvent accessible surface area term for the non-polar contribution, and rotamer counting for the sidechain entropy contribution. Comparison of experimental and calculated results indicate that because of pKa shifts in some of the mutated residues, primarily those involving Aspartate or Glutamate, proton uptake or release occurs in binding. When this effect was incorporated into the binding free energy calculations, the agreement with experiment improved significantly, and resulted in a mean error of about 1.9 kcal/mole. Thus these calculations predict that there should be a significant pH dependence to the change in binding caused by these mutations. The other major contributions to binding energy changes comes from solvation and charge charge interactions, which tend to oppose each other. Smaller contributions come from nonpolar interactions and sidechain entropy changes. The structures of the HyHel-10-HEL complexes with mutant HEL were obtained by modeling, and the effect of the modeled structure on the calculations was also examined. “Knowledge based” modeling and automatic generation of models using molecular mechanics produced comparable results. Proteins 33:39–48, 1998. © 1998 Wiley-Liss, Inc.     

Interaction of tetrasubstituted cationic aluminum phthalocyanine with artificial and natural membranes

A study of the properties of water-soluble tetrasubstituted cationic aluminum phthalocyanine (AlPcN₄) revealed efficient binding of this photosensitizer to phospholipid membranes as compared with tetrasulfonated aluminum and zinc phthalocyanine complexes. This also manifested itself in enhanced photodynamic activity of AlPcN₄ as measured by the photosensitized damage of gramicidin channels in a planar bilayer lipid membrane. The largest difference in the photodynamic activity of cationic and anionic phthalocyanines was observed in a membrane containing negatively charged lipids, thereby pointing to significant contribution of electrostatic interactions to the binding of photosensitizers to a membrane. Fluoride anions suppressed the photodynamic activity and binding to membrane of both tetraanionic and tetracationic aluminum phthalocyanines, which supports our hypothesis that interaction of charged metallophthalocyanines with phospholipid membranes is mostly determined by coordination of the central metal atom with the phosphate group of lipid.     

Electrostatic interactions in sperm whale myoglobin. Site specificity, roles in structural elements, and external electrostatic potential distributions

The electrostatic free energy contribution to the stability of sperm whale ferrimyoglobin was evaluated according to the static accessibility modified Tanford-Kirkwood model. The electrostatic free energy contribution of each distinct structural element was divided into one term arising from interactions between it and other elements (interelemental) and another from interactions within the particular element itself (intraelemental). At pH 7 the majority of the terms were found to be stabilizing. The interelemental terms are the dominant ones for most structural elements. The small interelemental terms of the C and D helices are compensated by large intraelemental interactions which stabilize these short helices. Perturbations in pH can be accommodated by the structural elements through a redistribution of stabilizing and destabilizing interactions. The electrostatic potentials calculated at the surface of the protein indicate that the internal compensation of local potentials achieved during folding results in a generally neutral protein-solvent interface save for two distinct areas of nonzero potential. The accessibility of each charged atom to solvent was analyzed in terms of the surface area lost to charged, polar and nonpolar atoms separately. The net solvent accessibility lost parallels closely that lost to nonpolar atoms alone, indicating a specific role for nonpolar atoms in defining dielectric shielding of charged atoms, aside from their participation in the well-known hydrophobic interactions.