Intrinsic protein electric fields: basic non-covalent interactions and relationship to protein-induced Stark effects

Rethinking Vibrational Stark Spectroscopy: Peak Shifts, Line Widths, and the Role of Non-Stark Solvent Coupling

A vibration's transition frequency is partly determined by the first-order Stark effect, which accounts for the electric field experienced by the mode. Using ultrafast infrared pump-probe and FT-IR spectroscopies, we characterized both the 0 → 1 and 1 → 2 vibrational transitions' field-dependent peak positions and line widths of the CN stretching mode of benzonitrile (BZN) and phenyl selenocyanate (PhSeCN) in ten solvents. We present a theoretical model that decomposes the observed line width into a field-dependent Stark contribution and a field-independent non-Stark solvent coupling contribution (NSC). The model demonstrates that the field-dependent peak position is independent of the line width, even when the NSC dominates the latter. Experiments show that when the Stark tuning rate is large compared to the NSC (PhSeCN), the line width has a field dependence, albeit with major NSC-induced excursions from linearity. When the Stark tuning rate is small relative to the NSC (BZN), the line width is field-independent. BZN's line widths are substantially larger for the 1 → 2 transition, indicating a 1 → 2 transition enhancement of the NSC. Additionally, we examine, theoretically and experimentally, the difference in the 0 → 1 and 1 → 2 transitions' Stark tuning rates. Second-order perturbation theory combined with density functional theory explain the difference and show that the 1 → 2 transition's Stark tuning rate is ∼10% larger. The Stark tuning rate of PhSeCN is larger than BZN's for both transitions, consistent with the theoretical calculations. This study provides new insights into vibrational line shape components and a more general understanding of the vibrational response to external electric fields.

Frequency-Dependent Polarizabilities of Amino Acids as Calculated by an Electrostatic Interaction Model

The frequency-dependent polarizability of the 20 essential amino acids has been calculated by an electrostatic interaction model where an Unsöld-type of model has been adopted for the frequency dependence. The interaction model has previously been parametrized from Hartree-Fock calculations on a set of molecules, and the model is in this work extended by sulfur parameters by including a set of 18 small sulfur compounds. The results for the amino acids by using the interaction model compare well with Hartree-Fock calculations with deviations of around 5% for the isotropic polarizability. Furthermore, the intrinsic (or optical) dielectric constant related to the polarizability has been calculated for three small proteins, ribonuclease inhibitor, lysozyme, and green fluorescent protein, adopting the interaction model. The results are consistent with the intrinsic dielectric constants found for proteins in the literature.

Critical role of RAGE in lung physiology and tumorigenesis: a potential target of therapeutic intervention?

Lung cancer is one of the most common malignancies in the world and one of the leading causes of death from cancer. In the search for molecules that may be involved in lung tumor induction and progression, the receptor of advanced glycation end products (RAGE) comes across as a critical regulator of lung physiology. RAGE is a multiligand receptor that presents a differential expression pattern in lung epithelial cells compared to other cell types being gradually increased from fetal to birth and adult life. Under stress conditions, RAGE expression and activation are rapidly elevated resulting in chronic inflammation, which, in turn, in many instances, promotes epithelial cell malignant transformation. RAGE overexpression in normal lung alveolar type I epithelial cells is followed by rapid downregulation upon malignant transformation, being associated with increased aggressiveness. This is a striking paradox, since in every other cell type the pattern of RAGE expression follows the opposite direction, suggesting the involvement of RAGE in the well-functioning of lung cells. Additionally, RAGE has been attributed with the role of adhesion molecule, since it can stabilize mature alveolar epithelial cells to their substrate (basal lamina) by interacting electrostatically with other molecules. However, the reduction of RAGE observed in lung tumorigenesis interrupts cell-to-cell and cell-to-substrate communication, which is a critical step for cancer cell induction, progression and migration. This review addresses the differential properties of RAGE in lung physiology and carcinogenesis, providing evidence of therapeutic possibilities.

Model for vaccine design by prediction of B-epitopes of IEDB given perturbations in peptide sequence, in vivo process, experimental techniques, and source or host organisms

