p53-induced DNA bending and twisting: p53 tetramer binds on the outer side of a DNA loop and increases DNA twisting

DNA binding activity of p53 is crucial for its tumor suppressor function. Our recent studies have shown that four molecules of the DNA binding domain of human p53 (p53DBD) bind the response elements with high cooperativity and bend the DNA. By using A-tract phasing experiments, we find significant differences between the bending and twisting of DNA by p53DBD and by full-length human wild-type (wt) p53. Our data show that four subunits of p53DBD bend the DNA by 32–36°, whereas wt p53 bends it by 51–57°. The directionality of bending is consistent with major groove bends at the two pentamer junctions in the consensus DNA response element. More sophisticated phasing analyses also demonstrate that p53DBD and wt p53 overtwist the DNA response element by ≈35° and ≈70°, respectively. These results are in accord with molecular modeling studies of the tetrameric complex. Within the constraints imposed by the protein subunits, the DNA can assume a range of conformations resulting from correlated changes in bend and twist angles such that the p53–DNA tetrameric complex is stabilized by DNA overtwisting and bending toward the major groove at the CATG tetramers. This bending is consistent with the inherent sequence-dependent anisotropy of the duplex. Overall, the four p53 moieties are placed laterally in a staggered array on the external side of the DNA loop and have numerous interprotein interactions that increase the stability and cooperativity of binding. The novel architecture of the p53 tetrameric complex has important functional implications including possible p53 interactions with chromatin.

Energy-based de novo protein folding by conformational space annealing and an off-lattice united-residue force field: Application to the 10-55 fragment of staphylococcal protein A and to apo calbindin D9K

The conformational space annealing (CSA) method for global optimization has been applied to the 10-55 fragment of the B-domain of staphylococcal protein A (protein A) and to a 75-residue protein, apo calbindin D9K (PDB ID code 1CLB), by using the UNRES off-lattice united-residue force field. Although the potential was not calibrated with these two proteins, the native-like structures were found among the low-energy conformations, without the use of threading or secondary-structure predictions. This is because the CSA method can find many distinct families of low-energy conformations. Starting from random conformations, the CSA method found that there are two families of low-energy conformations for each of the two proteins, the native-like fold and its mirror image. The CSA method converged to the same low-energy folds in all cases studied, as opposed to other optimization methods. It appears that the CSA method with the UNRES force field, which is based on the thermodynamic hypothesis, can be used in prediction of protein structures in real time.

Recognition of nonhybridizing base pairs during nucleotide excision repair of DNA

Nondistorting C4′ backbone adducts serve as molecular tools to analyze the strategy by which a limited number of human nucleotide excision repair (NER) factors recognize an infinite variety of DNA lesions. We have constructed composite DNA substrates containing a noncomplementary site adjacent to a nondistorting C4′ adduct to show that the loss of hydrogen bonding contacts between partner strands is an essential signal for the recruitment of NER enzymes. This specific conformational requirement for excision is mediated by the affinity of xeroderma pigmentosum group A (XPA) protein for nonhybridizing sites in duplex DNA. XPA recognizes defective Watson–Crick base pair conformations even in the absence of DNA adducts or other covalent modifications, apparently through detection of hydrophobic base components that are abnormally exposed to the double helical surface. This recognition function of XPA is enhanced by replication protein A (RPA) such that, in combination, XPA and RPA constitute a potent molecular sensor of denatured base pairs. Our results indicate that the XPA–RPA complex may promote damage recognition by monitoring Watson–Crick base pair integrity, thereby recruiting the human NER system preferentially to sites where hybridization between complementary strands is weakened or entirely disrupted.

NMR investigation of ferricytochrome c unfolding: Detection of an equilibrium unfolding intermediate and residual structure in the denatured state

Horse ferricytochrome c (cyt c) undergoes exchange of one of its axial heme ligands (Met-80) for one or more non-native ligands under denaturing conditions. We have used 1H NMR spectroscopy to detect two conformations of paramagnetic cyt c with non-native heme ligation through a range of urea concentrations. One non-native form is an equilibrium unfolding intermediate observed under partially denaturing conditions and is attributed to replacement of Met-80 with one or more Lys side chains. The second non-native form, in which the native Met ligand is replaced by a His, is observed under strongly denaturing conditions. Thermodynamic analysis of these data indicates a relatively small ΔG (17 kJ/mol) for the transition from native to the Lys-ligated intermediate and a significantly larger ΔG (47 kJ/mol) for the transition from native to the His-ligated species. Although CD and fluorescence data indicate that the equilibrium unfolding of cyt c is a two-state process, these NMR results implicate an intermediate with His-Lys ligation.

