Microwave dielectric relaxation in muscle. A second look.

The dielectric permittivity and conductivity of muscle fibers from the giant barnacle, Balanus nubilus, have been measured at 1, 25, and 37 degrees C, between 10 MHz and 17 GHz. The dominant microwave dielectric relaxation process in these fibers is due to dipolar relaxation of the tissue water, which shows a characteristic relaxation frequency equal to that of pure water, ranging from 9 GHz (1 degree C) to 25 GHz (37 degree C). The total permittivity decrease, epsilon 0 -- epsilon infinity, due to this process accounts for approximately 95% of the water content of the tissue; thus, the major fraction of tissue water is dielectrically identical to the pure fluid on a picosecond time scale. A second dielectric process contributes significantly to the tissue dielectric properties between 0.1 and 1--5 GHz, and arises in part form Maxwell-Wagner effects due to the electrolyte content of the tissue, and in part from dielectric relaxation of the tissue proteins themselves.     

The effect of protein relaxation on charge-charge interactions and dielectric constants of proteins.

The effect of the reorganization of the protein polar groups on charge-charge interaction and the corresponding effective dielectric constant (epsilon(eff)) is examined by the semimicroscopic version of the Protein Dipole Langevin Dipoles (PDLD/S) method within the framework of the Linear Response Approximation (LRA). This is done by evaluating the interactions between ionized residues in the reaction center of Rhodobacter sphaeroides, while taking into account the protein reorganization energy. It is found that an explicit consideration of the protein relaxation leads to a significant increase in epsilon(eff) and that semimicroscopic models that do not take this relaxation into account force one to use a large value for the so-called "protein dielectric constant," epsilon(p), of the Poisson-Boltzmann model or for the corresponding epsilon(in) in the PDLD/S model. An additional increase in epsilon(eff) is expected from the reorganization of ionized residues and from changes in the degree of water penetration. This finding provides further support for the idea that epsilon(in) (or epsilon(p)) represents contributions that are not considered explicitly. The present study also provides a systematic illustration of the nature of epsilon(eff), supporting our previously reported view that charge-charge interactions correspond to a large value of this "dielectric constant," even in protein interiors. It is also pointed out that epsilon(eff) for the interaction between ionizable groups in proteins is very different from the effective dielectric constant, epsilon'(eff), that determines the free energy of ion pairs in proteins (epsilon'(eff) reflects the effect of preoriented protein dipoles). Finally, the problems associated with the search for a general epsilon(in) are discussed. It is clarified that the epsilon(in) that reproduces the effect of protein relaxation on charge-charge interaction is not equal to the epsilon(in) that reproduces the corresponding effect upon formation of individual charges. This reflects fundamental inconsistencies in attempts to cast microscopic concepts in a macroscopic model. Thus one should either use a large epsilon(in) for charge-charge interactions and a small epsilon(in) for charge-dipole interactions or consider the protein relaxation microscopically.     

Comparison of dielectric response models for simulating electrostatic effects in proteins

Two recent approaches for calculating pK shifts in proteins are compared. The first of these uses Coulomb's law with a distance-dependent dielectric permittivity, epsilon (r), to model the screening effects of the environment, and the second uses a finite difference approach to solve Poisson's equation. It is shown that an explicit form of epsilon (r) which has been fitted to experimentally determined values of the dielectric permittivity in a range from 1 to 21 A can be approximated by a linear form in the functionally significant range of charge separations of approximately 3-10 A, but for distances greater than 10 A the effective permittivity is strongly nonlinear. A statistical analysis of the errors in calculated pK shifts due to electrostatic interactions between charges with separations greater than 10 A shows that there are only marginal differences in reliability between using Coulomb's law with an appropriate form of epsilon (r) or the finite difference approach for solving Poisson's equation. Thus it is concluded that pK shifts can be calculated just as well, and with considerably less effort, using Coulomb's law.     

Modeling implicit reorganization in continuum descriptions of protein-protein interactions

