Electrostatics and hydration at the homeodomain-DNA interface: chemical probes of an interfacial water cavity

What drives proteins into the major or minor grooves of DNA?

The energetic profiles of a significant number of protein-DNA systems at 20 degrees C reveal that, despite comparable Gibbs free energies, association with the major groove is primarily an enthalpy-driven process, whereas binding to the minor groove is characterized by an unfavorable enthalpy that is compensated by favorable entropic contributions. These distinct energetic signatures for major versus minor groove binding are irrespective of the magnitude of DNA bending and/or the extent of binding-induced protein refolding. The primary determinants of their different energetic profiles appear to be the distinct hydration properties of the major and minor grooves; namely, that the water in the A+T-rich minor groove is in a highly ordered state and its removal results in a substantial positive contribution to the binding entropy. Since the entropic forces driving protein binding into the minor groove are a consequence of displacing water ordered by the regular arrangement of polar contacts, they cannot be regarded as hydrophobic.

Regulation of sexual dimorphism: mutational and chemogenetic analysis of the doublesex DM domain

Doublesex (dsx) is a transcription factor in Drosophila that regulates somatic sexual differentiation. Male- and female-specific splicing isoforms of DSX share a novel DNA-binding domain, designated the DM motif. Broadly conserved among metazoan sex-determining factors, the DM domain contains a nonclassical zinc module and binds in the DNA minor groove. Here, we characterize the DM motif by site-directed and random mutagenesis using a yeast one-hybrid (Y1H) system and extend this analysis by chemogenetic complementation in vitro. The Y1H system is based on a sex-specific Drosophila enhancer element and validated through studies of intersexual dsx mutations. We demonstrate that the eight motif-specific histidines and cysteines engaged in zinc coordination are each critical and cannot be interchanged; folding also requires conserved aliphatic side chains in the hydrophobic core. Mutations that impair DNA binding tend to occur at conserved positions, whereas neutral substitutions occur at nonconserved sites. Evidence for a specific salt bridge between a conserved lysine and the DNA backbone is obtained through the synthesis of nonstandard protein and DNA analogs. Together, these results provide molecular links between the structure of the DM domain and its function in the regulation of sexual dimorphism.

SRY and human sex determination: the basic tail of the HMG box functions as a kinetic clamp to augment DNA bending

Human testis-determining factor SRY contains a high-mobility-group (HMG) box, an alpha-helical DNA-binding domain that binds within an expanded minor groove to induce DNA bending. This motif is flanked on the C-terminal end by a basic tail, which functions both as a nuclear localization signal and accessory DNA-binding element. Whereas the HMG box is broadly conserved among otherwise unrelated transcription factors, tails differ in sequence and mode of DNA binding. Contrasting examples are provided by SRY and lymphoid enhancer factor 1 (LEF-1): whereas the SRY tail remains in the minor groove distal to the HMG box, the LEF-1 tail binds back across the center of the bent DNA site. The LEF-1 tail relieves electrostatic repulsion that would otherwise be incurred within the compressed major groove to enable sharp DNA bending with high affinity. Here, we demonstrate that the analogous SRY tail functions as a "kinetic clamp" to regulate the lifetime of the bent DNA complex. As in LEF-1, partial truncation of the distal SRY tail reduces specific DNA affinity and DNA bending, but these perturbations are modest: binding is reduced by only 1.8-fold, and bending by only 7-10 degrees . "Tailed" and truncated SRY complexes exhibit similar structures (as probed by NMR) and distributions of long-range conformational substates (as probed by time-resolved fluorescence resonance energy transfer). Surprisingly, however, the SRY tail retards dissociation of the protein-DNA complex by 20-fold. The marked and compensating changes in rates of association and dissociation observed on tail truncation, disproportionate to perturbations in affinity or structure, suggest that this accessory element functions as a kinetic clamp to regulate the lifetime of the SRY-DNA complex. We speculate that the kinetic stability of a bent DNA complex is critical to the assembly and maintenance of a sex-specific transcriptional pre-initiation complex.

A semi-implicit solvent model for the simulation of peptides and proteins

