Unusual sequence effects on nucleotide excision repair of arylamine lesions: DNA bending/distortion as a primary recognition factor.

The environmental arylamine mutagens are implicated in the etiology of various sporadic human cancers. Arylamine-modified dG lesions were studied in two fully paired 11-mer duplexes with a -G*CN- sequence context, in which G* is a C8-substituted dG adduct derived from fluorinated analogs of 4-aminobiphenyl (FABP), 2-aminofluorene (FAF) or 2-acetylaminofluorene (FAAF), and N is either dA or dT. The FABP and FAF lesions exist in a simple mixture of 'stacked' (S) and 'B-type' (B) conformers, whereas the N-acetylated FAAF also samples a 'wedge' (W) conformer. FAAF is repaired three to four times more efficiently than FABP and FAF. A simple A- to -T polarity swap in the G*CA/G*CT transition produced a dramatic increase in syn-conformation and resulted in 2- to 3-fold lower nucleotide excision repair (NER) efficiencies in Escherichia coli. These results indicate that lesion-induced DNA bending/thermodynamic destabilization is an important DNA damage recognition factor, more so than the local S/B-conformational heterogeneity that was observed previously for FAF and FAAF in certain sequence contexts. This work represents a novel 3'-next flanking sequence effect as a unique NER factor for bulky arylamine lesions in E. coli.

Structural mechanism associated with domain opening in gain-of-function mutations in SHP2 phosphatase.

The SHP2 phosphatase plays a central role in a number of signaling pathways were it dephosphorylates various substrate proteins. Regulation of SHP2 activity is, in part, achieved by an intramolecular interaction between the PTP domain of the protein, which contains the catalytic site, and the N-SH2 domain leading to a "closed" protein conformation and autoinhibition. Accordingly, "opening" of the N-SH2 and PTP domains is required for the protein to become active. Binding of phosphopeptides to the N-SH2 domain is known to induce the opening event, while a number of gain-of-function (GOF) mutants, implicated in Noonan's Syndrome and childhood leukemias, are thought to facilitate opening. In the present study, a combination of computational and experimental methods are used to investigate the structural mechanism of opening of SHP2 and the impact of three GOF mutants, D61G, E76K, and N308D, on the opening mechanism. Calculated free energies of opening indicate that opening must be facilitated by effector molecules, possibly the protein substrates themselves, as the calculated free energies preclude spontaneous opening. Simulations of both wild type (WT) SHP2 and GOF mutants in the closed state indicate GOF activity to involve increased solvent exposure of selected residues, most notably Arg362, which in turn may enhance interactions of SHP2 with its substrate proteins and thereby aid opening. In addition, GOF mutations cause structural changes in the phosphopeptide-binding region of the N-SH2 domain leading to conformations that mimic the bound state. Such conformational changes are suggested to enhance binding of phosphopeptides and/or decrease interactions between the PTP and N-SH2 domains thereby facilitating opening. Experimental assays of the impact of effector molecules on SHP2 phosphatase activity against both small molecule and peptide substrates support the hypothesized mechanism of GOF mutant action. The present calculations also suggest a role for the C-SH2 domain of SHP2 in stabilizing the overall conformation of the protein in the open state, thereby aiding conformational switching between the open active and closed inactive states.

Nearest-neighbor nonparametric method for estimating the configurational entropy of complex molecules.

A method for estimating the configurational (i.e., non-kinetic) part of the entropy of internal motion in complex molecules is introduced that does not assume any particular parametric form for the underlying probability density function. It is based on the nearest-neighbor (NN) distances of the points of a sample of internal molecular coordinates obtained by a computer simulation of a given molecule. As the method does not make any assumptions about the underlying potential energy function, it accounts fully for any anharmonicity of internal molecular motion. It provides an asymptotically unbiased and consistent estimate of the configurational part of the entropy of the internal degrees of freedom of the molecule. The NN method is illustrated by estimating the configurational entropy of internal rotation of capsaicin and two stereoisomers of tartaric acid, and by providing a much closer upper bound on the configurational entropy of internal rotation of a pentapeptide molecule than that obtained by the standard quasi-harmonic method. As a measure of dependence between any two internal molecular coordinates, a general coefficient of association based on the information-theoretic quantity of mutual information is proposed. Using NN estimates of this measure, statistical clustering procedures can be employed to group the coordinates into clusters of manageable dimensions and characterized by minimal dependence between coordinates belonging to different clusters.

