Electrostatic interactions between rotor and stator in the bacterial flagellar motor

Bacterial flagellar motors rotate, obtaining power from the membrane gradient of protons or, in some species, sodium ions. Torque generation in the flagellar motor must involve interactions between components of the rotor and components of the stator. Sites of interaction between the rotor and stator have not been identified. Mutational studies of the rotor protein FliG and the stator protein MotA showed that both proteins contain charged residues essential for motor rotation. This suggests that functionally important electrostatic interactions might occur between the rotor and stator. To test this proposal, we examined double mutants with charged-residue substitutions in both the rotor protein FliG and the stator protein MotA. Several combinations of FliG mutations with MotA mutations exhibited strong synergism, whereas others showed strong suppression, in a pattern that indicates that the functionally important charged residues of FliG interact with those of MotA. These results identify a functionally important site of interaction between the rotor and stator and suggest a hypothesis for electrostatic interactions at the rotor–stator interface.

Potentiation of proton transfer function by electrostatic interactions in photosynthetic reaction centers from Rhodobacter sphaeroides: First results from site-directed mutation of the H subunit.

The x-ray crystallographic structure of the photosynthetic reaction center (RC) has proven critical in understanding biological electron transfer processes. By contrast, understanding of intraprotein proton transfer is easily lost in the immense richness of the details. In the RC of Rhodobacter (Rb.) sphaeroides, the secondary quinone (QB) is surrounded by amino acid residues of the L subunit and some buried water molecules, with M- and H-subunit residues also close by. The effects of site-directed mutagenesis upon RC turnover and quinone function have implicated several L-subunit residues in proton delivery to QB, although some species differences exist. In wild-type Rb. sphaeroides, Glu L212 and Asp L213 represent an inner shell of residues of particular importance in proton transfer to QB. Asp L213 is crucial for delivery of the first proton, coupled to transfer of the second electron, while Glu L212, possibly together with Asp L213, is necessary for delivery of the second proton, after the second electron transfer. We report here the first study, by site-directed mutagenesis, of the role of the H subunit in QB function. Glu H173, one of a cluster of strongly interacting residues near QB, including Asp L213, was altered to Gln. In isolated mutant RCs, the kinetics of the first electron transfer, leading to formation of the semiquinone, QB-, and the proton-linked second electron transfer, leading to the formation of fully reduced quinol, were both greatly retarded, as observed previously in the Asp L213 --> Asn mutant. However, the first electron transfer equilibrium, QA-QB <==> QAQB-, was decreased, which is opposite to the effect of the Asp L213 --> Asn mutation. These major disruptions of events coupled to proton delivery to QB were largely reversed by the addition of azide (N3-). The results support a major role for electrostatic interactions between charged groups in determining the protonation state of certain entities, thereby controlling the rate of the second electron transfer. It is suggested that the essential electrostatic effect may be to "potentiate" proton transfer activity by raising the pK of functional entities that actually transfer protons in a coupled fashion with the second electron transfer. Candidates include buried water (H3O+) and Ser L223 (serine-OH2+), which is very close to the O5 carbonyl of the quinone.

On the magnitude of the electrostatic contribution to ligand-DNA interactions.

A model based on the nonlinear Poisson-Boltzmann equation is used to study the electrostatic contribution to the binding free energy of a simple intercalating ligand, 3,8-diamino-6-phenylphenanthridine, to DNA. We find that the nonlinear Poisson-Boltzmann model accurately describes both the absolute magnitude of the pKa shift of 3,8-diamino-6-phenylphenanthridine observed upon intercalation and its variation with bulk salt concentration. Since the pKa shift is directly related to the total electrostatic binding free energy of the charged and neutral forms of the ligand, the accuracy of the calculations implies that the electrostatic contributions to binding are accurately predicted as well. Based on our results, we have developed a general physical description of the electrostatic contribution to ligand-DNA binding in which the electrostatic binding free energy is described as a balance between the coulombic attraction of a ligand to DNA and the disruption of solvent upon binding. Long-range coulombic forces associated with highly charged nucleic acids provide a strong driving force for the interaction of cationic ligands with DNA. These favorable electrostatic interactions are, however, largely compensated for by unfavorable changes in the solvation of both the ligand and the DNA upon binding. The formation of a ligand-DNA complex removes both charged and polar groups at the binding interface from pure solvent while it displaces salt from around the nucleic acid. As a result, the total electrostatic binding free energy is quite small. Consequently, nonpolar interactions, such as tight packing and hydrophobic forces, must play a significant role in ligand-DNA stability.

