An evaluation of molecular models of the cytochrome P450 Streptomyces griseolus enzymes P450SU1 and P450SU2

Summary P450SU1 and P450SU2 are herbicide-inducible bacterial cytochrome P450 enzymes from Streptomyces griseolus. They have two of the highest sequence identities to camphor hydroxylase (P450cam from Pseudomonas putida), the cytochrome P450 with the first known crystal structure. We have built several models of these two proteins to investigate the variability in the structures that can occur from using different modeling protocols. We looked at variability due to alignment methods, backbone loop conformations and refinement methods. We have constructed two models for each protein using two alignment algorithms, and then an additional model using an identical alignment but different loop conformations for both buried and surface loops. The alignments used to build the models were created using the Needleman-Wunsch method, adapted for multiple sequences, and a manual method that utilized both a dotmatrix search matrix and the Needleman-Wunsch method. After constructing the initial models, several energy minimization methods were used to explore the variability in the final models caused by the choice of minimization techniques. Features of cytochrome P450cam and the cytochrome P450 superfamily, such as the ferredoxin binding site, the heme binding site and the substrate binding site were used to evaluate the validity of the models. Although the final structures were very similar between the models with different alignments, active-site residues were found to be dependent on the conformations of buried loops and early stages of energy minimization. We show which regions of the active site are the most dependent on the particular methods used, and which parts of the structures seem to be independent of the methods.     

Distributed data parallel coupled-cluster algorithm: Application to the 2-hydroxypyridine/2-pyridone tautomerism

The recently developed parallel coupled-cluster algorithm of Rendell, Lee, and Lindh [Chem. Phys. Lett., 194 , 84 (1992)] is extended to allow four-indexed quantities containing one or two indices in the virtual orbital space to be stored across the global memory of distributed-memory parallel processors. Quantities such as the double-excitation amplitudes can now be distributed over multiple nodes, with blocks of data retrieved from remote nodes by the use of interrupt handlers. As an application of the new code, we have investigated the potential energy surface of the 2-hydroxypyridine/2-pyridone tautomers. Using large basis sets, the structure of each tautomer and the transition state connecting the two minima has been determined at the SCF level. The relative energy difference and the activation energy were then redetermined using the MP2, CCSD, and CCSD(T) methods. All calculations have been performed on Intel distributed-memory supercomputers. The largest coupled-cluster calculations contained over 2 million double-excitation amplitudes. © John Wiley & Sons, Inc.     

Flexibility in the partial reduction of 2,5-disubstituted pyrroles: application to the synthesis of DMDP

[reaction: see text] The partial reduction of electron-deficient 2,5-disubstituted pyrroles has been developed into a flexible procedure that gives control of relative stereochemistry by variation of the reduction conditions. After the reaction, the pyrroline products were dihydroxylated at C-3,4 to give either the cis or trans isomers; this flexibility means that a variety of polyhydroxylated pyrrolidines can be prepared in a short sequence. Finally, this method was applied to a synthesis of the naturally occurring glycosidase inhibitor DMDP.     

An efficient and highly flexible synthesis of (α),β-unsaturated γ-oxoesters

Siloxy substituted methyl cyclopropanecar?ylates are cleaved by bromine to afford α-bromo-γ-oxoester and, after treatment with triethylamine, good yields of α,β-unsaturated γ-oxoesters.     

Hydration free energies and entropies for water in protein interiors

Free energy calculations for the transfer of a water molecule from the pure liquid to an interior cavity site in a protein are presented. Two different protein cavities, in bovine pancreatic trypsin inhibitor (BPTI) and in the I76A mutant of barnase, represent very different environments for the water molecule: one which is polar, forming four water-protein hydrogen bonds, and one which is more hydrophobic, forming only one water-protein hydrogen bond. The calculations give very different free energies for the different cavities, with only the polar BPTI cavity predicted to be hydrated. The corresponding entropies for the transfer to the interior cavities are calculated as well and show that the transfer to the polar cavity is significantly entropically unfavorable while the transfer to the nonpolar cavity is entropically favorable. For both proteins an analysis of the fluctuations in the positions of the protein atoms shows that the addition of a water molecule makes the protein more flexible. This increased flexibility appears to be due to an increased length and weakened strength of protein-protein hydrogen bonds near the cavity.     

