Vibrational population relaxation of carbon monoxide in the heme pocket of photolyzed carbonmonoxy myoglobin: Comparison of time-resolved mid-IR absorbance experiments and molecular dynamics simulations

The vibrational energy relaxation of carbon monoxide in the heme pocket of sperm whale myoglobin was studied by using molecular dynamics simulation and normal mode analysis methods. Molecular dynamics trajectories of solvated myoglobin were run at 300 K for both the δ- and ɛ-tautomers of the distal His-64. Vibrational population relaxation times of 335 ± 115 ps for the δ-tautomer and 640 ± 185 ps for the ɛ-tautomer were estimated by using the Landau–Teller model. Normal mode analysis was used to identify those protein residues that act as the primary “doorway” modes in the vibrational relaxation of the oscillator. Although the CO relaxation rates in both the ɛ- and δ-tautomers are similar in magnitude, the simulations predict that the vibrational relaxation of the CO is faster in the δ-tautomer with the distal His playing an important role in the energy relaxation mechanism. Time-resolved mid-IR absorbance measurements were performed on photolyzed carbonmonoxy hemoglobin (Hb13CO). From these measurements, a T1 time of 600 ± 150 ps was determined. The simulation and experimental estimates are compared and discussed.

Molecular cloning, sequence, expression, and processing of the interleukin 16 precursor

Interleukin 16 (IL-16) has been shown to function as chemoattractant factor, as a modulator of T-cell activation, and as an inhibitor of immunodeficiency virus replication. The recent identification of inconsistencies in published IL-16 cDNA nucleotide sequences led to the proposal that IL-16 is synthesized in the form of a large precursor protein (pro-IL-16). To identify the true transcriptional start of the IL-16 mRNA rapid amplification of cDNA ends methods were applied. The complete pro-IL-16 cDNA was subsequently molecularly cloned, sequenced, and expressed in COS-7 cells. We report here that pro-IL-16 is most likely synthesized as a 67-kDa protein and is encoded from a major 2.6-kb transcript. Recombinant pro-IL-16 polypeptides are specifically cleaved in lysates of CD8(+) cells, suggesting that the naturally secreted bioactive form of IL-16 is smaller than the originally published 130 amino acids fragment. Moreover, in contrast to other interleukins such as IL-15, IL-16 mRNA expression is almost exclusively limited to lymphatic tissues underlining the potential of IL-16 as an immune regulatory molecule.

Molecular mechanisms underlying differential odor responses of a mouse olfactory receptor

The prevailing paradigm for G protein-coupled receptors is that each receptor is narrowly tuned to its ligand and closely related agonists. An outstanding problem is whether this paradigm applies to olfactory receptor (ORs), which is the largest gene family in the genome, in which each of 1,000 different G protein-coupled receptors is believed to interact with a range of different odor molecules from the many thousands that comprise “odor space.” Insights into how these interactions occur are essential for understanding the sense of smell. Key questions are: (i) Is there a binding pocket? (ii) Which amino acid residues in the binding pocket contribute to peak affinities? (iii) How do affinities change with changes in agonist structure? To approach these questions, we have combined single-cell PCR results [Malnic, B., Hirono, J., Sato, T. & Buck, L. B. (1999) Cell 96, 713–723] and well-established molecular dynamics methods to model the structure of a specific OR (OR S25) and its interactions with 24 odor compounds. This receptor structure not only points to a likely odor-binding site but also independently predicts the two compounds that experimentally best activate OR S25. The results provide a mechanistic model for olfactory transduction at the molecular level and show how the basic G protein-coupled receptor template is adapted for encoding the enormous odor space. This combined approach can significantly enhance the identification of ligands for the many members of the OR family and also may shed light on other protein families that exhibit broad specificities, such as chemokine receptors and P450 oxidases.

Simultaneous A8344G heteroplasmy and mitochondrial DNA copy number quantification in Myoclonus Epilepsy and Ragged-Red Fibers (MERRF) syndrome by a multiplex Molecular Beacon based real-time fluorescence PCR

