Structure-activity relationships derived by machine learning: the use of atoms and their bond connectivities to predict mutagenicity by inductive logic programming.

We present a general approach to forming structure-activity relationships (SARs). This approach is based on representing chemical structure by atoms and their bond connectivities in combination with the inductive logic programming (ILP) algorithm PROGOL. Existing SAR methods describe chemical structure by using attributes which are general properties of an object. It is not possible to map chemical structure directly to attribute-based descriptions, as such descriptions have no internal organization. A more natural and general way to describe chemical structure is to use a relational description, where the internal construction of the description maps that of the object described. Our atom and bond connectivities representation is a relational description. ILP algorithms can form SARs with relational descriptions. We have tested the relational approach by investigating the SARs of 230 aromatic and heteroaromatic nitro compounds. These compounds had been split previously into two subsets, 188 compounds that were amenable to regression and 42 that were not. For the 188 compounds, a SAR was found that was as accurate as the best statistical or neural network-generated SARs. The PROGOL SAR has the advantages that it did not need the use of any indicator variables handcrafted by an expert, and the generated rules were easily comprehensible. For the 42 compounds, PROGOL formed a SAR that was significantly (P < 0.025) more accurate than linear regression, quadratic regression, and back-propagation. This SAR is based on an automatically generated structural alert for mutagenicity.

High yield conversion of doxorubicin to 2-pyrrolinodoxorubicin, an analog 500-1000 times more potent: structure-activity relationship of daunosamine-modified derivatives of doxorubicin.

A convenient, high yield conversion of doxorubicin to 3'-deamino-3'-(2''-pyrroline-1''-yl)doxorubicin is described. This daunosamine-modified analog of doxorubicin is 500-1000 times more active in vitro than doxorubicin. The conversion is effected by using a 30-fold excess of 4-iodobutyraldehyde in anhydrous dimethylformamide. The yield is higher than 85%. A homolog of this compound, 3'-deamino-3'-(1'',3''-tetrahydropyridine-1''-yl)doxorubicin, was also synthesized by using 5-iodovaleraldehyde. In this homolog, the daunosamine nitrogen is incorporated into a six- instead of a five-membered ring. This analog was 30-50 times less active than its counterpart with a five-membered ring. A similar structure-activity relationship was found when 3'-deamino-3'-(3''-pyrrolidone-1''-yl)doxorubicin (containing a five-membered ring) and 3'-deamino-3'-(3''-piperidone-1''-yl)doxorubicin (with a six-membered ring) were tested in vitro, the former being 5 times more potent than the latter. To further elucidate structure-activity relationships, 3'-deamino-3'-(pyrrolidine-1''-yl)doxorubicin, 3'-deamino-3'-(isoindoline-2''-yl)doxorubicin, 3'-deamino-3'-(2''-methyl-2''-pyrroline-1''-yl)doxorubicin, and 3'-deamino-3'-(3''-pyrroline-1''-yl)doxorubicin were also synthesized and tested. All the analogs were prepared by using reactive halogen compounds for incorporating the daunosamine nitrogen of doxorubicin into a five- or six-membered ring. These highly active antineoplastic agents can be used for incorporation into targeted cytotoxic analogs of luteinizing hormone-releasing hormone intended for cancer therapy.

On the reaction mechanism of l-lactate oxidase: Quantitative structure-activity analysis of the reaction with para-substituted l-mandelates

The rate constants for reduction of the flavoenzyme, l-lactate oxidase, and a mutant (in which alanine 95 is replaced by glycine), by a series of para-substituted mandelates, in both the 2-1H- and 2-2H- forms, have been measured by rapid reaction spectrophotometry. In all cases, significant isotope effects (1H/2H = 3–7) on the rate constants of flavin reduction were found, indicating that flavin reduction is a direct measure of α-C-H bond breakage. The rate constants show only a small influence of the electronic characteristics of the substituents, but show a good correlation when combined with some substituent volume parameters. A surprisingly good correlation is found with the molecular mass of the substrate. The results are compatible with any mechanism in which there is little development of charge in the transition state. This could be a transfer of hydride to the flavin N(5) position or a synchronous mechanism in which the α-C-H is formally abstracted as a H+ while the resulting charge is simultaneously neutralized by another event.