The determination of free energies that govern protein-protein recognition is essential for a detailed molecular understanding of biological specificity. Continuum models of macromolecular interactions, in which the solvent is treated by an implicit representation and the proteins are treated semi-microscopically, are computationally tractable for estimating free energies, yet many questions remain concerning their accuracy. This article reports a continuum analysis of the free-energy changes underlying the binding of 31 interfacial alanine substitutions of two complexes of the antihen egg white lysozyme (HEL) antibody D1.3 bound with HEL or the antibody E5.2. Two implicit schemes for modeling the effects of protein and solvent relaxation were examined, in which the protein environment was treated as either homogeneous with a "protein dielectric constant" of epsilon(p) = 4 or inhomogeneous, with epsilon(p) = 4 for neutral residues and epsilon(p) = 25 for ionized residues. The results showed that the nonuniform dielectric model reproduced the experimental differences better, with an average absolute error of +/-1.1 kcal/mol, compared with +/-1.4 kcal/mol for the uniform model. More importantly, the error for charged residues in the nonuniform model is +/-0.8 kcal/mol and is nearly half of that corresponding to the uniform model. Several substitutions were clearly problematic in determining qualitative trends and probably required explicit structural reorganization at the protein-protein interface.     

Assembly and function of the T cell antigen receptor. Requirement of either the lysine or arginine residues in the transmembrane region of the alpha chain

The T cell receptor (TCR) for antigen consists, on the majority of peripheral lymphocytes, of an immunoglobulin-like, disulfide-linked heterodimeric glycoprotein: the alpha and beta chain. These proteins are noncovalently linked to at least four nonvariant proteins which comprise the CD3 complex: CD3 gamma, delta, epsilon, and zeta. Whereas the TCR alpha and beta proteins have positively charged residues in the transmembrane region, all the CD3 proteins have similarly placed negatively charged amino acid residues. It has been suggested that these basic and acidic amino acid residues may play an important role in TCR.CD3 complex assembly and/or function. In this paper, the structural and functional role of the lysine and arginine residues of the TCR alpha chain was addressed using oligonucleotide mediated site directed mutagenesis. The Arg256 and Lys261 residues of the TCR alpha cDNA of the HPB-ALL cell line were mutated to either Gly256 and/or Ile261. The altered cDNAs were transfected into a TCR alpha negative recipient mutant cell line of REX, clone 20A. Metabolic labeling of the T cell transfectants showed that mutation of either the Arg256 or Lys261 amino acid residues had no effect on the ability of the TCR alpha chain to form either a heterodimer with the TCR beta chain or a complex with the CD3 gamma, delta, and epsilon proteins. Consequently, the Arg256 to Gly256 and Lys261 to Ile261 mutations did not prevent the formation of a mature, functional TCR.CD3 complex on the cell surface as determined by immunofluorescence, cell surface radioiodination, and the ability of the transfectants to mobilize intracellular calcium after stimulation with a mitogenic anti-CD3 epsilon monoclonal antibody. In contrast, a mutant cDNA in which both the Arg256 and Lys261 residues were mutated to Gly256 and Ile261, respectively, failed to reconstitute the cell surface expression of the TCR.CD3 complex and, consequently, the ability to respond to mitogenic stimuli. In the absence of both the Arg256 and Lys261 residues, TCR alpha beta heterodimer formation was not observed. Cotransfection studies in COS cells showed that the failure of assembly of a heterodimer was likely due to an inability of the mutated TCR alpha chain to form a subcomplex with either the CD3 gamma, delta, epsilon, or zeta proteins.(ABSTRACT TRUNCATED AT 400 WORDS)     

Semi-continuum electrostatic calculations of redox potentials in photosystem I