The association of a particular mitochondrial DNA (mtDNA) mutation with different clinical phenotypes is a well-known feature of mitochondrial diseases. A simple genotype–phenotype correlation has not been found between mutation load and disease expression. Tissue and intercellular mosaicism as well as mtDNA copy number are thought to be responsible for the different clinical phenotypes. As disease expression of mitochondrial tRNA mutations is mostly in postmitotic tissues, studies to elucidate disease mechanisms need to be performed on patient material. Heteroplasmy quantitation and copy number estimation using small patient biopsy samples has not been reported before, mainly due to technical restrictions. In order to resolve this problem, we have developed a robust assay that utilizes Molecular Beacons to accurately quantify heteroplasmy levels and determine mtDNA copy number in small samples carrying the A8344G tRNALys mutation. It provides the methodological basis to investigate the role of heteroplasmy and mtDNA copy number in determining the clinical phenotypes.

Molecular Analysis of (R)-(+)-Mandelonitrile Lyase Microheterogeneity in Black Cherry1

The flavoprotein (R)-(+)-mandelonitrile lyase (MDL; EC 4.1.2.10), which plays a key role in cyanogenesis in rosaceous stone fruits, occurs in black cherry (Prunus serotina Ehrh.) homogenates as several closely related isoforms. Biochemical and molecular biological methods were used to investigate MDL microheterogeneity and function in this species. Three novel MDL cDNAs of high sequence identity (designated MDL2, MDL4, and MDL5) were isolated. Like MDL1 and MDL3 cDNAs (Z. Hu, J.E. Poulton [1997] Plant Physiol 115: 1359–1369), they had open reading frames that predicted a flavin adenine dinucleotide-binding site, multiple N-glycosylation sites, and an N-terminal signal sequence. The N terminus of an MDL isoform purified from seedlings matched the derived amino acid sequence of the MDL4 cDNA. Genomic sequences corresponding to the MDL1, MDL2, and MDL4 cDNAs were obtained by polymerase chain reaction amplification of genomic DNA. Like the previously reported mdl3 gene, these genes are interrupted at identical positions by three short, conserved introns. Given their overall similarity, we conclude that the genes mdl1, mdl2, mdl3, mdl4, and mdl5 are derived from a common ancestral gene and constitute members of a gene family. Genomic Southern-blot analysis showed that this family has approximately eight members. Northern-blot analysis using gene-specific probes revealed differential expression of the genes mdl1, mdl2, mdl3, mdl4, and mdl5.

Alternate protein frameworks for molecular recognition.

In an effort to determine whether proteins with structures other than the immunoglobulin fold can be used to mimic the ligand binding properties of antibodies, we generated a library from the four-helix bundle protein cytochrome b562 in which the two loops were randomized. Panning of this library against the bovine serum albumin (BSA) conjugate of N-methyl-p-nitrobenzylamine derivative 1 by phage display methods yielded cytochromes in which residues Trp-20, Arg-21, and Ser-22 in loop A and Arg-83 and Trp-84 in loop B were conserved. The individual mutants, which fold into native-like structure, bind selectively to the BSA-1 conjugate with micromolar dissociation constants (Kd), in comparison to a monoclonal antibody that binds selectively to 1 with a Kd of 290 nM. These and other antibody-like receptors may prove useful as therapeutic agents or as reagents for both intra- and extracellular studies.

Stabilization of APIs and nutraceuticals through cocrystallization and molecular confinement

