Computational characterization of the molecular structure and properties of Dye 7 for organic photovoltaics

Organic dyes have great potential for its use in solar cells. In this recent work, the molecular structure and properties of Dye 7 were obtained using density functional theory (DFT) and different levels of calculation. Upon comparing the molecular structure and the ultraviolet visible spectrum with experimental data reported in the literature, it was found that the M05-2X/6-31G(d) level of calculation gave the best approximation. Once the appropriate methodology had been obtained, the molecule was characterized by obtaining the infrared spectrum, dipole moment, total energy, isotropic polarizability, molecular orbital energies, free energy of solvation in different solvents, and the chemical reactivity sites using the condensed Fukui functions.     

Optimal reconstruction of a sequence from its probes.

An important combinatorial problem, motivated by DNA sequencing in molecular biology, is the reconstruction of a sequence over a small finite alphabet from the collection of its probes (the sequence spectrum), obtained by sliding a fixed sampling pattern over the sequence. Such construction is required for Sequencing-by-Hybridization (SBH), a novel DNA sequencing technique based on an array (SBH chip) of short nucleotide sequences (probes). Once the sequence spectrum is biochemically obtained, a combinatorial method is used to reconstruct the DNA sequence from its spectrum. Since technology limits the number of probes on the SBH chip, a challenging combinatorial question is the design of a smallest set of probes that can sequence an arbitrary DNA string of a given length. We present in this work a novel probe design, crucially based on the use of universal bases [bases that bind to any nucleotide (Loakes and Brown, 1994)] that drastically improves the performance of the SBH process and asymptotically approaches the information-theoretic bound up to a constant factor. Furthermore, the sequencing algorithm we propose is substantially simpler than the Eulerian path method used in previous solutions of this problem.     

Positional statistical significance in sequence alignment.

Beginning with the concept of near-optimal sequence alignments, we can assign a probability that each element in one sequence is paired in an alignment with each element in another sequence. This involves a sum over the set of all possible pairwise alignments. The method employs a designed hidden Markov model (HMM) and the rigorous forward and forward-backward algorithms of Rabiner. The approach can use any standard sequence-element-to-element probabilistic similarity measures and affine gap penalty functions. This allows the positional alignment statistical significance to be obtained as a function of such variables. A measure of the probabilistic relationship between any single sequence and a set of sequences can be directly obtained. In addition, the employed HMM with the Viterbi algorithm provides a simple link to the standard dynamic programming optimal alignment algorithms.     

Isolation and characterization of the nuclear matrix in Friend erythroleukemia cells: chromatin and hnRNA interactions with the nuclear matrix.

Nuclear matrices from undifferentiated and differentiated Friend erythroleukemia cells have been obtained by a method which removes DNA in a physiological buffer. These matrices preserved the characteristic topographical distribution of condensed and diffuse "chromatin" regions, as do nuclei in situ or isolated nuclei. Histone H1 was released from the nuclear matrix of undifferentiated cells by 0.3 M KCl; inner core histones were released by 1 M KCl. Nuclear matrix from differentiated cells did not maintain H1, and histone cores were fully released in 0.7 M KCl. KCl removed the core histones as an octameric structure with no evidence of preferential release of any single histone. Electron microscopy of KCl-treated matrix revealed no condensed regions but rather a network of fibrils in the whole DNA-depleted nuclei. When nuclear matrices from both types of cell were exposed to conditions of very low ionic strength, inner core histones and condensed regions remained. These observations support the contention that inner core histones are bound to matrix through natural ionic bonds or saline-labile elements, and that these interactions are implicated in chromatin condensation. hnRNA remained undegraded and tenaciously associated to the matrix fibrils, and was released only by chemical means which, by breaking hydrophobic and hydrogen bonds, produced matrix lysis. Very few nonhistone proteins were released upon complete digestion of DNA from either type of nuclei. The remaining nonhistone proteins represent a large number of species of which the majority may be matrix components. The molecular architecture in both condensed and diffuse regions of interphase nuclei appears to be constructed of two distinct kinds of fibers; the thicker chromatin fibers are interwoven with the thinner matrix fibers. The latter are formed by a heteropolymer of many different proteins.     

Structural bioinformatics of DNA: a web-based tool for the analysis of molecular dynamics results and structure prediction

