PubMed Health⌕ Search

Biomedical subjects

D S Goodsell

Publications and source records attributed to D S Goodsell.

At least 73 records · Page 4Linked to original sources

Design of stapled DNA-minor-groove-binding molecules with a mutable atom simulated annealing method.

We report the design of optimal linker geometries for the synthesis of stapled DNA-minor-groove-binding molecules. Netropsin, distamycin, and lexitropsins bind side-by-side to mixed-sequence DNA and offer an opportunity for the design of sequence-reading molecules. Stapled molecules, with two molecules covalently linked side-by-side, provide entropic gains and restrain the position of one molecule relative to its neighbor. Using a free-atom simulated annealing technique combined with a discrete mutable atom definition, optimal lengths and atomic composition for covalent linkages are determined, and a novel hydrogen bond 'zipper' is proposed to phase two molecules accurately side-by-side.

Base Sequence↗

Distributed automated docking of flexible ligands to proteins: parallel applications of AutoDock 2.4.

AutoDock 2.4 predicts the bound conformations of a small, flexible ligand to a nonflexible macromolecular target of known structure. The technique combines simulated annealing for conformation searching with a rapid grid-based method of energy evaluation based on the AMBER force field. AutoDock has been optimized in performance without sacrificing accuracy; it incorporates many enhancements and additions, including an intuitive interface. We have developed a set of tools for launching and analyzing many independent docking jobs in parallel on a heterogeneous network of UNIX-based workstations. This paper describes the current release, and the results of a suite of diverse test systems. We also present the results of a systematic investigation into the effects of varying simulated-annealing parameters on molecular docking. We show that even for ligands with a large number of degrees of freedom, root-mean-square deviations of less than 1 A from the crystallographic conformation are obtained for the lowest-energy dockings, although fewer dockings find the crystallographic conformation when there are more degrees of freedom.

Azepines↗

Structure of a dicationic monoimidazole lexitropsin bound to DNA.

An X-ray crystal structure has been solved of the complex of a dicationic lexitropsin with a B-DNA duplex of sequence CGCGAATTCGCG. The lexitropsin is identical to netropsin except for replacement of the first methylpyrrole ring by methylimidazole, converting a =CH- to =N-. Crystals are isomorphous with those of the DNA dodecamer in the absence of drug. Although the =N- for =CH- substitution was intended to make that locus on the drug molecule compatible with a G.C base pair, electrostatic attraction for the two cationic ends of the drug predominates, and this lexitropsin binds to the same central AATT site as does the parent netropsin. But unlike netropsin, this lexitropsin exhibits end-for-end disorder in the crystal. Both orientations were refined separately to completion. Final residual errors at 2.25 A resolution for the 2358 reflections above 2 sigma in F are R = 0.165 for one orientation (LexA) with 37 water molecules and 0.164 for the inverted drug orientation (LexB) with 40 water molecules. This molecular disorder is probably attributable to a weakening of binding to the AATT site occasioned by the imidazole-for-pyrrole substitution.

Antibiotics, Antineoplastic↗

Refinement of netropsin bound to DNA: bias and feedback in electron density map interpretation.

The X-ray crystal structure of the complex of the B-DNA dodecamer CGCGAATTCGCG with the antitumor drug netropsin has been reexamined to locate the drug accurately for computer-based drug design. The optimum solution is with the drug centered in the AATT region of the minor groove, making three good bifurcated hydrogen bonds with adenine N3 and thymine O2 atoms along the floor of the groove. Pyrrole rings of netropsin are packed against the C2 positions of adenines, leaving no room for the amine group of guanine and, hence, providing a structural rationale for the A.T specificity of netropsin. An alternative positioning in which the drug is shifted along the minor groove by ca. one-half base pair step is rejected on the basis of free R factor calculations and the appearance of the original drug-free difference maps. Final omit maps, although of more pleasing appearance, are not a dependable means of discriminating between right and wrong structures. The shifted alternative drug position ignores potential hydrogen bonding along the floor of the groove, provides no explanation for netropsin's observed A.T specificity, and is contradicted by NMR results [Patel, D. J. (1982) Proc. Natl. Acad. Sci. U.S.A. 79, 6424].

