Dynamical signature of abasic damage in DNA.

Time-dependent Stokes shift (TDSS) responses in proteins and DNA exhibit a broad range of long time scales (>10 ps) that are not present in bulk aqueous solution. The physical interpretation of the long TDSS time scales in biomolecular systems is a matter of considerable debate because of the many different components present in the sample (water, biomolecule, counterions), which have highly correlated motions and intrinsically different abilities to adapt to local perturbations. Here we use molecular dynamics (MD) simulations to show that the surprisingly slow (∼10 ns) TDSS response of coumarin 102 (C102), a base pair replacement, reflects a distinct dynamical signature for DNA damage. When the C102 molecule is covalently incorporated into DNA, an abasic site is created on the strand opposite the C102 probe. The abasic sugar exhibits a reversible interchange between intra- and extrahelical conformations that are kinetically stable on a nanosecond time scale. This conformational change, only possible in damaged DNA, was found to be responsible for the long time scales in the measured TDSS response. For the first time, a TDSS measurement has been attributed to a specific biomolecular motion. This finding directly contradicts the prevailing notion that the TDSS response in biomolecular contexts is dominated by hydration dynamics. It also suggests that TDSS experiments can be used to study ultrafast biomolecular dynamics that are inaccessible to other techniques.

Nitrile groups as vibrational probes of biomolecular structure and dynamics: an overview.

This paper presents an overview of recent experiments and theoretical developments aimed at using vibrational spectroscopy to understand the structure and dynamics of nitrile-labeled biomolecules. Nitrile groups are excellent vibrational probes of proteins and DNA because they absorb in a region of the spectrum that is relatively free of absorption due to the biomolecule, and they have high extinction coefficients. The vibrational frequency of nitrile groups is also extraordinarily sensitive to its local environment, and thus C[triple bond, length as m-dash]N bonds have been employed in both linear and 2-D infrared (IR) spectroscopy experiments and also as vibrational Stark probes of electric fields in proteins. The interpretation and design of these experiments would be enhanced by accurate calculations of IR spectra from molecular dynamics simulations. Recently, theoretical developments towards computing the vibrational spectrum of nitrile groups in the condensed-phase have been highly successful. A strong synergy between experiment and theory will further promote the use of vibrational spectroscopy of nitrile-labeled biomolecules to address fundamental questions of structure and dynamics that are elusive to other techniques.

Effects of Long-Range Electrostatics on Time-Dependent Stokes Shift Calculations.

Molecular dynamics simulations are essential to the correct interpretation of the response measured in time-dependent Stokes shift (TDSS) experiments of fluorescent probe molecules in biological environments. Within linear response theory, the TDSS response is the time correlation function of the fluctuations of ΔE(t), the difference between the solute environment interaction energy with the probe, modeled in both its electronically excited and ground states. ΔE(t) is dominated by electrostatic interactions between the environment and the ground- and excited-state charge distributions of the probe. The treatment of the long-ranged electrostatics in the calculation of the TDSS response in MD simulations is systematically investigated for three probes in aqueous solution: a model diatomic, coumarin 102, and Hoechst 33258. Nine different protocols for the treatment of the electrostatics were compared to particle mesh Ewald (PME), which was utilized as a reference standard. A computationally efficient pairwise alternative to PME, the damped shifted force method, was shown to reproduce the TDSS response calculated with PME for all three systems. In contrast, neglecting the role of the long-ranged electrostatics in the calculation of the TDSS response results in artifacts.

The dynamics of water at DNA interfaces: computational studies of Hoechst 33258 bound to DNA.

