Cation-pi interactions in aromatics of biological and medicinal interest: electrostatic potential surfaces as a useful qualitative guide.

The cation-pi interaction is an important, general force for molecular recognition in biological receptors. Through the sidechains of aromatic amino acids, novel binding sites for cationic ligands such as acetylcholine can be constructed. We report here a number of calculations on prototypical cation-pi systems, emphasizing structures of relevance to biological receptors and prototypical heterocycles of the type often of importance in medicinal chemistry. Trends in the data can be rationalized using a relatively simple model that emphasizes the electrostatic component of the cation-pi interaction. In particular, plots of the electrostatic potential surfaces of the relevant aromatics provide useful guidelines for predicting cation-pi interactions in new systems.

Direct molecular level measurements of the electrostatic properties of a protein surface

In this work, we used direct measurements with the surface force apparatus to determine the pH-dependent electrostatic charge density of a single binding face of streptavidin. Mean field calculations have been used with considerable success to model electrostatic potential fields near protein surfaces, but these models and their inherent assumptions have not been tested directly at the molecular level. Using the force apparatus and immobilized, oriented monolayers of streptavidin, we measured a pI of 5–5.5 for the biotin-binding face of the protein. This differs from the pI of 6.3 for the soluble protein and confirms that we probed the local electrostatic features of the macromolecule. With finite difference solutions of the linearized Poisson–Boltzmann equation, we then calculated the pH-dependent charge densities adjacent to the same face of the protein. These calculated values agreed quantitatively with those obtained by direct force measurements. Although our study focuses on the pH-dependence of surface electrostatics, this direct approach to probing the electrostatic features of proteins is applicable to investigations of any perturbations that alter the charge distribution of the surfaces of immobilized molecules.

Electrostatic steering and ionic tethering in enzyme–ligand binding: Insights from simulations

To bind at an enzyme’s active site, a ligand must diffuse or be transported to the enzyme’s surface, and, if the binding site is buried, the ligand must diffuse through the protein to reach it. Although the driving force for ligand binding is often ascribed to the hydrophobic effect, electrostatic interactions also influence the binding process of both charged and nonpolar ligands. First, electrostatic steering of charged substrates into enzyme active sites is discussed. This is of particular relevance for diffusion-influenced enzymes. By comparing the results of Brownian dynamics simulations and electrostatic potential similarity analysis for triose-phosphate isomerases, superoxide dismutases, and β-lactamases from different species, we identify the conserved features responsible for the electrostatic substrate-steering fields. The conserved potentials are localized at the active sites and are the primary determinants of the bimolecular association rates. Then we focus on a more subtle effect, which we will refer to as “ionic tethering.” We explore, by means of molecular and Brownian dynamics simulations and electrostatic continuum calculations, how salt links can act as tethers between structural elements of an enzyme that undergo conformational change upon substrate binding, and thereby regulate or modulate substrate binding. This is illustrated for the lipase and cytochrome P450 enzymes. Ionic tethering can provide a control mechanism for substrate binding that is sensitive to the electrostatic properties of the enzyme’s surroundings even when the substrate is nonpolar.

Theoretical investigation of the [1,2]-sigmatropic hydrogen migration in the mechanism of oxidation of 2-aminobenzoyl-CoA by 2-aminobenzoyl-CoA monooxygenase/reductase

The flavin hydroperoxide at the active site of the mixed-function oxidase 2-aminobenzoyl-CoA monooxygenase/reductase (Azoarcus evansii) transfers an oxygen to the 5-position of the 2-aminobenzoyl-CoA substrate to provide the alkoxide intermediate II−. Hydrogen migration from C5 to C6 follows this monooxygenation. The nature of the monooxygenation intermediate and plausible competing reactions leading to hydrogen migration have been considered. Ab initio molecular orbital theory has been used to calculate structures and electron distributions in intermediate and transition state structures. Electrostatic potential surface calculations establish that the transition state and product, associated with the C5 to C6 hydrogen transfer, are stabilized by electron distribution to the benzoyl-CoA thioester carbonyl oxygen. This is not so for the transition state and product associated with hydrogen transfer from C5 to C4. The activation energy for the 5,6-shift is 2.5 kcal/mol lower than that for the 5,4-shift. In addition, the product of the hydrogen 5,6-shift is more stable than is the product of the hydrogen 5,4-shift, by ≈6 kcal/mol. These results explain why only the shift of hydrogen from C5 to C6 is observed experimentally. Oxygen transfer and hydrogen migration almost coincide in the gas phase (activation energy of ≈0.6 kcal/mol, equivalent to a single bond vibration). Enzymatic formation of alkoxide II− requires its stabilization; thus, the rate constant for its breakdown would be slower than in the gas phase.