Base Sequence↗

Crystal structure of C-T-C-T-C-G-A-G-A-G. Implications for the structure of the Holliday junction.

The structure of the B-DNA decamer of sequence C-T-C-T-C-G-A-G-A-G shows a crossed arrangement of helices in the C2 crystal lattice. This is the fourth example of a crossed arrangement of B-DNA oligomers in a crystal, and in spite of the fact that each of these four crystallizes in a different lattice, all have nearly identical structures at the crossing contact. This ubiquitous crossing arrangement may be used to generate a structure for the Holliday junction that is fully consistent with the available physical data.

Base Sequence↗

Design of B-DNA cross-linking and sequence-reading molecules.

We report the design of hybrid molecules to bind in the minor groove of B-DNA, which combine DNA alkylating and cross-linking ability for increased chemotherapeutic efficacy, with sequence specificity, to minimize side effects. Optimal linkage geometries have been determined for the synthesis of bis-anthramycin and anthramycin-netropsin hybrid molecules. Earlier studies on linked drugs have typically been based on molecular mechanics calculations. This work, in contrast, uses the observed crystal structures of a netropsin/DNA complex and a new anthramycin/DNA complex to determine the exact spacing between two individual drugs when bound in the minor groove of B-DNA. Molecular linkers then are designed and tested between these two experimental positions, to form a chimeric or bis-linked compound molecule. A linked anthramycin-netropsin molecule has been designed specifically to target the polypurine tract second-strand primer site of the reverse transcriptase of HIV-1.

Base Sequence↗

1994 Molecular Graphics Art Show and Video Show.

The 1994 Molecular Graphics Art Show and Video Show were presented at the 13th annual international meeting of the Molecular Graphics and Modelling Society. The art show--shown in the Mary & Leigh Block Gallery on the campus of Northwestern University, Evanston, Illinois--included original artworks by eighteen artists and the video show included nine original animated works. All were chosen for their ability to present the complexity, diversity, and beauty of the molecular world in visual form. Works from a wide range of disciplines were represented, including work by scientists actively involved in structural research, by commercial illustrators presenting these results to students and physicians, and by fine artists exploring the meanings and implications of these molecules in our lives. Included in this issue of the Journal of Molecular Graphics are comments by the juror of the show, T.J. O'Donnell, a catalogue of the art show, and a catalogue of the video show.

Art↗

Bending and curvature calculations in B-DNA.

A simple program, BEND, has been written to calculate the magnitude of local bending and macroscopic curvature at each point along an arbitrary B-DNA sequence, using any desired bending model that specifies values of twist, roll and tilt as a function of sequence. The program has been used to evaluate six different DNA bending models in three categories. Two are bent non-A-tract models: (a) A new model based on the nucleosome positioning data of Satchwell et al 1986 (J. Mol. Biol. 191, 659-675), (b) The model of Calladine et al 1988 (J. Mol. Biol. 201, 127-137). Three are bent A-tract models: (c) The wedge model of Bolshoy et al 1991 (Proc. Natl. Acad. Sci. USA 88, 2312-2316), (d) The model of Cacchione et al 1989 (Biochem. 28, 8706-8713), (e) A reversed version of model (b). The last is a junction model: (f) The model of Koo & Crothers 1988 (Proc. Natl. Acad. Sci. USA 85, 1763-1767). Although they have widely different assumptions and values for twist, roll and tilt, all six models correctly predict experimental A-tract curvature as measured by gel retardation and cyclization kinetics, but only the new nucleosome positioning model is successful in predicting curvature in regions containing phased GGGCCC sequences. This model--showing local bending at mixed sequence DNA, strong bends at the sequence GGC, and straight, rigid A-tracts--is the only model consistent with both solution data from gel retardation and cyclization kinetics and structural data from x-ray crystallography.

Base Sequence↗

Crystal structure of a covalent DNA-drug adduct: anthramycin bound to C-C-A-A-C-G-T-T-G-G and a molecular explanation of specificity.

