Dynamics and thermodynamics of water around EcoRI bound to a minimally mutated DNA chain

Water plays an important role in protein-DNA interactions. Here, we examine using molecular dynamics simulations the differences in the dynamic and thermodynamic properties of water in the interfacial and intercalating regions of EcoRI bound to the cognate and to a minimally mutated noncognate DNA chain. The results show that the noncognate complex is not only more hydrated than the cognate complex, but the interfacial waters in the noncognate complex exhibit a faster dynamics, which in turn reduces the hydrogen-bond lifetimes. Thus, the higher hydration, faster reorientation dynamics and faster hydrogen-bond-relaxation times of water, taken together, indicate that, even with a minimal mutation of the DNA sequence, the interfacial regions of the noncognate complex are more poised to allowing the protein to diffuse away than to promoting the formation of a stable complex. Alternatively, the results imply that the slowed water dynamics in the interfacial regions when the protein chances upon a cognate sequence allow the formation of a stable specific protein-DNA complex leading to catalytic action.     

Dynamic properties of water around a protein-DNA complex from molecular dynamics simulations

Formation of protein-DNA complex is an important step in regulation of genes in living organisms. One important issue in this problem is the role played by water in mediating the protein-DNA interactions. In this work, we have carried out atomistic molecular dynamics simulations to explore the heterogeneous dynamics of water molecules present in different regions around a complex formed between the DNA binding domain of human TRF1 protein and a telomeric DNA. It is demonstrated that such heterogeneous water motions around the complex are correlated with the relaxation time scales of hydrogen bonds formed by those water molecules with the protein and DNA. The calculations reveal the existence of a fraction of extraordinarily restricted water molecules forming a highly rigid thin layer in between the binding motifs of the protein and DNA. It is further proved that higher rigidity of water layers around the complex originates from more frequent reformations of broken water-water hydrogen bonds. Importantly, it is found that the formation of the complex affects the transverse and longitudinal degrees of freedom of surrounding water molecules in a nonuniform manner.     

Nonspecific protein-DNA interactions: complexation of alpha-chymotrypsin with a genomic DNA

In this contribution, we report studies on nonspecific protein-DNA interactions of an enzyme protein bovine pancreatic alpha-chymotrypsin (CHT) with genomic DNA (from salmon testes) using two biologically common fluorescent probes: 1-anilinonaphthalene-8-sulfonate (ANS) and 2,6-p-toluidinonaphthalene sulfonate (TNS). TNS molecules that are nonspecifically bound to positively charged basic residues at the surface sites, not in the hydrophobic cavities of the protein, are preferentially displaced upon complexation of TNS-labeled CHT with DNA. The time-resolved fluorescence anisotropy of TNS molecules bound to hydrophobic cavities/clefts of CHT reveals that global tumbling motion of the protein is almost frozen in the protein-DNA complex. A control study on TNS-labeled human serum albumin (HSA) upon interaction with DNA clearly indicates that the ligands in the deep pockets of the protein cannot be displaced by interaction with DNA. We have also found that ANS, which binds to a specific surface site of CHT, is not displaced by DNA. The intactness of the ANS binding in CHT upon complexation with DNA offers the opportunity to measure the distance between the ANS binding site and the contact point of the ethidium bromide (EB)-labeled DNA using the F?rster resonance energy transfer (FRET) technique. Enzymatic activity studies on CHT on a substrate (Ala-Ala-Phe 7-amido-4-methyl coumarin) reveal that the active site of the enzyme remains open for the substrate even in the protein-DNA complex. Circular dichroism (CD) studies on CHT upon complexation with DNA confirm the structural integrity of CHT in the complex. Our studies have attempted to explore an application of nonspecific protein-DNA interactions in the characterization of ligand binding of a protein in solution.     

Direct observation of enhanced translocation of a homeodomain between DNA cognate sites by NMR exchange spectroscopy

