Induced fit DNA recognition by a minor groove binding analogue of Hoechst 33258: fluctuations in DNA A tract structure investigated by NMR and molecular dynamics simulations

NMR analysis and molecular dynamics simulations of d(GGTAATTACC)2 and its complex with a tetrahydropyrimidinium analogue of Hoechst 33258 suggest that DNA minor groove recognition in solution involves a combination of conformational selection and induced fit, rather than binding to a preorganised site. Analysis of structural fluctuations in the bound and unbound states suggests that the degree of induced fit observed is primarily a consequence of optimising van der Waals contacts with the walls of the minor groove resulting in groove narrowing through: (i) changes in base step parameters, including increased helical twist and propeller twist; (ii) changes to the sugar–phosphate backbone conformation to engulf the bound ligand; (iii) suppression of bending modes at the TpA steps. In contrast, the geometrical arrangement of hydrogen bond acceptors on the groove floor appears to be relatively insensitive to DNA conformation (helical twist and propeller twist). We suggest that effective recognition of DNA sequences (in this case an A tract structure) appears to depend to a significant extent on the sequence being flexible enough to be able to adopt the geometrically optimal conformation compatible with the various binding interactions, rather than involving ‘lock and key’ recognition.

Quantification of the hydrophobic interaction by simulations of the aggregation of small hydrophobic solutes in water

The hydrophobic interaction, the tendency for nonpolar molecules to aggregate in solution, is a major driving force in biology. In a direct approach to the physical basis of the hydrophobic effect, nanosecond molecular dynamics simulations were performed on increasing numbers of hydrocarbon solute molecules in water-filled boxes of different sizes. The intermittent formation of solute clusters gives a free energy that is proportional to the loss in exposed molecular surface area with a constant of proportionality of 45 ± 6 cal/mol⋅Å2. The molecular surface area is the envelope of the solute cluster that is impenetrable by solvent and is somewhat smaller than the more traditional solvent-accessible surface area, which is the area transcribed by the radius of a solvent molecule rolled over the surface of the cluster. When we apply a factor relating molecular surface area to solvent-accessible surface area, we obtain 24 cal/mol⋅Å2. Ours is the first direct calculation, to our knowledge, of the hydrophobic interaction from molecular dynamics simulations; the excellent qualitative and quantitative agreement with experiment proves that simple van der Waals interactions and atomic point-charge electrostatics account for the most important driving force in biology.

Electrostatic steering and ionic tethering in enzyme–ligand binding: Insights from simulations

To bind at an enzyme’s active site, a ligand must diffuse or be transported to the enzyme’s surface, and, if the binding site is buried, the ligand must diffuse through the protein to reach it. Although the driving force for ligand binding is often ascribed to the hydrophobic effect, electrostatic interactions also influence the binding process of both charged and nonpolar ligands. First, electrostatic steering of charged substrates into enzyme active sites is discussed. This is of particular relevance for diffusion-influenced enzymes. By comparing the results of Brownian dynamics simulations and electrostatic potential similarity analysis for triose-phosphate isomerases, superoxide dismutases, and β-lactamases from different species, we identify the conserved features responsible for the electrostatic substrate-steering fields. The conserved potentials are localized at the active sites and are the primary determinants of the bimolecular association rates. Then we focus on a more subtle effect, which we will refer to as “ionic tethering.” We explore, by means of molecular and Brownian dynamics simulations and electrostatic continuum calculations, how salt links can act as tethers between structural elements of an enzyme that undergo conformational change upon substrate binding, and thereby regulate or modulate substrate binding. This is illustrated for the lipase and cytochrome P450 enzymes. Ionic tethering can provide a control mechanism for substrate binding that is sensitive to the electrostatic properties of the enzyme’s surroundings even when the substrate is nonpolar.

Isolation and Characterization of Effector-Loop Mutants of CDC42 in Yeast