We present a new model of biomolecules hydration based on macroscopic electrostatic theory, that can both describe the microscopic details of solvent-solute interactions and allow for an efficient evaluation of the electrostatic hydration free energy. This semi-implicit model considers the solvent as an ensemble of polarizable pseudoparticles whose induced dipole describe both the electronic and orientational solvent polarization. In the presented version of the model, there is no mutual dipolar interaction between the particles, and they only interact through short-ranged Lennard-Jones interactions. The model has been integrated into a molecular dynamics code, and offers the possibility to simulate efficiently the conformational evolution of biomolecules. It is able to provide estimations of the electrostatic solvation free energy within short time windows during the simulation. It has been applied to the study of two small peptides, the octaalanine and the N-terminal helix of ribonuclease A, and two proteins, the bovine pancreatic trypsin inhibitor and the B1 immunoglobin-binding domain of streptococcal protein G. Molecular dynamics simulations of these biomolecules, using a slightly modified Amber force field, provide stable and meaningful trajectories in overall agreement with experiments and all-atom simulations. Correlations with respect to Poisson-Boltzmann electrostatic solvation free energies are also presented to discuss the parameterization of the model and its consequences.

Stepwise induced fit in the pico- to nanosecond time scale governs the complexation of the even-skipped transcriptional repressor homeodomain to DNA

Induced fit effects in the complex of a DNA decamer with two even-skipped transcriptional repressor homeodomain molecules were investigated by means of molecular dynamics simulations. Dynamics of these effects are found to be in the time scale from pico- to nanoseconds. First steps are made by the fast-moving DNA backbone phosphates, which upon binding change their B(I)/B(II) substate distribution. Further rearrangements in the DNA double helix induced upon complexation, like bending of the helix axis, changes of the minor groove width, and of different helical parameters, are slower and occur within a few nanoseconds. The flexibility of the DNA, especially of its backbone, seems thereby to play an important role for specific DNA ligand recognition.

Roles of hydration water molecules in molecular packing of the killer toxin from Pichia farinosa in its crystalline state investigated by cryogenic X-ray crystallography

The hydration structures around the killer toxin from Pichia farinosa were investigated by cryogenic X-ray crystallography. In particular, those contributing to the molecular association and the crystal contacts were analyzed with respect to the geometry and the networks of hydrogen bonds. The hydration water molecules attached on the surface so as to make up the surface shape in the contact complementary and mediated the intermolecular interactions through the networks of hydrogen bonds. Careful inspection of the contact area led to a proposal as to the molecular association mode of the toxin to determine the biological function in cells. In addition, the water-associated protein-protein interactions were approximated well by a simple theoretical equation on the solvation force expected in confined geometry. The present analysis may provide a way to analyze the crystal contact and molecular recognition in macromolecules in aqueous solution.

The tetrahymena ribozyme cleaves a 5'-methylene phosphonate monoester approximately 10(2)-fold faster than a normal phosphate diester: implications for enzyme catalysis of phosphoryl transfer reactions

