Dielectric relaxation and solvation dynamics of water in complex chemical and biological systems

Broadband 10-300 GHz stimulus-response sensing for chemical and biological entities

By illuminating the sample with a broadband 10-300 GHz stimulus and coherently detecting the response, we obtain reflection and transmission spectra of common powdered substances, and compare them as a starting point for distinguishing concealed threats in envelopes and on personnel. Because these samples are irregular and their dielectric properties cannot be modulated, however, the spectral information we obtain is largely qualitative. To show how to gain quantitative information on biological species at micro- and millimetre-wave frequencies, we introduce thermal modulation of a globular protein in solution, and show that changes in single-wavelength microwave reflections coincide with accepted visible absorption spectra, pointing the way towards gaining quantitative chemical and biological spectra from broadband terahertz systems.

Cys(x)His(y)-Zn2+ interactions: thiol vs. thiolate coordination

In zinc proteins, the Zn2+ cation frequently binds with a tetrahedral coordination to cysteine and histidine side chains, for example, in many DNA-binding proteins, where it plays primarily a structural role. We examine the possibility of thiolate protonation in Cys(x)His(y)-Zn2+ groups, both in proteins and in solution, through a combination of theoretical calculations and database analysis. Seventy-five percent of the thiolate-coordinated zincs in the Cambridge Structural Database are tetrahedral, while di-alkanethiol coordination always involves five or more ligands. Ab initio quantum calculations are performed on (ethanethiol/thiolate)(3)imidazole-Zn2+ complexes in vacuum, yielding geometries and gas phase basicities. Protonating one (respectively two) thiolates increases the Zn-S(thiol) distance by 0.4 A (respectively 0.3 A), providing a structural marker for protonation. The stabilities of the complexes in solution are compared by combining the gas phase basicities with continuum dielectric solvation calculations. In a continuum solvent with permittivity epsilon = 4, 20, or 80, one of three thiolates is predicted to be protonated at neutral pH. By extension, Cys4-Zn2+ groups are expected to be protonated in the same conditions. In contrast, most Cys3His and Cys4 geometries in the Protein Data Bank (PDB) appear consistent with all-thiolate Zn2+ coordination. This apparent discrepancy is resolved by two recent surveys of zinc protein structures, which suggest that these all-thiolate sites are stabilized by charged and polar groups nearby in the protein, thus overcoming their intrinsic instability. However, the experimental resolution is not sufficient in all the PDB structures to rule out a thiol/thiolate mixture, and protonated thiolates may occur in some proteins not solved at high resolution or not represented in the PDB, as suggested by recent mass spectrometry experiments; this possibility should be allowed for in X-ray structure refinement.

Hydrogen-bond dynamics near a micellar surface: origin of the universal slow relaxation at complex aqueous interfaces

The dynamics of hydrogen bonds among water molecules themselves and with the polar head groups (PHG) at a micellar surface have been investigated by long molecular dynamics simulations. The lifetime of the hydrogen bond between a PHG and a water molecule is found to be much longer than that between any two water molecules, and is likely to be a general feature of hydrophilic surfaces of organized assemblies. Analyses of individual water trajectories suggest that water molecules can remain bound to the micellar surface for more than 100 ps. The activation energy for such a transition from the bound to a free state for the water molecules is estimated to be about 3.5 kcal/mol.

On the competition of the purine bases, functionalities of peptide side chains, and protecting agents for the coordination sites of dicationic cisplatin derivatives