The highly conserved small GTPase Cdc42p is a key regulator of cell polarity and cytoskeletal organization in eukaryotic cells. Multiple effectors of Cdc42p have been identified, although it is unclear how their activities are coordinated to produce particular cell behaviors. One strategy used to address the contributions made by different effector pathways downstream of small GTPases has been the use of “effector-loop” mutants of the GTPase that selectively impair only a subset of effector pathways. We now report the generation and preliminary characterization of a set of effector-loop mutants of Saccharomyces cerevisiae CDC42. These mutants define genetically separable pathways influencing actin or septin organization. We have characterized the phenotypic defects of these mutants and the binding defects of the encoded proteins to known yeast Cdc42p effectors in vitro. The results suggest that these effectors cannot account for the observed phenotypes, and therefore that unknown effectors exist that affect both actin and septin organization. The availability of partial function alleles of CDC42 in a genetically tractable system serves as a useful starting point for genetic approaches to identify such novel effectors.

Timing of metamorphosis and the onset of the negative feedback loop between the thyroid gland and the pituitary is controlled by type II iodothyronine deiodinase in Xenopus laevis

Two important features of amphibian metamorphosis are the sequential response of tissues to different concentrations of thyroid hormone (TH) and the development of the negative feedback loop between the pituitary and the thyroid gland that regulates TH synthesis by the thyroid gland. At the climax of metamorphosis in Xenopus laevis (when the TH level is highest), the ratio of the circulating precursor thyroxine (T4) to the active form 3,5,3′-triiodothyronine (T3) in the blood is many times higher than it is in tissues. This difference is because of the conversion of T4 to T3 in target cells of the tadpole catalyzed by the enzyme type II iodothyronine deiodinase (D2) and the local effect (cell autonomy) of this activity. Limb buds and tails express D2 early and late in metamorphosis, respectively, correlating with the time that these organs undergo TH-induced change. T3 is required to complete metamorphosis because the peak concentration of T4 that is reached at metamorphic climax cannot induce the final morphological changes. At the climax of metamorphosis, D2 expression is activated specifically in the anterior pituitary cells that express the genes for thyroid-stimulating hormone but not in the cells that express proopiomelanocortin. Physiological concentrations of T3 but not T4 can suppress thyrotropin subunit β gene expression. The timing and the remarkable specificity of D2 expression in the thyrotrophs of the anterior pituitary coupled with the requirement for locally synthesized T3 strongly support a role for D2 in the onset of the negative feedback loop at the climax of metamorphosis.

Dependence of agonist activation on a conserved apolar residue in the third intracellular loop of the AT1 angiotensin receptor.

The coupling of agonist-activated seven transmembrane domain receptors to G proteins is known to involve the amino-terminal region of their third cytoplasmic loop. Analysis of the amino acids in this region of the rat type in angiotensin (AT1a) receptor identified Leu-222 as an essential residue in receptor activation by the physiological agonist, angiotensin II (Ang II). Nonpolar replacements for Leu-222 yielded functionally intact AT1 receptors, while polar or charged residues caused progressive impairment of Ang II-induced inositol phosphate generation. The decrease in agonist-induced signal generation was associated with a parallel reduction of receptor internalization, and was most pronounced for the Lys-222 mutant receptor. Although this mutant showed normal binding of the peptide antagonist, [Sar1,Ile6]Ang II, its affinity for Ang II was markedly reduced, consistent with its inability to adopt the high-affinity conformation. A search revealed that many Gq-coupled receptors contain an apolar amino acid (frequently leucine) in the position corresponding to Leu-222 of the AT1 receptor. These findings suggest that such a conserved apolar residue in the third intracellular loop is a crucial element in the agonist-induced activation of the AT1 and possibly many other G protein-coupled receptors.

Glass transition in DNA from molecular dynamics simulations.

Molecular dynamics simulations of the oligonucleotide duplex d(CGCGCG)2 in aqueous solution are used to investigate the glass transition phenomenon. The simulations were performed at temperatures in the 20 K to 340 K range. The mean square atomic fluctuations showed that the behavior of the oligonucleotide duplex was harmonic at low temperatures. A glass transition temperature at 223 K to 234 K was inferred for the oligonucleotide duplex, which is in agreement with experimental observations. The largest number of hydrogen bounds between the polar atoms of the oligonucleotide duplex and the water molecules was obtained at the glass transition temperature. With increasing temperature we observed a decrease in the average lifetime of the hydrogen bonds to water molecules.