We report here the release of a web-based tool (MDDNA) to study and model the fine structural details of DNA on the basis of data extracted from a set of molecular dynamics (MD) trajectories of DNA sequences involving all the unique tetranucleotides. The dynamic web interface can be employed to analyze the first neighbor sequence context effects on the 10 unique dinucleotide steps of DNA. Functionality is included to build all atom models of any user-defined sequence based on the MD results. The backend of this interface is a relational database storing the conformational details of DNA obtained in 39 different MD simulation trajectories comprising all the 136 unique tetranucleotide steps. Examples of the use of this data to predict DNA structures are included. Availability: http://humphry.chem.wesleyan.edu:8080/MDDNA. Supplementary information: Supplementary data including color figures are available at Bioinformatics online.     

Methane in Water: An Ab Initio Study

Abstract

In preceding publications we discussed some properties of pure water in condensed phases using an ab initio approach. Here this study is used as a basis of comparison for analysing the behaviour of water as a solvent in the presence of an apolar molecule. Our analysis is focused on the process of organization of the hydrogen bonding network around the solute. For this purpose we perform some ab initio calculations for a system of 32 water molecules and one methane molecule at 300 K; in particular, the average molecular dipole moment of water is determined and the result is compared with that of pure water. Next the attention is switched to the methane molecule; related properties such as excluded volume and sphericity of its shape are illustrated and discussed. A comparison with results obtained using classical approaches suggests that some classical models of water can be considered to be still valid when they are used to analyse the water-methane system.     

Imaging system for morphometric assessment of absorption or fluorescence in stained cells

An image acquisition and processing system has been developed for quantitative microscopy of absorption or fluorescence in stained cells. Three different light transducers are used in the system to exploit the best characteristics of these sensors for different biological measurements. A digital scanner, in the form of a linear array charge-coupled device (CCD), acquires data with high spatial and photometric resolution. A color (RGB) camera is employed when spectral information is required for the segmentation of cellular subcomponents. An image-intensified charged-injection device (CID) camera provides for very low light intensity measurements, primarily for fluorescence-labeled cells. Properties of these transducers, such as contrast transfer function, linearity, and photo-response nonuniformity, have been measured. Two dedicated image processing units were incorporated into the system. The front-end processor, based on a digital signal processor, provides functions such as object detection, raw image calibration, compression, artifact removal, and filtering. The second image processor is associated with the frame memory and includes a histogram processor, a dedicated arithmetic logic unit for image processing functions, and a graphics module for one-bit overlay functions. An interactive program was developed to acquire cell images and to experiment with a range of segmentation algorithms, feature extractions, and other image processing functions. The results of any image operation are displayed on the video monitor. Once a desired processing sequence is determined, the sequence may be stored to become part of a command library and can be executed thereafter as a single instruction.     

Electron Correlated Ab Initio Study of Amino Group Flexibility for Improvement of Molecular Mechanics Simulations on Nucleic Acid Conformations and Interactions

High level ab initio studies demonstrate substantial conformational flexibility of amino groups of nucleic acid bases. This flexibility is important for biological functions of DNA. Existing force field models of molecular mechanics do not describe this phenomenon due to a lack of quantitative experimental data necessary for an adjustment of empirical parameters. We have performed extensive calculations of nucleic acid bases at the MP2/6-31G(d,p) level of ab initio theory for broad set of amino group configurations. Two-dimensional maps of energy and geometrical characteristics as functions of two amino hydrogen torsions have been constructed. We approximate the maps by polynomial expressions, which can be used in molecular mechanics calculations. Detailed considerations of these maps enable us to propose a method for determination of numerical coefficients in the developed formulae using restricted sets of points obtained via higher-level calculations.     

Genetic Distances Based on Quantitative Traits

Morphological data showing continuous distributions, polygenically controlled, may be particularly useful in intergroup classification below the species level; an appropriate distance analysis based on these traits is an important tool in evolutionary biology and in plant and animal breeding.--The interpretation of morphological distances in genetic terms is not easy because simple phenotypic data may lead to biased estimates of genetic distances. Convenient estimates can be obtained whenever it is possible to breed populations according to a suitable crossing design and to derive information from genetic parameters.--A general method for determining genetic distances is proposed. The procedure of multivariate analysis of variance is extended to estimate appropriate genetic parameters (genetic effects). Not only are optimal statistical estimates of parameters obtained but also the procedure allows the measurement of genetic distances between populations as linear functions of the estimated parameters, providing an appropriate distance matrix that can be defined in terms of these parameters. The use of the T2 statistic, defined in terms of the vector of contrasts specifying the distance, permits the testing of the significance of any distance between any pair of populations that may be of interest from a genetic point of view.--A numerical example from maize diallel data is reported in order to illustrate the procedure. In particular, heterosis effects are used as the basis for estimates of genetic divergence between populations.     

用磁珠富集法从AFLP片段中分离微卫星DNA标记

