Thermal motions of the E. coli glucose-galactose binding protein studied using well-sampled, semi-atomistic simulations.

The E. coli glucose-galactose chemosensory receptor is a 309 residue, 32 kDa protein consisting of two distinct structural domains. We used two computational methods to examine the protein's thermal fluctuations, including both the large-scale interdomain movements that contribute to the receptor's mechanism of action, as well as smaller-scale motions. We primarily employ extremely fast, "semi-atomistic" Library-Based Monte Carlo (LBMC) simulations, which include all backbone atoms but "implicit" side chains. Our results were compared with previous experiments and all-atom molecular dynamics (MD) simulation. Both LBMC and MD simulations were performed using both the apo and glucose-bound form of the protein, with LBMC exhibiting significantly larger fluctuations. The LBMC simulations are in general agreement with the disulfide trapping experiments of Careaga & Falke (J. Mol. Biol., 1992, Vol. 226, 1219-35), which indicate that distant residues in the crystal structure (i.e. beta carbons separated by 10 to 20 angstroms) form spontaneous transient contacts in solution. Our simulations illustrate several possible "mechanisms" (configurational pathways) for these fluctuations. We also observe several discrepancies between our calculations and experimental rate constants. Nevertheless, we believe that our semi-atomistic approach could be used to study fluctuations in other proteins, perhaps for ensemble docking or other analyses of protein flexibility in virtual screening studies.

Binding of netropsin to several DNA constructs: evidence for at least two different 1:1 complexes formed from an -AATT-containing ds-DNA construct and a single minor groove binding ligand.

Isothermal titration calorimetry, ITC, has been used to determine the thermodynamics (DeltaG, DeltaH, and -TDeltaS) for binding netropsin to a number of DNA constructs. The DNA constructs included: six different 20-22mer hairpin forming sequences and an 8-mer DNA forming a duplex dimer. All DNA constructs had a single -AT-rich netropsin binding with one of the following sequences, (A(2)T(2))(2), (ATAT)(2), or (AAAA/TTTT). Binding energetics are less dependent on site sequence than on changes in the neighboring single stranded DNA (hairpin loop size and tail length). All of the 1:1 complexes exhibit an enthalpy change that is dependent on the fractional saturation of the binding site. Later binding ligands interact with a significantly more favorable enthalpy change (partial differential DeltaH(1-2) from 2 to 6 kcal/mol) and a significantly less favorable entropy change (partial differential (-TDeltaS(1-2))) from -4 to -9 kcal/mol). The ITC data could only be fit within expected experimental error by use of a thermodynamic model that includes two independent binding processes with a combined stoichiometry of 1 mol of ligand per 1 mol of oligonucleotide. Based on the biophysical evidence reported here, including theoretical calculations for the energetics of "trapping" or structuring of a single water molecule and molecular docking computations, it is proposed that there are two modes by which flexible ligands can bind in the minor groove of duplex DNA. The higher affinity binding mode is for netropsin to lay along the floor of the minor groove in a bent conformation and exclude all water from the groove. The slightly weaker binding mode is for the netropsin molecule to have a slightly more linear conformation and for the required curvature to be the result of a water molecule that bridges between the floor of the minor groove and two of the amidino nitrogens located at one end of the bound netropsin molecule.

Which aminoglycoside ring is most important for binding? A hydropathic analysis of gentamicin, paromomycin, and analogues.

The NMR structures of gentamicin and paromomycin in complex with the A-site of Escherichia coli 16S ribosomal RNA were modified with molecular modeling to 12 analogues. The intermolecular interactions between these molecules and RNA were examined using the HINT (Hydropathic INTeractions) computational model to obtain interaction scores that have been shown previously to be related to free energy. The calculations correlated well with experimental binding data, and the interaction scores were used to analyze the specific structural features of each aminoglycoside that contribute to the overall binding with the 16S rRNA. Our calculations indicate that, while ring I binds to the main binding pocket of the rRNA A-site, ring IV of paromomycin-based aminoglycosides contributes significantly to the overall binding.

Cloning and sequencing of the medium-chain S-acyl fatty acid synthetase thioester hydrolase cDNA from rat mammary gland.

cDNA clones coding for the medium-chain S-acyl fatty acid synthetase thioester hydrolase (thioesterase II) from rat mammary gland were identified in a bacteriophage lambda gt11 library and their nucleotide sequences were determined. The predicted coding region spans 263 amino acid residues and includes a sequence identical with that of a peptide derived from the enzyme active site. The rat thioesterase II cDNA sequence exhibits homology with that of a thioesterase found in duck uropygial glands.

Hemoglobin Long Island is caused by a single mutation (adenine to cytosine) resulting in a failure to cleave amino-terminal methionine.

Hemoglobin Long Island has two separate amino acid abnormalities of beta-globin structure: an extension of the NH2 terminus by a methionine residue and a histidine-to-proline substitution at the normal second position. The NH2-terminal methionine residue, the translation product of an AUG initiation codon, is present only transiently in nascent proteins. Because of the general biological implications of this abnormality, we investigated the nature of the genetic defect of this mutant. We determined the sequence of the relevant portion of the beta-globin mRNA by means of dideoxynucleotide chain termination of the complementary DNA (cDNA) in which an oligonucleotide complementary to codons 10-17 was used as a primer for reverse transcriptase. A histidine-to-proline substitution was confirmed in the mutant mRNA by identifying an adenine-to-cytosine transversion in the second codon. However, we were unable to find any other abnormality at either the AUG initiation codon or in the 56 bases upstream from the adenine-to-cytosine transversion (encompassing most of the 5' untranslated region of the mutant beta-globin mRNA). Thus, it appears that this single lesion probably interferes with the poorly understood methionine-cleaving mechanism that modulates most of prokaryotic and eukaryotic proteins.

Concerted action of beta-glucuronidase and beta-acetylglucosaminidase on hyaluronodextrins.

A kinetic analysis of the stepwise alternating action of beta-glucuronidase and beta-acetylglucosaminidase on oligosaccharides and dextrins derived from hyaluronic acid was undertaken, for better definition of the contribution of this process to hyaluronate catabolism. Production of monosaccharide from larger dextrins by action of either enzyme is powerfully inhibited by electrolyts. In the study, as in mammalian tissues, beta-glucuronidase is present in excess so that the concentration of beta-acetylglucosaminidase is rate controlling in the action on dextrin substrates. For this action, Vmax shows limited variation with ionic strength or molecular weight of substrate. At ionic strength 0.03, but not 0.18, Km decreases some 100-fold for increase of molecular weight from 2,000 to 15,000. It is specifically this decrease in Km that accounts for the prominent electrolyte inhibition observed with larger dextrins. The extremely low values of Km are attributed to multiple ionic enzyme-substrate interactions at sites remote from the catalytic center. The previously reported stimulation by electrolyte of the action of beta-glucuronidase and beta-acetylglucosaminidase on aryl glycosides, studied briefly, is apparently unrelated to the electrolyte effects seen with dextrins. The catabolic contribution of beta-glucuronidase and beta-acetylglucosaminidase appears to be restricted to hydrolysis of the smaller oligosaccharides produced by action of hyaluronidase, since, for any reasonable assumptions regarding cellular environment, the extent of their action on polymeric hyaluronate or larger dextrins must be limited.