A novel approach is presented for studying the kinetics of specific protein-DNA interactions by NMR exchange spectroscopy. The experimental design involves the direct observation of translocation of a homeodomain between cognate sites on two oligonucleotide duplexes, differing by only a single base pair at the edge of the DNA recognition sequence. The single base-pair change perturbs the 1H-15N correlation spectrum of a number of residues, while leaving the affinity for the DNA unchanged. The exchange process has apparent rate constants in the 5-20 s-1 range which are linearly dependent upon the concentration of free DNA. These rates are about 3 orders of magnitude larger than the dissociation rate constant determined by gel shift assays at nanomolar DNA concentrations. The complete NMR exchange data set, comprising auto- and cross-peak intensities as a function of mixing time at five concentrations of free DNA, can be fit simultaneously to a simple model in which protein translocation between DNA duplexes occurs via a second-order process (with rate constants of approximately 6 x 104 M-1 s-1) involving direct collision of a protein-DNA complex with free DNA. This is akin to intersegmental transfer, and a physical model for the process is discussed. Rapid translocation at high concentrations of free DNA observed directly by NMR exchange spectroscopy reconciles the long half-lives of protein-DNA complexes measured by biochemical analysis in vitro with the highly dynamic behavior of such complexes observed in vivo. The relevance of this mechanism to the kinetics of protein-DNA interactions within the cell is discussed.     

Quantitative distance information on protein-DNA complexes determined in polyacrylamide gels by fluorescence resonance energy transfer

In polyacrylamide gels, we have quantitatively determined Forster transfer (fluorescense resonance energy transfer, FRET) between two fluorescent dyes attached to DNA in protein-DNA complexes. The donor-dye fluorescein labeled to DNA retains its free mobility in the polyacrylamide gel, however, its fluorescence properties change. The different quantum yield of fluorescein in the gel is found to be independent of the gel concentration and can thus be quantitatively taken into account by a reduced Forster distance R0 of 46 A compared to 50 A in solution. We have determined global structural properties of two proteins binding to double-labeled DNA using a novel gel-based fluorescence resonance energy transfer assay. In polyacrylamide gels we have analyzed the binding of integration host factor (IHF) and the high mobility group protein NHP6a to their substrate DNA. The measured Forster transfer efficiency allows us to deduce information on the overall shape and the DNA bending angle in the complex.     

CF3 Rotation in 3-(Trifluoromethyl)phenanthrene: Solid State 19F and 1H NMR relaxation and Bloch-Wangsness-Redfield theory

We have observed and modeled the 1H and 19F solid-state nuclear spin relaxation process in polycrystalline 3-(trifluoromethyl)phenanthrene. The relaxation rates for the two spin species were observed from 85 to 300 K at the low NMR frequencies of omega/2pi = 22.5 and 53.0 MHz where CF3 rotation, characterized by a mean time tau between hops, is the only motion on the NMR time scale. All motional time scales (omegatau < 1, omegatau approximately 1, and omegatau > 1) are observed. The 1H spins are immobile on the NMR time scale but are coupled to the 19F spins via the unlike-spin dipole-dipole interaction. The temperature dependence of the observed relaxation rates (the relaxation is biexponential) shows considerable structure and a thorough analysis of Bloch-Wangsness-Redfield theory for this coupled spin system is provided. The activation energy for CF3 rotation is 11.5 +/- 0.7 kJ/mol, in excellent agreement with the calculation in a 13-molecule cluster provided in the companion paper where the crystal structure is reported and detailed ab initio electronic structure calculations are performed [Wang, X.; Mallory F. B.; Mallory, C. W; Beckmann, P. A.; Rheingold, A. L.; Francl, M. M J. Phys. Chem. A 2006, 110, 3954].     

Structure of the charge separated state P865(+)Q(A)- in the photosynthetic reaction centers of Rhodobacter sphaeroides by quantum beat oscillations and high-field electron paramagnetic resonance: evidence for light-induced Q(A)- reorientation