Single-atom substrate modifications have revealed an intricate network of transition state interactions in the Tetrahymena ribozyme reaction. So far, these studies have targeted virtually every oxygen atom near the reaction center, except one, the 5'-bridging oxygen atom of the scissile phosphate. To address whether interactions with this atom play any role in catalysis, we used a new type of DNA substrate in which the 5'-oxygen is replaced with a methylene (-CH2-) unit. Under (kcat/Km)S conditions, the methylene phosphonate monoester substrate dCCCUCUT(mp)TA4 (where mp indicates the position of the phosphonate linkage) unexpectedly reacts approximately 10(3)-fold faster than the analogous control substrates lacking the -CH2- modification. Experiments with DNA-RNA chimeric substrates reveal that the -CH2- modification enhances docking of the substrates into the catalytic core of the ribozyme by approximately 10-fold and stimulates the chemical cleavage by approximately 10(2)-fold. The docking effect apparently arises from the ability of the -CH2- unit to suppress inherently deleterious effects caused by the thymidine residue that immediately follows the cleavage site. To analyze the -O- to -CH2- modification in the absence of this thymidine residue, we prepared oligonucleotide substrates containing methyl phosphate or ethyl phosphonate at the reaction center, thereby eliminating the 3'-terminal TA4 nucleotidyl group. In this context, the -O- to -CH2-modification has no effect on docking but retains the approximately 10(2)-fold effect on the chemical step. To investigate further the stimulatory influence on the chemical step, we measured the "intrinsic" effect of the -O- to -CH2- modification in nonenzymatic reactions with nucleophiles. We found that in solution, the -CH2- modification stimulates chemical reactivity of the phosphorus center by <5-fold, substantially lower in magnitude than the stimulatory effect in the catalytic core of the ribozyme. The greater stimulatory effect of the -CH2- modification in the active site compared to in solution may arise from fortuitous changes in molecular geometry that allow the ribozyme to accommodate the phosphonate transition state better than the natural phosphodiester transition state. As the -CH2- unit lacks lone pair electrons, its effectiveness in the ribozyme reaction suggests that the 5'-oxygen of the scissile phosphate plays no role in catalysis via hydrogen bonding or metal ion coordination. Finally, we show by analysis of physical organic data that such interactions in general provide little catalytic advantage to RNA and protein phosphoryl transferases because the 5'-oxygen undergoes only a small buildup of negative charge during the reaction. In addition to its mechanistic significance for the Tetrahymena ribozyme reaction and phosphoryl transfer reactions in general, this work suggests that phosphonate monoesters may provide a novel molecular tool for determining whether the chemical step limits the rate of an enzymatic reaction. As methylene phosphonate monoesters react modestly faster than phosphate diesters in model reactions, a similarly modest stimulatory effect on an enzymatic reaction upon -CH2- substitution would suggest rate-limiting chemistry.

Role of salt bridges in homeodomains investigated by structural analyses and molecular dynamics simulations

Homeodomains are a class of helix-turn-helix DNA-binding protein motifs that play an important role in the control of cellular development in eukaryotes. They fold in a three alpha-helix structural module, where the third helix is the recognition helix that fits into the major groove of DNA. Structural analysis of the members of the homeodomain family led to the identification of interactions likely to stabilize the protein domains. Linking the helices pairwise, three salt bridges were found to be well preserved within the family. Also well conserved were two cation-pi interactions between aromatic and positively charged side chains. To analyze the structural role of the salt bridges, molecular dynamics simulations (MD) were carried out on the wild-type homeodomain from the Drosophila paired protein (1fjl) and on three mutants, which lack one or two salt bridges and mimic natural mutations in other homeodomains. Analysis of the trajectories revealed only small structural rearrangements of the three helices in all MD simulations, thereby suggesting that the salt bridges have no essential stabilizing role at room temperature, but rather might be important for improving thermostability. The latter hypothesis is supported by a good correlation between the melting midpoint temperatures of several homeodomains and the number of salt bridges and cation-pi interactions that connect secondary structures.

Exploring the role of glutamine 50 in the homeodomain-DNA interface: crystal structure of engrailed (Gln50 --> ala) complex at 2.0 A

We have determined the crystal structure of a complex containing the engrailed homeodomain Gln50 --> Ala variant (QA50) bound to the wild-type optimal DNA site (TAATTA) at 2.0 A resolution. Biochemical and genetic studies by other groups have suggested that residue 50 is an important determinant of differential DNA-binding specificity among homeodomains (distinguishing among various sites of the general form TAATNN). However, biochemical studies of the QA50 variant had revealed that it binds almost as tightly as the wild-type protein and with only modest changes in specificity. We have now determined the crystal structure of the QA50 variant to help understand the role of residue 50 in site-specific recognition. Our cocrystal structure shows some interesting changes in the water structure at the site of the substitution and shows some changes in the conformations of neighboring side chains. However, the structure, like the QA50 biochemical data, suggests that Gln50 plays a relatively modest role in determining the affinity and specificity of the engrailed homeodomain.

Sexual dimorphism in diverse metazoans is regulated by a novel class of intertwined zinc fingers