Perturbation methods add variation terms to a known experimental solution of one problem to approach a solution for a related problem without known exact solution. One problem of this type in immunology is the prediction of the possible action of epitope of one peptide after a perturbation or variation in the structure of a known peptide and/or other boundary conditions (host organism, biological process, and experimental assay). However, to the best of our knowledge, there are no reports of general-purpose perturbation models to solve this problem. In a recent work, we introduced a new quantitative structure-property relationship theory for the study of perturbations in complex biomolecular systems. In this work, we developed the first model able to classify more than 200,000 cases of perturbations with accuracy, sensitivity, and specificity >90% both in training and validation series. The perturbations include structural changes in >50000 peptides determined in experimental assays with boundary conditions involving >500 source organisms, >50 host organisms, >10 biological process, and >30 experimental techniques. The model may be useful for the prediction of new epitopes or the optimization of known peptides towards computational vaccine design.

Modeling molecular response in uniform and non-uniform electric fields

The response of a molecule to an electric field E, often a model of environment, can be expressed in terms of a sum of power series expansions. We investigate the accuracy and limits of applicability of this expression using one-, two-, and three-dimensional models of the hydrogen-bonded complex, ClH:NH(3). Energetic, structural, and vibrational spectroscopic characteristics are determined at first- and second-order in E and [nabla]E and compared with ab initio values for a range of uniform and non-uniform electric fields chosen to simulate molecular environments. It is found that even at field strengths large enough to cause dramatic structural change in the complex, energetic, structural, and vibrational spectroscopic characteristics are accurately calculated using only terms linear in E and [nabla]E. These results suggest that knowledge of the zero-field molecular potential energy, dipole, and quadrupole moment surfaces may be sufficient to accurately model the interaction of a molecule with a wide range of chemical environments.

Hole-burning spectroscopy as a probe of nano-environments and processes in biomolecules: a review

Hole-burning spectroscopy, a high-resolution spectroscopic technique, allows details of heterogeneous nano-environments in biological systems to be obtained from broad absorption bands. Recently, this technique has been applied to proteins, nucleic acids, cells, and substructures of water to probe the electrostatic conditions created by macromolecules and the surrounding solvent. Starting with the factors that obscure the homogeneous linewidth of a chromophore within an inhomogeneously broadened absorption or emission band, we describe properties and processes in biological systems that are reflected in the measured hole spectra. The technique also lends itself to the resolution of perturbation experiments, such as temperature cycling to elucidate energy landscape barriers, applied external electric fields (Stark effect) to measure net internal electric fields, and applied hydrostatic pressure to find the volume compressibility of proteins.

Heterotrimerization of heat-shock factors 1 and 2 provides a transcriptional switch in response to distinct stimuli

Organisms respond to circumstances threatening the cellular protein homeostasis by activation of heat-shock transcription factors (HSFs), which play important roles in stress resistance, development, and longevity. Of the four HSFs in vertebrates (HSF1-4), HSF1 is activated by stress, whereas HSF2 lacks intrinsic stress responsiveness. The mechanism by which HSF2 is recruited to stress-inducible promoters and how HSF2 is activated is not known. However, changes in the HSF2 expression occur, coinciding with the functions of HSF2 in development. Here, we demonstrate that HSF1 and HSF2 form heterotrimers when bound to satellite III DNA in nuclear stress bodies, subnuclear structures in which HSF1 induces transcription. By depleting HSF2, we show that HSF1-HSF2 heterotrimerization is a mechanism regulating transcription. Upon stress, HSF2 DNA binding is HSF1 dependent. Intriguingly, when the elevated expression of HSF2 during development is mimicked, HSF2 binds to DNA and becomes transcriptionally competent. HSF2 activation leads to activation of also HSF1, revealing a functional interdependency that is mediated through the conserved trimerization domains of these factors. We propose that heterotrimerization of HSF1 and HSF2 integrates transcriptional activation in response to distinct stress and developmental stimuli.

The importance of vibronic perturbations in ferrocytochrome c spectra: a reevaluation of spectral properties based on low-temperature optical absorption, resonance Raman, and molecular-dynamics simulations

