Mode matches and their locations in the hydrophobic free energy sequences of peptide ligands and their receptor eigenfunctions

Patterns in sequences of amino acid hydrophobic free energies predict secondary structures in proteins. In protein folding, matches in hydrophobic free energy statistical wavelengths appear to contribute to selective aggregation of secondary structures in “hydrophobic zippers.” In a similar setting, the use of Fourier analysis to characterize the dominant statistical wavelengths of peptide ligands’ and receptor proteins’ hydrophobic modes to predict such matches has been limited by the aliasing and end effects of short peptide lengths, as well as the broad-band, mode multiplicity of many of their frequency (power) spectra. In addition, the sequence locations of the matching modes are lost in this transformation. We make new use of three techniques to address these difficulties: (i) eigenfunction construction from the linear decomposition of the lagged covariance matrices of the ligands and receptors as hydrophobic free energy sequences; (ii) maximum entropy, complex poles power spectra, which select the dominant modes of the hydrophobic free energy sequences or their eigenfunctions; and (iii) discrete, best bases, trigonometric wavelet transformations, which confirm the dominant spectral frequencies of the eigenfunctions and locate them as (absolute valued) moduli in the peptide or receptor sequence. The leading eigenfunction of the covariance matrix of a transmembrane receptor sequence locates the same transmembrane segments seen in n-block-averaged hydropathy plots while leaving the remaining hydrophobic modes unsmoothed and available for further analyses as secondary eigenfunctions. In these receptor eigenfunctions, we find a set of statistical wavelength matches between peptide ligands and their G-protein and tyrosine kinase coupled receptors, ranging across examples from 13.10 amino acids in acid fibroblast growth factor to 2.18 residues in corticotropin releasing factor. We find that the wavelet-located receptor modes in the extracellular loops are compatible with studies of receptor chimeric exchanges and point mutations. A nonbinding corticotropin-releasing factor receptor mutant is shown to have lost the signatory mode common to the normal receptor and its ligand. Hydrophobic free energy eigenfunctions and their transformations offer new quantitative physical homologies in database searches for peptide-receptor matches.

Comparing protein–ligand interactions in solution and single crystals by Raman spectroscopy

By using a Raman microscope, we show that it is possible to probe the conformational states in protein crystals and crystal fragments under growth conditions (in hanging drops). The flavin cofactor in the enzyme para-hydroxybenzoate hydroxylase can assume two conformations: buried in the protein matrix (“in”) or essentially solvent-exposed (“out”). By using Raman difference spectroscopy, we previously have identified characteristic flavin marker bands for the in and out conformers in the solution phase. Now we show that the flavin Raman bands can be used to probe these conformational states in crystals, permitting a comparison between solution and crystal environments. The in or out marker bands are similar for the respective conformers in the crystal and in solution; however, significant differences do exist, showing that the environments for the flavin's isoalloxazine ring are not identical in the two phases. Moreover, the Raman-band widths of the flavin modes are narrower for both in and out conformers in the crystals, indicating that the flavin exists in a more limited range of closely related conformational states in the crystal than in solution. In general, the ability to compare detailed Raman data for complexes in crystals and solution provides a means of bridging crystallographic and solution studies.

Shedding light on the dark and weakly fluorescent states of green fluorescent proteins

Recent experiments on various similar green fluorescent protein (GFP) mutants at the single-molecule level and in solution provide evidence of previously unknown short- and long-lived “dark” states and of related excited-state decay channels. Here, we present quantum chemical calculations on cis-trans photoisomerization paths of neutral, anionic, and zwitterionic GFP chromophores in their ground and first singlet excited states that explain the observed behaviors from a common perspective. The results suggest that favorable radiationless decay channels can exist for the different protonation states along these isomerizations, which apparently proceed via conical intersections. These channels are suggested to rationalize the observed dramatic reduction of fluorescence in solution. The observed single-molecule fast blinking is attributed to conversions between the fluorescent anionic and the dark zwitterionic forms whereas slow switching is attributed to conversions between the anionic and the neutral forms. The predicted nonadiabatic crossings are seen to rationalize the origins of a variety of experimental observations on a common basis and may have broad implications for photobiophysical mechanisms in GFP.

CC-1065 and the duocarmycins: unraveling the keys to a new class of naturally derived DNA alkylating agents.