将花生(Arachis hypogaea L．)基因组DNA经酶切后转变成AFLP DNA片段,然后用生物素标记的简单重复序列(SSR)作探针与其杂交,杂交复合物固定到包被有链亲和素的磁珠上,经过一系列的洗涤过程,含有SSR的AFLP片段被吸附在磁珠表面。这些片段经洗脱下来后,先用对应的AFLP引物扩增,再进行克隆和测序,根据SSR两端的保守序列设计引物,经过多态性分析后,便可得到微卫星DNA标记。整个实验过程操作简单、消耗少,可在一周内完成,可作为从植物中分离SSR的一种简单有效的方法。     

Assessing the significance of local sequence homologies

The extent of homology between two sequences is generally expressed quantitatively by using a set of rules to assign to aligned residue pairs, numerical values that depend on some measure of the similarity between the pairs. A homology score for the alignment is obtained by summing the numerical value for each pair, and possibly also adding some appropriate penalty for insertions and deletions in the alignment. Whatever the set of rules used in arriving at the best homology score, the fundamental question that remains is whether the score implies that the sequences are in some way related. In the absence of biological criteria, such a question is generally given a preliminary answer by a Monte Carlo evaluation of the statistical significance of the score. By randomizing the two sequences a large number of times and noting the homology on each throw, the probability of a given score can be obtained for each sequence, and hence the likelihood that the observed homology could have occurred by chance for sequences with the given composition, is easily assessed. In addition to this statistical assessment, however, there is a second statistical question that must also be evaluated—particularly for short homologous stretches (<15 residues). When homology searches are performed against databases containing several hundred thousand residues, an important number to know is the probability that the homologous sequence would have occurred simply as the result of the large number of arrangements of short sequences that must be present in any large collection of disparate sequences. In this note we evaluate different versions of this question using different sets of rules. Our results indicate that homologies of roughly six or more residues out of ten are statistically significant and that the best procedure for finding homologies and determining their significance is with the use of the Dayhoff mutation matrix, assigning a weight of at least —8 as a penalty for gaps.     

A space-variant differential operator for visual sensitivity

All the elements of a Fourier analysis can be derived from the experiments of Graham and Robson on contrast sensitivity. Once their experiment is posed as an eigenvalue problem, a complete orthonormal set of eigenfunctions results from solving the associated differential equation. Neither sine and cosine nor Gabor functions result. Instead, the Hermite functions arise as the eigenfunctions of a space-variant differential operator used to model the contrast sensitivity of human observers. These functions, up to a constant, are their own Fourier transforms, and in principle can be used to exactly represent the Fourier transform of naturally occuring visual images.     

The translation of DNA primary base sequence into three-dimensional structure

A procedure is outlined to obtain a reliable computer-generatedrepresentation of the DNA duplex from its primary sequence ofbase pairs. The calculations are based on the potential energiesof interaction of adjacent side groups. The methods are, however,completely general and can be adapted to any set of base sequencedependent conformational rules. Static representations of theDNA are compared with the distributions of conformations obtainedfrom Monte Carlo simulation studies. Direct matrix generatorcalculations of the average (equilibrium) extension and orientationof various sequences and numerical estimates of the flexibilityof the chains as a whole are also reported. The methods areapplied to three short fragments of kinetoplast DNA from Crithidiafasciculata which exhibit dramatically different behavior onnon-denaturing poly-acrylamide gels. Received on August 17, 1987; accepted on December 19, 1987     

Virtual and solution conformations of oligosaccharides

The possibility that observed nuclear Overhauser enhancements and bulk longitudinal relaxation times, parameters measured by 1H NMR and often employed in determining the preferred solution conformation of biologically important molecules, are the result of averaging over many conformational states is quantitatively evaluated. Of particular interest was to ascertain whether certain 1H NMR determined conformations are "virtual" in nature; i.e., the fraction of the population of molecules actually found at any time within the subset of conformational space defined as the "solution conformation" is vanishingly small. A statistical mechanics approach was utilized to calculate an ensemble average relaxation matrix from which (NOE)'s and (T1)'s are calculated. Model glycosidic linkages in four oligosaccharides were studied. The solution conformation at any glycosidic linkage is properly represented by a normalized, Boltzmann distribution of conformers generated from an appropriate potential energy surface. The nature of the resultant population distributions is such that 50% of the molecular population is found within 1% of available microstates, while 99% of the molecular population occupies about 10% of the ensemble microstates, a number roughly equal to that sterically allowed. From this analysis we conclude that in many cases quantitative interpretation of NMR relaxation data, which attempts to define a single set of allowable torsion angle values consistent with the observed data, will lead to solution conformations that are either virtual or reflect torsion angle values possessed by a minority of the molecular population. On the other hand, calculation of ensemble average NMR relaxation data yields values in agreement with experimental results. Observed values of NMR relaxation data are the result of the complex interdependence of the population distribution and NOE (or T1) surfaces in conformational space. In conformational analyses, NMR data can therefore be used to test different population distributions calculated from empirical potential energy functions.     

