Estimating equations for parameters in means and covariances of multivariate discrete and continuous responses.

Generalized estimating equations are introduced in an ad hoc fashion for the covariance matrix of a multivariate response. These equations are to be solved jointly with score equations from a generalized linear model for mean parameters. A class of quadratic exponential models is used to develop joint estimating equations for mean and covariance parameters in a more systematic fashion, and proposals for the use of such equations are developed. Comments on the relative merits of the ad hoc and model-based approaches to estimation are given and a regression illustration with a bivariate response is provided.     

Large scale genotype-phenotype correlation analysis based on phylogenetic trees

We provide two methods for identifying changes in genotype that are correlated with changes in a phenotype implied by phylogenetic trees. The first method, VENN, works when the number of branches over which the change occurred are modest. VENN looks for genetic changes that are completely penetrant with phenotype changes on a tree. The second method, CCTSWEEP, allows for a partial matching between changes in phenotypes and genotypes and provides a score for each change using Maddison's concentrated changes test. The mutations that are highly correlated with phenotypic change can be ranked by score. We use these methods to find SNPs correlated with resistance to Bacillus anthracis in inbred mouse strains. Our findings are consistent with the current biological literature, and also suggest potential novel candidate genes.     

A theoretical treatment of interval mapping of a disease gene using transmission disequilibrium tests

The genetic basis of the transmission disequilibrium test (TDT) for two-marker loci is explored from first principles. In this case, parents doubly heterozygous for a given haplotype at the pair of marker loci that are each in linkage disequilibrium with the disease gene with the further possibility of a second-order linkage disequilibrium are considered. The number of times such parents transmit the given haplotype to their affected offspring is counted and compared with the frequencies of haplotypes that are not transmitted. This is done separately for the coupling and repulsion phases of doubly heterozygous genotypes. Expectations of the counts for each of the sixteen cells possible with four-marker gametic types (transmitted vs not transmitted) are derived. Based on a test of symmetry in a square 4 x 4 contingency table, chi-square tests are proposed for the null hypothesis of no linkage between the markers and the disease gene. The power of the tests is discussed in terms of the corresponding non-centrality parameters for the alternative hypothesis that both the markers are linked with the disease locus. The results indicate that the power increases with the decrease in recombination probability and that it is higher for a lower frequency of the disease gene. Taking a pair of markers in an interval for exploring the linkage with the disease gene seems to be more informative than the single-marker case since the values of the non-centrality parameters tend to be consistently higher than their counterparts in the single-marker case. Limitations of the proposed test are also discussed.     

Differential elution of histones from gel-trapped nuclei

A method has been developed for characterizing at, submicrogram levels the heterogeneity of histones from purified nuclei. The histones are eluted with a smooth concentration gradient from nuclei trapped in polyacrylamide–gel threads and are collected in a micro fraction collector suitable for volumes in the 10–100 μl range. The gradient and fraction collection systems are governed by cam driven syringes. Samples obtained are subjected to electrophoresis in a starch-gel system and the gels are stained with a highly sensitive stain specific for guanidinium groups. Seven major and a similar number of minor components are demonstrated in the histones. The method of differential elution of trapped macromolecules is suitable for use with systems other than nuclei and histones.     