Solution structure of the antiapoptotic protein bcl-2

The structures of two isoforms of Bcl-2 that differ by two amino acids have been determined by NMR spectroscopy. Because wild-type Bcl-2 behaved poorly in solution, the structures were determined by using Bcl-2/Bcl-xL chimeras in which part of the putative unstructured loop of Bcl-2 was replaced with a shortened loop from Bcl-xL. These chimeric proteins have a low pI compared with the wild-type protein and are soluble. The structures of the two Bcl-2 isoforms consist of 6 α-helices with a hydrophobic groove on the surface similar to that observed for the homologous protein, Bcl-xL. Comparison of the Bcl-2 structures to that of Bcl-xL shows that although the overall fold is the same, there are differences in the structural topology and electrostatic potential of the binding groove. Although the structures of the two isoforms of Bcl-2 are virtually identical, differences were observed in the ability of the proteins to bind to a 25-residue peptide from the proapoptotic Bad protein and a 16-residue peptide from the proapoptotic Bak protein. These results suggest that there are subtle differences in the hydrophobic binding groove in Bcl-2 that may translate into differences in antiapoptotic activity for the two isoforms.

Crystal structure at 2.6-A resolution of human macrophage migration inhibitory factor.

Macrophage migration inhibitory factor (MIF) was the first cytokine to be described, but for 30 years its role in the immune response remained enigmatic. In recent studies, MIF has been found to be a novel pituitary hormone and the first protein identified to be released from immune cells on glucocorticoid stimulation. Once secreted, MIF counterregulates the immunosuppressive effects of steroids and thus acts as a critical component of the immune system to control both local and systemic immune responses. We report herein the x-ray crystal structure of human MIF to 2.6 angstrom resolution. The protein is a trimer of identical subunits. Each monomer contains two antiparallel alpha-helices that pack against a four-stranded beta-sheet. The monomer has an additional two beta-strands that interact with the beta-sheets of adjacent subunits to form the interface between monomers. The three beta-sheets are arranged to form a barrel containing a solvent-accessible channel that runs through the center of the protein along a molecular 3-fold axis. Electrostatic potential maps reveal that the channel has a positive potential, suggesting that it binds negatively charged molecules. The elucidated structure for MIF is unique among cytokines or hormonal mediators, and suggests that this counterregulator of glucocorticoid action participates in novel ligand-receptor interactions.

Negative electrostatic surface potential of protein sites specific for anionic ligands.

Determination of the crystal structure of an "open" unliganded active mutant (T141D) form of the Escherichia coli phosphate receptor for active transport has allowed calculation of the electrostatic surface potential for it and two other comparably modeled receptor structures (wild type and D137N). A discovery of considerable implication is the intensely negative potential of the phosphate-binding cleft. We report similar findings for a sulfate transport receptor, a DNA-binding protein, and, even more dramatically, redox proteins. Evidently, for proteins such as these, which rely almost exclusively on hydrogen bonding for anion interactions and electrostatic balance, a noncomplementary surface potential is not a barrier to binding. Moreover, experimental results show that the exquisite specificity and high affinity of the phosphate and sulfate receptors for unions are insensitive to modulations of charge potential, but extremely sensitive to conditions that leave a hydrogen bond donor or acceptor unpaired.

Design of fast enzymes by optimizing interaction potential in active site

The diffusional encounter between substrate and enzyme, and hence catalytic efficiency, can be enhanced by mutating charged residues on the surface of the enzyme. In this paper we present a simple method for screening such mutations. This is based on our earlier result that electrostatic enhancement of the enzyme-substrate binding rate constant can be accounted for just by the interaction potential within the active site. Assuming that catalytic and structural integrity is maintained, the catalytic efficiency can be optimized by surface charge mutations which lead to stronger interaction potential within the active site. Application of the screening method on superoxide dismutase shows that only charge mutations close to the active site will have practical effect on the catalytic efficiency. This rationalizes a large number of findings obtained in previous simulation and experimental studies.

Relationship between the oxidation potential and electron spin density of the primary electron donor in reaction centers from Rhodobacter sphaeroides