LC/MS/MS analysis of the endogenous dimethyltryptamine hallucinogens,their precursors,and major metabolites in rat pineal gland microdialysate

We report a qualitative liquid chromatography–tandem mass spectrometry (LC/MS/MS) method for the simultaneous analysis of the three known N,N‐dimethyltryptamine endogenous hallucinogens, their precursors and metabolites, as well as melatonin and its metabolic precursors. The method was characterized using artificial cerebrospinal fluid (aCSF) as the matrix and was subsequently applied to the analysis of rat brain pineal gland‐aCSF microdialysate. The method describes the simultaneous analysis of 23 chemically diverse compounds plus a deuterated internal standard by direct injection, requiring no dilution or extraction of the samples. The results demonstrate that this is a simple, sensitive, specific and direct approach to the qualitative analysis of these compounds in this matrix. The protocol also employs stringent MS confirmatory criteria for the detection and confirmation of the compounds examined, including exact mass measurements. The excellent limits of detection and broad scope make it a valuable research tool for examining the endogenous hallucinogen pathways in the central nervous system. We report here, for the first time, the presence of N,N‐dimethyltryptamine in pineal gland microdialysate obtained from the rat. Copyright © 2013 John Wiley & Sons, Ltd.     

Inferring epidemiological parameters from phylogenetic information for the HIV-1 epidemic among MSM

The HIV-1 epidemic in Europe is primarily sustained by a dynamic topology of sexual interactions among MSM who have individual immune systems and behavior. This epidemiological process shapes the phylogeny of the virus population. Both fields of epidemic modeling and phylogenetics have a long history, however it remains difficult to use phylogenetic data to infer epidemiological parameters such as the structure of the sexual network and the per-act infectiousness. This is because phylogenetic data is necessarily incomplete and ambiguous. Here we show that the cluster-size distribution indeed contains information about epidemiological parameters using detailed numberical experiments. We simulate the HIV epidemic among MSM many times using the Monte Carlo method with all parameter values and their ranges taken from literature. For each simulation and the corresponding set of parameter values we calculate the likelihood of reproducing an observed cluster-size distribution. The result is an estimated likelihood distribution of all parameters from the phylogenetic data, in particular the structure of the sexual network, the per-act infectiousness, and the risk behavior reduction upon diagnosis. These likelihood distributions encode the knowledge provided by the observed cluster-size distrbution, which we quantify using information theory. Our work suggests that the growing body of genetic data of patients can be exploited to understand the underlying epidemiological process.     

A velocity program using the Kanade–Lucas–Tomasi feature-tracking algorithm with demonstration for pressure and electroosmosis conditions

Investigating microfluidic flow profiles is of interest in the microfluidics field for the determination of various characteristics of a lab-on-a-chip system. Microparticle tracking velocimetry uses computational methods upon recording video footage of microfluidic flow to ultimately visualize motion within a microfluidic system across all frames of a video. Current methods are computationally expensive or require extensive instrumentation. A computational method suited to microparticle tracking applications is the robust Kanade–Lucas–Tomasi (KLT) feature-tracking algorithm. This work explores a microparticle tracking velocimetry program using the KLT feature-tracking algorithm. The developed program is demonstrated using pressure-driven and EOF and compared with the respective mathematical fluid flow models. An electrostatics analysis of EOF conditions is performed in the development of the mathematical using a Poisson's Equation solver. This analysis is used to quantify the zeta potential of the electroosmotic system. Overall, the KLT feature-tracking algorithm presented in this work proved to be highly reliable and computationally efficient for investigations of pressure-driven and EOF in a microfluidic system.     

Molecular organization in two-dimensional films of liquid II. Langmuir and Langmuir-Blodgett films of a perylene-like dye mixed with liquid crystals having a terminal cyano group crystalline mixtures

Langmuir and Langmuir-Blodgett (LB) films of a perylene-like compound and its binary mixtures with 4-octyl-4'-cyanobiphenyl (8CB) and 4-pentyl-4"-cyano-p-terphenyl (5CT) have been studied. On the basis of the surface pressure-area isotherms, the molecular organization on the air-water interface has been estimated. Information about the miscibility or the phase separation of components in the binary mixtures has been obtained. The spectroscopic study of the LB films has allowed conclusions to be drawn about the arrangement of the molecules on the quartz slides. The fluorescence spectra of the LB films of the perylene-like compound have revealed the formation of self-aggregates.