BACKGROUND: Most soluble proteins are active as low-number oligomers. Statistical surveys of oligomeric proteins have defined the roles of hydrophobicity and complementarity in the stability of protein interfaces, but tend to average structural features over a diverse set of protein-protein interfaces, blurring information on how individual interfaces are stabilized. RESULTS: We report a visual survey of 136 homodimeric proteins from the Brookhaven Protein Data Bank, with images that highlight the major structural features of each protein-protein interaction surface. Nearly all of these proteins have interfaces formed between two globular subunits. Surprisingly, the pattern of hydrophilicity over the surface of these interfaces is quite variable. Approximately one-third of the interfaces show a recognizable hydrophobic core, with a single large, contiguous, hydrophobic patch surrounded by a ring of intersubunit polar interactions. The remaining two-thirds of the proteins show a varied mixture of small hydrophobic patches, polar interactions and water molecules scattered over the entire interfacial area. Ten proteins in the survey have intertwined interfaces formed by extensive interdigitation of the two subunit chains. These interfaces are very hydrophobic and are associated with proteins that require both stability and internal symmetry. CONCLUSIONS: The archetypal protein interface, with a defined hydrophobic core, is present in only a minority of the surveyed homodimeric proteins. Most homodimeric proteins are stabilized by a combination of small hydrophobic patches, polar interactions and a considerable number of bridging water molecules. The presence or absence of a hydrophobic core within these interfaces does not correlate with specific protein functions.

Adenylosuccinate Synthase↗

PubMed1998

The theoretical limits of DNA sequence discrimination by linked polyamides.

Linked polyamides bind in the minor groove of double-stranded DNA in a partially sequence-specific manner. This report analyzes the theoretical limits of DNA sequence discrimination by linked polyamides composed of two to four different types of heterocyclic rings, determining (i) the optimal choice of base-binding specificity for each ring and (ii) the optimal design for a polyamide composed of these rings to target a given DNA sequence and designed to maximize the fraction of the total polyamide binding to the specified target sequence relative to all other sequences. The results show that, fortuitously, polyamides composed of pyrrole, a naturally occurring G-excluding element, and imidazole, a rationally designed G-favoring element, have features similar to the theoretical optimum design for polyamides composed of two different rings. The results also show that, in polyamides composed of two or three types of heterocyclic rings, choosing a nonspecific "placeholder" ring, which binds equally strongly to each of the four bases, along with one or two base-specific rings will often enhance sequence specificity over a polyamide composed entirely of base-specific rings.

Amides↗

PubMed1998

Interactive modeling of supramolecular assemblies.

The modeling of supramolecular structure presents two major challenges: (1) managing the large amount of sequence, structural and biochemical data, and (2) presenting the data to the user in a flexible and comprehensible manner that addresses these problems. We describe a visualization environment for the creation and analysis of supramolecular models. A set of modular symmetry tools, collectively called SymGen, has been created, providing a flexible platform for the creation of complex assemblies, with interactive control of all symmetry elements and their parameters. A second tool, SymSearch, allows a range of parameters defined within SymGen to be sampled and the resulting conformations to be evaluated. The environment avoids information overload, caused by the large number of atoms in supramolecular complexes, by using a multiresolution spherical harmonic representation that allows the user to display only essential features. Spherical harmonics also enables control of the triangulation level, allowing the user to reduce the complexity of the geometric description to retain interactive speed. The visual fidelity of the surface data is retained by using texture maps that are independent of the resolution of the underlying triangulation. We describe the design and implementation of this environment, and three illustrative examples of its utility.

Capsid↗

PubMed1998

Automated docking and the search for HIV protease inhibitors.

This article will discuss the motivations, technologies, and future directions of computational automated docking in the context of the structure-based rational design of HIV-1 protease inhibitors. Docking simulations are widely used for screening of compound libraries to identify new drug leads, employing a simple model for rapid testing of thousands of compounds. Docking simulations are also useful for lead enhancement, using more detailed models to analyze the atomic interactions between inhibitors and target macromolecules. Major advances have been reported in the development of empirical force fields, which now allow assessment of relative binding strength and drug specificity, and extensions of automated docking techniques allow de novo drug design.

Computer Simulation↗

PubMed1998

An analysis of a class of DNA sequence reading molecules.

Linked polyamides are a class of designed molecules that bind in the minor groove of double-stranded DNA in a partially sequence-specific manner but have limited sequence discriminatory abilities. This suggests a need for design alternatives to create molecules with enhanced sequence specificity. In this report we present formal proofs of the theoretical limits of the DNA sequence specificity of hypothetical sequence reading molecules as a function of their base recognition properties and sequence content and length of their target sequence. We prove that molecules containing nonspecific readers at critical positions within the molecule may have enhanced sequence specificity over molecules composed entirely of base specific reading elements. We also determine optimal patterns of base recognition for molecules in order to optimize their target sequence specificity. We also examine the effect of the length of a polyamide (i.e., the number of base pairs it binds) on its sequence discriminatory ability and determine necessary concentration dependent constraints on the binding free energies in order for longer polyamides to have greater sequence specificity than shorter ones. We show that unless the discriminatory ability of a ring for its preferred base is very strong, longer polyamides do not necessarily have greater sequence specificity over shorter ones when compared at the same molar concentration.