We have measured and analyzed the low-temperature (T=10 K) absorption spectrum of reduced horse heart and yeast cytochrome c. Both spectra show split and asymmetric Q(0) and Q(upsilon) bands. The spectra were first decomposed into the individual split vibronic sidebands assignable to B(1g) (nu15) and A(2g) (nu19, nu21, and nu22) Herzberg-Teller active modes due to their strong intensity in resonance Raman spectra acquired with Q(0) and Q(upsilon) excitations. The measured band splittings and asymmetries cannot be rationalized solely in terms of electronic perturbations of the heme macrocycle. On the contrary, they clearly point to the importance of considering not only electronic perturbations but vibronic perturbations as well. The former are most likely due to the heterogeneity of the electric field produced by charged side chains in the protein environment, whereas the latter reflect a perturbation potential due to multiple heme-protein interactions, which deform the heme structure in the ground and excited states. Additional information about vibronic perturbations and the associated ground-state deformations are inferred from the depolarization ratios of resonance Raman bands. The results of our analysis indicate that the heme group in yeast cytochrome c is more nonplanar and more distorted along a B(2g) coordinate than in horse heart cytochrome c. This conclusion is supported by normal structural decomposition calculations performed on the heme extracted from molecular-dynamic simulations of the two investigated proteins. Interestingly, the latter are somewhat different from the respective deformations obtained from the x-ray structures.

Optical band splitting and electronic perturbations of the heme chromophore in cytochrome C at room temperature probed by visible electronic circular dichroism spectroscopy

We have measured the electronic circular dichroism (ECD) of the ferri- and ferro-states of several natural cytochrome c derivatives (horse heart, chicken, bovine, and yeast) and the Y67F mutant of yeast in the region between 300 and 750 nm. Thus, we recorded the ECD of the B- and Q-band region as well as the charge-transfer band at approximately 695 nm. The B-band region of the ferri-state displays a nearly symmetric couplet at the B0-position that overlaps with a couplet 790 cm-1 higher in energy, which we assigned to a vibronic side-band transition. For the ferro-state, the couplet is greatly reduced, but still detectable. The B-band region is dominated by a positive Cotton effect at energies lower than B0 that is attributed to a magnetically allowed iron-->heme charge-transfer transition as earlier observed for nitrosyl myoglobin and hemoglobin. The Q-band region of the ferri-state is poorly resolved, but displays a pronounced positive signal at higher wavenumbers. This must result from a magnetically allowed transition, possibly from the methionine ligand to the dxy-hole of Fe3+. For the ferro-state, the spectra resolve the vibronic structure of the Qv-band. A more detailed spectral analysis reveals that the positively biased spectrum can be understood as a superposition of asymmetric couplets of split Q0 and Qv-states. Substantial qualitative and quantitative differences between the respective B-state and Q-state ECD spectra of yeast and horse heart cytochrome c can clearly be attributed to the reduced band splitting in the former, which results from a less heterogeneous internal electric field. Finally, we investigated the charge-transfer band at 695 nm in the ferri-state spectrum and found that it is composed of at least three bands, which are assignable to different taxonomic substates. The respective subbands differ somewhat with respect to their Kuhn anisotropy ratio and their intensity ratios are different for horse and yeast cytochrome c. Our data therefore suggests different substate populations for these proteins, which is most likely assignable to a structural heterogeneity of the distal Fe-M80 coordination of the heme chromophore.

Functionally Relevant Electric-Field Induced Perturbations of the Prosthetic Group of Yeast Ferrocytochrome c Mutants Obtained from a Vibronic Analysis of Low-Temperature Absorption Spectra