The structure of the secondary radical pair, P865(+)Q(A)-, in fully deuterated and Zn-substituted reaction centers (RCs) of the purple bacterium Rhodobacter sphaeroides R-26 has been determined by high-time resolution and high-field electron paramagnetic resonance (EPR). A computer analysis of quantum beat oscillations, observed in a two-dimensional Q-band (34 GHz) EPR experiment, provides the orientation of the various magnetic tensors of P865(+)Q(A)- with respect to a magnetic reference frame. The orientation of the g-tensor of P865(+) in an external reference system is adapted from a single-crystal W-band (95 GHz) EPR study [Klette, R.; T?rring, J. T.; Plato, M.; M?bius, K.; B?nigk, B.; Lubitz, W. J. Phys. Chem. 1993, 97, 2015-2020]. Thus, we obtain the three-dimensional structure of the charge separated state P865(+)Q(A)- on a nanosecond time scale after light-induced charge separation. Comparison with crystallographic data reveals that the position of the quinone is essentially the same as that in the X-ray structure. However, the head group of Q(A)- has undergone a 60 degrees rotation in the ring plane relative to its orientation in the crystal structure. Analysis suggests that the two different QA conformations are functionally relevant states which control the electron-transfer kinetics from Q(A)- to the secondary quinone acceptor QB. It appears that the rate-limiting step of this reaction is a reorientation of Q(A)- in its binding pocket upon light-induced reduction. The new kinetic model accounts for striking observations by Kleinfeld et al. who reported that electron transfer from Q(A)- to QB proceeds in RCs cooled to cryogenic temperature under illumination but does not proceed in RCs cooled in the dark [Kleinfeld, D.; Okamura, M. Y.; Feher, G. Biochemistry 1984, 23, 5780-5786].     

Advantages and drawbacks of nanospray for studying noncovalent protein-DNA complexes by mass spectrometry

The noncovalent complexes between the BlaI protein dimer (wild-type and GM2 mutant) and its double-stranded DNA operator were studied by nanospray mass spectrometry and tandem mass spectrometry (MS/MS). Reproducibility problems in the nanospray single-stage mass spectra are emphasized. The relative intensities depend greatly on the shape of the capillary tip and on the capillary-cone distance. This results in difficulties in assessing the relative stabilities of the complexes simply from MS(1) spectra of protein-DNA mixtures. Competition experiments using MS/MS are a better approach to determine relative binding affinities. A competition between histidine-tagged BlaIWT (BlaIWTHis) and the GM2 mutant revealed that the two proteins have similar affinities for the DNA operator, and that they co-dimerize to form heterocomplexes. The low sample consumption of nanospray allows MS/MS spectra to be recorded at different collision energies for different charge states with 1 microL of sample. The MS/MS experiments on the dimers reveal that the GM2 dimer is more kinetically stable in the gas phase than the wild-type dimer. The MS/MS experiments on the complexes shows that the two proteins require the same collision energy to dissociate from the complex. This indicates that the rate-limiting step in the monomer loss from the protein-DNA complex arises from the breaking of the protein-DNA interface rather than the protein-protein interface. The dissociation of the protein-DNA complex proceeds by the loss of a highly charged monomer (carrying about two-thirds of the total charge and one-third of the total mass). MS/MS experiments on a heterocomplex also show that the two proteins BlaIWTHis and BlaIGM2 have slightly different charge distributions in the fragments. This emphasizes the need for better understanding the dissociation mechanisms of biomolecular complexes.     

Slow manifold for a bimolecular association mechanism

Finding the slow manifold for two-variable ordinary differential equation (ODE) models of chemical reactions with a single equilibrium is generally simple. In such planar ODEs the slow manifold is the unique trajectory corresponding to the slow relaxation of the system as it moves towards the equilibrium point. One method of finding the slow manifold is to use direct iteration of a functional equation; another method is to obtain a series solution of the trajectory differential equation of the system. In some cases these two methods agree order-by-order in the singular perturbation parameter controlling the fast relaxation of the intermediate (complex). However, de la Llave has found a model ODE where the series method always diverges. Bimolecular association is an example of a chemical reaction where the series method for finding the slow manifold diverges but the iterative method converges. In this mechanism a complex is formed which can then undergo unimolecular decay, i.e., [reaction: see text]. The kinetics of this reaction are investigated and its properties compared with two other two-step mechanisms where series expansion and iteration methods are equivalent: the Michaelis-Menten mechanism for enzyme kinetics, and the Lindemann-Christiansen mechanism of unimolecular decay in gas kinetics.     

An 15N NMR spin relaxation dispersion study of the folding of a pair of engineered mutants of apocytochrome b562