Sex determination is regulated by diverse pathways. Although upstream signals vary, a cysteine-rich DNA-binding domain (the DM motif) is conserved within downstream transcription factors ofDrosophila melanogaster (Doublesex) and Caenorhabditis elegans (MAB-3). Vertebrate DM genes have likewise been identified and, remarkably, are associated with human sex reversal (46, XY gonadal dysgenesis). Here we demonstrate that the structure of the Doublesex domain contains a novel zinc module and disordered tail. The module consists of intertwined CCHC and HCCC Zn2+-binding sites; the tail functions as a nascent recognition α-helix. Mutations in either Zn2+-binding site or tail can lead to an intersex phenotype. The motif binds in the DNA minor groove without sharp DNA bending. These molecular features, unusual among zinc fingers and zinc modules, underlie the organization of a Drosophila enhancer that integrates sex- and tissue-specific signals. The structure provides a foundation for analysis of DM mutations affecting sexual dimorphism and courtship behavior.

Phosphorothioate substitution can substantially alter RNA conformation

Phosphorothioate substitution-interference experiments, routinely used to stereospecifically identify phosphoryl oxygen sites that participate in RNA-ligand binding and RNA-directed catalysis, rest in their interpretation on the untested assumption that substitution does not alter the conformation of the modified molecule from its biologically active state. Using NMR spectroscopy, we have tested this assumption by determining the structural effect of stereospecific phosphorothioate substitution at five positions in an RNA hairpin containing the binding site for bacteriophage MS2 capsid protein. At most sites, substitution has little or no effect, causing minor perturbations in the phosphate backbone and increasing the stacking among nucleotides in the hairpin loop. At one site, however, phosphorothioate substitution causes an unpaired adenine necessary for formation of the capsid protein-RNA complex to loop out of the RNA helix into the major groove. These results indicate that phosphorothioate substitution can substantially alter the conformation of RNA at positions of irregular secondary structure, complicating the use of substitution-interference experiments to study RNA structure and function.

Theory of biomolecular recognition

Specific, noncovalent binding of biomolecules can only be understood by considering structural, thermodynamic, and kinetic issues. The theoretical foundations for such analyses have been clarified in the past year. Computational techniques for both particle-based and continuum models continue to improve and to yield useful insights into an ever wider range of biomolecular systems.

Crystal structure of the yeast MATalpha2/MCM1/DNA ternary complex

The structure of a complex containing the homeodomain repressor protein MATalpha2 and the MADS-box transcription factor MCM1 bound to DNA has been determined by X-ray crystallography at 2.25 A resolution. It reveals the protein-protein interactions responsible for cooperative binding of MATalpha2 and MCM1 to DNA. The otherwise flexible amino-terminal extension of the MATalpha2 homeodomain forms a beta-hairpin that grips the MCM1 surface through parallel beta-strand hydrogen bonds and close-packed, predominantly hydrophobic, side chains. DNA bending induced by MCM1 brings the two proteins closer together, facilitating their interaction. An unusual feature of the complex is that an eight-amino-acid sequence adopts an alpha-helical conformation in one of two copies of the MATalpha2 monomer and a beta-strand conformation in the other. This 'chameleon' sequence of MATalpha2 may be important for recognizing natural operator sites.

Structure of the human NF-kappaB p52 homodimer-DNA complex at 2.1 A resolution

The crystal structure of human NF-kappaB p52 in its specific complex with the natural kappaB DNA binding site MHC H-2 has been solved at 2.1 A resolution. Whereas the overall structure resembles that of the NF-kappaB p50-DNA complex, pronounced differences are observed within the 'insert region'. This sequence segment differs in length between different Rel proteins. Compared with NF-kappaB p50, the compact alpha-helical insert region element is rotated away from the core of the N-terminal domain, opening up a mainly polar cleft. The insert region presents potential interaction surfaces to other proteins. The high resolution of the structure reveals many water molecules which mediate interactions in the protein-DNA interface. Additional complexity in Rel protein-DNA interaction comes from an extended interfacial water cavity that connects residues at the edge of the dimer interface to the central DNA bases. The observed water network might acount for differences in binding specificity between NF-kappaB p52 and NF-kappaB p50 homodimers.