We have measured the low temperature (T = 20 K) absorption spectra of the N52A, N52V, N52I, Y67F, and N52AY67F mutants of ferrous Saccharomyces cerevisiae (baker's yeast) cytochrome c. All the bands in the Q0- and Q(v)-band region are split, and the intensity distributions among the split bands are highly asymmetric. The spectra were analyzed by a decomposition into Voigtian profiles. The spectral parameters thus obtained were further analyzed in terms of the vibronic coupling model of Schweitzer-Stenner and Bigman (Schweitzer-Stenner, R.; Bigman, D. J. Phys. Chem. B 2001, 7064-7073) to identify parameters related to electronic and vibronic perturbations of the heme macrocycle. We report that the electronic perturbation is of B(1g) symmetry and reflects the heterogeneity of the electric field at the heme, that is, the difference between the gradients along the perpendicular N-Fe-N axis of the heme core. We found that all the investigated mutations substantially increase this electronic perturbation, so that the spectral properties become similar to those of horse heart cytochrome c. Moreover, the electronic perturbation was found to correlate nonlinearly with the enthalpy changes associated with the reduction of the heme iron. Group theoretical arguments are invoked to propose a simple model which explains how a perturbation of the obtained symmetry can stabilize the reduced state of the heme iron. Finally, vibronic coupling parameters obtained from the analysis of the Q(v)-band region suggest that the investigated mutations decrease the nonplanar deformations of the heme group. This finding was reproduced by a normal mode structural decomposition (NSD) analysis of the N52V and N52VY67F heme conformations obtained from a 1 ns molecular dynamics simulation. We argue that the reduced nonplanarity contributes to the stabilization of the reduced state.

Stability of proteins: Temperature, pressure and the role of the solvent

We focus on the various aspects of the physics related to the stability of proteins. We review the pure thermodynamic aspects of the response of a protein to pressure and temperature variations and discuss the respective stability phase diagram. We relate the experimentally observed shape of this diagram to the low degree of correlation between the fluctuations of enthalpy and volume changes associated with the folding-denaturing transition and draw attention to the fact that one order parameter is not enough to characterize the transition. We discuss in detail microscopic aspects of the various contributions to the free energy gap of proteins and put emphasis on how a cosolvent may either enlarge or diminish this gap. We review briefly the various experimental approaches to measure changes in protein stability induced by cosolvents, denaturants, but also by pressure and temperature. Finally, we discuss in detail our own molecular dynamics simulations on cytochrome c and show what happens under high pressure, how glycerol influences structure and volume fluctuations, and how all this compares with experiments.

Temperature excursion infrared (TEIR) spectroscopy used to study hydrogen bonding between water and biomolecules

Water is a highly polar molecule that is capable of making four H-bonding linkages. Stability and specificity of folding of water-soluble protein macromolecules are determined by the interplay between water and functional groups of the protein. Yet, under some conditions, water can be replaced with sugar or other polar protic molecules with retention of protein structure. Infrared (IR) spectroscopy allows one to probe groups on the protein that interact with solvent, whether the solvent is water, sugar or glycerol. The basis of the measurement is that IR spectral lines of functional groups involved in H-bonding show characteristic spectral shifts with temperature excursion, reflecting the dipolar nature of the group and its ability to H-bond. For groups involved in H-bonding to water, the stretching mode absorption bands shift to lower frequency, whereas bending mode absorption bands shift to higher frequency as temperature decreases. The results indicate increasing H-bonding and decreasing entropy occurring as a function of temperature, even at cryogenic temperatures. The frequencies of the amide group modes are temperature dependent, showing that as temperature decreases, the amide group H-bonds to water strengthen. These results are relevant to protein stability as a function of temperature. The influence of solvent relaxation is demonstrated for tryptophan fluorescence over the same temperature range where the solvent was examined by infrared spectroscopy.

R-state hemoglobin bound to heterotropic effectors: models of the DPG, IHP and RSR13 binding sites

We performed a docking study followed by a 500-ps molecular dynamics simulation of R-state human adult hemoglobin (HbA) complexed to different heterotropic effectors [2,3-diphosphoglycerate (DPG), inositol hexaphosphate (IHP), and 2-[4-[(3,5-dichlorophenylcarbamoyl)-]methyl]-phenoxy]-2-methylpropionic acid (RSR13)) to propose a molecular basis for recently reported interactions of effectors with oxygenated hemoglobin. The simulations were carried out with counterions and explicit solvation. As reported for T-state HbA, the effector binding sites are also located in the central cavity of the R-state and differ depending on effector anionic character. DPG and IHP bind between the alpha-subunits and the RSR13 site spans the alpha1-, alpha2- and beta2-subunits. The generated models provide the first report of the molecular details of R-state HbA bound to heterotropic effectors.

Compressibility and uncoupling of cytochrome P450cam: high pressure FTIR and activity studies

The effect of the hydrostatic pressure on the CO ligand stretch vibration in cytochrome P450cam-CO bound with various substrates is studied by FTIR. The vibration frequency is linearily shifted to lower values with increasing pressure. The slope of the shift gives the isothermal compressibility of the heme pocket and is found to be related to the high-spin state content in an opposite direction to that previously observed from the pressure-induced shift of the Soret band. This opposite behaviour is explained by the dual effect of heme pocket water molecules both on the CO ligand and on electrostatic potentials produced by the protein at the distal side. The latter effect disturbs ligand-distal side contacts which are needed for a specific proton transfer in oxygen activation when dioxygen is the ligand. Their loss results in uncoupled H(2)O(2) formation.