Key studies defining the DNA alkylation properties and selectivity of a new class of exceptionally potent, naturally occurring antitumor antibiotics including CC-1065, duocarmycin A, and duocarmycin SA are reviewed. Recent studies conducted with synthetic agents containing deep-seated structural changes and the unnatural enantiomers of the natural products and related analogs have defined the structural basis for the sequence-selective alkylation of duplex DNA and fundamental relationships between chemical structure, functional reactivity, and biological properties. The agents undergo a reversible, stereoelectronically controlled adenine-N3 addition to the least substituted carbon of the activated cyclopropane within selected AT-rich sites. The preferential AT-rich non-covalent binding selectivity of the agents within the narrower, deeper AT-rich minor groove and the steric accessibility to the alkylation site that accompanies deep AT-rich minor groove penetration control the sequence-selective DNA alkylation reaction and stabilize the resulting adduct. For the agents that possess sufficient reactivity to alkylate DNA, a direct relationship between chemical or functional stability and biological potency has been defined.

Subtle atomic group discrimination in the RNA minor groove

As a problem in molecular recognition and for drug discovery, great interest has developed around the possibility that RNA structures could be discriminated by peptides and other small molecules. Although small peptides have been shown to have the capacity to discriminate specific bulges and loops in RNA molecules, discrimination of double helical regions by a peptide binder has not been reported. Indeed, the most accessible part of an RNA helix is the minor groove, and fundamental stereochemical considerations have suggested that discrimination of at least some base pairs would be difficult in the minor groove. Here we report the design and isolation of a peptide binder that manifests the most subtle kind of discrimination of base pair differences in the RNA minor groove. Functional discrimination of a single atomic group is demonstrated as well as the difference between two different angular orientations of the same group. This report of RNA helix discrimination by a peptide binder suggests a richer potential for RNA minor groove recognition than previously thought.

Carbonic anhydrase is an ancient enzyme widespread in prokaryotes

Carbonic anhydrases catalyze the reversible hydration of CO2 and are ubiquitous in highly evolved eukaryotes. The recent identification of a third class of carbonic anhydrase (γ class) in a methanoarchaeon and our present finding that the β class also extends into thermophilic species from the Archaea domain led us to initiate a systematic search for these enzymes in metabolically and phylogenetically diverse prokaryotes. Here we show that carbonic anhydrase is widespread in the Archaea and Bacteria domains, and is an ancient enzyme. The occurrence in chemolithoautotrophic species occupying deep branches of the universal phylogenetic tree suggests a role for this enzyme in the proposed autotrophic origin of life. The presence of the β and γ classes in metabolically diverse species spanning the Archaea and Bacteria domains demonstrates that carbonic anhydrases have a far more extensive and fundamental role in prokaryotic biology than previously recognized.

Toward functional genomics in bacteria: Analysis of gene expression in Escherichia coli from a bacterial artificial chromosome library of Bacillus cereus

As the study of microbes moves into the era of functional genomics, there is an increasing need for molecular tools for analysis of a wide diversity of microorganisms. Currently, biological study of many prokaryotes of agricultural, medical, and fundamental scientific interest is limited by the lack of adequate genetic tools. We report the application of the bacterial artificial chromosome (BAC) vector to prokaryotic biology as a powerful approach to address this need. We constructed a BAC library in Escherichia coli from genomic DNA of the Gram-positive bacterium Bacillus cereus. This library provides 5.75-fold coverage of the B. cereus genome, with an average insert size of 98 kb. To determine the extent of heterologous expression of B. cereus genes in the library, we screened it for expression of several B. cereus activities in the E. coli host. Clones expressing 6 of 10 activities tested were identified in the library, namely, ampicillin resistance, zwittermicin A resistance, esculin hydrolysis, hemolysis, orange pigment production, and lecithinase activity. We analyzed selected BAC clones genetically to identify rapidly specific B. cereus loci. These results suggest that BAC libraries will provide a powerful approach for studying gene expression from diverse prokaryotes.

Computer simulations of enzyme catalysis: Finding out what has been optimized by evolution

