Mutational analysis and molecular modeling of the nonapeptide hormone binding domains of the [Arg8]vasotocin receptor.

To identify determinants that form nonapeptide hormone binding domains of the white sucker Catostomus commersoni [Arg8]vasotocin receptor, chimeric constructs encoding parts of the vasotocin receptor and parts of the isotocin receptor have been analyzed by [(3,5-3H)Tyr2, Arg8]vasotocin binding to membranes of human embryonic kidney cells previously transfected with the different cDNA constructs and by functional expression studies in Xenopus laevis oocytes injected with mutant cRNAs. The results indicate that the N terminus and a region spanning the second extracellular loop and its flanking transmembrane segments, which contains a number of amino acid residues that are conserved throughout the nonapeptide receptor family, contribute to the affinity of the receptor for its ligand. Nonapeptide selectivity, however, is mainly defined by transmembrane region VI and the third extracellular loop. These results are complemented by a molecular model of the vasotocin receptor obtained by aligning its sequence with those of other G-protein coupled receptors as well as that of bacteriorhodopsin. The model indicates that amino acid residues of transmembrane regions II-VII that are located close to the extracellular surface also contribute to the binding of vasotocin.

Flows of knowledge from universities and federal laboratories: Modeling the flow of patent citations over time and across institutional and geographic boundaries

The extent to which new technological knowledge flows across institutional and national boundaries is a question of great importance for public policy and the modeling of economic growth. In this paper we develop a model of the process generating subsequent citations to patents as a lens for viewing knowledge diffusion. We find that the probability of patent citation over time after a patent is granted fits well to a double-exponential function that can be interpreted as the mixture of diffusion and obsolescense functions. The results indicate that diffusion is geographically localized. Controlling for other factors, within-country citations are more numerous and come more quickly than those that cross country boundaries.

A combined kinetic and modeling study of the catalytic center subsites of human angiogenin.

Kinetic analysis and molecular modeling have been used to map the ribonucleolytic center of angiogenin (Ang). Pyrimidine nucleotides were found to interact very weakly with Ang, consistent with the inaccessible B1 pyrimidine binding site revealed by x-ray crystallography. Ang also lacks an effective phosphate binding site on the 5' side of B1. Although the B2 site that preferentially binds purines on the 3' side of B1 is also weak, its associated phosphate subsites make substantial contributions: both 3',5'-ADP and 5'-ADP have Ki values 6-fold lower than for 5'-AMP, and adding a 3'-phosphate to the substrate CpA increases Kcat/Km by 9-fold. Thus Ang has a functional P2 site on the 3' side of B2 and a site for a second phosphate on the 5' side of B2. Modeling of an Ang-d(ApTpApA) complex suggested that Arg-5 forms part of the P2 site and that a 2'-phosphate might bind more tightly than a 3'-phosphate. Both predictions were confirmed kinetically. The subsite map obtained by this combined approach indicated that 5'-diphosphoadenosine 2'-phosphate might be a more potent inhibitor than any of the nucleotides tested thus far. Indeed, its Ki value of 150 microM is 50-fold lower than that for the best nucleotide previously reported and 400-fold lower than the Km for the best dinucleotide substrate. This compound may serve as a suitable starting point for the eventual design of tight-binding inhibitors of Ang as antiangiogenic agents for human therapy.

Analysis of human cytochrome P450 3A4 cooperativity: Construction and characterization of a site-directed mutant that displays hyperbolic steroid hydroxylation kinetics