The endogenous calcium ions of horseradish peroxidase C are required to maintain the functional nonplanarity of the heme

Horseradish peroxidase C (HRPC) binds 2 mol calcium per mol of enzyme with binding sites located distal and proximal to the heme group. The effect of calcium depletion on the conformation of the heme was investigated by combining polarized resonance Raman dispersion spectroscopy with normal coordinate structural decomposition analysis of the hemes extracted from models of Ca(2+)-bound and Ca(2+)-depleted HRPC generated and equilibrated using molecular dynamics simulations. Results show that calcium removal causes reorientation of heme pocket residues. We propose that these rearrangements significantly affect both the in-plane and out-of-plane deformations of the heme. Analysis of the experimental depolarization ratios are clearly consistent with increased B(1g)- and B(2g)-type distortions in the Ca(2+)-depleted species while the normal coordinate structural decomposition results are indicative of increased planarity for the heme of Ca(2+)-depleted HRPC and of significant changes in the relative contributions of three of the six lowest frequency deformations. Most noteworthy is the decrease of the strong saddling deformation that is typical of all peroxidases, and an increase in ruffling. Our results confirm previous work proposing that calcium is required to maintain the structural integrity of the heme in that we show that the preferred geometry for catalysis is lost upon calcium depletion.

Electric susceptibility of unsolvated glycine-based peptides

The DC electric susceptibilities of unsolvated glycine-based peptides, WGn (W = tryptophan and G = glycine) with n = 1-5, have been measured by deflection of a molecular beam in an electric field. These are the first electric deflection measurements performed on peptides. At 300 K the susceptibilities are in the range of 200-400 A(3). By far the largest contribution to the susceptibilities is from the permanent dipole moment of the peptides. The results indicate that the peptides do not have rigid conformations with fixed dipoles. Instead the dipole is averaged as the peptides explore their energy landscape. For a given WGn peptide, all molecules have almost the same average dipole, which suggests that they all explore a similar energy landscape on the microsecond time scale of the measurement. The measured susceptibilities are in good overall agreement with values calculated from the average dipole moment deduced from Monte Carlo simulations.

Cytochrome P-450-CO and substrates: lessons from ligand binding under high pressure

An overview of the application of high-pressure studies on the carbon monoxide complex of cytochrome P-450 is given. Different approaches to characterize ligand binding steps, the conformational states and substates and the compressibility of the ligand-bound complex are reviewed. A particular focus is the effect of substrates on these properties. It is shown that substrate mobility, compressibility and water accessibility are interrelated and may have functional meaning.

Aromatic substrate specificity of horseradish peroxidase C studied by a combined fluorescence line narrowing/energy minimization approach: the effect of localized side-chain reorganization

Horseradish peroxidase C binds a wide variety of small H-donor compounds such as benzohydroxamic acid (BHA) and 2-naphthohydroxamic acid (NHA). In this work, we use the Mg(II)-mesoporphyrin prosthetic group derivative as a spectroscopic probe of the active site and of the interaction with the substrates. We report on high-resolution fluorescence line-narrowed spectra which show that the effects of substrate binding on the electronic transitions are similar for both substrates and present data on the normal vibrational modes that are active in the vibronic spectra. Analysis of the vibrational frequencies shows that the Mg(II) ion is 5-coordinate in all cases, thus ruling out a solvent water as sixth ligand. The frequency shifts observed as a result of substrate binding are also indicative of a more rigid prosthetic group upon substrate binding. We present models for MgMP-HRP and its complexes with both substrates and compare the resulting structures on the basis of a modeling approach combining energy minimization to finite difference Poisson--Boltzmann calculations which partitions the various relative protein contributions to substrate binding. We show that the electrostatic potential of the prosthetic group is modified by the binding event. Analysis of the unbound and bound energy-minimized structures shows that the enzyme modulates substrate binding by subtle charge reorganization in the vicinity of the catalytic site and that this rearrangement is not attributable to significant secondary structure conformational changes but to side-chain reorganization.

Activation of signal-transduction mechanisms may underlie the therapeutic effects of an applied electric field