The midpoint redox potentials (E(m)) of all cofactors in photosystem I from Synechococcus elongatus as well as of the iron-sulfur (Fe(4)S(4)) clusters in two soluble ferredoxins from Azotobacter vinelandii and Clostridium acidiurici were calculated within the framework of a semi-continuum dielectric approach. The widely used treatment of proteins as uniform media with single dielectric permittivity is oversimplified, particularly, because permanent charges are considered both as a source for intraprotein electric field and as a part of dielectric polarizability. Our approach overcomes this inconsistency by using two dielectric constants: optical epsilon(o)=2.5 for permanent charges pre-existing in crystal structure, and static epsilon(s) for newly formed charges. We also take into account a substantial dielectric heterogeneity of photosystem I revealed by photoelectric measurements and a liquid junction potential correction for E(m) values of relevant redox cofactors measured in aprotic solvents. We show that calculations based on a single permittivity have the discrepancy with experimental data larger than 0.7 V, whereas E(m) values calculated within our approach fall in the range of experimental estimates. The electrostatic analysis combined with quantum chemistry calculations shows that (i) the energy decrease upon chlorophyll dimerization is essential for the downhill mode of primary charge separation between the special pair P(700) and the primary acceptor A(0); (ii) the primary donor is apparently P(700) but not a pair of accessory chlorophylls; (iii) the electron transfer from the A branch quinone Q(A) to the iron-sulfur cluster F(X) is most probably downhill, whereas that from the B branch quinone Q(B) to F(X) is essentially downhill.     

Correlation of electron transfer rate in photosynthetic reaction centers with intraprotein dielectric properties

A number of the electrogenic reactions in photosystem I, photosystem II, and bacterial reaction centers (RC) were comparatively analyzed, and the variation of the dielectric permittivity (epsilon) in the vicinity of electron carriers along the membrane normal was calculated. The value of epsilon was minimal at the core of the complexes and gradually increased towards the periphery. We found that the rate of electron transfer (ET) correlated with the value of the dielectric permittivity: the fastest primary ET reactions occur in the low-polarity core of the complexes within the picosecond time range, whereas slower secondary reactions take place at the high-polarity periphery of the complexes within micro- to millisecond time range. The observed correlation was quantitatively interpreted in the framework of the Marcus theory. We calculated the reorganization energy of ET carriers using their van der Waals volumes and experimentally determined epsilon values. The electronic coupling was calculated by the empirical Moser-Dutton rule for the distance-dependent electron tunneling rate in nonadiabatic ET reactions. We concluded that the local dielectric permittivity inferred from the electrometric measurements could be quantitatively used to estimate the rate constant of ET reactions in membrane proteins with resolved atomic structure with the accuracy of less than one order of magnitude.     

Comparison of calculation and experiment implicates significant electrostatic contributions to the binding stability of barnase and barstar

The contributions of electrostatic interactions to the binding stability of barnase and barstar were studied by the Poisson-Boltzmann model with three different protocols: a), the dielectric boundary specified as the van der Waals (vdW) surface of the protein along with a protein dielectric constant (epsilon (p)) of 4; b), the dielectric boundary specified as the molecular (i.e., solvent-exclusion (SE)) surface along with epsilon (p) = 4; and c), "SE + epsilon (p) = 20." The "vdW + epsilon (p) = 4" and "SE + epsilon (p) = 20" protocols predicted an overall electrostatic stabilization whereas the "SE + epsilon (p) = 4" protocol predicted an overall electrostatic destabilization. The "vdW + epsilon (p) = 4" protocol was most consistent with experiment. It quantitatively reproduced the observed effects of 17 mutations neutralizing charged residues lining the binding interface and the measured coupling energies of six charge pairs across the interface and reasonably rationalized the experimental ionic strength and pH dependences of the binding constant. In contrast, the "SE + epsilon (p) = 4" protocol predicted significantly larger coupling energies of charge pairs whereas the "SE + epsilon (p) = 20" protocol did not predict any pH dependence. This study calls for further scrutiny of the different Poisson-Boltzmann protocols and demonstrates potential danger in drawing conclusions on electrostatic contributions based on a particular calculation protocol.     

Effect of different low electric field conditions on the dielectric properties of Ehrlich tumor