Vascular imprints of neuronal activity: Relationships between the dynamics of cortical blood flow, oxygenation, and volume changes following sensory stimulation

Modern functional neuroimaging methods, such as positron-emission tomography (PET), optical imaging of intrinsic signals, and functional MRI (fMRI) utilize activity-dependent hemodynamic changes to obtain indirect maps of the evoked electrical activity in the brain. Whereas PET and flow-sensitive MRI map cerebral blood flow (CBF) changes, optical imaging and blood oxygenation level-dependent MRI map areas with changes in the concentration of deoxygenated hemoglobin (HbR). However, the relationship between CBF and HbR during functional activation has never been tested experimentally. Therefore, we investigated this relationship by using imaging spectroscopy and laser-Doppler flowmetry techniques, simultaneously, in the visual cortex of anesthetized cats during sensory stimulation. We found that the earliest microcirculatory change was indeed an increase in HbR, whereas the CBF increase lagged by more than a second after the increase in HbR. The increased HbR was accompanied by a simultaneous increase in total hemoglobin concentration (Hbt), presumably reflecting an early blood volume increase. We found that the CBF changes lagged after Hbt changes by 1 to 2 sec throughout the response. These results support the notion of active neurovascular regulation of blood volume in the capillary bed and the existence of a delayed, passive process of capillary filling.

ACTIVITY: a database on DNA/RNA sites activity adapted to apply sequence-activity relationships from one system to another

ACTIVITY is a database on DNA/RNA site sequences with known activity magnitudes, measurement systems, sequence-activity relationships under fixed experimental conditions and procedures to adapt these relationships from one measurement system to another. This database deposits information on DNA/RNA affinities to proteins and cell nuclear extracts, cutting efficiencies, gene transcription activity, mRNA translation efficiencies, mutability and other biological activities of natural sites occurring within promoters, mRNA leaders, and other regulatory regions in pro- and eukaryotic genomes, their mutant forms and synthetic analogues. Since activity magnitudes are heavily system-dependent, the current version of ACTIVITY is supplemented by three novel sub-databases: (i) SYSTEM, measurement systems; (ii) KNOWLEDGE, sequence-activity relationships under fixed experimental conditions; and (iii) CROSS_TEST, procedures adapting a relationship from one measurement system to another. These databases are useful in molecular biology, pharmacogenetics, metabolic engineering, drug design and biotechnology. The databases can be queried using SRS and are available through the Web, http://wwwmgs.bionet.nsc.ru/systems/Activity/.

Structure-activity analysis of thanatin, a 21-residue inducible insect defense peptide with sequence homology to frog skin antimicrobial peptides.

Immune challenge to the insect Podisus maculiventris induces synthesis of a 21-residue peptide with sequence homology to frog skin antimicrobial peptides of the brevinin family. The insect and frog peptides have in common a C-terminally located disulfide bridge delineating a cationic loop. The peptide is bactericidal and fungicidal, exhibiting the largest antimicrobial spectrum observed so far for an insect defense peptide. An all-D-enantiomer is nearly inactive against Gram-negative bacteria and some Gram-positive strains but is fully active against fungi and other Gram-positive bacteria, suggesting that more than one mechanism accounts for the antimicrobial activity of this peptide. Studies with truncated synthetic isoforms underline the role of the C-terminal loop and flanking residues for the activity of this molecule for which we propose the name thanatin.

Structural studies of p21Waf1/Cip1/Sdi1 in the free and Cdk2-bound state: conformational disorder mediates binding diversity.