The glycine binding site of the N-methyl-D-aspartate receptor subunit NR1: identification of novel determinants of co-agonist potentiation in the extracellular M3-M4 loop region.

The N-methyl-D-aspartate (NMDA) subtype of ionotropic glutamate receptors is a heterooligomeric membrane protein composed of homologous subunits. Here, the contribution of the M3-M4 loop of the NR1 subunit to the binding of glutamate and the co-agonist glycine was investigated by site-directed mutagenesis. Substitution of the phenylalanine residues at positions 735 or 736 of the M3-M4 loop produced a 15- to 30-fold reduction in apparent glycine affinity without affecting the binding of glutamate and the competitive glycine antagonist 7-chlorokynurenic acid; mutation of both residues caused a >100-fold decrease in glycine affinity. These residues are found in a C-terminal region of the M3-M4 loop that shows significant sequence similarity to bacterial amino acid-binding proteins. Epitope tagging revealed both the N-terminus and the M3-M4 loop to be exposed extracellularly, whereas a C-terminal epitope was localized intracellularly. These results indicate that the M3-M4 loop is part of the ligand-binding pocket of the NR1 subunit and provide the basis for a refined model of the glycine-binding site of the NMDA receptor.

A single residue in the M2-M3 loop is a major determinant of coupling between binding and gating in neuronal nicotinic receptors.

Binding of agonists to nicotinic acetylcholine receptors generates a sequence of changes that activate a cation-selective conductance. By measuring electrophysiological responses in chimeric alpha7/alpha3 receptors expressed in Xenopus oocytes, we have showed the involvement of the M2-M3 loop in coupling agonist binding to the channel gate. An aspartate residue therein, Asp-266 in the alpha7 subunit, was identified by site-directed mutagenesis as crucial, since mutants at this position exhibited very poor functional responses to three different nicotinic agonists. We have extended this investigation to another neuronal nicotinic receptor (alpha3/beta4), and found that a homologous residue in the beta4 subunit, Asp-268, played a similar role in coupling. These findings are consistent with a hypothesis that the aspartate residue in the M2-M3 loop, which is conserved in all homomer-forming alpha-type subunits and all neuronal beta-type subunits that combine to form functional receptors, is a major determinant of information transmission from binding site to channel gate in all neuronal nicotinic receptors.

Recombinational repair of gaps in DNA is asymmetric in Ustilago maydis and can be explained by a migrating D-loop model.

Recombinational repair of double-stranded DNA gaps was investigated in Ustilago maydis. The experimental system was designed for analysis of repair of an autonomously replicating plasmid containing a cloned gene disabled by an internal deletion. It was discovered that crossing over rarely accompanied gap repair. The strong bias against crossing over was observed in three different genes regardless of gap size. These results indicate that gap repair in U. maydis is unlikely to proceed by the mechanism envisioned in the double-stranded break repair model of recombination, which was developed to account for recombination in Saccharomyces cerevisiae. Experiments aimed at exploring processing of DNA ends were performed to gain understanding of the mechanism responsible for the observed bias. A heterologous insert placed within a gap in the coding sequence of two different marker genes strongly inhibited repair if the DNA was cleaved at the promoter-proximal junction joining the insert and coding sequence but had little effect on repair if the DNA was cleaved at the promoter-distal junction. Gene conversion of plasmid restriction fragment length polymorphism markers engineered in sequences flanking both sides of a gap accompanied repair but was directionally biased. These results are interpreted to mean that the DNA ends flanking a gap are subject to different types of processing. A model featuring a single migrating D-loop is proposed to explain the bias in gap repair outcome based on the observed asymmetry in processing the DNA ends.

Enhanced conformational sampling in Monte Carlo simulations of proteins: application to a constrained peptide.