NMR structure of the bovine prion protein

The NMR structures of the recombinant 217-residue polypeptide chain of the mature bovine prion protein, bPrP(23–230), and a C-terminal fragment, bPrP(121–230), include a globular domain extending from residue 125 to residue 227, a short flexible chain end of residues 228–230, and an N-terminal flexibly disordered “tail” comprising 108 residues for the intact protein and 4 residues for bPrP(121–230), respectively. The globular domain contains three α-helices comprising the residues 144–154, 173–194, and 200–226, and a short antiparallel β-sheet comprising the residues 128–131 and 161–164. The best-defined parts of the globular domain are the central portions of the helices 2 and 3, which are linked by the only disulfide bond in bPrP. Significantly increased disorder and mobility is observed for helix 1, the loop 166–172 leading from the β-strand 2 to helix 2, the end of helix 2 and the following loop, and the last turn of helix 3. Although there are characteristic local differences relative to the conformations of the murine and Syrian hamster prion proteins, the bPrP structure is essentially identical to that of the human prion protein. On the other hand, there are differences between bovine and human PrP in the surface distribution of electrostatic charges, which then appears to be the principal structural feature of the “healthy” PrP form that might affect the stringency of the species barrier for transmission of prion diseases between humans and cattle.

Molecular basis for CD40 signaling mediated by TRAF3

Tumor necrosis factor receptors (TNFR) are single transmembrane-spanning glycoproteins that bind cytokines and trigger multiple signal transduction pathways. Many of these TNFRs rely on interactions with TRAF proteins that bind to the intracellular domain of the receptors. CD40 is a member of the TNFR family that binds to several different TRAF proteins. We have determined the crystal structure of a 20-residue fragment from the cytoplasmic domain of CD40 in complex with the TRAF domain of TRAF3. The CD40 fragment binds as a hairpin loop across the surface of the TRAF domain. Residues shown by mutagenesis and deletion analysis to be critical for TRAF3 binding are involved either in direct contact with TRAF3 or in intramolecular interactions that stabilize the hairpin. Comparison of the interactions of CD40 with TRAF3 vs. TRAF2 suggests that CD40 may assume different conformations when bound to different TRAF family members. This molecular adaptation may influence binding affinity and specific cellular triggers.

Position-specific trapping of topoisomerase I–DNA cleavage complexes by intercalated benzo[a]- pyrene diol epoxide adducts at the 6-amino group of adenine

DNA topoisomerase I (top1) is the target of potent anticancer agents, including camptothecins and DNA intercalators, which reversibly stabilize (trap) top1 catalytic intermediates (cleavage complexes). The aim of the present study was to define the structural relationship between the site(s) of covalently bound intercalating agents, whose solution conformations in DNA are known, and the site(s) of top1 cleavage. Two diastereomeric pairs of oligonucleotide 22-mers, derived from a sequence used to determine the crystal structure of top1–DNA complexes, were synthesized. One pair contained either a trans-opened 10R- or 10S-benzo[a]pyrene 7,8-diol 9,10-epoxide adduct at the N6-amino group of a central 2′-deoxyadenosine residue in the scissile strand, and the other pair contained the same two adducts in the nonscissile strand. These adducts were derived from the (+)-(7R,8S,9S,10R)- and (−)-(7S,8R,9R,10S)-7,8-diol 9,10-epoxides in which the benzylic 7-hydroxyl group and the epoxide oxygen are trans. On the basis of analogy with known solution conformations of duplex oligonucleotides containing these adducts, we conclude that top1 cleavage complexes are trapped when the hydrocarbon adduct is intercalated between the base pairs flanking a preexisting top1 cleavage site, or between the base pairs immediately downstream (3′ relative to the scissile strand) from this site. We propose a model with the +1 base rotated out of the duplex, and in which the intercalated adduct prevents religation of the corresponding nucleotide at the 5′ end of the cleaved DNA. These results suggest mechanisms whereby intercalating agents interfere with the normal function of human top1.

Energetics of the interaction between water and the helical peptide group and its role in determining helix propensities