The Pt-L bond energies of simple triammineplatinum(II) complexes, [Pt(NH(3))(3)L](2+), with oxygen-, nitrogen-, and sulfur-containing donor ligands L have been predicted and rationalized using density functional theory. The ligands L have been chosen as models for functionalities of peptide side chains, for sulfur-containing protecting agents, and for adenine and guanine sites of the DNA as the ultimate target of platinum anticancer drugs. Calculation of the Pt-L bond energy in [Pt(NH(3))(3)L](2+) reveals that the soft metal center of triammineplatinum(II) prefers N ligands over S ligands. This remarkable result has been discussed in light of several interpretations of the hard and soft acids and bases principle. The concept of orbital-symmetry-based energy decomposition has been employed for the determination of the contributions from sigma and pi orbital interactions, electrostatics, and intramolecular hydrogen bonding to the Pt-L bond energy. The calculations show that considerable differences in the bond energies of the triammineplatinum(II) complexes with N-heterocycles such as 1-methylimidazole, 9-methyladenine, and 9-methylguanine arise from electrostatics rather than from orbital interactions. Surprisingly, the net stabilization by hydrogen bonding between the (Pt)N-H group and the oxygen of 9-methylguanine is as weak as the intramolecular hydrogen bond in the aqua complex [Pt(NH(3))(3)(H(2)O)](2+), challenging the common hypothesis that DNA-active anticancer drugs require carrier ligands with N-H functionalities because of their hydrogen-bonding ability. The influence of a polarizable environment on the stability of the complexes has been investigated systematically with the dependence of the dielectric constant epsilon. With increasing epsilon, the complexes with S-containing ligands are more strongly stabilized than the complexes of the N-containing heterocycles. At epsilon = 78.4, the dielectric constant of water, 9-methylguanine remains the only purine derivative investigated which is competitive to neutral sulfur ligands. These findings are particularly important for a rationalization of the results from recent experimental studies on the competition of biological donor ligands L for coordination with the metal center of [Pt(dien)L](2+) (dien = 1,5-diamino 3-azapentane).

Gaussian fluctuations and linear response in an electron transfer protein

In response to charge separation or transfer, polar liquids respond in a simple linear fashion. A similar linear response for proteins might be expected from the central limit theorem and is postulated in widely used theories of protein electrostatics, including the Marcus electron transfer theory and dielectric continuum theories. Although these theories are supported by a variety of experimental data, the exact validity of a linear protein dielectric response has been difficult to determine. Molecular dynamics simulations are presented that establish a linear dielectric response of both protein and surrounding solvent over the course of a biologically relevant electron transfer reaction: oxido-reduction of yeast cytochrome c in solution. Using an umbrella-sampling free energy approach with long simulations, an accurate treatment of long-range electrostatics and both classical and quantum models of the heme, good agreement is obtained with experiment for the redox potential relative to a heme-octapeptide complex. We obtain a reorganization free energy that is only half that for heme-octapeptide and is reproduced with a dielectric continuum model where the heme vicinity has a dielectric constant of only 1.1. This value implies that the contribution of protein reorganization to the electron transfer free energy barrier is reduced almost to the theoretical limit (a dielectric of one), and that the fluctuations of the electrostatic potential on the heme have a simple harmonic form, in accord with Marcus theory, even though the fluctuations of many individual protein groups (especially at the protein surface) are anharmonic.

Femtosecond dynamics of rubredoxin: tryptophan solvation and resonance energy transfer in the protein

We report here studies of tryptophan (Trp) solvation dynamics in water and in the Pyrococcus furiosus rubredoxin protein, including the native and its apo and denatured forms. We also report results on energy transfer from Trp to the iron-sulfur [Fe-S] cluster. Trp fluorescence decay with the onset of solvation dynamics of the chromophore in water was observed with femtosecond resolution ( approximately 160 fs; 65% component), but the emission extended to the picosecond range (1.1 ps; 35% component). In contrast, the decay is much slower in the native rubredoxin; the Trp fluorescence decay extends to 10 ps and longer, reflecting the local rigidity imposed by residues and by the surface water layer. The dynamics of resonance energy transfer from the two Trps to the [Fe-S] cluster in the protein was observed to follow a temporal behavior characterized by a single exponential (15-20 ps) decay. This unusual observation in a protein indicates that the resonance transfer is to an acceptor of a well-defined orientation and separation. From studies of the mutant protein, we show that the two Trp residues have similar energy-transfer rates. The critical distance for transfer (R(0)) was determined, by using the known x-ray data, to be 19.5 A for Trp-36 and 25.2 A for Trp-3, respectively. The orientation factor (kappa(2)) was deduced to be 0.13 for Trp-36, clearly indicating that molecular orientation of chromophores in the protein cannot be isotropic with kappa(2) value of 2/3. These studies of solvation and energy-transfer dynamics, and of the rotational anisotropy, of the wild-type protein, the (W3Y, I23V, L32I) mutant, and the fmetPfRd variant at various pH values reveal a dynamically rigid protein structure, which is probably related to the hyperthermophilicity of the protein.

Ultrafast dynamics in aqueous polyacrylamide solutions