DAPI (4',6-diamidino-2-phenylindole) binds differently to DNA and RNA: minor-groove binding at AT sites and intercalation at AU sites.

The interaction of DAPI and propidium with RNA (polyA.polyU) and corresponding DNA (polydA.polydT) sequences has been compared by spectroscopic, kinetic, viscometric, Tm, and molecular modeling methods. Spectral changes of propidium are similar on binding to the AT and AU sequences but are significantly different for binding of DAPI. Spectral changes for DAPI with the DNA sequence are consistent with the expected groove-binding mode. All spectral changes for complexes of propidium with RNA and DNA and for DAPI with RNA, however, are consistent with an intercalation binding mode. When complexed with RNA, for example, DAPI aromatic protons signals shift significantly upfield, and the DAPI UV-visible spectrum shows significantly larger changes than when complexed with DNA. Slopes of log kd (dissociation rate constants) versus-log [Na+] plots are similar for complexes of propidium with RNA and DNA and for the DAPI-RNA complex and are in the range expected for an intercalation complex. The slope for the DAPI-DNA complex, however, is much larger and is in the range expected for a groove-binding complex. Association kinetics results also support an intercalation binding mode for the DAPI-RNA complex. The viscosity of polyA.polyU solutions increases significantly on addition of both propidium and DAPI, again in agreement with an intercalation binding mode for both molecules with RNA. Molecular modeling studies completely support the experimental findings and indicate that DAPI forms a very favorable intercalation complex with RNA. DAPI also forms a very stable complex in the minor groove of AT sequences of DNA, but the stabilizing interactions are considerably reduced in the wide, shallow minor groove of RNA. Modeling studies,thus,indicate that DAPI interaction energetics are more favorable for minor-groove binding in AT sequences but are more favorable for interaction in RNA.     

Optical characteristics of flowing blood

The transmitted light strength (TS) through a thin blood layer changes with variation in blood flow, such as positive streaming transparency for low hematocrits and negative streaming transparency for high hematocrits. These phenomena are examined theoretically and experimentally. Maxwell’s equations are solved assuming that erythrocytes are oblate spheroids to investigate these phenomena due to flowing blood. The theoretical results reveal that the scattering and absorption cross sections for flowing blood are larger than those for stagnant blood. Experimental results indicate that the TS for both oxygenated and deoxygenated flowing blood, with a hematocrit of up to approximately 20%, was stronger than that for stagnant blood. The TS decreased for flowing blood with a hematocrit of approximately 20% or greater. Applying the theoretical scattering and absorption cross sections to the absorption and multiple scattering theory of Victor Twersky, the changes in the TS due to flowing blood are obtained theoretically. From the theoretical and experimental results, the positive streaming transparency phenomenon of flowing blood with a low hematocrit and the negative streaming transparency phenomenon with a high hematocrit are found to result from increased scattering and absorption cross sections because of the orientation of flowing erythrocytes.     

POPULATION GENETICS OF INTRAGAMETOPHYTIC SELFING

Intragametophytic selfing is a mode of reproduction occurring in homosporous ferns where two gametes from the same haploid gametophyte form a completely homozygous sporophyte. The inbreeding equilibrium is derived for a population with partial intragametophytic selfing, selfing, and outcrossing. Procedures for directly estimating the extent of intragametophytic selfing and selfing using parent-offspring data are given. The conditions for a stable polymorphism from a heterozygous-advantage fitness model are more restrictive for partial intragametophytic selfing than for selfing. The rate of decay of gametic disequilibrium is slower for partial intragametophytic selfing than for selfing. Based on these findings, one would predict that plants with intragametophytic selfing would have less polymorphism for loci with a heterozygous advantage and more gametic disequilibrium between neutral loci than is expected for populations with an equivalent amount of selfing. Data from several studies are consistent with these predictions.     

Inferring qualitative relations in genetic networks and metabolic pathways

MOTIVATION: Inferring genetic network architecture from time series data of gene expression patterns is an important topic in bioinformatics. Although inference algorithms based on the Boolean network were proposed, the Boolean network was not sufficient as a model of a genetic network. RESULTS: First, a Boolean network model with noise is proposed, together with an inference algorithm for it. Next, a qualitative network model is proposed, in which regulation rules are represented as qualitative rules and embedded in the network structure. Algorithms are also presented for inferring qualitative relations from time series data. Then, an algorithm for inferring S-systems (synergistic and saturable systems) from time series data is presented, where S-systems are based on a particular kind of nonlinear differential equation and have been applied to the analysis of various biological systems. Theoretical results are shown for Boolean networks with noises and simple qualitative networks. Computational results are shown for Boolean networks with noises and S-systems, where real data are not used because the proposed models are still conceptual and the quantity and quality of currently available data are not enough for the application of the proposed methods.     

Application of particle swarm optimization and proximal support vector machines for fault detection

This paper presents a novel application of particle swarm optimization (PSO) in combination with another computational intelligence (CI) technique, namely, proximal support vector machine (PSVM) for machinery fault detection. Both real-valued and binary PSO algorithms have been considered along with linear and nonlinear versions of PSVM. The time domain vibration signals of a rotating machine with normal and defective bearings are processed for feature extraction. The features extracted from original and preprocessed signals are used as inputs to the classifiers (PSVM) for detection of machine condition. Input features are selected using a PSO algorithm. The classifiers are trained with a subset of experimental data for known machine conditions and are tested using the remaining data. The procedure is illustrated using the experimental vibration data of a rotating machine. The influences of the number of features, PSO algorithms and type of classifiers (linear or nonlinear PSVM) on the detection success are investigated. Results are compared with a genetic algorithm (GA) and principal component analysis (PCA). The PSO based approach gave test classification success above 90% which were comparable with the GA and much better than PCA. The results show the effectiveness of the selected features and classifiers in detection of machine condition.     

Kinetic properties of the binding of alpha-lytic protease to peptide boronic acids

The kinetic parameters for peptide boronic acids in their interaction with alpha-lytic protease were determined and found to be similar to those of other serine proteases [Kettner, C., & Shenvi, A. B. (1984) J. Biol. Chem. 259, 15106-15114]. alpha-Lytic protease hydrolyzes substrates with either alanine or valine in the P1 site and has a preference for substrate with a P1 alanine. The most effective inhibitors are tri- and tetrapeptide analogues that have a -boroVal-OH residue in the P1 site. At pH 7.5, MeOSuc-Ala-Ala-Pro-boroVal-OH has a Ki of 6.4 nM and Boc-Ala-Pro-boroVal-OH has a Ki of 0.35 nM. Ac-boroVal-OH and Ac-Pro-boroVal-OH are 220,000- and 500-fold less effective, respectively, than the tetrapeptide analogue. The kinetic properties of the tri- and tetrapeptide analogues are consistent with the mechanism for slow-binding inhibition, E + I in equilibrium EI in equilibrium EI*, while the less effective inhibitors are simple competitive inhibitors. MeO-Suc-Ala-Ala-Pro-boroAla-OH is a simple competitive inhibitor with a Ki of 67 nM at pH 7.5. Other peptide boronic acids, which are analogues of nonsubstrates, are less effective than substrate analogues but still are effective competitive inhibitors. For example, MeOSuc-Ala-Ala-Pro-boroPhe-OH has a Ki of 0.54 microM although substrates with a phenylalanine in the P1 position are not hydrolyzed. Binding for boronic acid analogues of both substrate and nonsubstrate analogues is pH dependent with higher affinity near pH 7.5. Similar binding properties have been observed for pancreatic elastase. Both enzymes have almost identical requirements for an extended peptide inhibitor sequence in order to exhibit highly effective binding and slow-binding characteristics.(ABSTRACT TRUNCATED AT 250 WORDS)     

Efficient methods for filtering and ranking fragments for the prediction of structurally variable regions in proteins

The prediction of protein 3D structures close to insertions and deletions or, more generally, loop prediction, is still one of the major challenges in homology modeling projects. In this article, we developed ranking criteria and selection filters to improve knowledge-based loop predictions. These criteria were developed and optimized for a test data set containing 678 insertions and deletions. The examples are, in principle, predictable from the used loop database with an RMSD < 1 A and represent realistic modeling situations. Four noncorrelated criteria for the selection of fragments are evaluated. A fast prefilter compares the distance between the anchor groups in the template protein with the stems of the fragments. The RMSD of the anchor groups is used for fitting and ranking of the selected loop candidates. After fitting, repulsive close contacts of loop candidates with the template protein are used for filtering, and fragments with backbone torsion angles, which are unfavorable according to a knowledge-based potential, are eliminated. By the combined application of these filter criteria to the test set, it was possible to increase the percentage of predictions with a global RMSD < 1 A to over 50% among the first five ranks, with average global RMSD values for the first rank candidate that are between 1.3 and 2.2 A for different loop lengths. Compared to other examples described in the literature, our large numbers of test cases are not self-predictions, where loops are placed in a protein after a peptide loop has been cut out, but are attempts to predict structural changes that occur in evolution when a protein is affected by insertions and deletions.     

Bacteriophage and Bacteria in a Flow Reactor

The Levin-Stewart model of bacteriophage predation of bacteria in a chemostat is modified for a flow reactor in which bacteria are motile, phage diffuse, and advection brings fresh nutrient and removes medium, cells and phage. A fixed latent period for phage results in a system of delayed reaction-diffusion equations with non-local nonlinearities. Basic reproductive numbers are obtained for bacteria and for phage which predict survival of each in the bio-reactor. These are expressed in terms of physical and biological parameters. Persistence and extinction results are obtained for both bacteria and phage. Numerical simulations are in general agreement with those for the chemostat model.     

Detection of mutations in human DNA

Efficient methods for the detection of mutations are of fundamental importance in research and in diagnostics. By detection of a DNA sequence alteration that cosegregates with a clinical phenotype in an affected family, the gene at fault may be identified and assigned a function. Mutation detection methods are also a rate-limiting factor for the clinical application of DNA diagnostics. Currently a large number of techniques are in use to scan for new mutations and to distinguish among previously established sequence variants. Here, some of the problems connected with mutation detection are discussed together with principles on which current and future mutation detection assays can be based.     

Study of the surface barrier discharge generated by electrodes in the form of a series of parallel metal strips

A surface discharge in a system where metal electrodes in the form of a series of parallel strips are positioned on the dielectric surface is studied. Analytical formulas for calculating the spatial distribution of the potential and the electric field in a discharge cell are derived. It is shown that the geometry of the metal electrodes should be taken into account (along with physical and chemical characteristics of the dielectric, the voltage applied to the electrodes, and other parameters of the system) for generation of the electric field with optimal configuration in the discharge cell. The obtained results are also applicable for analysis of discharge cells with a coplanar barrier discharge where metal electrodes are positioned in the dielectric at small depths. The results are of interest since a barrier discharge is one of the efficient methods for generating non-equilibrium plasma at high pressures for a variety of technological applications.     

Cephalic morphology of Hymenopus coronatus (Insecta: Mantodea) and its phylogenetic implications

External and internal head structures of the mantodean Hymenopus coronatus are examined and described in detail. The results are elaborately compared with the literature. Strong crests on the anterior tentorial arms that articulate with the subantennal suture, a parietal suture and glossae and paraglossae with anteriorly bent tips are proposed as new potential apomorphies for Mantodea while a head capsule being wider than long, enlarged compound eyes, the presence of a frontal shield or scutellum, lateral lobes in the anterior tentorial arms, the presence of a transverse and an interantennal suture and the reduction of the mentum are confirmed as apomorphies, As potential apomorphies for Dictyoptera the reduction of Musculus tentoriobuccalis lateralis (M. 49) is newly presented and a “perforate” tentorium, lacinial incisivi that are located in a galeal pouch and the presence of a postmola are confirmed. The present study shows the value of cephalic morphology for phylogenetic analysis but also points out that further studies including evolutionary key taxa are essential for resolving the evolutionary adaptations among dictyopterans.     

Bull breeding soundness evaluation in Southern Africa

The motivation for and process leading up to the publication of a new bull breeding soundness certification standard endorsed by the South African Veterinary Association is described. The veterinary certificate of bull breeding soundness and explanatory notes and minimum standards are shown. The first component of the certificate is a declaration by the veterinarian that the bull complies with the minimum standards set for examinations for the selected purpose, these being for use as a natural service sire, as a donor of semen for distribution, and for insurance purposes. This is followed by the details of the bull and owner, and a list of the recommended examinations and tests for the bull with provision for which were performed. Certificates are available in book form with the explanatory notes and minimum standards on the reverse, and a carbon copy which remains in the book.The clarity and ease of completion of the document are regarded as being positive features. Bulls are either classified as breeding sound or not, with no actual parameters indicated on the document and no certificate issued for those which do not meet the set criteria. Contact details of the parties involved are shown on the certificate to allow for communication as a means of avoiding disputes.     

Process technology for production and recovery of heterologous proteins with Pichia pastoris

Developments in process techniques for production and recovery of heterologous proteins with Pichia pastoris are presented. Limitations for the standard techniques are described, and alternative techniques that solve the limitations problems are reviewed together with the methods that resulted in higher productivity of the P. pastoris processes. The main limitations are proteolysis of the secreted products and cell death in the high cell density bioreactor cultures. As a consequence, both low productivity and lower quality of the feedstock for downstream processing are achieved in processes hampered with these problems. Methods for exploring proteolysis and cell death are also presented. Solving the problems makes the conditions for downstream processing superior for the P. pastoris expression systems compared to other systems, which either need complex media or rely on intracellular production. These improved conditions allow for interfacing of cultivation with downstream processing in an integrated fashion.     

The importance of being exhaustive. Optimization of bridging structural water molecules and water networks in models of biological systems

Algorithms and protocols are described for the optimization for H-bonding of isolated singular H2O molecules and entire networks of H2O molecules. Unlike other approaches that are prone to being trapped in local energy minima, these methods rely on exhaustive searches of orientation space for the H2O molecules. The results are scored with the HINT hydropathic interaction model, but the algorithms should be general for any energy-scoring computation. Two examples are provided: 1) the tightly-bound H2O molecule 301 of HIV-1 protease is shown to be more reasonably oriented in terms of forming H-bonds with this method than with a molecular mechanics energy minimization method; and 2) the H2O network surrounding carbonmonoxymyoglobin is constructed and analyzed for a 1.80-A neutron-diffraction structure. The H-atom positions calculated with this method show a somewhat better agreement with the experimental results than do the H-atom positions calculated with molecular mechanics, and both are considerably better than random.     