15N relaxation dispersion NMR spectroscopy has been used to study exchange dynamics in a pair of mutants of Rd-apocyt b562, a redesigned four-helix-bundle protein. An analysis of the relaxation data over a range of temperatures establishes that exchange in both proteins is best modeled as two-state and that it derives from the folding/unfolding transition. These results are in accord with predictions based on the reaction coordinate for the folding of the protein determined from native-state hydrogen exchange data [Chu, R.; Pei, W.; Takei, J.; Bai, Y. Biochemistry 2002, 41, 7998-8003]. The kinetics and thermodynamics of the folding transition have been characterized in detail. Although only a narrow range of temperatures could be examined, it is clear that the folding rate temperature profile is distinctly non-Arrhenius for both mutants, with the folding barrier for at least one of them entropic.     

Catalysis by methyltrioxorhenium(VII): reduction of hydronium ions by europium(II) and reduction of perchlorate ions by europium(II) and chromium(II)

The title reactions occur stepwise, the first and fastest being MeReO3 + Eu2+ --> Re(VI) + Eu3+ (k298 = 2.7 x 10(4) L mol(-1) s(-1)), followed by rapid reduction of Re(VI) by Eu2+ to MeReO2. The latter species is reduced by a third Eu2+ to Re(IV), a metastable species characterized by an intense charge transfer band, epsilon410 = 910 L mol(-1) cm(-1) at pH 1; the rate constant for its formation is 61.3 L mol(-1) s(-1), independent of [H+]. Yet another reduction step occurs, during which hydrogen is evolved at a rate v = k[Re(IV)][Eu2+][H+](-1), with k = 2.56 s(-1) at mu = 0.33 mol L(-1). The 410 nm Re(IV) species bears no ionic charge on the basis of the kinetic salt effect. We attribute hydrogen evolution to a reaction between H-ReVO and H3O+, where the hydrido complex arises from the unimolecular rearrangement of Re(III)-OH in a reaction that cannot be detected directly. Chromium(II) ions do not evolve H2, despite E(Cr) degrees approximately E(EU) degrees. We attribute this lack of reactivity to the Re(IV) intermediate being captured as [Re(IV)-O-Cr(III)]2+, with both metals having substitutionally inert d3 electronic configurations. Hydrogen evolution occurs in chloride or triflate media; with perchlorate present, MeReO2 reduces perchlorate to chloride, as reported previously [Abu-Omar, M. M.; Espenson, J. H. Inorg. Chem. 1995, 34, 6239-6240].     

Ruthenium-catalyzed intermolecular coupling reactions of arylamines with ethylene and 1,3-dienes: mechanistic insight on hydroamination vs ortho-C-H bond activation

[reaction: see text]. The cationic ruthenium complex [(PCy3)2(CO)(Cl)Ru=CHCH=C(CH3)2]+BF4- was found to be an effective catalyst for the coupling reaction of aniline and ethylene to form a approximately 1:1 ratio of N-ethylaniline and 2-methylquinoline products. The analogous reaction with 1,3-dienes resulted in the preferential formation of Markovnikov addition products. The normal isotope effect of k(NH)/k(ND) = 2.2 (aniline and aniline-d7 at 80 degrees C) and the Hammett rho = -0.43 (correlation of para-substituted p-X-C6H4NH2) suggest an N-H bond activation rate-limiting step for the catalytic reaction.     

Polarized transient absorption to resolve electron transfer between tryptophans in DNA photolyase

Transient absorption spectroscopy is a powerful tool for studying biological electron-transfer chains, provided that their members give rise to distinct changes of their absorption spectra. There are, however, chains that contain identical molecules, so that electron transfer between them does not change net absorption. An example is the chain flavin adenine dinucleotide (FAD)-W382-W359-W306 in DNA photolyase from E. coli. Upon absorption of a photon, the excited state of FADH* (neutral FAD radical) abstracts an electron from the tryptophan residue W382 in approximately 30 ps (monitored by transient absorption). The cation radical W382*+ is presumably reduced by W359 and W359*+ by W306. The latter two reactions could not be monitored directly so far because the absorption changes of the partners compensate in each step. To overcome this difficulty, we used linearly polarized flashes for excitation of FADH*, thus inducing a preferential axis in the a priori unoriented sample (photoselection). Because W359 and W306 are very differently oriented within the protein, detection with polarized light should allow us to distinguish them. To demonstrate this, W306 was mutated to redox-inert phenylalanine. We show that the resulting anisotropy spectrum of the initial absorption changes (measured at 10 ns time resolution) is in line with W359 being oxidized. The corresponding spectrum in wildtype photolyase is clearly different and identifies W306 as the oxidized species. These findings set an upper limit of 10 ns for electron transfer from W306 to W359*+ in wildtype DNA photolyase, consistent with previous, more indirect evidence [Aubert, C.; Vos, M. H.; Mathis, P.; Eker, A. P. M.; Brettel, K. Nature 2000, 405, 586-590].     