The alanine helix provides a model system for studying the energetics of interaction between water and the helical peptide group, a possible major factor in the energetics of protein folding. Helix formation is enthalpy-driven (−1.0 kcal/mol per residue). Experimental transfer data (vapor phase to aqueous) for amides give the enthalpy of interaction with water of the amide group as ≈−11.5 kcal/mol. The enthalpy of the helical peptide hydrogen bond, computed for the gas phase by quantum mechanics, is −4.9 kcal/mol. These numbers give an enthalpy deficit for helix formation of −7.6 kcal/mol. To study this problem, we calculate the electrostatic solvation free energy (ESF) of the peptide groups in the helical and β-strand conformations, by using the delphi program and parse parameter set. Experimental data show that the ESF values of amides are almost entirely enthalpic. Two key results are: in the β-strand conformation, the ESF value of an interior alanine peptide group is −7.9 kcal/mol, substantially less than that of N-methylacetamide (−12.2 kcal/mol), and the helical peptide group is solvated with an ESF of −2.5 kcal/mol. These results reduce the enthalpy deficit to −1.5 kcal/mol, and desolvation of peptide groups through partial burial in the random coil may account for the remainder. Mutant peptides in the helical conformation show ESF differences among nonpolar amino acids that are comparable to observed helix propensity differences, but the ESF differences in the random coil conformation still must be subtracted.

Molecular determinants of tissue selectivity in estrogen receptor modulators

Interaction of the estrogen receptor/ligand complex with a DNA estrogen response element is known to regulate gene transcription. In turn, specific conformations of the receptor-ligand complex have been postulated to influence unique subsets of estrogen-responsive genes resulting in differential modulation and, ultimately, tissue-selective outcomes. The estrogen receptor ligands raloxifene and tamoxifen have demonstrated such tissue-specific estrogen agonist/antagonist effects. Both agents antagonize the effects of estrogen on mammary tissue while mimicking the actions of estrogen on bone. However, tamoxifen induces significant stimulation of uterine tissue whereas raloxifene does not. We postulate that structural differences between raloxifene and tamoxifen may influence the conformations of their respective receptor/ligand complexes, thereby affecting which estrogen-responsive genes are modulated in various tissues. These structural differences are 4-fold: (A) the presence of phenolic hydroxyls, (B) different substituents on the basic amine, (C) incorporation of the stilbene moiety into a cyclic benzothiophene framework, and (D) the imposition of a carbonyl “hinge” between the basic amine-containing side chain and the olefin. A series of raloxifene analogs that separately exemplify each of these differences have been prepared and evaluated in a series of in vitro and in vivo assays. This strategy has resulted in the development of a pharmacophore model that attributes the differences in effects on the uterus between raloxifene and tamoxifen to a low-energy conformational preference imparting an orthogonal orientation of the basic side chain with respect to the stilbene plane. This three-dimensional array is dictated by a single carbon atom in the hinge region of raloxifene. These data indicate that differences in tissue selective actions among benzothiophene and triarylethylene estrogen receptor modulators can be ascribed to discrete ligand conformations.

Forced unfolding modulated by disulfide bonds in the Ig domains of a cell adhesion molecule

Cell adhesion molecules (CAMs) mediate cell attachment and stress transfer through extracellular domains. Here we forcibly unfold the Ig domains of a prototypical Ig superfamily CAM that contains intradomain disulfide bonds. The Ig domains of all such CAMs have conformations homologous to cadherin extracellular domains, titin Ig-type domains, and fibronectin type-III (FNIII) domains. Atomic force microscopy has been used to extend the five Ig domains of Mel-CAM (melanoma CAM)—a protein that is overexpressed in metastatic melanomas—under conditions where the disulfide bonds were either left intact or disrupted through reduction. Under physiological conditions where intradomain disulfide bonds are intact, partial unfolding was observed at forces far smaller than those reported previously for either titin's Ig-type domains or tenascin's FNIII domains. This partial unfolding under low force may be an important mechanism for imparting elasticity to cell–cell contacts, as well as a regulatory mechanism for adhesive interactions. Under reducing conditions, Mel-CAM's Ig domains were found to fully unfold through a partially folded state and at slightly higher forces. The results suggest that, in divergent evolution of all such domains, stabilization imparted by disulfide bonds relaxes requirements for strong, noncovalent, folded-state interactions.

RegulonDB (version 3.2): transcriptional regulation and operon organization in Escherichia coli K-12

RegulonDB is a database on mechanisms of transcription regulation and operon organization in Escherichia coli K-12. The current version has considerably increased numbers of regulatory elements such as promoters, binding sites and terminators. The complete repertoire of known and predicted DNA-binding transcriptional regulators can be considered to be included in this version. The database now distinguishes different allosteric conformations of regulatory proteins indicating the one active in binding and regulating the different promoters. A new set of operon predictions has been incorporated. The relational design has been modified accordingly. Furthermore, a major improvement is a graphic display enabling browsing of the database with a Java-based graphic user interface with three zoom-levels connected to properties of each chromo­somal element. The purpose of these modifications is to make RegulonDB a useful tool and control set for tran­scriptome experiments. RegulonDB can be accessed on the web at the URL: http://www.cifn.unam.mx/Computational_Biology/regulondb/