The role of London forces in defining noncentrosymmetric order of high dipole moment–high hyperpolarizability chromophores in electrically poled polymeric thin films

Graphs of second harmonic generation coefficients and electro-optic coefficients (measured by ellipsometry, attenuated total reflection, and two-slit interference modulation) as a function of chromophore number density (chromophore loading) are experimentally observed to exhibit maxima for polymers containing chromophores characterized by large dipole moments and polarizabilities. Modified London theory is used to demonstrated that this behavior can be attributed to the competition of chromophore-applied electric field and chromophore–chromophore electrostatic interactions. The comparison of theoretical and experimental data explains why the promise of exceptional macroscopic second-order optical nonlinearity predicted for organic materials has not been realized and suggests routes for circumventing current limitations to large optical nonlinearity. The results also suggest extensions of measurement and theoretical methods to achieve an improved understanding of intermolecular interactions in condensed phase materials including materials prepared by sequential synthesis and block copolymer methods.

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.

The N-terminal domain of PsaF: Precise recognition site for binding and fast electron transfer from cytochrome c6 and plastocyanin to photosystem I of Chlamydomonas reinhardtii

The PsaF-deficient mutant 3bF of Chlamydomonas reinhardtii was used to modify PsaF by nuclear transformation and site-directed mutagenesis. Four lysine residues in the N-terminal domain of PsaF, which have been postulated to form the positively charged face of a putative amphipathic α-helical structure were altered to K12P, K16Q, K23Q, and K30Q. The interactions between plastocyanin (pc) or cytochrome c6 (cyt c6) and photosystem I (PSI) isolated from wild type and the different mutants were analyzed using crosslinking techniques and flash absorption spectroscopy. The K23Q change drastically affected crosslinking of pc to PSI and electron transfer from pc and cyt c6 to PSI. The corresponding second order rate constants for binding of pc and cyt c6 were reduced by a factor of 13 and 7, respectively. Smaller effects were observed for mutations K16Q and K30Q, whereas in K12P the binding was not changed relative to wild type. None of the mutations affected the half-life of the microsecond electron transfer performed within the intermolecular complex between the donors and PSI. The fact that these single amino acid changes within the N-terminal domain of PsaF have different effects on the electron transfer rate constants and dissociation constants for both electron donors suggests the existence of a rather precise recognition site for pc and cyt c6 that leads to the stabilization of the final electron transfer complex through electrostatic interactions.

On the electronic nature of low-barrier hydrogen bonds in enzymatic reactions

The electronic nature of low-barrier hydrogen bonds (LBHBs) in enzymatic reactions is discussed based on combined low temperature neutron and x-ray diffraction experiments and on high level ab initio calculations by using the model substrate benzoylacetone. This molecule has a LBHB, as the intramolecular hydrogen bond is described by a double-well potential with a small barrier for hydrogen transfer. From an “atoms in molecules” analysis of the electron density, it is found that the hydrogen atom is stabilized by covalent bonds to both oxygens. Large atomic partial charges on the hydrogen-bonded atoms are found experimentally and theoretically. Therefore, the hydrogen bond gains stabilization from both covalency and from the normal electrostatic interactions found for long, weak hydrogen bonds. Based on comparisons with other systems having short-strong hydrogen bonds or LBHBs, it is proposed that all short-strong and LBHB systems possess similar electronic features of the hydrogen-bonded region, namely polar covalent bonds between the hydrogen atom and both heteroatoms in question.

The mechanism of proton pumping by cytochrome c oxidase

Cytochrome c oxidase catalyzes the reduction of oxygen to water that is accompanied by pumping of four protons across the mitochondrial or bacterial membrane. Triggered by the results of recent x-ray crystallographic analyses, published data concerning the coupling of individual electron transfer steps to proton pumping are reanalyzed: Conversion of the conventional oxoferryl intermediate F to the fully oxidized form O is connected to pumping of only one proton. Most likely one proton is already pumped during the double reduction of O, and only three protons during conversion of the “peroxy” forms P to O via the oxoferryl form F. Based on the available structural, spectroscopic, and mutagenesis data, a detailed mechanistic model, carefully considering electrostatic interactions, is presented. In this model, each of the four reductions of heme a during the catalytic cycle is coupled to the uptake of one proton via the D-pathway. These protons, but never more than two, are temporarily stored in the regions of the heme a and a3 propionates and are driven to the outside (“pumped”) by electrostatic repulsion from protons entering the active site during turnover. The first proton is pumped by uptake of one proton via the K-pathway during reduction, the second and third proton during the P → F transition when the D-pathway and the active site become directly connected, and the fourth one upon conversion of F to O. Atomic structures are assigned to each intermediate including F′ with an alternative route to O.