A 2.3-A X-ray crystal structure analysis has been carried out on the antitumor drug anthramycin, covalently bound to a ten base pair DNA double helix of sequence C-C-A-A-C-G-T-T-G-G. One drug molecule sits within the minor groove at each end of the helix, covalently bound through its C11 position to the N2 amine of the penultimate guanine of the chain. The stereochemical conformation is C11S, C11aS. The natural twist of the anthramycin molecule in the C11aS conformation matches the twist of the minor groove as it winds along the helix; a C11aR drug would only fit into a left-handed helix. The C11S attachment is roughly equatorial to the overall plane of the molecule, whereas a C11R attachment would be axial and would obstruct the fitting of the drug into the groove. The six-membered ring of anthramycin points toward the 3' end of the chain to which it is covalently attached or toward the end of the helix. The acrylamide tail attached to the five-membered ring extends back along the minor groove toward the center of the helix, binding in a manner reminiscent of netropsin or distamycin. The drug-DNA complex is stabilized by hydrogen bonds from C9-OH, N10, and the end of the acrylamide tail to base pair edges on the floor of the minor groove. The origin of anthramycin specificity for three successive purines arises not from specific hydrogen bonds but from the low twist angles adopted by purine-purine steps in a B-DNA helix. Binding of anthramycin induces a low twist at T-G in the T-G-G sequence of this DNA-drug complex, by comparison with the structure of the free DNA. The origin of anthramycin's preference for adenines flanking the alkylated guanine arises from a netropsin-like fitting of the acrylamide tail into the minor groove.

Anthramycin↗

The crystal structure of C-C-A-T-T-A-A-T-G-G. Implications for bending of B-DNA at T-A steps.

The single-crystal X-ray analysis of trigonal C-C-A-T-T-A-A-T-G-G, and its comparison with orthorhombic C-G-A-T-T-A-A-T-C-G, have shown that the A-T-T-A-A-T sequence has limited polymorphism under the influence of packing forces from neighboring molecules in the crystal. The T-A step is intrinsically variable. It is not inconsistent with a large propeller twist, a narrow minor groove, and a single spine of hydration, as has sometimes been claimed on theoretical grounds. The T-A step does show a persistent positive roll, in a direction that compresses the major groove, and this may be a significant factor in macroscopic DNA curvature induced by phased A-tracts. A-tracts, as understood in this paper, include A-A and A-T steps, but not the T-A step, which is disruptive. Three conclusions regarding A-tract-induced curvature can be drawn from this and other X-ray crystal structure analyses, and from key gel retardation experiments: (1) The A-tract bending model is disqualified on two grounds: (i) tilt-wedge bending within A-tracts is incompatible with the observed direction of curvature; (ii) roll-wedge bending within A-tracts is contradicted by every crystal structure analysis, and is inconsistent with gel retardation results for (G-C-A-A-A-A-T-T-T-T)n and for (A-A-A-A-A-T-T-T-T-T)n. (2) The junction bend model is contradicted by crystallography because: (i) the inclination of base-pairs does not change between A-tract and non-A-tract regions of helix; and (ii) the observed bends at GC/AT junctions are roll-wedge bends, not tilt-wedge as the junction bend model demands. (3) The non-A-tract bending model is consistent with both gel retardation data and with X-ray crystallography, and must be regarded as the only consistent model for A-tract bending.

Base Sequence↗

"...the tyranny of the lattice...".

A systematic comparison of crystal structures of nine different B-DNA dodecamers, in three different space groups, with and without A-tracts, shows that crystal packing or lattice forces are of secondary importance for helix axis bending, minor-groove width, and propeller twist. While other local helix parameters may be influenced or even established by crystal packing, the properties just enumerated are determined primarily by base sequence. One and the same crystal packing scheme can accommodate a bend in one of two different directions, or no bend at all. A-tract regions of B-DNA are inherently straight and unbent, with base-pair inclination no different from that of general-sequence B-DNA. Where bends are observed at junctions between G.C and A.T regions, they always involve a roll about base-pair long axes in a direction that compresses the wide major groove and, hence, are 90 degrees away from that necessary for the correctness of the junction model of A-tract bending. The G.C/A.T junction appears to be a flexible hinge, capable of adopting either a straight or a bent conformation under the local influence of weak crystal packing forces. Such forces therefore are a source of information about DNA deformability and not a curse to be deplored. But as an indication of the weakness of crystal packing forces, introduction of a single bromine atom in the major groove is sufficient to eliminate a bend, although brominated and unbrominated crystals are isomorphous.