Dynamic characteristic of amitriptyline in water by ultrasonic relaxation method and molecular orbital calculation

Ultrasonic absorption coefficients in the frequency range of 0.8-220 MHz have been measured in aqueous solution of amitriptyline (3-(10,11-dihydro-5H-dibenzo[a,d]cycloheptene-5-ylidene)-N,N-dimethyl-1-propanamine) in the concentration range from 0.20 to 0.60 mol dm(-3) at 25 °C. A single relaxational phenomenon has been observed, and the relaxation frequency is independent of the concentration. It has been also observed that the amplitude of the relaxational absorption increases linearly with the analytical concentration. From these ultrasonic relaxation data, it has been concluded that the relaxation is associated with a unimolecular reaction due to a conformational change of the solute molecule, such as a structural change due to a rotational motion of a group in the solute molecule. Molecular orbital semiempirical methods using AM1 (Austin model 1) and PM3 (modified neglect of diatomic overlap parametric method 3) have been applied to obtain the standard enthalpy of formation for amitriptyline molecule at various dihedral angles around one of the bonds in alkylamine side chain. The results have shown the two clear minimum standard enthalpies of formation for amitriptyline. From the difference of the two values, the standard enthalpy change between the two stable conformers has been calculated be 2.9 kJ mol(-1). On a rough assumption that the standard enthalpy change reflects the standard free energy change, the equilibrium constant for the rotational isomers has been estimated to be 0.31. Combining this value with the experimental ultrasonic relaxation frequency, the backward and forward rate constants have been evaluated. The standard enthalpy change of the reaction has been also estimated from the concentration dependence of the maximum absorption per wavelength, and it has been close to that calculated by the semiempirical methods. The ultrasonic absorption measurements have been also carried out in amitriptyline solution in the presence of β-cyclodextrin. However, the ultrasonic relaxation has not been found in the above frequency range. The result has been discussed in relation to the host-guest complex formation between β-cyclodextrin and amitriptyline.     

E(2)-Elimination in the Decomposition of N-Bromoamino Acid Anions

The kinetics of oxidation of eight structurally different amino acids by hypobromite ion (BrO(-)) in the presence of hydroxide ion has been studied. The reactions proceed through the rapid formation of N-bromoamino acid anion which then decomposes in the rate-limiting step. The decomposition of N-bromoamino acid anions is found to proceed through an unimolecular and a base-catalyzed reaction. The large negative values of DeltaS() and the products formed suggest that the hydroxide ion-catalyzed reactions proceed through an E(2) mechanism. N-Bromoamino acid anion gives an absorption spectrum with lambda(max) at approximately 290 nm.     

2.4 A crystal structure of an oxaliplatin 1,2-d(GpG) intrastrand cross-link in a DNA dodecamer duplex