Cytochrome P450 3A4 is generally considered to be the most important human drug-metabolizing enzyme and is known to catalyze the oxidation of a number of substrates in a cooperative manner. An allosteric mechanism is usually invoked to explain the cooperativity. Based on a structure–activity study from another laboratory using various effector–substrate combinations and on our own studies using site-directed mutagenesis and computer modeling of P450 3A4, the most likely location of effector binding is in the active site along with the substrate. Our study was designed to test this hypothesis by replacing residues Leu-211 and Asp-214 with the larger Phe and Glu, respectively. These residues were predicted to constitute a portion of the effector binding site, and the substitutions were designed to mimic the action of the effector by reducing the size of the active site. The L211F/D214E double mutant displayed an increased rate of testosterone and progesterone 6β-hydroxylation at low substrate concentrations and a decreased level of heterotropic stimulation elicited by α-naphthoflavone. Kinetic analyses of the double mutant revealed the absence of homotropic cooperativity with either steroid substrate. At low substrate concentrations the steroid 6β-hydroxylase activity of the wild-type enzyme was stimulated by a second steroid, whereas L211F/D214E displayed simple substrate inhibition. To analyze L211F/D214E at a more mechanistic level, spectral binding studies were carried out. Testosterone binding by the wild-type enzyme displayed homotropic cooperativity, whereas substrate binding by L211F/D214E displayed hyperbolic behavior.

Genetic definition of a protein-splicing domain: Functional mini-inteins support structure predictions and a model for intein evolution

Inteins are protein-splicing elements, most of which contain conserved sequence blocks that define a family of homing endonucleases. Like group I introns that encode such endonucleases, inteins are mobile genetic elements. Recent crystallography and computer modeling studies suggest that inteins consist of two structural domains that correspond to the endonuclease and the protein-splicing elements. To determine whether the bipartite structure of inteins is mirrored by the functional independence of the protein-splicing domain, the entire endonuclease component was deleted from the Mycobacterium tuberculosis recA intein. Guided by computer modeling studies, and taking advantage of genetic systems designed to monitor intein function, the 440-aa Mtu recA intein was reduced to a functional mini-intein of 137 aa. The accuracy of splicing of several mini-inteins was verified. This work not only substantiates structure predictions for intein function but also supports the hypothesis that, like group I introns, mobile inteins arose by an endonuclease gene invading a sequence encoding a small, functional splicing element.

RNA-controlled polymorphism in the in vivo assembly of 180-subunit and 120-subunit virions from a single capsid protein

Repeated, specific interactions between capsid protein (CP) subunits direct virus capsid assembly and exemplify regulated protein–protein interactions. The results presented here reveal a striking in vivo switch in CP assembly. Using cryoelectron microscopy, three-dimensional image reconstruction, and molecular modeling, we show that brome mosaic virus (BMV) CP can assemble in vivo two remarkably distinct capsids that selectively package BMV-derived RNAs in the absence of BMV RNA replication: a 180-subunit capsid indistinguishable from virions produced in natural infections and a previously unobserved BMV capsid type with 120 subunits arranged as 60 CP dimers. Each such dimer contains two CPs in distinct, nonequivalent environments, in contrast to the quasi-equivalent CP environments throughout the 180-subunit capsid. This 120-subunit capsid utilizes most of the CP interactions of the 180-subunit capsid plus nonequivalent CP–CP interactions. Thus, the CP of BMV, and perhaps other viruses, can encode CP–CP interactions that are not apparent from mature virions and may function in assembly or disassembly. Shared structural features suggest that the 120- and 180-subunit capsids share assembly steps and that a common pentamer of CP dimers may be an important assembly intermediate. The ability of a single CP to switch between distinct capsids by means of alternate interactions also implies reduced evolutionary barriers between different capsid structures. The in vivo switch between alternate BMV capsids is controlled by the RNA packaged: a natural BMV genomic RNA was packaged in 180-subunit capsids, whereas an engineered mRNA containing only the BMV CP gene was packaged in 120-subunit capsids. RNA features can thus direct the assembly of a ribonucleoprotein complex between alternate structural pathways.

Stretched and overwound DNA forms a Pauling-like structure with exposed bases

We investigate structural transitions within a single stretched and supercoiled DNA molecule. With negative supercoiling, for a stretching force >0.3 pN, we observe the coexistence of B-DNA and denatured DNA from σ ≈ −0.015 down to σ = −1. Surprisingly, for positively supercoiled DNA (σ > +0.037) stretched by 3 pN, we observe a similar coexistence of B-DNA and a new, highly twisted structure. Experimental data and molecular modeling suggest that this structure has ≈2.62 bases per turn and an extension 75% larger than B-DNA. This structure has tightly interwound phosphate backbones and exposed bases in common with Pauling’s early DNA structure [Pauling, L. & Corey, R. B. (1953), Proc. Natl. Acad. Sci. USA 39, 84–97] and an unusual structure proposed for the Pf1 bacteriophage [Liu, D. J. & Day, L. A. (1994) Science 265, 671–674].

