Interaction between like-charged colloidal spheres in electrolyte solutions

How colloidal particles interact with each other is one of the key issues that determines our ability to interpret experimental results for phase transitions in colloidal dispersions and our ability to apply colloid science to various industrial processes. The long-accepted theories for answering this question have been challenged by results from recent experiments. Herein we show from Monte-Carlo simulations that there is a short-range attractive force between identical macroions in electrolyte solutions containing divalent counterions. Complementing some recent and related results by others, we present strong evidence of attraction between a pair of spherical macroions in the presence of added salt ions for the conditions where the interacting macroion pair is not affected by any other macroions that may be in the solution. This attractive force follows from the internal-energy contribution of counterion mediation. Contrary to conventional expectations, for charged macroions in an electrolyte solution, the entropic force is repulsive at most solution conditions because of localization of small ions in the vicinity of macroions. Both Derjaguin–Landau–Verwey–Overbeek theory and Sogami–Ise theory fail to describe the attractive interactions found in our simulations; the former predicts only repulsive interaction and the latter predicts a long-range attraction that is too weak and occurs at macroion separations that are too large. Our simulations provide fundamental “data” toward an improved theory for the potential of mean force as required for optimum design of new materials including those containing nanoparticles.

Intrinsic compressibility and volume compression in solvated proteins by molecular dynamics simulation at high pressure.