An isolated, surface-expressed I domain of the integrin αLβ2 is sufficient for strong adhesive function when locked in the open conformation with a disulfide bond

We introduced disulfide bonds to lock the integrin αLβ2 I domain in predicted open, ligand binding or closed, nonbinding conformations. Transfectants expressing αLβ2 heterodimers containing locked-open but not locked-closed or wild-type I domains constitutively adhered to intercellular adhesion molecule-1 (ICAM-1) substrates. Locking the I domain closed abolished constitutive and activatable adhesion. The isolated locked-open I domain bound as well as the activated αLβ2 heterodimer, and binding was abolished by reduction of the disulfide. Lovastatin, which binds under the conformationally mobile C-terminal α-helix of the I domain, inhibited binding to ICAM-1 by αLβ2 with wild-type, but not locked-open I domains. These data establish the importance of conformational change in the αL I domain for adhesive function and show that this domain is sufficient for full adhesive activity.

Comparing protein–ligand interactions in solution and single crystals by Raman spectroscopy

By using a Raman microscope, we show that it is possible to probe the conformational states in protein crystals and crystal fragments under growth conditions (in hanging drops). The flavin cofactor in the enzyme para-hydroxybenzoate hydroxylase can assume two conformations: buried in the protein matrix (“in”) or essentially solvent-exposed (“out”). By using Raman difference spectroscopy, we previously have identified characteristic flavin marker bands for the in and out conformers in the solution phase. Now we show that the flavin Raman bands can be used to probe these conformational states in crystals, permitting a comparison between solution and crystal environments. The in or out marker bands are similar for the respective conformers in the crystal and in solution; however, significant differences do exist, showing that the environments for the flavin's isoalloxazine ring are not identical in the two phases. Moreover, the Raman-band widths of the flavin modes are narrower for both in and out conformers in the crystals, indicating that the flavin exists in a more limited range of closely related conformational states in the crystal than in solution. In general, the ability to compare detailed Raman data for complexes in crystals and solution provides a means of bridging crystallographic and solution studies.

Herpes simplex virus DNA packaging sequences adopt novel structures that are specifically recognized by a component of the cleavage and packaging machinery

The product of the herpes simplex virus type 1 UL28 gene is essential for cleavage of concatemeric viral DNA into genome-length units and packaging of this DNA into viral procapsids. To address the role of UL28 in this process, purified UL28 protein was assayed for the ability to recognize conserved herpesvirus DNA packaging sequences. We report that DNA fragments containing the pac1 DNA packaging motif can be induced by heat treatment to adopt novel DNA conformations that migrate faster than the corresponding duplex in nondenaturing gels. Surprisingly, these novel DNA structures are high-affinity substrates for UL28 protein binding, whereas double-stranded DNA of identical sequence composition is not recognized by UL28 protein. We demonstrate that only one strand of the pac1 motif is responsible for the formation of novel DNA structures that are bound tightly and specifically by UL28 protein. To determine the relevance of the observed UL28 protein–pac1 interaction to the cleavage and packaging process, we have analyzed the binding affinity of UL28 protein for pac1 mutants previously shown to be deficient in cleavage and packaging in vivo. Each of the pac1 mutants exhibited a decrease in DNA binding by UL28 protein that correlated directly with the reported reduction in cleavage and packaging efficiency, thereby supporting a role for the UL28 protein–pac1 interaction in vivo. These data therefore suggest that the formation of novel DNA structures by the pac1 motif confers added specificity on recognition of DNA packaging sequences by the UL28-encoded component of the herpesvirus cleavage and packaging machinery.

Point mutations and sequence variability in proteins: Redistributions of preexisting populations

Here we study the effect of point mutations in proteins on the redistributions of the conformational substates. We show that regardless of the location of a mutation in the protein structure and of its type, the observed movements of the backbone recur largely at the same positions in the structures. Despite the different interactions that are disrupted and formed by the residue substitution, not only are the conformations very similar, but the regions that move are also the same, regardless of their sequential or spatial distance from the mutation. This observation leads us to conclude that, apart from some extreme cases, the details of the interactions are not critically important in determining the protein conformation or in specifying which parts of the protein would be more prone to take on different local conformations in response to changes in the sequence. This finding further illustrates why proteins manifest a robustness toward many mutational events. This nonuniform distribution of the conformer population is consistently observed in a variety of protein structural types. Topology is critically important in determining folding pathways, kinetics, building block cutting, and anatomy trees. Here we show that topology is also very important in determining which regions of the protein structure will respond to sequence changes, regardless of the sequential or spatial location of the mutation.