Base Sequence↗

PubMed1998

Computational coevolution of antiviral drug resistance.

An understanding of antiviral drug resistance is important in the design of effective drugs. Comprehensive features of the interaction between drug designs and resistance mutations are difficult to study experimentally because of the very large numbers of drugs and mutants involved. We describe a computational framework for studying antiviral drug resistance. Data on HIV-1 protease are used to derive an approximate model that predicts interaction of a wide range of mutant forms of the protease with a broad class of protease inhibitors. An algorithm based on competitive coevolution is used to find highly resistant mutant forms of the protease, and effective inhibitors against such mutants, in the context of the model. We use this method to characterize general features of inhibitors that are effective in overcoming resistance, and to study related issues of selection pathways, cross-resistance, and combination therapies.

Algorithms↗

PubMed1997

Defining GC-specificity in the minor groove: side-by-side binding of the di-imidazole lexitropsin to C-A-T-G-G-C-C-A-T-G.

BACKGROUND: Polyamide drugs, such as netropsin, distamycin and their lexitropsin derivatives, can be inserted into a narrow B-DNA minor groove to form 1:1 complexes that can distinguish AT base pairs from GC, but cannot detect end-for-end base-pair reversals such as TA for AT. In contrast, 2:1 side-by-side polyamide drug complexes potentially are capable of such discrimination. Imidazole (Im) and pyrrole (Py) rings side-by-side read a GC base pair with the Im ring recognizing the guanine side. But the reason for this specific G-Im association is unclear because the guanine NH2 group sits in the center of the groove. A 2:1 drug:DNA complex that presents Im at both ends of a GC base pair should help unscramble the issue of imidazole reading specificity. RESULTS: We have determined the crystal structure of a 2:1 complex of a di-imidazole lexitropsin (DIM), an analogue of distamycin, and a DNA decamer with the sequence C-A-T-G-G-C-C-A-T-G. The two DIM molecules sit antiparallel to one another in a broad minor groove, with their cationic tails widely separated. Im rings of one drug molecule stack against amide groups of the other. DIM1 rests against nucleotides C7A8T9G10 of strand 1 of the helix, whereas DIM2 rests against G14G15C16C17 on strand 2. All DIM amide nitrogens donate hydrogen bonds to N and O atoms on the floor of the DNA groove and, in addition, the two Im rings on DIM2 accept hydrogen bonds from guanine N2 amines, thereby providing specific reading. The guanine N2 amine can bond to Im on its own side of the groove, but not on the cytosine side, because of limits on close approach of the two Im rings and the geometry of sp2 hybridization about the amide nitrogen. CONCLUSIONS: Im and Py rings distinguish AT from GC base pairs because of steric factors involving the bulk of the guanine amine, and the ability of Im to form a hydrogen bond with the amine. Side-by-side Im and Py rings differentiate GC from CG base pairs because of tight steric contacts and sp2 hybridization at the amine nitrogen atom, with the favored conformations being G/Im,Py/C and C/Py,Im/G. Discrimination between AT and TA base pairs may be possible using bulkier rings, such as thiazole to select the A end of the base pair.

Binding Sites↗

PubMed1997

Estimation of the DNA sequence discriminatory ability of hairpin-linked lexitropsins.

Three- and four-ring polyamides containing N-methylimidazole and N-methylpyrrole, and their hairpin-linked derivatives, bind side-by-side in the minor groove of DNA in a sequence-specific manner. The sequences recognized by side-by-side molecules are dependent on the pairings of the polyamide rings to the bases. In this study we report a mathematical model for estimating the free energies of binding for gamma-aminobutyric acid-linked polyamides to 5- and 6-bp DNA sequences. The model parameters are calibrated by a least-squares fit to 35 experimental binding constants. The model performs well in cross-validation experiments and the parameters are consistent with previously proposed empirical rules of polyamide-DNA binding. We apply the model to the design of targeted polyamides, evaluating the ability of the proposed polyamides to bind to a DNA sequence of interest while minimizing binding to the remaining DNA sequences.

Antineoplastic Agents↗

PubMed1997

Progress in the design of DNA sequence-specific lexitropsins.

Sequence-specific polyamides that bind in the minor groove of DNA are attractive candidates for antibiotics, cancer chemotherapeutics, and transcriptional antagonists. This paper reviews the progress of structure-based design of minor-groove-binding polyamides, from the first structure of netropsin with DNA, to the effective linked polyamides currently under study. A theory of polyamide specificity is also reviewed, introducing methods to determine the optimal strategies for targeting a given DNA sequence within a genome of competing sequences.

Antineoplastic Agents↗