Constant pressure and temperature molecular dynamics techniques have been employed to investigate the changes in structure and volumes of two globular proteins, superoxide dismutase and lysozyme, under pressure. Compression (the relative changes in the proteins' volumes), computed with the Voronoi technique, is closely related with the so-called protein intrinsic compressibility, estimated by sound velocity measurements. In particular, compression computed with Voronoi volumes predicts, in agreement with experimental estimates, a negative bound water contribution to the apparent protein compression. While the use of van der Waals and molecular volumes underestimates the intrinsic compressibilities of proteins, Voronoi volumes produce results closer to experimental estimates. Remarkably, for two globular proteins of very different secondary structures, we compute identical (within statistical error) protein intrinsic compressions, as predicted by recent experimental studies. Changes in the protein interatomic distances under compression are also investigated. It is found that, on average, short distances compress less than longer ones. This nonuniform contraction underlines the peculiar nature of the structural changes due to pressure in contrast with temperature effects, which instead produce spatially uniform changes in proteins. The structural effects observed in the simulations at high pressure can explain protein compressibility measurements carried out by fluorimetric and hole burning techniques. Finally, the calculation of the proteins static structure factor shows significant shifts in the peaks at short wavenumber as pressure changes. These effects might provide an alternative way to obtain information concerning compressibilities of selected protein regions.

Direct interaction of Gβγ with a C-terminal Gβγ-binding domain of the Ca2+ channel α1 subunit is responsible for channel inhibition by G protein-coupled receptors

Several classes of voltage-gated Ca2+ channels (VGCCs) are inhibited by G proteins activated by receptors for neurotransmitters and neuromodulatory peptides. Evidence has accumulated to indicate that for non-L-type Ca2+ channels the executing arm of the activated G protein is its βγ dimer (Gβγ). We report below the existence of two Gβγ-binding sites on the A-, B-, and E-type α1 subunits that form non-L-type Ca2+ channels. One, reported previously, is in loop 1 connecting transmembrane domains I and II. The second is located approximately in the middle of the ca. 600-aa-long C-terminal tails. Both Gβγ-binding regions also bind the Ca2+ channel β subunit (CCβ), which, when overexpressed, interferes with inhibition by activated G proteins. Replacement in α1E of loop 1 with that of the G protein-insensitive and Gβγ-binding-negative loop 1 of α1C did not abolish inhibition by G proteins, but the exchange of the α1E C terminus with that of α1C did. This and properties of α1E C-terminal truncations indicated that the Gβγ-binding site mediating the inhibition of Ca2+ channel activity is the one in the C terminus. Binding of Gβγ to this site was inhibited by an α1-binding domain of CCβ, thus providing an explanation for the functional antagonism existing between CCβ and G protein inhibition. The data do not support proposals that Gβγ inhibits α1 function by interacting with the site located in the loop I–II linker. These results define the molecular mechanism by which presynaptic G protein-coupled receptors inhibit neurotransmission.

Folding and activity of circularly permuted forms of a polytopic membrane protein

The transmembrane subunit of the Glc transporter (IICBGlc), which mediates uptake and concomitant phosphorylation of glucose, spans the membrane eight times. Variants of IICBGlc with the native N and C termini joined and new N and C termini in the periplasmic and cytoplasmic surface loops were expressed in Escherichia coli. In vivo transport/in vitro phosphotransferase activities of the circularly permuted variants with the termini in the periplasmic loops 1 to 4 were 35/58, 32/37, 0/3, and 0/0% of wild type, respectively. The activities of the variants with the termini in the cytoplasmic loops 1 to 3 were 0/25, 0/4 and 24/70, respectively. Fusion of alkaline phosphatase to the periplasmic C termini stabilized membrane integration and increased uptake and/or phosphorylation activities. These results suggest that internal signal anchor and stop transfer sequences can function as N-terminal signal sequences in a circularly permuted α-helical bundle protein and that the orientation of transmembrane segments is determined by the amino acid sequence and not by the sequential appearance during translation. Of the four IICBGlc variants with new termini in periplasmic loops, only the one with the discontinuity in loop 4 is inactive. The sequences of loop 4 and of the adjacent TM7 and TM8 are conserved in all phosphoenolpyruvate-dependent carbohydrate:phosphotransferase system transporters of the glucose family.

Alternate protein frameworks for molecular recognition.

In an effort to determine whether proteins with structures other than the immunoglobulin fold can be used to mimic the ligand binding properties of antibodies, we generated a library from the four-helix bundle protein cytochrome b562 in which the two loops were randomized. Panning of this library against the bovine serum albumin (BSA) conjugate of N-methyl-p-nitrobenzylamine derivative 1 by phage display methods yielded cytochromes in which residues Trp-20, Arg-21, and Ser-22 in loop A and Arg-83 and Trp-84 in loop B were conserved. The individual mutants, which fold into native-like structure, bind selectively to the BSA-1 conjugate with micromolar dissociation constants (Kd), in comparison to a monoclonal antibody that binds selectively to 1 with a Kd of 290 nM. These and other antibody-like receptors may prove useful as therapeutic agents or as reagents for both intra- and extracellular studies.

The role of enzyme distortion in the single displacement mechanism of family 19 chitinases

By using molecular dynamics simulations, we have examined the binding of a hexaNAG substrate and two potential hydrolysis intermediates (an oxazoline ion and an oxocarbenium ion) to a family 19 barley chitinase. We find the hexaNAG substrate binds with all sugars in a chair conformation, unlike the family 18 chitinase which causes substrate distortion. Glu 67 is in a position to protonate the anomeric oxygen linking sugar residues D and E whereas Asn 199 serves to hydrogen bond with the C2′ N-acetyl group of sugar D, thus preventing the formation of an oxazoline ion intermediate. In addition, Glu 89 is part of a flexible loop region allowing a conformational change to occur within the active site to bring the oxocarbenium ion intermediate and Glu 89 closer by 4–5 Å. A hydrolysis product with inversion of the anomeric configuration occurs because of nucleophilic attack by a water molecule that is coordinated by Glu 89 and Ser 120. Issues important for the design of inhibitors specific to family 19 chitinases over family 18 chitinases also are discussed.

The ability of HIV type 1 to use CCR-3 as a coreceptor is controlled by envelope V1/V2 sequences acting in conjunction with a CCR-5 tropic V3 loop

Although infection by primary HIV type 1 (HIV-1) isolates normally requires the functional interaction of the viral envelope protein with both CD4 and the CCR-5 coreceptor, a subset of such isolates also are able to use the distinct CCR-3 receptor. By analyzing the ability of a series of wild-type and chimeric HIV-1 envelope proteins to mediate CCR-3-dependent infection, we have determined that CCR-3 tropism maps to the V1 and V2 variable region of envelope. Although substitution of the V1/V2 region of a CCR-3 tropic envelope into the context of a CCR-5 tropic envelope is both necessary and sufficient to confer CCR-3 tropism, this same substitution has no phenotypic effect when inserted into a CXCR-4 tropic HIV-1 envelope context. However, this latter chimera acquires both CCR-3 and CCR-5 tropism when a CCR-5 tropic V3 loop sequence also is introduced. These data demonstrate that the V1/2 region of envelope can, like the V3 loop region, encode a particular coreceptor requirement and suggest that a functional envelope:CCR-3 interaction may depend on the cooperative interaction of CCR-3 with both the V1/V2 and the V3 region of envelope.

AtGRP7, a nuclear RNA-binding protein as a component of a circadian-regulated negative feedback loop in Arabidopsis thaliana

The endogenous clock that drives circadian rhythms is thought to communicate temporal information within the cell via cycling downstream transcripts. A transcript encoding a glycine-rich RNA-binding protein, Atgrp7, in Arabidopsis thaliana undergoes circadian oscillations with peak levels in the evening. The AtGRP7 protein also cycles with a time delay so that Atgrp7 transcript levels decline when the AtGRP7 protein accumulates to high levels. After AtGRP7 protein concentration has fallen to trough levels, Atgrp7 transcript starts to reaccumulate. Overexpression of AtGRP7 in transgenic Arabidopsis plants severely depresses cycling of the endogenous Atgrp7 transcript. These data establish both transcript and protein as components of a negative feedback circuit capable of generating a stable oscillation. AtGRP7 overexpression also depresses the oscillation of the circadian-regulated transcript encoding the related RNA-binding protein AtGRP8 but does not affect the oscillation of transcripts such as cab or catalase mRNAs. We propose that the AtGRP7 autoregulatory loop represents a “slave” oscillator in Arabidopsis that receives temporal information from a central “master” oscillator, conserves the rhythmicity by negative feedback, and transduces it to the output pathway by regulating a subset of clock-controlled transcripts.

Multiple loop structures critical for ligand binding of the integrin α4 subunit in the upper face of the β-propeller mode 1

A non-I-domain integrin, α4β1, recognizes vascular cell adhesion molecule 1 (VCAM-1) and the IIICS portion of fibronectin. To localize regions of α4 critical for ligand binding, we swapped several predicted loops within or near the putative ligand-binding site of α4 (which spans repeats 2–5 of the seven N-terminal repeats) with the corresponding regions of α5. Swapping residues 112–131 in repeat 2, or residues 237–247 in repeat 4, completely blocked adhesion to immobilized VCAM-1 and connecting segment 1 (CS-1) peptide. However, swapping residues 40–52 in repeat 1, residues 151–164 in repeat 3, or residues 282–288 (which contain a putative cation binding motif) in repeat 5 did not affect or only slightly reduced adhesion to these ligands. The binding of several function-blocking antibodies is blocked by swapping residues 112–131, 151–164, and 186–191 (which contain previously identified residues critical for ligand binding, Tyr-187 and Gly-190). These results are consistent with the recently published β-propeller folding model of the integrin α4 subunit [Springer, T. A. (1997) Proc. Natl. Acad. Sci. USA 94, 65–72], in which seven four-stranded β-sheets are arranged in a torus around a pseudosymmetric axis. The regions of α4 critical for ligand binding are adjacent to each other and are located in the upper face, the predicted ligand-binding site, of the β-propeller model, although they are not adjacent in the primary structure.

Identification of specific Rp-phosphate oxygens in the tRNA anticodon loop required for ribosomal P-site binding

tRNA binding to the ribosomal P site is dependent not only on correct codon–anticodon interaction but also involves identification of structural elements of tRNA by the ribosome. By using a phosphorothioate substitution–interference approach, we identified specific nonbridging Rp-phosphate oxygens in the anticodon loop of tRNAPhe from Escherichia coli which are required for P-site binding. Stereo-specific involvement of phosphate oxygens at these positions was confirmed by using synthetic anticodon arm analogues at which single Rp- or Sp-phosphorothioates were incorporated. Identical interference results with yeast tRNAPhe and E. coli tRNAfMet indicate a common backbone conformation or common recognition elements in the anticodon loop of tRNAs. N-ethyl-N-nitrosourea modification–interference experiments with natural tRNAs point to the importance of the same phosphates in the loop. Guided by the crystal structure of tRNAPhe, we propose that specific Rp-phosphate oxygens are required for anticodon loop (“U-turn”) stabilization or are involved in interactions with the ribosome on correct tRNA–mRNA complex formation.

A conserved serine-rich stretch in the glutamate transporter family forms a substrate-sensitive reentrant loop

Neuronal and glial glutamate transporters remove the excitatory neurotransmitter glutamate from the synaptic cleft. The proteins belong to a large family of secondary transporters, which includes bacterial glutamate transporters. The C-terminal half of the glutamate transporters is well conserved and thought to contain the translocation path and the binding sites for substrate and coupling ions. A serine-rich sequence motif in this part of the proteins is located in a putative intracellular loop. Cysteine-scanning mutagenesis was applied to this loop in the glutamate transporter GltT of Bacillus stearothermophilus. The loop was found to be largely intracellular, but three consecutive positions in the conserved serine-rich motif (S269, S270, and E271) are accessible from both sides of the membrane. Single-cysteine mutants in the serine-rich motif were still capable of glutamate transport, but modification with N-ethylmaleimide blocked the transport activity in six mutants (T267C, A268C, S269C, S270C, E271C, and T272C). Two milimolars l-glutamate effectively protected against the modification of the cysteines at position 269–271 from the periplasmic side of the membrane but was unable to protect cysteine modification from the cytoplasmic side of the membrane. The results indicate that the conserved serine-rich motif in the glutamate transporter forms a reentrant loop, a structure that is found in several ion channels but is unusual for transporter proteins. The reentrant loop is of crucial importance for the function of the glutamate transporter.

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.