A theoretical treatment of interval mapping of a disease gene using transmission disequilibrium tests

The genetic basis of the transmission disequilibrium test (TDT) for two-marker loci is explored from first principles. In this case, parents doubly heterozygous for a given haplotype at the pair of marker loci that are each in linkage disequilibrium with the disease gene with the further possibility of a second-order linkage disequilibrium are considered. The number of times such parents transmit the given haplotype to their affected offspring is counted and compared with the frequencies of haplotypes that are not transmitted. This is done separately for the coupling and repulsion phases of doubly heterozygous genotypes. Expectations of the counts for each of the sixteen cells possible with four-marker gametic types (transmitted vs not transmitted) are derived. Based on a test of symmetry in a square 4 × 4 contingency table, chi-square tests are proposed for the null hypothesis of no linkage between the markers and the disease gene. The power of the tests is discussed in terms of the corresponding non-centrality parameters for the alternative hypothesis that both the markers are linked with the disease locus. The results indicate that the power increases with the decrease in recombination probability and that it is higher for a lower frequency of the disease gene. Taking a pair of markers in an interval for exploring the linkage with the disease gene seems to be more informative than the single-marker case since the values of the non-centrality parameters tend to be consistently higher than their counterparts in the single-marker case. Limitations of the proposed test are also discussed.     

Construction of Selection Indexes with Restrictions

In this paper the problem of unistage selection with inequality constraints is formulated. If the predictor and criterion variables are all normally distributed, this problem can be written as a convex programming problem, with a linear objective function and with linear constraints and a quadratic constraint. Using the duality theory, for convex nonlinear programming it is proved, that its dual problem can be transformed into a convex minimization problem with non-negativity conditions. Good computational methods are known for solving this problem. By the help of the dual problem sufficient conditions for a solution of the original primal problem are derived and illustrated by an example of practical interest.