The origin of the catalytic power of enzymes is discussed, paying attention to evolutionary constraints. It is pointed out that enzyme catalysis reflects energy contributions that cannot be determined uniquely by current experimental approaches without augmenting the analysis by computer simulation studies. The use of energy considerations and computer simulations allows one to exclude many of the popular proposals for the way enzymes work. It appears that the standard approaches used by organic chemists to catalyze reactions in solutions are not used by enzymes. This point is illustrated by considering the desolvation hypothesis and showing that it cannot account for a large increase in kcat relative to the corresponding kcage for the reference reaction in a solvent cage. The problems associated with other frequently invoked mechanisms also are outlined. Furthermore, it is pointed out that mutation studies are inconsistent with ground state destabilization mechanisms. After considering factors that were not optimized by evolution, we review computer simulation studies that reproduced the overall catalytic effect of different enzymes. These studies pointed toward electrostatic effects as the most important catalytic contributions. The nature of this electrostatic stabilization mechanism is far from being obvious because the electrostatic interaction between the reacting system and the surrounding area is similar in enzymes and in solution. However, the difference is that enzymes have a preorganized dipolar environment that does not have to pay the reorganization energy for stabilizing the relevant transition states. Apparently, the catalytic power of enzymes is stored in their folding energy in the form of the preorganized polar environment.

Homophilic adhesion mediated by the neural cell adhesion molecule involves multiple immunoglobulin domains.

The neural cell adhesion molecule (N-CAM) mediates homophilic binding between a variety of cell types including neurons, neurons and glia, and neurons and muscle cells. The mechanism by which N-CAM on one cell interacts with N-CAM on another, however, is unknown. Attempts to identify which of the five immunoglobulin-like domains (Ig I-V) and the two fibronectin type III repeats (FnIII 1-2) in the extracellular region of N-CAM are involved in this process have led to ambiguous results. We have generated soluble recombinant proteins corresponding to each of the individual immunoglobulin domains and the combined FnIII 1-2 and prepared polyclonal antibodies specific for each. The purified proteins and antibodies were used in aggregation experiments with fluorescent microspheres and chicken embryo brain cells to determine possible contributions of each domain to homophilic adhesion. The recombinant domains were tested for their ability to bind to purified native N-CAM, to bind to each other, and to inhibit the aggregation of N-CAM on microspheres and the aggregation of neuronal cells. Each of the immunoglobulin domains bound to N-CAM, and in solution all of the immunoglobulin domains inhibited the aggregation of N-CAM-coated microspheres. Soluble Ig II, Ig III, and Ig IV inhibited neuronal aggregation; antibodies against whole N-CAM, the Ig III domain, and the Ig I domain all strongly inhibited neuronal aggregation, as well as the aggregation of N-CAM-coated microspheres. Of all the domains, the third immunoglobulin domain alone demonstrated the ability to self-aggregate, whereas Ig I bound to Ig V and Ig II bound to Ig IV. The combined FnIII 1-2 exhibited a slight ability to self-aggregate but did not bind to any of the immunoglobulin-like domains. These results suggest that N-CAM-N-CAM binding involves all five immunoglobulin domains and prompt the hypothesis that in homophilic cell-cell binding mediated by N-CAM these domains may interact pairwise in an antiparallel orientation.

Linguistic aspects of speech synthesis.

The conversion of text to speech is seen as an analysis of the input text to obtain a common underlying linguistic description, followed by a synthesis of the output speech waveform from this fundamental specification. Hence, the comprehensive linguistic structure serving as the substrate for an utterance must be discovered by analysis from the text. The pronunciation of individual words in unrestricted text is determined by morphological analysis or letter-to-sound conversion, followed by specification of the word-level stress contour. In addition, many text character strings, such as titles, numbers, and acronyms, are abbreviations for normal words, which must be derived. To further refine these pronunciations and to discover the prosodic structure of the utterance, word part of speech must be computed, followed by a phrase-level parsing. From this structure the prosodic structure of the utterance can be determined, which is needed in order to specify the durational framework and fundamental frequency contour of the utterance. In discourse contexts, several factors such as the specification of new and old information, contrast, and pronominal reference can be used to further modify the prosodic specification. When the prosodic correlates have been computed and the segmental sequence is assembled, a complete input suitable for speech synthesis has been determined. Lastly, multilingual systems utilizing rule frameworks are mentioned, and future directions are characterized.