Successful treatment of various medical complaints with an applied electric field has been reported over the years. The identities of the cellular mechanisms that are influenced by this type of treatment and facilitate the positive effects, remain elusive. A study of many in vitro and in vivo reports revealed that the beneficial effects can be attributed to the activation of membrane proteins, and specifically proteins involved in signal-transduction mechanisms. Not only may the proteins be affected but it is now well established that enhanced Ca(2+)influx, observed to follow electric stimulation of cells, also contributes to many calcium-dependent cellular processes which can be linked to the therapeutic effects discussed in this paper. An hypothesis of the physical changes caused by an applied, relatively small (10(3)to 10(4)V m(-1)range), electric field with low to moderate frequency (below 150 Hz), is postulated.

Engineering the pH-optimum of a triglyceride lipase: from predictions based on electrostatic computations to experimental results

The optimisation of enzymes for particular purposes or conditions remains an important target in virtually all protein engineering endeavours. Here, we present a successful strategy for altering the pH-optimum of the triglyceride lipase cutinase from Fusarium solani pisi. The computed electrostatic pH-dependent potentials in the active site environment are correlated with the experimentally observed enzymatic activities. At pH-optimum a distinct negative potential is present in all the lipases and esterases that we studied so far. This has prompted us to propose the "The Electrostatic Catapult Model" as a model for product release after cleavage of the ester bond. The origin of the negative potential is associated with the titration status of specific residues in the vicinity of the active site cleft. In the case of cutinase, the role of Glu44 was systematically investigated by mutations into Ala and Lys. Also, the neighbouring Thr45 was mutated into Proline, with the aim of shifting the spatial location of Glu44. All the charge mutants displayed altered titration behaviour of active site electrostatic potentials. Typically, the substitution of the residue Glu44 pushes the onset of the active site negative potential towards more alkaline conditions. We, therefore, predicted more alkaline pH optima, and this was indeed the experimentally observed. Finally, it was found that the pH-dependent computed Coulombic energy displayed a strong correlation with the observed melting temperatures of native cutinase.

Macromolecular electrostatics: continuum models and their growing pains

Theoretical understanding of macromolecular electrostatics has advanced substantially over the past year. Continuum models have given promising results for calculating protein-ligand binding free energy differences, as well as pK(a)s and redox properties, particularly with explicit treatment of multiple conformers. Generalized Born and other techniques have led to the first molecular dynamics simulations of proteins and RNA with continuum solvent. Continuum and microscopic descriptions of dielectric relaxation have been critically compared.

A study on the mechanism of the proton transport in bacteriorhodopsin: the importance of the water molecule

The mechanism of proton transport around the Schiff base in bacteriorhodopsin was investigated by ab initio molecular orbital (MO) calculations. Computations were performed for the case where there is a water molecule between the Schiff base and the Asp residue and for the case where there is no water molecule. Changes in the atomic configuration and potential energy through the proton transport process were compared between two cases. In the absence of water, the protonated Schiff base was not stable, and a proton was spontaneously detached from the Schiff base. On the other hand, a stable structure of the protonated Schiff base was obtained in the presence of water. This suggests that the presence of a water molecule is required for stability in the formation of a protonated Schiff base.

Sampling field heterogeneity at the heme of c-type cytochromes by spectral hole burning spectroscopy and electrostatic calculations

We report on a comparative investigation of the heme pocket fields of two Zn-substituted c-type cytochromes-namely yeast and horse heart cytochromes c-using a combination of hole burning Stark spectroscopy and electrostatic calculations. The spectral hole burning experiments are consistent with different pocket fields experienced at the hemes of the respective cytochromes. In the case of horse heart Zn-cytochrome c, two distinguishable electronic origins with different electrostatic properties are observed. The yeast species, on the other hand, displays a single electronic origin. Electrostatic calculations and graphics modeling using the linearized finite-difference Poisson-Boltzmann equation performed at selected time intervals on nanosecond-molecular dynamics trajectories show that the hemes of the respective cytochromes sample different potentials as they explore conformational space. The electrostatic potentials generated by the protein matrix at the heme show different patterns in both cytochromes, and we suggest that the cytochromes differ by the number of "electrostatic substates" that they can sample, thus accounting for the different spectral populations observed in the two cytochromes.

Protonation of porphyrin in iron-free cytochrome c: spectral properties of monocation free base porphyrin, a charge analogue of ferric heme

