c-Abl is normally regulated by an autoinhibitory mechanism, the disruption of which leads to chronic myelogenous leukemia. The details of this mechanism have been elusive because c-Abl lacks a phosphotyrosine residue that triggers the assembly of the autoinhibited form of the closely related Src kinases by internally engaging the SH2 domain. Crystal structures of c-Abl show that the N-terminal myristoyl modification of c-Abl 1b binds to the kinase domain and induces conformational changes that allow the SH2 and SH3 domains to dock onto it. Autoinhibited c-Abl forms an assembly that is strikingly similar to that of inactive Src kinases but with specific differences that explain the differential ability of the drug STI-571/Gleevec/imatinib (STI-571) to inhibit the catalytic activity of Abl, but not that of c-Src.

A sequence of ordered solvent peaks in the electron density map of the minor groove region of ApT-rich tracts of the double helix is a characteristic of B-form DNA well established from crystallography. This feature, termed the “spine of hydration”, has been discussed as a stabilizing feature of B-DNA, the structure of which is known to be sensitive to environmental effects. Nanosecond-range molecular dynamics simulations on the DNA duplex of sequence d(CGCGAATTCGCG) have been carried out, including explicit consideration of ∼4000 water molecules and 22 Na+ counterions, and based on the new AMBER 4.1 force field with the particle mesh Ewald summation used in the treatment of long-range interactions. The calculations support a dynamical model of B-DNA closer to the B form than any previously reported. Analysis of the dynamical structure of the solvent revealed that, in over half of the trajectory, a Na+ ion is found in the minor groove localized at the ApT step. This position, termed herein the “ApT pocket”, was noted previously (Lavery, R.; Pullman, B. J. Biomol. Struct. Dyn. 1985, 5, 1021) to be of uniquely low negative electrostatic potential relative to other positions of the groove, a result supported by the location of a Na+ ion in the crystal structure of the dApU miniduplex [Seeman, N.; et al. J. Mol. Biol. 1976, 104, 109) and by additional calculations described herein based on continuum electrostatics. The Na+ ion in the ApT pocket interacts favorably with the thymine O2 atom on opposite strands of the duplex and is well articulated with the water molecules which constitute the remainder of the minor groove spine. This result indicates that counterions may intrude on the minor groove spine of hydration on B-form DNA and subsequently influence the environmental structure and thermodynamics in a sequence-dependent manner. The observed narrowing of the minor groove in the AATT region of the d(CGCGAATTCGCG) structure may be due to direct binding effects and also to indirect modulation of the electrostatic repulsions that occur when a counterion resides in the minor groove “AT pocket”. The idea of localized complexation of otherwise mobile counterions in electronegative pockets in the grooves of DNA helices introduces a heretofore mostly unappreciated source of sequence-dependent effects on local conformational, helicoidal, and morphological structure and may have important implications in understanding the functional energetics and specificity of the interactions of DNA and RNA with regulatory proteins, pharmaceutical agents, and other ligands.

The improper activation of the Abl tyrosine kinase results in chronic myeloid leukemia (CML). The recognition of an inactive conformation of Abl, in which a catalytically important Asp-Phe-Gly (DFG) motif is flipped by approximately 180° with respect to the active conformation, underlies the specificity of the cancer drug imatinib, which is used to treat CML. The DFG motif is not flipped in crystal structures of inactive forms of the closely related Src kinases, and imatinib does not inhibit c-Src. We present a structure of the kinase domain of Abl, determined in complex with an ATP–peptide conjugate, in which the protein adopts an inactive conformation that resembles closely that of the Src kinases. An interesting aspect of the Src-like inactive structure, suggested by molecular dynamics simulations and additional crystal structures, is the presence of features that might facilitate the flip of the DFG motif by providing room for the phenylalanine to move and by coordinating the aspartate side chain as it leaves the active site. One class of mutations in BCR–Abl that confers resistance to imatinib appears more likely to destabilize the inactive Src-like conformation than the active or imatinib-bound conformations. Our results suggest that interconversion between distinctly different inactive conformations is a characteristic feature of the Abl kinase domain.

The effect of C-terminal tyrosine phosphorylation on molecular motions in the Src kinases Hck and c-Src is investigated by molecular dynamics simulations. The SH2 and SH3 domains of the inactive kinases are seen to be tightly coupled by the connector between them, impeding activation. Dephosphorylation of the tail reduces the coupling between the SH2 and SH3 domains in the simulations, as does replacement of connector residues with glycine. A mutational analysis of c-Src expressed in Schizosaccharomyces pombe demonstrates that replacement of residues in the SH2-SH3 connector with glycine activates c-Src. The SH2-SH3 connector appears to be an inducible "snap lock" that clamps the SH2 and SH3 domains upon tail phosphorylation, but which allows flexibility when the tail is released.