We have investigated the ultrafast dynamics of aqueous polyacrylamide ([-CH(2)CH(CONH(2))-](n), or PAAm) solutions using femtosecond optical heterodyne-detected Raman-induced Kerr effect spectroscopy (OHD-RIKES). The observed aqueous PAAm dynamics are nearly identical for both M(w) = 1500 and 10 000. Aqueous propionamide (CH(3)CH(2)CONH(2), or PrAm) solutions were also studied, because PrAm is an exact model for the PAAm constitutional repeat unit (CRU). The longest time scale dynamics observed for both aqueous PAAm and PrAm solutions occur in the 4-10 ps range. Over the range of concentrations from 0 to 40 wt %, the picosecond reorientation time constants for the aqueous PAAm and PrAm solutions scale linearly with the solution concentration, despite the fact that the solution shear viscosities vary exponentially from 1 to 264 cP. For a given value of solution concentration in weight percent, constant ratios of measured reorientation time constants for PAAm to PrAm are obtained. This ratio of PAAm to PrAm reorientation time constants is equal to the ratio of the volume for the PAAm constitutional repeat unit (-CH(2)CHCONH(2)-) to the molecular volume of PrAm. For these reasons, we assign the polymer reorientation dynamics to motions of the entire constitutional repeat unit, not only side group motions. Simple molecular dynamics simulations of H[-CH(2)CH(CONH(2))-](7)H in a periodic box with 180 water molecules support this assignment. Amide-amide and amide-water hydrogen-bonding interactions lead to strongly oscillatory femtosecond dynamics in the Kerr transients, peaking at 80, 410, and 750 fs.

Dielectric relaxation in an enzyme active site: molecular dynamics simulations interpreted with a macroscopic continuum model

Dielectric relaxation plays an important role in many chemical processes in proteins, including acid-base titration, ligand binding, and charge transfer reactions. Its complexity makes experimental characterization difficult, and so, theoretical approaches are valuable. The comparison of molecular dynamics free energy simulations with simpler models such as a dielectric continuum model is especially useful for obtaining qualitative insights. We have analyzed a charge insertion process that models deprotonation or mutation of an important side chain in the active site of the enzyme aspartyl-tRNA synthetase. Complexes with the substrate aspartate and the analogue asparagine were studied. The resulting dielectric relaxation was found to involve both ligand and side chain rearrangements in the active site and to account for a large part of the overall charging free energy. With the continuum model, charge insertion is performed along a two-step pathway: insertion into a static environment, followed by relaxation of the environment. These correspond to different physical processes and require different protein dielectric constants. A low value of approximately 1 is needed for the static step, consistent with the parametrization of the molecular mechanics charge set used. A value of 3-6 (depending on the exact insertion site and the nature of the ligand) is needed to describe the dielectric relaxation step. This moderate value indicates that, for this system, the local protein polarizability in the active site is within at most a factor of 2 of that expected at nonspecific positions in a protein interior.

Femtosecond dynamics of flavoproteins: charge separation and recombination in riboflavine (vitamin B2)-binding protein and in glucose oxidase enzyme

Flavoproteins can function as hydrophobic sites for vitamin B(2) (riboflavin) or, in other structures, with cofactors for catalytic reactions such as glucose oxidation. In this contribution, we report direct observation of charge separation and recombination in two flavoproteins: riboflavin-binding protein and glucose oxidase. With femtosecond resolution, we observed the ultrafast electron transfer from tryptophan(s) to riboflavin in the riboflavin-binding protein, with two reaction times: approximately 100 fs (86% component) and 700 fs (14%). The charge recombination was observed to take place in 8 ps, as probed by the decay of the charge-separated state and the recovery of the ground state. The time scale for charge separation and recombination indicates the local structural tightness for the dynamics to occur that fast and with efficiency of more than 99%. In contrast, in glucose oxidase, electron transfer between flavin-adenine-dinucleotide and tryptophan(s)/tyrosine(s) takes much longer times, 1.8 ps (75%) and 10 ps (25%); the corresponding charge recombination occurs on two time scales, 30 ps and nanoseconds, and the efficiency is still more than 97%. The contrast in time scales for the two structurally different proteins (of the same family) correlates with the distinction in function: hydrophobic recognition of the vitamin in the former requires a tightly bound structure (ultrafast dynamics), and oxidation-reduction reactions in the latter prefer the formation of a charge-separated state that lives long enough for chemistry to occur efficiently. Finally, we also studied the influence on the dynamics of protein conformations at different ionic strengths and denaturant concentrations and observed the sharp collapse of the hydrophobic cleft and, in contrast, the gradual change of glucose oxidase.