Charged groups reside mainly on protein surfaces, but for proteins that incorporate redox centers, a charge typically exists at the prosthetic group within the interior. How a protein accommodates a buried charge and the effect of redox changes on protein stability are thermodynamically related problems. To examine these problems in cytochrome c, the metal-free protein was used as a model. When pH is lowered, the neutral, monocation, and dication forms of the porphyrin are progressively formed as indicated by their characteristic absorption spectra. Infrared studies of the protein over this pH range show that the protein remains in a predominately alpha-helical structure, although the carboxyl groups of the dicarboxylic amino acids become protonated at lower pH. The monocation porphyrin form (which has not been previously reported in a protein and is a charge analogue of ferric heme) has a fluorescence maximum at 609 nm. The pKs for the respective one and two protonation of the porphyrin pyrrole Ns are 3.2 and 1.6 for the folded protein, and 4.4 and 3.1 for the unfolded protein. These values indicate that the protection of the polypeptide chain for protonation is approximately 3 kcal.

The individual tyrosines of proteins: their spectra may or may not differ from those in water or other solvents

The overall derivative spectrum of a protein is the sum of the individual derivative spectra just as the overall ultraviolet spectrum of a protein is the sum of its component parts. The RNase and DNA binding protein Sso7d has two tyrosines and one tryptophan. We used two mutant forms of the protein to show that the individual aromatics contribute derivative spectra that can be explained on the basis of their environments. We used mutant forms of iso-1-cytochrome c to estimate the contributions of the single tryptophan and three of the five tyrosines to the overall derivative spectrum. The tryptophan spectrum is not exceptional. The comparable tyrosine spectra are more complex. The derivative spectrum of individual tyrosines does not correspond to that expected on the basis of concentration. This is a reflection of two factors: (1) the extent to which mutations are sensed distally through the introduction and compression of packing defects; and (2) the extent to which electronic transitions of tyrosine are influenced by nearby atoms. This influence could take the form of tyrosine residing in an area where the dielectric coefficient is not uniform; it could also result from tyrosine bumping into neighboring atoms with lower frequency than it does in solution.

A mutational analysis of the acetylcholine receptor channel transmitter binding site

Mutagenesis and single-channel kinetic analysis were used to investigate the roles of four acetylcholine receptor channel (AChR) residues that are candidates for interacting directly with the agonist. The EC50 of the ACh dose-response curve was increased following alpha-subunit mutations Y93F and Y198F and epsilon-subunit mutations D175N and E184Q. Single-channel kinetic modeling indicates that the increase was caused mainly by a reduced gating equilibrium constant (Theta) in alphaY198F and epsilonD175N, by an increase in the equilibrium dissociation constant for ACh (KD) and a reduction in Theta in alphaY93F, and only by a reduction in KD in epsilonE184Q. This mutation altered the affinity of only one of the two binding sites and was the only mutation that reduced competition by extracellular K+. Additional mutations of epsilonE184 showed that K+ competition was unaltered in epsilonE184D and was virtually eliminated in epsilonE184K, but that neither of these mutations altered the intrinsic affinity for ACh. Thus there is an apparent electrostatic interaction between the epsilonE184 side chain and K+ ( approximately 1.7kBT), but not ACh+. The results are discussed in terms of multisite and induced-fit models of ligand binding to the AChR.

Fluorescence line narrowing applied to the study of proteins

Fluorescence line narrowing is a high resolution spectroscopic technique that uses low temperature and laser excitation to optically select specific subpopulations from the inhomogeneously broadened absorption band of the sample. When applied to the study of fluorescent groups in proteins one can obtain vibronically resolved spectra, which can be analyzed to give information on spectral line shapes, vibrational energies of both the ground and excited state molecule, and the inhomogeneous distribution function of the electronic transitions. These parameters reveal information about the chromophoric prosthetic group and the protein matrix and are functions of geometric strains and local electric fields imposed by the protein. Examples of the use of fluorescence line narrowing are discussed in investigations of heme proteins, photosynthetic systems and tryptophan-containing proteins.

The protein state of matter

Three ways are generally used to visualize proteins: (1) a static model in which the atomic positions are defined, (2) a dynamic model taking into account fluctuations, and (3) a reactive model that reflects the internal and external electric fields of the molecule. The properties of chromophoric prosthetic groups can be probed by optical spectroscopy, and when high resolution techniques are used, the results reveal information about the local electric fields in proteins, as influenced and determined by atomic positions and dynamics.