Together, spectroscopy combined with computational studies that relate directly to the experimental measurements have the potential to provide unprecedented insight into the dynamics of important biological processes. Recent time-resolved fluorescence experiments have shown that the time scales for collective reorganization at the interface of proteins and DNA with water are more than an order of magnitude slower than in bulk aqueous solution. The molecular interpretation of this change in the collective response is somewhat controversial some attribute the slower reorganization to dramatically retarded water motion, while others describe rapid water dynamics combined with a slower biomolecular response. To connect directly to solvation dynamics experiments of the fluorescent probe Hoechst 33258 (H33258) bound to DNA, we have generated 770 ns of molecular dynamics (MD) simulations and calculated the equilibrium and nonequilibrium solvation response to excitation of the probe. The calculated time scales for the solvation response of H33258 free in solution (0.17 and 1.4 ps) and bound to DNA (1.5 and 20 ps) are highly consistent with experiment (0.2 and 1.2 ps, 1.4 and 19 ps, respectively). Decomposition of the calculated response revealed that water solvating the probe bound to DNA was still relatively mobile, only slowing by a factor of 2-3, while DNA motion was responsible for the long-time component (approximately 20 ps).

Solvation dynamics of Hoechst 33258 in water: an equilibrium and nonequilibrium molecular dynamics study.

Integrated within an appropriate theoretical framework, molecular dynamics (MD) simulations are a powerful tool to complement experimental studies of solvation dynamics. Together, experiment, theory, and simulation have provided substantial insight into the dynamic behavior of polar solvents. MD investigations of solvation dynamics are especially valuable when applied to the heterogeneous environments found in biological systems, where the calculated response of the environment to the electrostatic perturbation of the probe molecule can easily be decomposed by component (e.g., aqueous solvent, biomolecule, ions), greatly aiding the molecular-level interpretation of experiments. A comprehensive equilibrium and nonequilibrium MD study of the solvation dynamics of the fluorescent dye Hoechst 33258 (H33258) in aqueous solution is presented. Many fluorescent probes employed in experimental studies of solvation dynamics in biological systems, such as the DNA minor groove binder H33258, have inherently more conformational flexibility than prototypical fused-ring chromophores. The role of solute flexibility was investigated by developing a fully flexible force-field for the H33258 molecule and by simulating its solvation response. While the timescales for the total solvation response calculated using both rigid (0.16 and 1.3 ps) and flexible (0.17 and 1.4 ps) models of the probe closely matched the experimentally measured solvation response (0.2 and 1.2 ps), there were subtle differences in the response profiles, including the presence of significant oscillations for the flexible probe. A decomposition of the total response of the flexible probe revealed that the aqueous solvent was responsible for the overall decay, while the oscillations result from fluctuations in the electrostatic terms in the solute intramolecular potential energy. A comparison of equilibrium and nonequilibrium approaches for the calculation of the solvation response confirmed that the solvation dynamics of H33258 in water is well-described by linear response theory for both rigid and flexible models of the probe.

Molecular dynamics simulations of arachidonic acid complexes with COX-1 and COX-2: insights into equilibrium behavior.

The cyclooxygenase (COX) enzymes are responsible for the committed step in prostaglandin biosynthesis, the generation of prostaglandin H(2). As a result, these enzymes are pharmacologically important targets for nonsteroidal antiinflammatory drugs, such as aspirin and newer COX-2 selective inhibitors. The cyclooxygenases are functional homodimers, and each subunit contains both a cyclooxygenase and a peroxidase active site. These enzymes are quite interesting mechanistically, as the conversion of arachidonic acid to prostaglandin H(2) requires two oxygenation and two cyclization reactions, resulting in the formation of five new chiral centers with nearly absolute regio- and stereochemical fidelity. We have used molecular dynamics (MD) simulations to investigate the equilibrium behavior of both COX-1 and COX-2 enzyme isoforms with bound arachidonate. These simulations were compared with reference simulations of arachidonate in solution to explore the effect of enzyme on substrate conformation and positioning in the active site. The simulations suggest that the substrate has greater conformational freedom in the COX-2 active site, consistent with the larger COX-2 active site volume observed in X-ray crystal structures. The simulations reveal different conformational behavior for arachidonate in each subunit over the course of extended equilibrium MD simulations. The simulations also provide detailed information for several protein channels that might be important for oxygen and water transport to or from active sites or for intermediate trafficking between the cyclooxygenase and peroxidase active sites. The detailed comparisons for COX-1 versus COX-2 active site structural fluctuations may also provide useful information for design of new isozyme-selective inhibitors.