A Monte Carlo simulation method for globular proteins, called extended-scaled-collective-variable (ESCV) Monte Carlo, is proposed. This method combines two Monte Carlo algorithms known as entropy-sampling and scaled-collective-variable algorithms. Entropy-sampling Monte Carlo is able to sample a large configurational space even in a disordered system that has a large number of potential barriers. In contrast, scaled-collective-variable Monte Carlo provides an efficient sampling for a system whose dynamics is highly cooperative. Because a globular protein is a disordered system whose dynamics is characterized by collective motions, a combination of these two algorithms could provide an optimal Monte Carlo simulation for a globular protein. As a test case, we have carried out an ESCV Monte Carlo simulation for a cell adhesive Arg-Gly-Asp-containing peptide, Lys-Arg-Cys-Arg-Gly-Asp-Cys-Met-Asp, and determined the conformational distribution at 300 K. The peptide contains a disulfide bridge between the two cysteine residues. This bond mimics the strong geometrical constraints that result from a protein's globular nature and give rise to highly cooperative dynamics. Computation results show that the ESCV Monte Carlo was not trapped at any local minimum and that the canonical distribution was correctly determined.

Mathematical analysis of activation thresholds in enzyme-catalyzed positive feedbacks: application to the feedbacks of blood coagulation.

A hierarchy of enzyme-catalyzed positive feedback loops is examined by mathematical and numerical analysis. Four systems are described, from the simplest, in which an enzyme catalyzes its own formation from an inactive precursor, to the most complex, in which two sequential feedback loops act in a cascade. In the latter we also examine the function of a long-range feedback, in which the final enzyme produced in the second loop activates the initial step in the first loop. When the enzymes generated are subject to inhibition or inactivation, all four systems exhibit threshold properties akin to excitable systems like neuron firing. For those that are amenable to mathematical analysis, expressions are derived that relate the excitation threshold to the kinetics of enzyme generation and inhibition and the initial conditions. For the most complex system, it was expedient to employ numerical simulation to demonstrate threshold behavior, and in this case long-range feedback was seen to have two distinct effects. At sufficiently high catalytic rates, this feedback is capable of exciting an otherwise subthreshold system. At lower catalytic rates, where the long-range feedback does not significantly affect the threshold, it nonetheless has a major effect in potentiating the response above the threshold. In particular, oscillatory behavior observed in simulations of sequential feedback loops is abolished when a long-range feedback is present.

Meso1, a basic-helix-loop-helix protein involved in mammalian presomitic mesoderm development.

To identify genes involved in the regulation of early mammalian development, we have developed a dominant-negative mutant basic-helix-loop-helix (bHLH) protein probe for interaction cloning and have isolated a member of the bHLH family of transcription factors, Meso1. Meso1-E2A heterodimers are capable of binding to oligonucleotide probes that contain a bHLH DNA recognition motif. In mouse embryos, Meso1 is expressed prior to MyoD1 family members. Meso1 expression is first detected at the neural plate stage of development in the paraxial mesoderm of the head and in presomitic mesodermal cells prior to their condensation into somites. Our findings suggest that Meso1 may be a key regulatory gene involved in the early events of vertebrate mesoderm differentiation.

Stable loop in the crystal structure of the intercalated four-stranded cytosine-rich metazoan telomere.

In most metazoans, the telomeric cytosine-rich strand repeating sequence is d(TAACCC). The crystal structure of this sequence was solved to 1.9-A resolution. Four strands associate via the cytosine-containing parts to form a four-stranded intercalated structure held together by C.C+ hydrogen bonds. The base-paired strands are parallel to each other, and the two duplexes are intercalated into each other in opposite orientations. One TAA end forms a highly stabilized loop with the 5' thymine Hoogsteen-base-paired to the third adenine. The 5' end of this loop is in close proximity to the 3' end of one of the other intercalated cytosine strands. Instead of being entirely in a DNA duplex, this structure suggests the possibility of an alternative conformation for the cytosine-rich telomere strands.

A sampling problem in molecular dynamics simulations of macromolecules.

Correlations in low-frequency atomic displacements predicted by molecular dynamics simulations on the order of 1 ns are undersampled for the time scales currently accessible by the technique. This is shown with three different representations of the fluctuations in a macromolecule: the reciprocal space of crystallography using diffuse x-ray scattering data, real three-dimensional Cartesian space using covariance matrices of the atomic displacements, and the 3N-dimensional configuration space of the protein using dimensionally reduced projections to visualize the extent to which phase space is sampled.