The cyclin-dependent kinase (Cdk) inhibitor p21Waf1/Cip1/Sdi1, important for p53-dependent cell cycle control, mediates G1/S arrest through inhibition of Cdks and possibly through inhibition of DNA replication. Cdk inhibition requires a sequence of approximately 60 amino acids within the p21 NH2 terminus. We show, using proteolytic mapping, circular dichroism spectropolarimetry, and nuclear magnetic resonance spectroscopy, that p21 and NH2-terminal fragments that are active as Cdk inhibitors lack stable secondary or tertiary structure in the free solution state. In sharp contrast to the disordered free state, however, the p21 NH2 terminus adopts an ordered stable conformation when bound to Cdk2, as shown directly by NMR spectroscopy. We have, thus, identified a striking disorder-order transition for p21 upon binding to one of its biological targets, Cdk2. This structural transition has profound implications in light of the ability of p21 to bind and inhibit a diverse family of cyclin-Cdk complexes, including cyclin A-Cdk2, cyclin E-Cdk2, and cyclin D-Cdk4. Our findings suggest that the flexibility, or disorder, of free p21 is associated with binding diversity and offer insights into the role for structural disorder in mediating binding specificity in biological systems. Further, these observations challenge the generally accepted view of proteins that stable secondary and tertiary structure are prerequisites for biological activity and suggest that a broader view of protein structure should be considered in the context of structure-activity relationships.

Rational interleukin 2 therapy for HIV positive individuals: daily low doses enhance immune function without toxicity.

When administered in high doses to HIV positive (HIV+) individuals, interleukin 2 (IL-2) causes extreme toxicity and markedly increases plasma HIV levels. Integration of the information from the structure-activity relationships of the IL-2 receptor interaction, the cellular distribution of the different classes of IL-2 receptors, and the pharmacokinetics of IL-2 provides for the rationale that low IL-2 doses should circumvent toxicity. Therefore, to identify a nontoxic, but effective and safe IL-2 treatment regimen that does not stimulate viral replication, doses of IL-2 from 62,500 to 250,000 IU/m2/day were administered subcutaneously for 6 months to 16 HIV+ individuals with 200-500 CD4+ T cells/mm3. IL-2 was already detectable in the plasma of most HIV+ individuals even before therapy. Peak plasma IL-2 levels were near saturating for high affinity IL-2 receptors in 10 individuals who received the maximum nontoxic dose, which ranged from 187,500 to 250,000 IU/m2/day. During the 6 months of treatment at this dose range, plasma levels of proinflammatory cytokines remained undetectable, and plasma HIV RNA levels did not change significantly. However, delayed type hypersensitivity responses to common recall antigens were markedly augmented, and there were IL-2 dose-dependent increases in circulating Natural Killer cells, eosinophils, monocytes, and CD4+ T cells. Expanded clinical trials of low dose IL-2 are now warranted, especially in combination with effective antivirals to test for the prevention of immunodeficiency and the emergence of drug-resistant mutants and for the eradication of residual virions.

Expression of human inducible nitric oxide synthase in a tetrahydrobiopterin (H4B)-deficient cell line: H4B promotes assembly of enzyme subunits into an active dimer.

Murine inducible nitric oxide (NO) synthase (iNOS) is catalytically active only in dimeric form. Assembly of its purified subunits into a dimer requires H4B. To understand the structure-activity relationships of human iNOS, we constitutively expressed recombinant human iNOS in NIH 3T3 cells by using a retroviral vector. These cells are deficient in de novo H4B biosynthesis and the role of H4B in the expression and assembly of active iNOS in an intact cell system could be studied. In the absence of added H4B, NO synthesis by the cells was minimal, whereas cells grown with supplemental H4B or the H4B precursor sepiapterin generated NO (74.1 and 63.3 nmol of nitrite per 10(6) cells per 24 h, respectively). NO synthesis correlated with an increase in intracellular H4B but no increase in iNOS protein. Instead, an increased percentage of dimeric iNOS was observed, rising from 20% in cytosols from unsupplemented cells to 66% in H4B-supplemented cell cytosols. In all cases, only dimeric iNOS displayed catalytic activity. Cytosols prepared from H4B-deficient cells exhibited little iNOS activity but acquired activity during a 60- to 120-min incubation with H4B, reaching final activities of 60-72 pmol of citrulline per mg of protein per min. Reconstitution of cytosolic NO synthesis activity was associated with conversion of monomers into dimeric iNOS during the incubation. Thus, human iNOS subunits dimerize to form an active enzyme, and H4B plays a critical role in promoting dimerization in intact cells. This reveals a post-translational mechanism by which intracellular H4B can regulate iNOS expression.