Base Composition↗

Crystallographic analysis of C-C-A-A-G-C-T-T-G-G and its implications for bending in B-DNA.

Stacked B-DNA double helices of sequence C-C-A-A-G-C-T-T-G-G exhibit the same 23 degrees bend at -T-G-G C-C-A- across the nonbonded junction between helices that is observed in the middle of the decamer helix of sequence C-A-T-G-G-C-C-A-T-G, even though the space group (hexagonal vs orthorhombic), crystal packing, and connectedness at the center of the bent segment are quite different. An identical bend occurs across the interhelix junction of every monoclinic crystal structure of sequence C-C-A-x-x-x-x-T-G-G, suggesting that T-G-G-C-C-A constitutes a natural bending element in B-DNA. The bend occurs by rolling stacked base pairs about their long axes; there is no "tilt" component. Of the three possible models for A-tract bending--bent-A-tract, junction bends, or bent-non-A--which cannot be distinguished by solution measurements, all crystallographic evidence over the past 10 years unanimously supports the non-A regions as the actual bending loci.

Base Sequence↗

Crystal structure of CATGGCCATG and its implications for A-tract bending models.

The single-crystal x-ray analysis of orthorhombic CATGGCCATG has revealed a previously unrecognized mode of intrinsic bending in DNA. The decamer shows a smooth bend of 23 degrees over the central four base pairs, caused by preferential stacking interactions at guanine bases. The bend is produced by a roll of base pairs along their long axes, in a direction that compresses the wide major groove of the double helix. This major-groove-compressing bend at GGC, plus the abundant crystallographic evidence that runs of successive adenine bases (A-tracts) are straight and unbent, requires rethinking of the models most commonly invoked to explain A-tract bending. A decade of excellent experimental work involving gel migration experiments, cyclization kinetics, and nucleosome phasing has clearly established that introduction of short A-tracts into a general DNA sequence in synchrony with the natural repeat of the helix leads to bending. But it does not logically and inevitably follow that the actual bending is to be found within these introduced A-tracts or even at junctions with general-sequence B-DNA.

Base Composition↗

Automated docking in crystallography: analysis of the substrates of aconitase.

Automated docking of substrates to proteins of known structure aids the process of crystallographic analysis in two ways. First, automated docking can be used to generate a small number of starting models for substrates using only protein coordinates from an early stage of refinement. Second, automated docking provides a method for exploring aspects of catalysis that are inaccessible to crystallography by postulating binding modes of catalytic intermediates. This paper describes the use of automated docking to explore the binding of substrates to aconitase. The technique starts with a substrate molecule in an arbitrary configuration and position and finds favorable docked configurations in a (static) protein active site based on a molecular mechanics type force field. Using protein coordinates from an early stage of refinement of an aconitase-isocitrate complex, we successfully predicted the binding configuration of isocitrate. Four configurations were found, the energetically most favorable of which fit the observed electron density well and was used as a starting model for further refinement. Two configurations were found in citrate docking experiments, the second of which approximates the mode of substrate binding in an aconitase-nitrocitrate complex. We were also able to propose two binding modes of the catalytic intermediate cis-aconitate. These correspond closely to the isocitrate and the citrate binding modes. The relation of these new results to the proposed reaction mechanism is discussed.

Aconitate Hydratase↗

Soluble proteins: size, shape and function.

Why are proteins so big? Why do cells build oligomeric proteins? A visual survey of the protein structures available in the Protein Data Bank sheds new light on these questions.

Proteins↗

Molecular illustration in black and white.

Two-dimensional image processing techniques are used to create black and white illustrations from conventional z-buffer images. The illustrations appear to be composed of lines--outlines, contour lines, and hatched shading--but are calculated as grey-scale raster images. They have the advantage of faithful reproduction in publication and ease of display on workstations and laser printers. A FORTRAN source code is provided.

Computer Graphics↗