We report the results of four new molecular dynamics (MD) simulations on the DNA duplex of sequence d(CGCGAATTCGCG)2, including explicit consideration of solvent water, and a sufficient number of Na+ counterions to provide electroneutrality to the system. Our simulations are configured particularly to characterize the latest MD models of DNA, and to provide a basis for examining the sensitivity of MD results to the treatment of boundary conditions, electrostatics, initial placement of solvent, and run lengths. The trajectories employ the AMBER 4.1 force field. The simulations use particle mesh Ewald summation for boundary conditions, and range in length from 500 ps to 5.0 ns. Analysis of the results is carried out by means of time series for conformationalm, helicoidal parameters, newly developed indices of DNA axis bending, and groove widths. The results support a dynamically stable model of B-DNA for d(CGCGAATTCGCG)2 over the entire length of the trajectory. The MD results are compared with corresponding crystallographic and NMR studies on the d(CGCGAATTCGCG)2 duplex, and placed in the context of observed behavior of B-DNA by comparisons with the complete crystallographic data base of B-form structures. The calculated distributions of mobile solvent molecules, both water and counterions, are displayed. The calculated solvent structure of the primary solvation shell is compared with the location of ordered solvent positions in the corresponding crystal structure. The results indicate that ordered solvent positions in crystals are roughly twice as structured as bulk water. Detailed analysis of the solvent dynamics reveals evidence of the incorporation of ions in the primary solvation of the minor groove B-form DNA. The idea of localized complexation of otherwise mobile counterions in electronegative pockets in the grooves of DNA helices introduces an additional source of sequence-dependent effects on local conformational, helicoidal, and morphological structure, and may have important implications for understanding the functional energetics and specificity of the interactions of DNA and RNA with regulatory proteins, pharmaceutical agents, and other ligands.

The dimeric ring-shaped sliding clamp of E. coli DNA polymerase III (beta subunit, homolog of eukaryotic PCNA) is loaded onto DNA by the clamp loader gamma complex (homolog of eukaryotic Replication Factor C, RFC). The delta subunit of the gamma complex binds to the beta ring and opens it. The crystal structure of a beta:delta complex shows that delta, which is structurally related to the delta' and gamma subunits of the gamma complex, is a molecular wrench that induces or traps a conformational change in beta such that one of its dimer interfaces is destabilized. Structural comparisons and molecular dynamics simulations suggest a spring-loaded mechanism in which the beta ring opens spontaneously once a dimer interface is perturbed by the delta wrench.

We describe herein a computationally intensive project aimed at carrying out molecular dynamics (MD) simulations including water and counterions on B-DNA oligomers containing all 136 unique tetranucleotide base sequences. This initiative was undertaken by an international collaborative effort involving nine research groups, the "Ascona B-DNA Consortium" (ABC). Calculations were carried out on the 136 cases imbedded in 39 DNA oligomers with repeating tetranucleotide sequences, capped on both ends by GC pairs and each having a total length of 15 nucleotide pairs. All MD simulations were carried out using a well-defined protocol, the AMBER suite of programs, and the parm94 force field. Phase I of the ABC project involves a total of approximately 0.6 mus of simulation for systems containing approximately 24,000 atoms. The resulting trajectories involve 600,000 coordinate sets and represent approximately 400 gigabytes of data. In this article, the research design, details of the simulation protocol, informatics issues, and the organization of the results into a web-accessible database are described. Preliminary results from 15-ns MD trajectories are presented for the d(CpG) step in its 10 unique sequence contexts, and issues of stability and convergence, the extent of quasiergodic problems, and the possibility of long-lived conformational substates are discussed.

The Abl and Src tyrosine kinases are key signaling proteins that are of considerable interest as drug targets in cancer and many other diseases. The regulatory mechanisms that control the activity of these proteins are complex, and involve large-scale conformational changes in response to phosphorylation and other modulatory signals. The success of the Abl inhibitor imatinib in the treatment of chronic myelogenous leukemia has shown the potential of kinase inhibitors, but the rise of drug resistance in patients has also shown that drugs with alternative modes of binding to the kinase are needed. The detailed understanding of mechanisms of protein-drug interaction and drug resistance through biophysical methods demands a method for the production of active protein on the milligram scale. We have developed a bacterial expression system for the kinase domains of c-Abl and c-Src, which allows for the quick expression and purification of active wild-type and mutant kinase domains by coexpression with the YopH tyrosine phosphatase. This method makes practical the use of isotopic labeling of c-Abl and c-Src for NMR studies, and is also applicable for constructs containing the SH2 and SH3 domains of the kinases.