Studies of low electric fields (LEFs) effects on the permeability of the cell membrane are of great interest in molecular medicine. Electroendocytosis is a novel technique depends on using LEFs to incorporate macromolecules as anticancer drugs or genes into the cells. There are wide debates about the optimum electric conditions for electroendocytosis. In this article, Ehrlich tumor tissues were exposed to different LEFs voltages and frequencies in vitro. Dielectric properties before and after the exposure were determined. The results indicated that the exposed groups have significant high permittivity and conductivity compared to unexposed group, as well as having significant low impedance. The results indicated that dielectric measurements can be used to indicate the efficiency of electroendocytosis that as permittivity and conductivity of cell membranes increase, more molecules can passed into the cells. It was also indicated that, as the pulse amplitude increases, the LEFs influence increases, while changing pulse frequency has no obvious effect on dielectric properties of Ehrlich tumor.     

An electrostatic basis for the stability of thermophilic proteins

Two factors provide key contributions to the stability of thermophilic proteins relative to their mesophilic homologues: electrostatic interactions of charged residues in the folded state and the dielectric response of the folded protein. The dielectric response for proteins in a "thermophilic series" globally modulates the thermal stability of its members, with the calculated dielectric constant for the protein increasing from mesophiles to hyperthermophiles. This variability results from differences in the distribution of charged residues on the surface of the protein, in agreement with structural and genetic observations. Furthermore, the contribution of electrostatic interactions to the stability of the folded state is more favorable for thermophilic proteins than for their mesophilic homologues. This leads to the conclusion that electrostatic interactions play an important role in determining the stability of proteins at high temperatures. The interplay between electrostatic interactions and dielectric response also provides further rationalization for the enhanced stability of thermophilic proteins with respect to cold-denaturation. Taken together, the distribution of charged residues and their fluctuations have been shown to be factors in modulating protein stability over the entire range of biologically relevant temperatures.     

Spatial distribution of dielectric shielding in the interior of Pyrococcus furiosus rubredoxin as sampled in the subnanosecond timeframe by hydrogen exchange

Experimental pK values of ionizable sidechains provide the most direct test for models representing dielectric shielding within the interior of a protein. However, only the strongly shifted pK values are particularly useful for discriminating among models. NMR titration studies have usually found only one or two such shifted pK values in each protein, so that the fitting of the experimental data to a uniform internal dielectric (epsilon(int)) model is not well constrained. The observed variation among proteins for such epsilon(int) estimates may reflect nonuniformity of dielectric shielding within each protein interior or qualitative differences between individual proteins. The differential amide kinetic acidities for a series of metal-substituted rubredoxins are shown to be consistent with Poisson-Boltzmann predictions of dielectric shielding that is relatively uniform for all of the amides that are sensitive to the metal charge, a region which corresponds to roughly 1/3 of the internal volume. The effective epsilon(int) values near 6 that are found in this study are significantly lower than many such estimates derived from sidechain pK measurements. The differing timeframes in which dielectric relaxation can respond to the highly transient peptide anion as compared to the longer lived states of the charged sidechains offers an explanation for the lower apparent dielectric constant deduced from these measurements.     

NMR studies of myelin basic protein. XIII. Assignment of histidine residues in rabbit, bovine and porcine proteins

Myelin basic protein from three species (rabbit, cow and pig) and peptides from enzymatic digests or cleavage of the proteins have been examined in aqueous solutions by proton nuclear magnetic resonance (NMR) at 400 MHz. The epsilon 1-CH and delta 2-CH resonances of all the histidine residues in the three proteins have been assigned and the pK values have been measured. The heterogeneity of chemical shifts among these resonances can be variously ascribed to persistent localized secondary structures and to effects arising from charged side-chains, particularly those of aspartic acid residues, and from side-chains of aromatic moieties.     

Dielectric properties of E. coli cell as simulated by the three-shell spheroidal model

Dielectric properties of E. coli cell have been re-studied by means of the three-shell spheroidal model, where the three shells correspond to the outer membrane, the periplasmic space and the inner membrane, respectively. With the model, a curve-fitting procedure has been developed to analyze the dielectric spectra. Although E. coli cell has been studied before, its special morphological structure was taken into account more comprehensively than any previous model in the present work. Dielectric properties of various cell components have been estimated from the observed dielectric spectra, especially the permittivity of the outer membrane, which was evaluated quantitatively for the first time. The values of epsilon(om) were 12 for kappa(om) of 0 to 10(-4) S/m and 34 for kappa(om) of 10(-3) S/m. The specific capacitance of the inner membrane was 0.6-0.70 microF/cm(2). The relative permittivity and the conductivity of the cytoplasm were about 100 and 0.22 S/m, respectively, and the conductivity of the periplasmic space was 2.2-3.2 S/m.     