From patterns to processes: Phase and density dependencies in the Canadian lynx cycle

Across the boreal forest of North America, lynx populations undergo 10-year cycles. Analysis of 21 time series from 1821 to the present demonstrates that these fluctuations are generated by nonlinear processes with regulatory delays. Trophic interactions between lynx and hares cause delayed density-dependent regulation of lynx population growth. The nonlinearity, in contrast, appears to arise from phase dependencies in hunting success by lynx through the cycle. Using a combined approach of empirical, statistical, and mathematical modeling, we highlight how shifts in trophic interactions between the lynx and the hare generate the nonlinear process primarily by shifting functional response curves during the increase and the decrease phases.

Selenoprotein oxidoreductase with specificity for thioredoxin and glutathione systems

Thioredoxin (Trx) and glutathione (GSH) systems are considered to be two major redox systems in animal cells. They are reduced by NADPH via Trx reductase (TR) or oxidized GSH (GSSG) reductase and further supply electrons for deoxyribonucleotide synthesis, antioxidant defense, and redox regulation of signal transduction, transcription, cell growth, and apoptosis. We cloned and characterized a pyridine nucleotide disulfide oxidoreductase, Trx and GSSG reductase (TGR), that exhibits specificity for both redox systems. This enzyme contains a selenocysteine residue encoded by the TGA codon. TGR can reduce Trx, GSSG, and a GSH-linked disulfide in in vitro assays. This unusual substrate specificity is achieved by an evolutionary conserved fusion of the TR and glutaredoxin domains. These observations, together with the biochemical probing and molecular modeling of the TGR structure, suggest a mechanism whereby the C-terminal selenotetrapeptide serves a role of a protein-linked GSSG and shuttles electrons from the disulfide center within the TR domain to either the glutaredoxin domain or Trx.

A family of lysozyme-like virulence factors in bacterial pathogens of plants and animals.

We describe a conserved family of bacterial gene products that includes the VirB1 virulence factor encoded by tumor-inducing plasmids of Agrobacterium spp., proteins involved in conjugative DNA transfer of broad-host-range bacterial plasmids, and gene products that may be involved in invasion by Shigella spp. and Salmonella enterica. Sequence analysis and structural modeling show that the proteins in this group are related to chicken egg white lysozyme and are likely to adopt a lysozyme-like structural fold. Based on their similarity to lysozyme, we predict that these proteins have glycosidase activity. Iterative data base searches with three conserved sequence motifs from this protein family detect a more distant relationship to bacterial and bacteriophage lytic transglycosylases, and goose egg white lysozyme. Two acidic residues in the VirB1 protein of Agrobacterium tumefaciens form a putative catalytic dyad, Each of these residues was changed into the corresponding amide by site-directed mutagenesis. Strains of A. tumefaciens that express mutated VirB1 proteins have a significantly reduced virulence. We hypothesize that many bacterial proteins involved in export of macromolecules belong to a widespread class of hydrolases and cleave beta-1,4-glycosidic bonds as part of their function.

Conversion of cucumber linoleate 13-lipoxygenase to a 9-lipoxygenating species by site-directed mutagenesis

Multiple lipoxygenase sequence alignments and structural modeling of the enzyme/substrate interaction of the cucumber lipid body lipoxygenase suggested histidine 608 as the primary determinant of positional specificity. Replacement of this amino acid by a less-space-filling valine altered the positional specificity of this linoleate 13-lipoxygenase in favor of 9-lipoxygenation. These alterations may be explained by the fact that H608V mutation may demask the positively charged guanidino group of R758, which, in turn, may force an inverse head-to-tail orientation of the fatty acid substrate. The R758L+H608V double mutant exhibited a strongly reduced reaction rate and a random positional specificity. Trilinolein, which lacks free carboxylic groups, was oxygenated to the corresponding (13S)-hydro(pero)xy derivatives by both the wild-type enzyme and the linoleate 9-lipoxygenating H608V mutant. These data indicate the complete conversion of a linoleate 13-lipoxygenase to a 9-lipoxygenating species by a single point mutation. It is hypothesized that H608V exchange may alter the orientation of the substrate at the active site and/or its steric configuration in such a way that a stereospecific dioxygen insertion at C-9 may exclusively take place.