Cocrystallization of the molecule of interest could be a smart and dainty way to tune solubility properties of solid phases leaving the molecule chemically unchanged, hence it is widely investigated by companies and by solid state scientists. Despite of this extremely high interest towards cocrystallization no particular emphasis has been paid to using it as a means to stabilize liquid molecules. In this work we define a benchmark of relevant molecules for human health that have been combined with suitable partners according to crystal engineering methods in order to obtain cocrystals. Solubility properties in different solvents of cocrystals new solid phases have been tested and compared to the properties of the drugs. A further approach to deal with volatile compounds is molecular confinement inside molecular scaffold. Nowadays metal organic frameworks (MOFs) are studied in many fields ranging from catalysis to trapping or storage of gases, such as hydrogen, methane, CO2 thanks to their extremely high porosity. Our goal is to confine liquid guests of biological relevance inside MOF pores, monitoring via X-ray diffraction, spectroscopy and thermal analysis the stabilization of the molecule of interest inside the cavities.

QSAR e dinâmica molecular no estudo de sistemas biomoleculares: predição da atividade biológica de antagonistas do receptor sigma-1 e simulações de bicamadas lipídicas

Diferentes abordagens teóricas têm sido utilizadas em estudos de sistemas biomoleculares com o objetivo de contribuir com o tratamento de diversas doenças. Para a dor neuropática, por exemplo, o estudo de compostos que interagem com o receptor sigma-1 (Sig-1R) pode elucidar os principais fatores associados à atividade biológica dos mesmos. Nesse propósito, estudos de Relações Quantitativas Estrutura-Atividade (QSAR) utilizando os métodos de regressão por Mínimos Quadrados Parciais (PLS) e Rede Neural Artificial (ANN) foram aplicados a 64 antagonistas do Sig-1R pertencentes à classe de 1-arilpirazóis. Modelos PLS e ANN foram utilizados com o objetivo de descrever comportamentos lineares e não lineares, respectivamente, entre um conjunto de descritores e a atividade biológica dos compostos selecionados. O modelo PLS foi obtido com 51 compostos no conjunto treinamento e 13 compostos no conjunto teste (r² = 0,768, q² = 0,684 e r²teste = 0,785). Testes de leave-N-out, randomização da atividade biológica e detecção de outliers confirmaram a robustez e estabilidade dos modelos e mostraram que os mesmos não foram obtidos por correlações ao acaso. Modelos também foram gerados a partir da Rede Neural Artificial Perceptron de Multicamadas (MLP-ANN), sendo que a arquitetura 6-12-1, treinada com as funções de transferência tansig-tansig, apresentou a melhor resposta para a predição da atividade biológica dos compostos (r²treinamento = 0,891, r²validação = 0,852 e r²teste = 0,793). Outra abordagem foi utilizada para simular o ambiente de membranas sinápticas utilizando bicamadas lipídicas compostas por POPC, DOPE, POPS e colesterol. Os estudos de dinâmica molecular desenvolvidos mostraram que altas concentrações de colesterol induzem redução da área por lipídeo e difusão lateral e aumento na espessura da membrana e nos valores de parâmetro de ordem causados pelo ordenamento das cadeias acil dos fosfolipídeos. As bicamadas lipídicas obtidas podem ser usadas para simular interações entre lipídeos e pequenas moléculas ou proteínas contribuindo para as pesquisas associadas a doenças como Alzheimer e Parkinson. As abordagens usadas nessa tese são essenciais para o desenvolvimento de novas pesquisas em Química Medicinal Computacional.

Molecular Genetic Basis of Opposite Gene by Environment Interactions in Reading Disability and Attention Deficit/Hyperactivity Disorder