Electrostatic effects on the alpha-helix and beta-strand formation of BPTI(16-36) studied by Monte Carlo simulated annealing.

The electrostatic effects on the secondary structure forming tendencies of a peptide fragment with residues 16-36 of bovine pancreatic trypsin inhibitor, BPTI(16-36), are studied using Monte Carlo simulated annealing simulations. We consider three dielectric functions epsilon(r) of distance r: constant dielectric function (epsilon = 2; strong electrostatic interactions) and sigmoidal functions varying from epsilon(0) = 2 to epsilon(infinity) = 47 (intermediate) and to epsilon(infinity) = 78 (weak). Simulations with epsilon = 2 suggest that this peptide exhibits a significant propensity for beta-strand formations in accordance with a beta-sheet structure of the relevant segment in native BPTI. The tendency for alpha-helix formations becomes almost comparable with that of beta-strands in the simulation with epsilon(infinity) = 47, and there appears no appreciable conformational propensity for this case. Finally, the results with epsilon(infinity) = 78 generate low-energy conformations with conspicuous alpha-helices. These findings suggest the possibility that the change in electrostatic interactions can be the key factor for the conformational transitions of peptides between alpha-helix and beta-sheet that have recently been observed in experiments. These changes in electrostatic interactions can arise from those in various environmental factors such as conformations of the rest of the protein molecule and solvent conditions.     

A frequency-dependent finite-difference time-domain formulation for induced current calculations in human beings.

The finite-difference time-domain (FDTD) method has been used to calculate SARs and induced currents involving whole-body or partial-body exposures of models to spatially uniform or nonuniform (far-field or near-field), to sinusoidally varying EM fields, or to transient fields such as those associated with electromagnetic pulses. However, a weakness of the FDTD algorithm is that the dispersion of the tissue's dielectric properties is ignored and frequency-independent properties are assumed. Although this is permissible for continuous-wave or narrow-band irradiation, the results may be highly erroneous for short pulses, in which ultra-wide bandwidths are involved. In some recent publications, procedures are described for one- and two-dimensional problems for media in which the complex permittivity epsilon * (omega) may be described by a single-order Debye relaxation equation or a modified version thereof. These procedures based on a convolution integral describing D(t) in terms of E(t) cannot be extended to human tissues for which multiterm Debye relaxation equations must generally be used. We describe here a new differential-equation approach that can be used for general dispersive media. We illustrate the use of this approach by one- and three-dimensional examples of media for which epsilon * (omega) is given by a multiterm Debye equation, and for an approximate two-thirds muscle-equivalent model of the human body. Based on a single run involving a Gaussian pulse, the frequency-dependent FDTD [(FD)2TD] method allows calculations of SARs and induced currents at various frequencies by taking the Fourier components of the induced E fields. The (FD)2TD method can also be used to calculate coupling of the short (ultra-wideband) pulses to the human body.     

Combining conformational flexibility and continuum electrostatics for calculating pK(a)s in proteins

Protein stability and function relies on residues being in their appropriate ionization states at physiological pH. In situ residue pK(a)s also provides a sensitive measure of the local protein environment. Multiconformation continuum electrostatics (MCCE) combines continuum electrostatics and molecular mechanics force fields in Monte Carlo sampling to simultaneously calculate side chain ionization and conformation. The response of protein to charges is incorporated both in the protein dielectric constant (epsilon(prot)) of four and by explicit conformational changes. The pK(a) of 166 residues in 12 proteins was determined. The root mean square error is 0.83 pH units, and >90% have errors of <1 pH units whereas only 3% have errors >2 pH units. Similar results are found with crystal and solution structures, showing that the method's explicit conformational sampling reduces sensitivity to the initial structure. The outcome also changes little with protein dielectric constant (epsilon(prot) 4-20). Multiconformation continuum electrostatics titrations show coupling of conformational flexibility and changes in ionization state. Examples are provided where ionizable side chain position (protein G), Asn orientation (lysozyme), His tautomer distribution (RNase A), and phosphate ion binding (RNase A and H) change with pH. Disallowing these motions changes the calculated pK(a).     