Molecular dynamics simulations of arachidonic acid-derived pentadienyl radical intermediate complexes with COX-1 and COX-2: insights into oxygenation regio- and stereoselectivity.

The two cyclooxygenase enzymes, COX-1 and COX-2, are responsible for the committed step in prostaglandin biosynthesis and are the targets of the nonsteroidal antiinflammatory drugs aspirin and ibuprofen and the COX-2 selective inhibitors, Celebrex, Vioxx, and Bextra. The enzymes are remarkable in that they catalyze two dioxygenations and two cyclizations of the native substrate, arachidonic acid, with near absolute regio- and stereoselectivity. Several theories have been advanced to explain the nature of enzymatic control over this series of reactions, including suggestions of steric shielding and oxygen channeling. As proposed here, selective radical trapping and spin localization in the substrate-derived pentadienyl radical intermediate can also be envisioned. Herein we describe the results of explicit, 10 ns molecular dynamics simulations of both COX-1 and COX-2 with the substrate-derived pentadienyl radical intermediate bound in the active site. The enzymes' influence on the conformation of the pentadienyl radical was investigated, along with the accessible space above and below the radical plane and the width of several channels to the active site that could function as access routes for molecular oxygen. Additional simulations demonstrated the extent of molecular oxygen mobility within the active site. The results suggest that spin localization is unlikely to play a role in enzymatic control of this reaction. Instead, a combination of oxygen channeling, steric shielding, and selective radical trapping appears to be responsible. This work adds a dynamic perspective to the strong foundation of static structural data available for these enzymes.

Three-dimensional models for beta-adrenergic receptor complexes with agonists and antagonists.

Molecular modeling methods have been used to construct three-dimensional models for agonist and antagonist complexes with beta-adrenergic receptors. The recent rhodopsin crystal structure was used as a template in standard homology modeling methods. The rhodopsin-based homology models were assessed for agreement with experimental results for beta-adrenergic receptors, and compared with receptor models developed using de novo modeling techniques. While the de novo and homology-derived receptor models are generally quite similar, there are some localized structural differences that impact the putative ligand-binding site significantly. The de novo receptor models appear to provide much better agreement with experimental data, particularly for receptor complexes with agonist ligands. The de novo receptor models also yield some interesting and testable hypotheses for the structural basis of beta-adrenergic receptor subtype ligand selectivity.

Refinement of the conformation of a critical region of charge-charge interaction between cholecystokinin and its receptor.

Insight into the molecular basis of cholecystokinin (CCK) binding to its receptor has come from receptor mutagenesis and photoaffinity labeling studies, with both contributing to the current hypothesis that the acidic Tyr-sulfate-27 residue within the peptide is situated adjacent to basic Arg(197) in the second loop of the receptor. Here, we refine our understanding of this region of interaction by examining a structure-activity series of these positions within both ligand and receptor and by performing three-dimensional molecular modeling of key pairs of modified ligand and receptor constructs. The important roles of Arg(197) and Tyr-sulfate-27 were supported by the marked negative impact on binding and biological response with their natural partner molecule when the receptor residue was replaced by acidic Asp or Glu and when the peptide residue was replaced by basic Arg, Lys, p-amino-Phe, p-guanidino-Phe, or p-methylamino-Phe. Complementary ligand-receptor charge-exchange experiments were unable to regain the lost function. This was supported by the molecular modeling, which demonstrated that the charge-reversed double mutants could not form a good interaction without extensive rearrangement of receptor conformation. The models further predicted that R197D and R197E mutations would lead to conformational changes in the extracellular domain, and this was experimentally supported by data showing that these mutations decreased peptide agonist and antagonist binding and increased nonpeptidyl antagonist binding. These receptor constructs also had increased susceptibility to trypsin degradation relative to the wild-type receptor. In contrast, the relatively conservative R197K mutation had modest negative impact on peptide agonist binding, again consistent with the modeling demonstration of loss of a series of stabilizing inter- and intramolecular bonds. The strong correlation between predicted and experimental results support the reported refinement in the three-dimensional structure of the CCK-occupied receptor.