From cell death to embryo arrest: Mathematical models of human preimplantation embryo development

Human preimplantation embryos exhibit high levels of apoptotic cells and high rates of developmental arrest during the first week in vitro. The relation between the two is unclear and difficult to determine by conventional experimental approaches, partly because of limited numbers of embryos. We apply a mixture of experiment and mathematical modeling to show that observed levels of cell death can be reconciled with the high levels of embryo arrest seen in the human only if the developmental competence of embryos is already established at the zygote stage, and environmental factors merely modulate this. This suggests that research on improving in vitro fertilization success rates should move from its current concentration on optimizing culture media to focus more on the generation of a healthy zygote and on understanding the mechanisms that cause chromosomal and other abnormalities during early cleavage stages.

Redesigning secondary structure to invert coenzyme specificity in isopropylmalate dehydrogenase.

Rational engineering of enzymes involves introducing key amino acids guided by a knowledge of protein structure to effect a desirable change in function. To date, all successful attempts to change specificity have been limited to substituting individual amino acids within a protein fold. However, the infant field of protein engineering will only reach maturity when changes in function can be generated by rationally engineering secondary structures. Guided by x-ray crystal structures and molecular modeling, site-directed mutagenesis has been used to systematically invert the coenzyme specificity of Thermus thermophilus isopropylmalate dehydrogenase from a 100-fold preference for NAD to a 1000-fold preference for NADP. The engineered mutant, which is twice as active as wild type, contains four amino acid substitutions and an alpha-helix and loop that replaces the original beta-turn. These results demonstrate that rational engineering of secondary structures to produce enzymes with novel properties is feasible.

A highly active decarboxylating dehydrogenase with rationally inverted coenzyme specificity.

The isocitrate dehydrogenase of Escherichia coli, which lacks the Rossmann fold common to other dehydrogenases, displays a 7000-fold preference for NADP over NAD (calculated as the ratio of kcat/Km). Guided by x-ray crystal structures and molecular modeling, site-directed mutagenesis has been used to introduce six substitutions in the adenosine binding pocket that systematically shift coenzyme preference toward NAD. The engineered enzyme displays an 850-fold preference for NAD over NADP, which exceeds the 140-fold preference displayed by a homologous NAD-dependent enzyme. Of the six mutations introduced, only one is identical in all related NAD-dependent enzyme sequences--strict adherence to homology as a criterion for replacing these amino acids impairs function. Two additional mutations at remote sites improve performance further, resulting in a final mutant enzyme with kinetic characteristics and coenzyme preference comparable to naturally occurring homologous NAD-dependent enzymes.

Toward the ultimate synthesis/recognition system.

This paper predicts speech synthesis, speech recognition, and speaker recognition technology for the year 2001, and it describes the most important research problems to be solved in order to arrive at these ultimate synthesis and recognition systems. The problems for speech synthesis include natural and intelligible voice production, prosody control based on meaning, capability of controlling synthesized voice quality and choosing individual speaking style, multilingual and multidialectal synthesis, choice of application-oriented speaking styles, capability of adding emotion, and synthesis from concepts. The problems for speech recognition include robust recognition against speech variations, adaptation/normalization to variations due to environmental conditions and speakers, automatic knowledge acquisition for acoustic and linguistic modeling, spontaneous speech recognition, naturalness and ease of human-machine interaction, and recognition of emotion. The problems for speaker recognition are similar to those for speech recognition. The research topics related to all these techniques include the use of articulatory and perceptual constraints and evaluation methods for measuring the quality of technology and systems.