Steric selectivity in Na channels arising from protein polarization and mobile side chains

Monte Carlo simulations of equilibrium selectivity of Na channels with a DEKA locus are performed over a range of radius R and protein dielectric coefficient epsilon(p). Selectivity arises from the balance of electrostatic forces and steric repulsion by excluded volume of ions and side chains of the channel protein in the highly concentrated and charged (approximately 30 M) selectivity filter resembling an ionic liquid. Ions and structural side chains are described as mobile charged hard spheres that assume positions of minimal free energy. Water is a dielectric continuum. Size selectivity (ratio of Na+ occupancy to K+ occupancy) and charge selectivity (Na+ to Ca2+) are computed in concentrations as low as 10(-5) M Ca2+. In general, small R reduces ion occupancy and favors Na+ over K+ because of steric repulsion. Small epsilon(p) increases occupancy and favors Na+ over Ca2+ because protein polarization amplifies the pore's net charge. Size selectivity depends on R and is independent of epsilon(p); charge selectivity depends on both R and epsilon(p). Thus, small R and epsilon(p) make an efficient Na channel that excludes K+ and Ca2+ while maximizing Na+ occupancy. Selectivity properties depend on interactions that cannot be described by qualitative or verbal models or by quantitative models with a fixed free energy landscape.     

The origin and impact of bound water around intrinsically disordered proteins

Proteins and water couple dynamically over a wide range of time scales. Motivated by their central role in protein function, protein-water dynamics and thermodynamics have been extensively studied for structured proteins, where correspondence to structural features has been made. However, properties controlling intrinsically disordered protein (IDP)-water dynamics are not yet known. We report results of megahertz-to-terahertz dielectric spectroscopy and molecular dynamics simulations of a group of IDPs with varying charge content along with structured proteins of similar size. Hydration water around IDPs is found to exhibit more heterogeneous rotational and translational dynamics compared with water around structured proteins of similar size, yielding on average more restricted dynamics around individual residues of IDPs, charged or neutral, compared with structured proteins. The on-average slower water dynamics is found to arise from excess tightly bound water in the first hydration layer, which is related to greater exposure to charged groups. The more tightly bound water to IDPs correlates with the smaller hydration shell found experimentally, and affects entropy associated with protein-water interactions, the contribution of which we estimate based on the dielectric measurements and simulations. Water-IDP dynamic coupling at terahertz frequencies is characterized by the dielectric measurements and simulations.     

Dielectric increment and dielectric dispersion of solutions containing simple charged linear macromolecules. I. Theory

A theory is derived for the static and frequency dependent value of the electric permittivity for model systems representing a solution of a macromolecule bearing a large number of identical charges. The polyion is represented either as a charged rigid rod (A) or as a sequence of charged rodlike subunits in an arbitrary but fixed configuration (B) and it is assumed that a certain fraction of the counterions is closely associated to the macromolecule. The dielectric properties are described in terms of fluctuations in the distribution of the associated counterions along the polyion. These fluctuations can occur locally between potential barriers marking the ends of the subunits (if considered) but can also extend over the whole molecule. Neglecting correlations between different associated counterions expressions for the static value of the dielectric increment are obtained which reveal its dependence on the fraction of bound ions, on the charge of the counterions and on the length of the molecule for model A or the radius of gyration for model B. The dynamic behaviour of A is distinguishable from that of B as the former will present one single dispersion curve of the frequency dependent electric permittivity while the latter may give rise to two different dispersion regions. This will be the case if both the exchange between bound and free ions and the rotation of the complete molecule are relatively slow in comparison to the local bound counterion density fluctuations and if these fluctuations occur on a much shorter time scale than the ion density fluctuations extending over the complete macromolecule.