Calmodulin Target Database

The intracellular calcium sensor protein calmodulin (CaM) interacts with a large number of proteins to regulate their biological functions in response to calcium stimulus. This molecular recognition process is diverse in its mechanism, but can be grouped into several classes based on structural and sequence information. We have developed a web-based database (http://calcium.uhnres.utoronto.ca/ctdb) for this family of proteins containing CaM binding sites or, as we propose to call it herein, CaM recruitment signaling (CRS) motifs. At present the CRS motif found in approximately 180 protein sequences in the databases can be divided into four subclasses, each subclass representing a distinct structural mode of molecular recognition involving CaM. The database can predict a putative CRS location within a given protein sequence, identify the subclass to which it may belong, and structural and biophysical parameters such as hydrophobicity, hydrophobic moment, and propensity for a -helix formation. This revised version was published online in July 2006 with corrections to the Cover Date.     

Deviations from Michaelis-Menten kinetics. Computation of the probabilities of obtaining complex curves from simple kinetic schemes.

1. It is possible to calculate the intrinsic probability associated with any curve shape that is allowed for rational functions of given degree when the coefficients are independent or dependent random variables with known probability distributions. 2. Computations of such probabilities are described when the coefficients of the rational function are generated according to several probability distribution functions and in particular when rate constants are varied randomly for several simple model mechanisms. 3. It is concluded that each molecular mechanism is associated with a specific set of curve-shape probabilities, and this could be of value in discriminating between model mechanisms. 4. It is shown how a computer program can be used to estimate the probability of new complexities such as extra inflexions and turning points as the degree of rate equations increases. 5. The probability of 3 : 3 rate equations giving 2 : 2 curve shapes is discussed for unrestricted coefficients and also for the substrate-modifier mechanisms. 6. The probability associated with the numerical values of coefficients in rate equations is also calculated for this mechanism, and a possible method for determining the approximate magnitude of product-release steps is given. 7. The computer programs used in the computations have been deposited as Supplement SUP 50113 (21 pages) with the British Library Lending Division, Boston Spa, Wetherby, West Yorkshire LS23 7BQ, U.K., from whom copies can be obtained on the terms indicated in Biochem, J. (1978) 169, 5.     

Application of discrete Fourier inter-coefficient difference for assessing genetic sequence similarity

Digital signal processing (DSP) techniques for biological sequence analysis continue to grow in popularity due to the inherent digital nature of these sequences. DSP methods have demonstrated early success for detection of coding regions in a gene. Recently, these methods are being used to establish DNA gene similarity. We present the inter-coefficient difference (ICD) transformation, a novel extension of the discrete Fourier transformation, which can be applied to any DNA sequence. The ICD method is a mathematical, alignment-free DNA comparison method that generates a genetic signature for any DNA sequence that is used to generate relative measures of similarity among DNA sequences. We demonstrate our method on a set of insulin genes obtained from an evolutionarily wide range of species, and on a set of avian influenza viral sequences, which represents a set of highly similar sequences. We compare phylogenetic trees generated using our technique against trees generated using traditional alignment techniques for similarity and demonstrate that the ICD method produces a highly accurate tree without requiring an alignment prior to establishing sequence similarity.     

An oligodeoxynucleotide affinity column for the isolation of sequence specific DNA binding proteins.

A nucleic acid affinity matrix containing a short oligodeoxynucleotide ligand has been prepared as an example of a material which can be used for the rapid and effective isolation of sequence specific DNA binding proteins. Two complementary oligodeoxynucleotides have been employed, one of which contains a small 5'-spacer arm with a terminal thiol group. Using this terminal thiol group, the ligand can be covalently coupled to Tresyl-activated Sepharose 4B or Epoxy-activated Sepharose 6B via a thioether linkage. This approach allows the specific attachment of the nucleic acid ligand via its 5'-terminus to the insoluble matrix. The double stranded affinity material was obtained by annealing of the complementary DNA fragment. As an example, we have used an eicosomer affinity column containing the sequence d(GAATTC) for the isolation of the Eco RI restriction endonuclease. Using a single column, the enzyme could be isolated by eluting the column with a single step or multistep gradient of increasing salt concentration. The enzyme was purified to 75%-85% homogeneity with yields of 0.1 mg to 0.2 mg from 0.5 g of cell paste.