(1R,2R-Diaminocyclohexane)oxalatoplatinum(II) (oxaliplatin) is a third-generation platinum anticancer compound that produces the same type of inter- and intrastrand DNA cross-links as cisplatin. In combination with 5-fluorouracil, oxaliplatin has been recently approved in Europe, Asia, and Latin America for the treatment of metastatic colorectal cancer. We present here the crystal structure of an oxaliplatin adduct of a DNA dodecanucleotide duplex having the same sequence as that previously reported for cisplatin (Takahara, P. M.; Rosenzweig, A. C.; Frederick, C. A.; Lippard, S. J. Nature 1995, 377, 649-652). Pt-MAD data were used to solve this first X-ray structure of a platinated DNA duplex derived from an active platinum anticancer drug other than cisplatin. The overall geometry and crystal packing of the complex, refined to 2.4 A resolution, are similar to those of the cisplatin structure, despite the fact that the two molecules crystallize in different space groups. The platinum atom of the [Pt(R,R-DACH)](2+) moiety forms a 1,2-intrastrand cross-link between two adjacent guanosine residues in the sequence 5'-d(CCTCTGGTCTCC), bending the double helix by approximately 30 degrees toward the major groove. Both end-to-end and end-to-groove packing interactions occur in the crystal lattice. The latter is positioned in the minor groove opposite the platinum cross-link. A novel feature of the present structure is the presence of a hydrogen bond between the pseudoequatorial NH hydrogen atom of the (R,R)-DACH ligand and the O6 atom of the 3'-G of the platinated d(GpG) lesion. This finding provides structural evidence for the importance of chirality in mediating the interaction between oxaliplatin and duplex DNA, calibrating previously published models used to explain the reactivity of enantiomerically pure vicinal diamine platinum complexes with DNA in solution. It also provides a new kind of chiral recognition between an enantiomerically pure metal complex and the DNA double helix.     

Kinetics of unfolding the human telomeric DNA quadruplex using a PNA trap

The kinetics of opening of the DNA quadruplex formed by the human telomeric repeat have been investigated using real-time fluorescence resonance energy transfer (FRET) measurements with a peptide nucleic acid (PNA) trap. It has been found that this opening is zero-order with respect to PNA, indicating that the initial step is a rate-limiting internal rearrangement of the quadruplex. A study of the temperature dependence of the rate of quadruplex opening was performed and the activation energy of the process estimated to be 98 +/- 8 kJ mol(-1).     

A combined ONIOM quantum chemical-molecular dynamics study of zinc-uracil bond breaking in yeast cytosine deaminase

A QM/MM method that combines ONIOM quantum chemistry and molecular dynamics is developed and applied to a step in the deamination of cytosine to uracil in yeast cytosine deaminase (yCD). A two-layer ONIOM calculation is used for the reaction complex, with an inner part treated at a high level for the chemical reaction (bond breaking) and a middle part treated at a lower level for relevant protein residues that are frozen in the quantum optimization. An outer layer (protein and solvent) is treated using MD. Configurations for the entire system are generated using MD and optimized with ONIOM. The method permits the use of high-level quantum calculations along with sufficient configurational sampling to approximate the potential of mean force for certain bond-breaking reactions. A previously proposed reaction mechanism for deamination (Sklenak, S.; Yao, L. S.; Cukier, R. I.; Yan, H. G. J. Am. Chem. Soc. 2004, 126, 14879) requires breaking the bond between a catalytic zinc and the O4 of uracil in order to permit product release. Using an ONIOM approach, direct bond cleavage was found to be energetically unfavorable. In the work presented here, the combined ONIOM MD method is used to show that the barrier for bond cleavage is small, approximately 3 kcal/mol, and, consequently, should not be the rate-limiting step in the reaction.     

Thermodynamic profiles for CO photodissociation from heme model compounds: effect of proximal ligands

Here we present a comprehensive study of the thermodynamic parameters (enthalpy, entropy, and volume changes) associated with carbon monoxide photodissociation and rebinding to Fe(II) microperoxidase-11 (Fe(II)MP11) and Fe(ll) tetrakis(4-sulfonatophenyl)porphine complex (FeII4SP) with water and 2-methylimidazole as proximal ligands. CO photodissociation from FeII4SP complexes is accompanied by a positive volume change of approximately 17 mL mol(-1). A smaller volume change of approximately 12 mL mol(-1) was observed for CO dissociation from Fe(II)MP-11. We attribute the positive volume change to cleavage of the Fe-CO covalent bond and to solvent reorganization due to the low-spin to high-spin transition. CO binding is an exothermic reaction with an enthalpy change of -17 kcal mol(-1) for the CO-FeII4SP complexes and -13 kcal mol(-1) for the CO-Fe(II)MP11 complex. In all cases, the ligand recombination occurs as a single-exponential process indicating that CO dissociation is followed by direct CO rebinding to a high-spin five-coordinate complex without concomitant dissociation of the proximal base. In addition, observed negative activation entropies and volumes for ligand binding to (2-Melm)FeII4SP and MP-11, respectively, suggest that CO rebinding can be described by an associative mechanism with bond formation being the rate-limiting step.