The goal of this study is to better understand the genetic basis of Reading Disability (RD) and Attention Deficit Hyperactivity Disorder (ADHD) by examining molecular G x E interactions with parental education for each disorder. Research indicates that despite sharing genetic risk factors, RD and ADHD are influenced by different types of G x E interactions with parental education - a diathesis stress interaction in the case of ADHD and a bioecological interaction in RD. In order to resolve this apparent paradox, we conducted a preliminary study using behavioral genetic methods to test for G x E interactions in RD and the inattentive subtype of ADHD (ADHD-I) in the same sample of monozygotic and dizygotic Colorado Learning Disabilities Research Center same-sex twin pairs (DeFries et al., 1997), and our findings were consistent with the literature. We posited a genetic hypothesis for this opposite pattern of interactions, which suggests that only genes specific to each disorder enter into these opposite interactions, not the shared genes underlying their comorbidity. This study sought to further investigate this paradox using molecular genetics methods. We examined multiple candidate genes identified for RD or related language phenotypes and those identified for ADHD for G x E interactions with parental education. The specific aims of this study were as follows: 1) partition known risk alleles for RD and/or related language phenotypes and ADHD-I into those which are pleiotropic and non-pleiotropic by testing each risk allele for association with both RD and ADHD-I, 2) explore the main effects of parental education on both RD and ADHD-I, 3) address G-E correlations, and 4) conduct exploratory G x E interaction analyses in order to test the genetic hypothesis. Analyses suggested a number of pleiotropic genes that influence both RD and ADHD; however, results did not remain after correcting for multiple comparisons. Although exploratory G x E interaction findings were not significant after multiple comparison correction, results suggested a G x E interaction in the bioecological direction with KIAA0319, parental education, and ADHD-I. Given the limited power in the current study, replication of these findings with larger samples is necessary.

Comportamento tribo-mecânico e desgaste adesivo de materiais em nanoescala: análises por dinâmica molecular e mecânica do contínuo.

As formulações baseadas na mecânica do contínuo, embora precisas até certo ponto, por vezes não podem ser utilizadas, ou não são conceitualmente corretas para o entendimento de fenômenos em escalas reduzidas. Estas limitações podem aparecer no estudo dos fenômenos tribológicos em escala nanométrica, que passam a necessitar de novos métodos experimentais, teóricos e computacionais que permitam explorar estes fenômenos com a resolução necessária. Simulações atomísticas são capazes de descrever fenômenos em pequena escala, porém, o número necessário de átomos modelados e, portanto, o custo computacional - geralmente torna-se bastante elevado. Por outro lado, os métodos de simulação associados à mecânica do contínuo são mais interessantes em relação ao custo computacional, mas não são precisos na escala atômica. A combinação entre essas duas abordagens pode, então, permitir uma compreensão mais realista dos fenômenos da tribologia. Neste trabalho, discutem-se os conceitos básicos e modelos de atrito em escala atômica e apresentam-se estudos, por meio de simulação numérica, para a análise e compreensão dos mecanismos de atrito e desgaste no contato entre materiais. O problema é abordado em diferentes escalas, e propõe-se uma abordagem conjunta entre a Mecânica do Contínuo e a Dinâmica Molecular. Para tanto, foram executadas simulações numéricas, com complexidade crescente, do contato entre superfícies, partindo-se de um primeiro modelo que simula o efeito de defeitos cristalinos no fenômeno de escorregamento puro, considerando a Dinâmica Molecular. Posteriormente, inseriu-se, nos modelos da mecânica do contínuo, considerações sobre o fenômeno de adesão. A validação dos resultados é feita pela comparação entre as duas abordagens e com a literatura.

Transport across two interacting quantum dots: Bulk Kondo, Kondo box, and molecular regimes

We analyze the transport properties of a double quantum dot device with both dots coupled to perfect conducting leads and to a finite chain of N noninteracting sites connecting both of them. The interdot chain strongly influences the transport across the system and the local density of states of the dots. We study the case of a small number of sites, so that Kondo box effects are present, varying the coupling between the dots and the chain. For odd N and small coupling between the interdot chain and the dots, a state with two coexisting Kondo regimes develops: the bulk Kondo due to the quantum dots connected to leads and the one produced by the screening of the quantum dot spins by the spin in the finite chain at the Fermi level. As the coupling to the interdot chain increases, there is a crossover to a molecular Kondo effect, due to the screening of the molecule (formed by the finite chain and the quantum dots) spin by the leads. For even N the two Kondo temperatures regime does not develop and the physics is dominated by the usual competition between Kondo and antiferromagnetism between the quantum dots. We finally study how the transport properties are affected as N is increased. For the study we used exact multiconfigurational Lanczos calculations and finite-U slave-boson mean-field theory at T=0. The results obtained with both methods describe qualitatively and also quantitatively the same physics.