Estimation of the absolute internal-rotation entropy of molecules with two torsional degrees of freedom from stochastic simulations.

A method of statistical estimation is applied to the problem of evaluating the absolute entropy of internal rotation in a molecule with two torsional degrees of freedom. The configurational part of the entropy is obtained as that of the joint probability density of an arbitrary form represented by a two-dimensional Fourier series, the coefficients of which are statistically estimated using a sample of the torsional angles of the molecule obtained by a stochastic simulation. The internal rotors in the molecule are assumed to be attached to a common frame, and their reduced moments of inertia are initially calculated as functions of the two torsional angles, but averaged over all the remaining internal degrees of freedom using the stochastic-simulation sample of the atomic configurations of the molecule. The torsional-angle dependence of the reduced moments of inertia can be also averaged out, and the absolute internal-rotation entropy of the molecule is obtained in a good approximation as the sum of the configurational entropy and a kinetic contribution fully determined by the averaged reduced moments of inertia. The method is illustrated using Monte Carlo simulations of isomers of stilbene and halogenated derivatives of propane. The two torsional angles in cis-stilbene are found to be much more strongly correlated than those in trans-stilbene, while the degree of the angular correlation in propane increases strongly on substitution of hydrogen atoms with chlorine.

Application of molecular dynamics simulations to spin-labeled oligonucleotides.

The EPR study of spin labeled macromolecules has provided insight into structural and dynamical properties of DNA, proteins, and related systems. While spin labeling has been useful, it is experimentally difficult to determine if the spin label significantly alters the conformation of the macromolecule to which it is attached. Molecular modeling has proven to be a powerful tool for studying structure and dynamics of biologically important molecules. Here, we have conducted molecular dynamics (MD) studies of spin labeled oligonucleotides (ONs) bearing a five (5sp) or six (6sp) membered ring nitroxide, and the corresponding unmodified ON using the suite of programs contained in Amber 5.0 with the Cornell et al. 94 force field (Cornell, W. D., Cieplak, P., Bayly, C. I., Gould, I. R., Merz, Jr., K. M. Ferguson, D. M., Spellmeyer, D. C., Fox, T., Caldwell, J. W., and Kollman, P. A. A Second Generation Force Field for the Simulation of Proteins and Nucleic Acids. J. Am. Chem. Soc. 117, 5179-5197 (1995)). Quantum mechanical calculations employing the B3LYP method with the standard 6-31G* basis set using Gaussian98 were performed and, together with available crystallographic data for analogous nitroxides, new parameters for the nitrogen, oxygen, nitroxide alpha-carbon, and sp-hybridized carbon atoms have been developed suitable for the Cornell et al. 94 force field. MD simulations on the double-stranded (ds) spin labeled ONs, along with the corresponding unmodified analogues, have been studied over the course of 4 ns and conformational properties of all ONs are described based on the analysis of the trajectories. The spin labels were found to alter the global conformation of the ONs to which they were attached to accommodate the spin labels. The major changes include widening the major groove, decreasing helical twist, and more negative X-displacement of the base pairs. The magnitude of the effect was dependent on the specific structure of the spin label. Average and 'most representative' structures derived from the molecular dynamics simulations correlate with the experimental data on the spin labeled ONs.

Probing triplex formation by EPR spectroscopy using a newly synthesized spin label for oligonucleotides.

Spin labels have been extensively used to study the dynamics of oligonucleotides. Spin labels that are more rigidly attached to a base in an oligonucleotide experience much larger changes in their range of motion than those that are loosely tethered. Thus, their electron paramagnetic resonance spectra show larger changes in response to differences in the mobility of the oligonucleotides to which they are attached. An example of this is 5-(2,2,5,5-tetramethyl-3-ethynylpyrrolidine-1-oxyl)-uridine (1). How ever, the synthesis of this modified DNA base is quite involved and, here, we report the synthesis of a new spin-labeled DNA base, 5-(2,2,6,6-tetramethyl-4-ethynylpiperidyl-3-ene-1-oxyl)-uridine (2). This spin label is readily prepared in half the number of steps required for 1, and yet behaves in a spectroscopically analogous manner to 1 in oligonucleotides. Finally, it is shown here that both spin labels 1 and 2 can be used to detect the formation of both double-stranded and triplex DNA.