The primary fibrin polymerization pocket: Three-dimensional structure of a 30-kDa C-terminal γ chain fragment complexed with the peptide Gly-Pro-Arg-Pro

After vascular injury, a cascade of serine protease activations leads to the conversion of the soluble fibrinogen molecule into fibrin. The fibrin monomers then polymerize spontaneously and noncovalently to form a fibrin gel. The primary interaction of this polymerization reaction is between the newly exposed N-terminal Gly-Pro-Arg sequence of the α chain of one fibrin molecule and the C-terminal region of a γ chain of an adjacent fibrin(ogen) molecule. In this report, the polymerization pocket has been identified by determining the crystal structure of a 30-kDa C-terminal fragment of the fibrin(ogen) γ chain complexed with the peptide Gly-Pro-Arg-Pro. This peptide mimics the N terminus of the α chain of fibrin. The conformational change in the protein upon binding the peptide is subtle, with electrostatic interactions primarily mediating the association. This is consistent with biophysical experiments carried out over the last 50 years on this fundamental polymerization reaction.

Molecular dynamics and free-energy calculations applied to affinity maturation in antibody 48G7

We investigated the relative free energies of hapten binding to the germ line and mature forms of the 48G7 antibody Fab fragments by applying a continuum model to structures sampled from molecular dynamics simulations in explicit solvent. Reasonable absolute and very good relative free energies were obtained. As a result of nine somatic mutations that do not contact the hapten, the affinity-matured antibody binds the hapten >104 tighter than the germ line antibody. Energetic analysis reveals that van der Waals interactions and nonpolar contributions to solvation are similar and drive the formations of both the germ line and mature antibody–hapten complexes. Affinity maturation of the 48G7 antibody therefore appears to occur through reorganization of the combining site geometry in a manner that optimizes the balance of gaining favorable electrostatic interactions with the hapten and losing those with solvent during the binding process. As reflected by lower rms fluctuations in the antibody–hapten complex, the mature complex undergoes more restricted fluctuations than the germ line complex. The dramatically increased affinity of the 48G7 antibody over its germ line precursor is thus made possible by electrostatic optimization.

Proton uptake by bacterial reaction centers: The protein complex responds in a similar manner to the reduction of either quinone acceptor

In bacterial photosynthetic reaction centers, the protonation events associated with the different reduction states of the two quinone molecules constitute intrinsic probes of both the electrostatic interactions and the different kinetic events occurring within the protein in response to the light-generated introduction of a charge. The kinetics and stoichiometries of proton uptake on formation of the primary semiquinone QA− and the secondary acceptor QB− after the first and second flashes have been measured, at pH 7.5, in reaction centers from genetically modified strains and from the wild type. The modified strains are mutated at the L212Glu and/or at the L213Asp sites near QB; some of them carry additional mutations distant from the quinone sites (M231Arg → Leu, M43Asn → Asp, M5Asn → Asp) that compensate for the loss of L213Asp. Our data show that the mutations perturb the response of the protein system to the formation of a semiquinone, how distant compensatory mutations can restore the normal response, and the activity of a tyrosine residue (M247Ala → Tyr) in increasing and accelerating proton uptake. The data demonstrate a direct correlation between the kinetic events of proton uptake that are observed with the formation of either QA− or QB−, suggesting that the same residues respond to the generation of either semiquinone species. Therefore, the efficiency of transferring the first proton to QB is evident from examination of the pattern of H+/QA− proton uptake. This delocalized response of the protein complex to the introduction of a charge is coordinated by an interactive network that links the Q− species, polarizable residues, and numerous water molecules that are located in this region of the reaction center structure. This could be a general property of transmembrane redox proteins that couple electron transfer to proton uptake/release reactions.

Structural basis of DNA bending and oriented heterodimer binding by the basic leucine zipper domains of Fos and Jun