Electronic structure methods applied to gas-carbon reactions

A review is given on the fundamental studies of gas-carbon reactions using electronic structure methods in the last several decades. The three types of electronic structure methods including semi-empirical, ab initio and density functional theory, methods are briefly introduced first, followed by the studies on carbon reactions with hydrogen and oxygen-containing gases (non-catalysed and catalysed). The problems yet to solve and possible promising directions are discussed. (C) 2002 Elsevier Science Ltd. All rights reserved.

Adsorption of Supercritical fluids on graphitised thermal carbon black: Molecular layer structure theory versus grand canonical Monte Carlo simulation

This paper presents a detailed analysis of adsorption of supercritical fluids on nonporous graphitized thermal carbon black. Two methods are employed in the analysis. One is the molecular layer structure theory (MLST), proposed recently by our group, and the other is the grand canonical Monte Carlo (GCMC) simulation. They were applied to describe the adsorption of argon, krypton, methane, ethylene, and sulfur hexafluoride on graphitized thermal carbon black. It was found that the MLST describes all the experimental data at various temperatures well. Results from GCMC simulations describe well the data at low pressure but show some deviations at higher pressures for all the adsorbates tested. The question of negative surface excess is also discussed in this paper.

Molecular characterization of the microbial species that colonize human ileal and colonic mucosa by using 16S rDNA sequence analysis

Aim: The aim of this study was to characterize the bacterial community adhering to the mucosa of the terminal ileum, and proximal and distal colon of the human digestive tract. Methods and Results: Pinch samples of the terminal ileum, proximal and distal colon were taken from a healthy 35-year-old, and a 68-year-old subject with mild diverticulosis. The 16S rDNA genes were amplified using a low number of PCR cycles, cloned, and sequenced. In total, 361 sequences were obtained comprising 70 operational taxonomic units (OTU), with a calculated coverage of 82.6%. Twenty-three per cent of OTU were common to the terminal ileum, proximal colon and distal colon, but 14% OTU were only found in the terminal ileum, and 43% were only associated with the proximal or distal colon. The most frequently represented clones were from the Clostridium group XIVa (24.7%), and the Bacteroidetes (Cytophaga-Flavobacteria-Bacteroides ) cluster (27.7%). Conclusion: Comparison of 16S rDNA clone libraries of the hindgut across mammalian species confirms that the distribution of phylogenetic groups is similar irrespective of the host species. Lesser site-related differences within groups or clusters of organisms, are probable. Significance and Impact: This study provides further evidence of the distribution of the bacteria on the mucosal surfaces of the human hindgut. Data contribute to the benchmarking of the microbial composition of the human digestive tract.

Genetic assignment methods for the direct, real-time estimation of migration rate: a simulation-based exploration of accuracy and power

Genetic assignment methods use genotype likelihoods to draw inference about where individuals were or were not born, potentially allowing direct, real-time estimates of dispersal. We used simulated data sets to test the power and accuracy of Monte Carlo resampling methods in generating statistical thresholds for identifying F-0 immigrants in populations with ongoing gene flow, and hence for providing direct, real-time estimates of migration rates. The identification of accurate critical values required that resampling methods preserved the linkage disequilibrium deriving from recent generations of immigrants and reflected the sampling variance present in the data set being analysed. A novel Monte Carlo resampling method taking into account these aspects was proposed and its efficiency was evaluated. Power and error were relatively insensitive to the frequency assumed for missing alleles. Power to identify F-0 immigrants was improved by using large sample size (up to about 50 individuals) and by sampling all populations from which migrants may have originated. A combination of plotting genotype likelihoods and calculating mean genotype likelihood ratios (D-LR) appeared to be an effective way to predict whether F-0 immigrants could be identified for a particular pair of populations using a given set of markers.