Vibrational population relaxation of carbon monoxide in the heme pocket of photolyzed carbonmonoxy myoglobin: Comparison of time-resolved mid-IR absorbance experiments and molecular dynamics simulations

The vibrational energy relaxation of carbon monoxide in the heme pocket of sperm whale myoglobin was studied by using molecular dynamics simulation and normal mode analysis methods. Molecular dynamics trajectories of solvated myoglobin were run at 300 K for both the δ- and ɛ-tautomers of the distal His-64. Vibrational population relaxation times of 335 ± 115 ps for the δ-tautomer and 640 ± 185 ps for the ɛ-tautomer were estimated by using the Landau–Teller model. Normal mode analysis was used to identify those protein residues that act as the primary “doorway” modes in the vibrational relaxation of the oscillator. Although the CO relaxation rates in both the ɛ- and δ-tautomers are similar in magnitude, the simulations predict that the vibrational relaxation of the CO is faster in the δ-tautomer with the distal His playing an important role in the energy relaxation mechanism. Time-resolved mid-IR absorbance measurements were performed on photolyzed carbonmonoxy hemoglobin (Hb13CO). From these measurements, a T1 time of 600 ± 150 ps was determined. The simulation and experimental estimates are compared and discussed.

Yeast chorismate mutase in the R state: Simulations of the active site

The isomerization of chorismate to prephenate by chorismate mutase in the biosynthetic pathway that forms Tyr and Phe involves C5—O (ether) bond cleavage and C1—C9 bond formation in a Claisen rearrangement. Development of negative charge on the ether oxygen, stabilized by Lys-168 and Glu-246, is inferred from the structure of a complex with a transition state analogue (TSA) and from the pH-rate profile of the enzyme and the E246Q mutant. These studies imply a protonated Glu-246 well above pH 7. Here, several 500-ps molecular dynamics simulations test the stability of enzyme–TSA complexes by using a solvated system with stochastic boundary conditions. The simulated systems are (i) protonated Glu-246 (stable), (ii) deprotonated Glu-246 (unstable), (iii) deprotonated Glu-246 plus one H2O between Glu-246 and the ether oxygen (unstable), (iv) the E246Q mutant (stable), and (v) addition of OH− between protonated Glu-246 and the ether oxygen. In (v), a local conformational change of Lys-168 displaced the OH− into the solvent region, suggesting a possible rate-determining step that precedes the catalytic step. In a 500-ps simulation of the enzyme complexed with the reactant chorismate or the product prephenate, no water molecule remained near the oxygen of the ligand. Calculations using the linearized Poisson–Boltzmann equation show that the effective pKa of Glu-246 is shifted from 5.8 to 8.1 as the negative charge on the ether oxygen of the TSA is changed from −0.56 electron to −0.9 electron. Altogether, these results support retention of a proton on Glu-246 to high pH and the absence of a water molecule in the catalytic steps.

The Stem-Loop Binding Protein (SLBP1) Is Present in Coiled Bodies of the Xenopus Germinal Vesicle

The stem-loop binding protein (SLBP1) binds the 3′ stem-loop of histone pre-mRNA and is required for efficient processing of histone transcripts in the nucleus. We examined the localization of SLBP1 in the germinal vesicle of Xenopus laevis oocytes. In spread preparations of germinal vesicle contents, an anti-SLBP1 antibody stained coiled bodies and specific chromosomal loci, including terminal granules, axial granules, and some loops. After injection of myc-tagged SLBP1 transcripts into the oocyte cytoplasm, newly translated myc-SLBP1 protein was detectable in coiled bodies within 4 h and in terminal and axial granules by 8 h. To identify the region(s) of SLBP1 necessary for subnuclear localization, we subcloned various parts of the SLBP1 cDNA and injected transcripts of these into the cytoplasm of oocytes. We determined that 113 amino acids at the carboxy terminus of SLBP1 are sufficient for coiled body localization and that disruption of a previously defined RNA-binding domain did not alter this localization. Coiled bodies also contain the U7 small nuclear ribonucleoprotein particle (snRNP), which participates in cleavage of the 3′ end of histone pre-mRNA. The colocalization of SLBP1 and the U7 snRNP in the coiled body suggests coordinated control of their functions, perhaps through a larger histone-processing particle. Some coiled bodies are attached to the lampbrush chromosomes at the histone gene loci, consistent with the view that coiled bodies in the oocyte recruit histone-processing factors to the sites of histone pre-mRNA transcription. The non-histone chromosomal sites at which SLBP1 is found include the genes coding for 5 S rRNA, U1 snRNA, and U2 snRNA, suggesting a wider role for SLBP1 in the biosynthesis of small non-spliced RNAs.

Oocytes are a source of catecholamines in the primate ovary: Evidence for a cell–cell regulatory loop

Catecholamines, thought to derive from the extrinsic innervation of the ovary, participate in the regulation of ovarian development and mature gonadal function. Recently, intraovarian neurons containing tyrosine hydroxylase (TH), the rate-limiting enzyme in catecholamine biosynthesis, were described in the ovary of nonhuman primates. We now show that the primate ovary expresses both the genes encoding TH and dopamine β-hydroxylase (DBH), the key enzymes in norepinephrine (NE) biosynthesis. Ovarian neurons were identified as a site of TH and DBH gene expression, and surprisingly, oocytes were identified as an exclusive site of DBH synthesis. Oocytes contain neither TH mRNA nor protein, indicating that they are unable to synthesize dopamine (DA). They did, however, express a DA transporter gene identical to that found in human brain. The physiological relevance of this transporter system and DBH in oocytes was indicated by the ability of isolated oocytes to metabolize exogenous DA into NE. Isolated follicles containing oocytes—but not those from which the oocytes had been removed—responded to DA with an elevation in cAMP levels; this elevation was prevented by propranolol, a β-adrenoreceptor antagonist. The results suggest that oocytes and somatic cells are linked by a neuroendocrine loop consisting of NE synthesized in oocytes from actively transported DA and cAMP produced by somatic follicular cells in response to NE-induced β-adrenoreceptor activation.

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.