The primary electron donor in bacterial reaction centers is a dimer of bacteriochlorophyll a molecules, labeled L or M based on their proximity to the symmetry-related protein subunits. The electronic structure of the bacteriochlorophyll dimer was probed by introducing small systematic variations in the bacteriochlorophyll–protein interactions by a series of site-directed mutations that replaced residue Leu M160 with histidine, tyrosine, glutamic acid, glutamine, aspartic acid, asparagine, lysine, and serine. The midpoint potentials for oxidation of the dimer in the mutants showed an almost continuous increase up to ≈60 mV compared with wild type. The spin density distribution of the unpaired electron in the cation radical state of the dimer was determined by electron–nuclear–nuclear triple resonance spectroscopy in solution. The ratio of the spin density on the L side of the dimer to the M side varied from ≈2:1 to ≈5:1 in the mutants compared with ≈2:1 for wild type. The correlation between the midpoint potential and spin density distribution was described using a simple molecular orbital model, in which the major effect of the mutations is assumed to be a change in the energy of the M half of the dimer, providing estimates for the coupling and energy levels of the orbitals in the dimer. These results demonstrate that the midpoint potential can be fine-tuned by electrostatic interactions with amino acids near the dimer and show that the properties of the electronic structure of a donor or acceptor in a protein complex can be directly related to functional properties such as the oxidation–reduction midpoint potential.

Electrostatic surface plasmon resonance: Direct electric field-induced hybridization and denaturation in monolayer nucleic acid films and label-free discrimination of base mismatches

We demonstrate that in situ optical surface plasmon resonance spectroscopy can be used to monitor hybridization kinetics for unlabeled DNA in tethered monolayer nucleic acid films on gold in the presence of an applied electrostatic field. The dc field can enhance or retard hybridization and can also denature surface-immobilized DNA duplexes. Discrimination between matched and mismatched hybrids is achieved by simple adjustment of the electrode potential. Although the electric field at the interface is extremely large, the tethered single-stranded DNA thiol probes remain bound and can be reused for subsequent hybridization reactions without loss of efficiency. Only capacitive charging currents are drawn; redox reactions are avoided by maintaining the gold electrode potential within the ideally polarizable region. Because of potential-induced changes in the shape of the surface plasmon resonance curve, we account for the full curve rather than simply the shift in the resonance minimum.

Location of a potential transport binding site in a sigma class glutathione transferase by x-ray crystallography.

The crystal structure of the sigma class glutathione transferase from squid digestive gland in complex with S-(3-iodobenzyl)glutathione reveals a third binding site for the glutathione conjugate besides the two in the active sites of the dimer. The additional binding site is near the crystallographic two-fold axis between the two alpha 4-turn-alpha 5 motifs. The principal binding interactions with the conjugate include specific electrostatic interactions between the peptide and the two subunits and a hydrophobic cavity found across the two-fold axis that accommodates the 3-iodobenzyl group. Thus, two identical, symmetry-related but mutually exclusive binding modes for the third conjugate are observed. The hydrophobic pocket is about 14 A from the hydroxyl group of Tyr-7 in the active site. This site is a potential transport binding site for hydrophobic molecules or their glutathione conjugates.

Kinetics of photo-induced electron transfer from high-potential iron-sulfur protein to the photosynthetic reaction center of the purple phototroph Rhodoferax fermentans.

The kinetics of photo-induced electrontransfer from high-potential iron-sulfur protein (HiPIP) to the photosynthetic reaction center (RC) of the purple phototroph Rhodoferarfermentans were studied. The rapid photooxidation of heme c-556 belonging to RC is followed, in the presence of HiPIP, by a slower reduction having a second-order rate constant of 4.8 x 10(7) M(-1) x s(-1). The limiting value of kobs at high HiPIP concentration is 95 s(-1). The amplitude of this slow process decreases with increasing HiPIP concentration. The amplitude of a faster phase, observed at 556 and 425 nm and involving heme c-556 reduction, increases proportionately. The rate constant of this fast phase, determined at 425 and 556 nm, is approximately 3 x 10(5) s(-1). This value is not dependent on HiPIP concentration, indicating that it is related to a first-order process. These observations are interpreted as evidence for the formation of a HiPIP-RC complex prior to the excitation flash, having a dissociation constant of -2.5 microM. The fast phase is absent at high ionic strength, indicating that the complex involves mainly electrostatic interactions. The ionic strength dependence of kobs for the slow phase yields a second-order rate constant at infinite ionic strength of 5.4 x 10(6) M(-1) x s(-1) and an electrostatic interaction energy of -2.1 kcal/mol (1 cal = 4.184 J). We conclude that Rhodoferar fermentans HiPIP is a very effective electron donor to the photosynthetic RC.