A paradigm for drug discovery employing encoded combinatorial libraries.

Very large combinatorial libraries of small molecules on solid supports can now be synthesized and each library element can be identified after synthesis by using chemical tags. These tag-encoded libraries are potentially useful in drug discovery, and, to test this utility directly, we have targeted carbonic anhydrase (carbonate dehydratase; carbonate hydro-lyase, EC 4.2.1.1) as a model. Two libraries consisting of a total of 7870 members were synthesized, and structure-activity relationships based on the structures predicted by the tags were derived. Subsequently, an active representative of each library was resynthesized (2-[N-(4-sulfamoylbenzoyl)-4'-aminocyclohexanespiro]-4-oxo-7 -hydroxy- 2,3-dihydrobenzopyran and [N-(4-sulfamoylbenzoyl)-L-leucyl]piperidine-3-carboxylic acid) and these compounds were shown to have nanomolar dissociation constants (15 and 4 nM, respectively). In addition, a focused sublibrary of 217 sulfamoylbenzamides was synthesized and revealed a clear, testable structure-activity relationship describing isozyme-selective carbonic anhydrase inhibitors.

Heterodimer formation and activity in the human enzyme galactose-1-phosphate uridylyltransferase.

One of the fundamental questions concerning expression and function of dimeric enzymes involves the impact of naturally occurring mutations on subunit assembly and heterodimer activity. This question is of particular interest for the human enzyme galactose-l-phosphate uridylyl-transferase (GALT), impairment of which results in the inherited metabolic disorder galactosemia, because many if not most patients studied to date are compound heterozygotes rather than true molecular homozygotes. Furthermore, the broad range of phenotypic severity observed in these patients raises the possibility that allelic combination, not just allelic constitution, may play some role in determining outcome. In the work described herein, we have selected two distinct naturally occurring null mutations of GALT, Q188R and R333W, and asked the questions (i) what are the impacts of these mutations on subunit assembly, and (ii) if heterodimers do form, are they active? To answer these questions, we have established a yeast system for the coexpression of epitope-tagged alleles of human GALT and investigated both the extent of specific GALT subunit interactions and the activity of defined heterodimer pools. We have found that both homodimers and heterodimers do form involving each of the mutant subunits tested and that both heterodimer pools retain substantial enzymatic activity. These results are significant not only in terms of their implications for furthering our understanding of galactosemia and GALT holoenzyme structure-function relationships but also because the system described may serve as a model for similar studies of other complexes composed of multiple subunits.

Sequencing of 3-O sulfate containing heparin decasaccharides with a partial antithrombin III binding site

Heparin- and heparan sulfate-like glycosaminoglycans (HLGAGs) represent an important class of molecules that interact with and modulate the activity of growth factors, enzymes, and morphogens. Of the many biological functions for this class of molecules, one of its most important functions is its interaction with antithrombin III (AT-III). AT-III binding to a specific heparin pentasaccharide sequence, containing an unusual 3-O sulfate on a N-sulfated, 6-O sulfated glucosamine, increases 1,000-fold AT-III's ability to inhibit specific proteases in the coagulation cascade. In this manner, HLGAGs play an important biological and pharmacological role in the modulation of blood clotting. Recently, a sequencing methodology was developed to further structure-function relationships of this important class of molecules. This methodology combines a property-encoded nomenclature scheme to handle the large information content (properties) of HLGAGs, with matrix-assisted laser desorption ionization MS and enzymatic and chemical degradation as experimental constraints to rapidly sequence picomole quantities of HLGAG oligosaccharides. Using the above property-encoded nomenclature-matrix-assisted laser desorption ionization approach, we found that the sequence of the decasaccharide used in this study is ΔU2SHNS,6SI2SHNS,6SI2SHNS,6SIHNAc,6SGHNS,3S,6S (±DDD4–7). We confirmed our results by using integral glycan sequencing and one-dimensional proton NMR. Furthermore, we show that this approach is flexible and is able to derive sequence information on an oligosaccharide mixture. Thus, this methodology will make possible both the analysis of other unusual sequences in HLGAGs with important biological activity as well as provide the basis for the structural analysis of these pharamacologically important group of heparin/heparan sulfates.

A High-Affinity Ca2+ Pump, ECA1, from the Endoplasmic Reticulum Is Inhibited by Cyclopiazonic Acid but Not by Thapsigargin1

To identify and characterize individual Ca2+ pumps, we have expressed an Arabidopsis ECA1 gene encoding an endoplasmic reticulum-type Ca2+-ATPase homolog in the yeast (Saccharomyces cerevisiae) mutant K616. The mutant (pmc1pmr1cnb1) lacks a Golgi and a vacuolar membrane Ca2+ pump and grows very poorly on Ca2+-depleted medium. Membranes isolated from the mutant showed high H+/Ca2+-antiport but no Ca2+-pump activity. Expression of ECA1 in endomembranes increased mutant growth by 10- to 20-fold in Ca2+-depleted medium. 45Ca2+ pumping into vesicles from ECA1 transformants was detected after the H+/Ca2+-antiport activity was eliminated with bafilomycin A1 and gramicidin D. The pump had a high affinity for Ca2+ (Km = 30 nm) and displayed two affinities for ATP (Km of 20 and 235 μm). Cyclopiazonic acid, a specific blocker of animal sarcoplasmic/endoplasmic reticulum Ca2+-ATPase, inhibited Ca2+ transport (50% inhibition dose = 3 nmol/mg protein), but thapsigargin (3 μm) did not. Transport was insensitive to calmodulin. These results suggest that this endoplasmic reticulum-type Ca2+-ATPase could support cell growth in plants as in yeast by maintaining submicromolar levels of cytosolic Ca2+ and replenishing Ca2+ in endomembrane compartments. This study demonstrates that the yeast K616 mutant provides a powerful expression system to study the structure/function relationships of Ca2+ pumps from eukaryotes.

Solution structure of the 2-amino-1- methyl-6-phenylimidazo[4,5-b]pyridine C8-deoxyguanosine adduct in duplex DNA

The carcinogenic heterocyclic amine (HA) 2-amino-1-methyl-6-phenylimidazo[4,5-b]pyridine (PhIP) is formed during the cooking of various meats. To enable structure/activity studies aimed at understanding how DNA damaged by a member of the HA class of compounds can ultimately lead to cancer, we have determined the first solution structure of an 11-mer duplex containing the C8-dG adduct formed by reaction with N-acetoxy-PhIP. A slow conformational exchange is observed in which the PhIP ligand either intercalates into the DNA helix by denaturing and displacing the modified base pair (main form) or is located outside the helix in a minimally perturbed B-DNA duplex (minor form). In the main base-displaced intercalation structure, the minor groove is widened, and the major groove is compressed at the lesion site because of the location of the bulky PhIP-N-methyl and phenyl ring in the minor groove; this distortion causes significant bending of the helix. The PhIP phenyl ring interacts with the phosphodiester-sugar ring backbone of the complementary strand and its fast rotation with respect to the intercalated imidazopyridine ring causes substantial distortions at this site, such as unwinding and bulging-out of the strand. The glycosidic torsion angle of the [PhIP]dG residue is syn, and the displaced guanine base is directed toward the 3′ end of the modified strand. This study contributes, to our knowledge, the first structural information on the biologically relevant HA class to a growing body of knowledge about how conformational similarities and differences for a variety of types of lesions can influence protein interactions and ultimately biological outcome.