Interactions among transcription factors that bind to separate sequence elements require bending of the intervening DNA and juxtaposition of interacting molecular surfaces in an appropriate orientation. Here, we examine the effects of single amino acid substitutions adjacent to the basic regions of Fos and Jun as well as changes in sequences flanking the AP-1 site on DNA bending. Substitution of charged amino acid residues at positions adjacent to the basic DNA-binding domains of Fos and Jun altered DNA bending. The change in DNA bending was directly proportional to the change in net charge for all heterodimeric combinations between these proteins. Fos and Jun induced distinct DNA bends at different binding sites. Exchange of a single base pair outside of the region contacted in the x-ray crystal structure altered DNA bending. Substitution of base pairs flanking the AP-1 site had converse effects on the opposite directions of DNA bending induced by homodimers and heterodimers. These results suggest that Fos and Jun induce DNA bending in part through electrostatic interactions between amino acid residues adjacent to the basic region and base pairs flanking the AP-1 site. DNA bending by Fos and Jun at inverted binding sites indicated that heterodimers bind to the AP-1 site in a preferred orientation. Mutation of a conserved arginine within the basic regions of Fos and transversion of the central C:G base pair in the AP-1 site to G:C had complementary effects on the orientation of heterodimer binding and DNA bending. The conformational variability of the Fos–Jun–AP-1 complex may contribute to its functional versatility at different promoters.

pH gating of ROMK (Kir1.1) channels: Control by an Arg-Lys-Arg triad disrupted in antenatal Bartter syndrome

Inward-rectifier K+ channels of the ROMK (Kir1.1) subtype are responsible for K+ secretion and control of NaCl absorption in the kidney. A hallmark of these channels is their gating by intracellular pH in the neutral range. Here we show that a lysine residue close to TM1, identified previously as a structural element required for pH-induced gating, is protonated at neutral pH and that this protonation drives pH gating in ROMK and other Kir channels. Such anomalous titration of this lysine residue (Lys-80 in Kir1.1) is accomplished by the tertiary structure of the Kir protein: two arginines in the distant N and C termini of the same subunit (Arg-41 and Arg-311 in Kir1.1) are located in close spatial proximity to the lysine allowing for electrostatic interactions that shift its pKa into the neutral pH range. Structural disturbance of this triad as a result from a number of point mutations found in patients with antenatal Bartter syndrome shifts the pKa of the lysine residue off the neutral pH range and results in channels permanently inactivated under physiological conditions. Thus, the results provide molecular understanding for normal pH gating of Kir channels as well as for the channel defects found in patients with antenatal Bartter syndrome.

Role of the J-domain in the cooperation of Hsp40 with Hsp70

The Escherichia coli Hsp40 DnaJ and Hsp70 DnaK cooperate in the binding of proteins at intermediate stages of folding, assembly, and translocation across membranes. Binding of protein substrates to the DnaK C-terminal domain is controlled by ATP binding and hydrolysis in the N-terminal ATPase domain. The interaction of DnaJ with DnaK is mediated at least in part by the highly conserved N-terminal J-domain of DnaJ that includes residues 2–75. Heteronuclear NMR experiments with uniformly 15N-enriched DnaJ2–75 indicate that the chemical environment of residues located in helix II and the flanking loops is perturbed on interaction with DnaK or a truncated DnaK molecule, DnaK2–388. NMR signals corresponding to these residues broaden and exhibit changes in chemical shifts in the presence of DnaK(MgADP). Addition of MgATP largely reversed the broadening, indicating that NMR signals of DnaJ2–75 respond to ATP-dependent changes in DnaK. The J-domain interaction is localized to the ATPase domain of DnaK and is likely to be dominated by electrostatic interactions. The results suggest that the J-domain tethers DnaK to DnaJ-bound substrates, which DnaK then binds with its C-terminal peptide-binding domain.

Electric field-induced reorganization of two-component supported bilayer membranes

Application of electric fields tangent to the plane of a confined patch of fluid bilayer membrane can create lateral concentration gradients of the lipids. A thermodynamic model of this steady-state behavior is developed for binary systems and tested with experiments in supported lipid bilayers. The model uses Flory’s approximation for the entropy of mixing and allows for effects arising when the components have different molecular areas. In the special case of equal area molecules the concentration gradient reduces to a Fermi–Dirac distribution. The theory is extended to include effects from charged molecules in the membrane. Calculations show that surface charge on the supporting substrate substantially screens electrostatic interactions within the membrane. It also is shown that concentration profiles can be affected by other intermolecular interactions such as clustering. Qualitative agreement with this prediction is provided by comparing phosphatidylserine- and cardiolipin-containing membranes.