Precise and global identification of phospholipid molecular species by an Orbitrap mass spectrometer and automated search engine Lipid Search

In the present research, we have established a new lipidomics approach for the comprehensive and precise identification of molecular species in a crude lipid mixture using a LTQ Orbitrap mass spectrometer (MS) and reverse-phase liquid chromatography (RPLC) combination with our newly developed lipid search engine “Lipid Search”. LTQ Orbitrap provides high mass accuracy MS spectra by Fourier-transform (FT) mass spectrometer mode and can perform rapid MSⁿ by ion trap (IT) mass spectrometer mode. In this study, the negative ion mode was selected to detect fragment ions from phospholipids, such as fatty acid anions, by MS2 or MS3. We selected the specific detection approach by neutral loss survey-dependent MS3, for the identification of molecular species of phosphatidylcholine, sphingomyelin and phosphatidylserine. Identification of molecular species was performed by using both the high mass accuracy of the mass spectrometric data obtained from FT mode and structural data obtained from fragments in IT mode. Some alkylacyl and alkenylacyl species have the same m/z value as molecular-related ions and fragment ions, thus, direct acid hydrolysis analysis was performed to identify alkylacyl and alkenylacyl species, and then the RPLC–LTQ Orbitrap method was applied. As a result, 290 species from mouse liver and 248 species from mouse brain were identified within six different classes of phospholipid, only those in manually detected and confirmed. Most of all manually detected mass peaks were also automatically detected by “Lipid Search”. Adding to differences in molecular species in different classes of phospholipids, many characteristic differences in molecular species were detected in mouse liver and brain. More variable number of saturated and monounsaturated fatty acid-containing molecular species were detected in mouse brain than liver.     

Conformational analysis of methylphenidate: comparison of molecular orbital and molecular mechanics methods

Summary Methylphenidate (MP) binds to the cocaine binding site on the dopamine transporter and inhibits reuptake of dopamine, but does not appear to have the same abuse potential as cocaine. This study, part of a comprehensive effort to identify a drug treatment for cocaine abuse, investigates the effect of choice of calculation technique and of solvent model on the conformational potential energy surface (PES) of MP and a rigid methylphenidate (RMP) analogue which exhibits the same dopamine transporter binding affinity as MP. Conformational analysis was carried out by the AM1 and AM1/SM5.4 semiempirical molecular orbital methods, a molecular mechanics method (Tripos force field with the dielectric set equal to that of vacuum or water) and the HF/6-31G* molecular orbital method in vacuum phase. Although all three methods differ somewhat in the local details of the PES, the general trends are the same for neutral and protonated MP. In vacuum phase, protonation has a distinctive effect in decreasing the regions of space available to the local conformational minima. Solvent has little effect on the PES of the neutral molecule and tends to stabilize the protonated species. The random search (RS) conformational analysis technique using the Tripos force field was found to be capable of locating the minima found by the molecular orbital methods using systematic grid search. This suggests that the RS/Tripos force field/vacuum phase protocol is a reasonable choice for locating the local minima of MP. However, the Tripos force field gave significantly larger phenyl ring rotational barriers than the molecular orbital methods for MP and RMP. For both the neutral and protonated cases, all three methods found the phenyl ring rotational barriers for the RMP conformers/invertamers (denoted as cte, tte, and cta) to be: cte, tte> MP > cta. Solvation has negligible effect on the phenyl ring rotational barrier of RMP. The B3LYP/6-31G* density functional method was used to calculate the phenyl ring rotational barrier for neutral MP and gave results very similar to those of the HF/6-31G* method.     

Batch tautomer generation with MolTPC

Besides all their conformational degrees of freedom, drug‐like molecules and natural products often also undergo tautomeric interconversions. Compared to the huge efforts made in experimental investigation of tautomerism, open and free algorithmic solutions for prototropic tautomer generation are surprisingly rare. The few freely available software packages limit their output to a subset of the possible configurational space by sometimes unwanted prior assumptions and complete neglection of ring‐chain tautomerism. Here, we describe an adjustable fully automatic tautomer enumeration approach, which is freely available and also incorporates the detection of ring‐chain variants. The algorithm is implemented in the MolTPC framework and accessible on SourceForge. © 2013 Wiley Periodicals, Inc.     

Testing an alternative search algorithm for compound identification with the ‘Wiley Registry of Tandem Mass Spectral Data,MSforID’

A tandem mass spectral database system consists of a library of reference spectra and a search program. State‐of‐the‐art search programs show a high tolerance for variability in compound‐specific fragmentation patterns produced by collision‐induced decomposition and enable sensitive and specific ‘identity search’. In this communication, performance characteristics of two search algorithms combined with the ‘Wiley Registry of Tandem Mass Spectral Data, MSforID’ (Wiley Registry MSMS, John Wiley and Sons, Hoboken, NJ, USA) were evaluated. The search algorithms tested were the MSMS search algorithm implemented in the NIST MS Search program 2.0g (NIST, Gaithersburg, MD, USA) and the MSforID algorithm (John Wiley and Sons, Hoboken, NJ, USA). Sample spectra were acquired on different instruments and, thus, covered a broad range of possible experimental conditions or were generated in silico. For each algorithm, more than 30 000 matches were performed. Statistical evaluation of the library search results revealed that principally both search algorithms can be combined with the Wiley Registry MSMS to create a reliable identification tool. It appears, however, that a higher degree of spectral similarity is necessary to obtain a correct match with the NIST MS Search program. This characteristic of the NIST MS Search program has a positive effect on specificity as it helps to avoid false positive matches (type I errors), but reduces sensitivity. Thus, particularly with sample spectra acquired on instruments differing in their setup from tandem‐in‐space type fragmentation, a comparably higher number of false negative matches (type II errors) were observed by searching the Wiley Registry MSMS. Copyright © 2013 John Wiley & Sons, Ltd.     

Efficiency of random search procedures along the silicon cluster series: Si(n) (n=5-10, 15, and 20)

The efficiency of the simplest isomeric search procedure consisting in random generation of sets of atomic coordinates followed by density functional theory geometry optimization is tested on the silicon cluster series (Si(5-10, 15, 20)). Criteria such as yield, isomer distributions and recurrences are used to clearly establish the performance of the approach with respect to increasing cluster size. The elimination of unphysical candidate structures and the use of distinct box shapes and theoretical levels are also investigated. For the smaller Si(n) (n=5-10) clusters, the generation of random coordinates within a spherical box is found to offer a reasonable alternative to more complex algorithms by allowing straightforward identification of every known low-lying local minima. The simple stochastic search of larger clusters (i.e. Si(15) and Si(20)) is however complicated by the exponentially increasing number of both low- and high-lying minima leading to rather arbitrary and non-comprehensive results.     

Exploring sequential quantum adiabatic switching across supersymmetric partners for finding the eigenstates of a system

We demonstrate that one can exhaustively determine the n‐bound eigenstates of a Hamiltonian H by constructing a sequence of supersymmetric (SUSY) partner Hamiltonians and invoking a time‐dependent quantum adiabatic switching algorithm for passage from the ground state of one to the other. The ground states of the initial pair H⁽⁰⁾ and H⁽¹⁾ are constructed by solving the Riccati equation for the superpotential ?⁽⁰⁾ for H⁽⁰⁾ and adiabatically switching from the ground state Ψ of H⁽⁰⁾ to the ground state Ψ of H⁽¹⁾. The charge operator Q is then used to recover the first excited state Ψ of H⁽⁰⁾. The procedure is repeated for the ground states of SUSY pairs H^{(n + 1)} and H^{(n + 2)}, and appropriate charge operators lead to the excited states Ψ of H⁽⁰⁾ with , thereby exhausting the full eigenspectrum of H⁽⁰⁾. The workability of the proposed method is shown with several well‐known examples. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Molecular conformation search by distance matrix perturbations

Three-dimensional molecular structure is fundamental in chemical function identification and computer-aided drug design. The enumeration of a small number of feasible conformations provides a rigorous way to determine the optimal or a few acceptable conformations. Our contribution concerns a heuristic enhancement of a method based on distance geometry, typically in relation with experiments of the NMR type. Distance geometry has been approached by different viewpoints; ours is expected to help in several subtasks arising in the process that determines 3D structure from distance information. More precisely, the input to our algorithm consists of a set of approximate distances of varying precision; some are specified by the covalent structure and others by Nuclear Magnetic Resonance (NMR) experiments (or X-ray crystallography which, however, requires crystallization). The output is a valid tertiary structure in a specified neighborhood of the input. Our approach should help in detecting outliers of the NMR experiments, and handles inputs with partial information. Moreover, our technique is able to bound the number of degrees of freedom of the conformation manifold. We have used numerical linear algebra algorithms for reasons of speed, and because they are well-implemented, fully documented and widely available. Our main tools include, besides distance matrices, structure-preserving matrix perturbations for minimizing singular values. Our MATLAB (or SCILAB) implementation is described and illustrated.AMS subject Classification: 92E10 Molecular structure, 92C40 Biochemistry, molecular biology, 65F15 Eigenvalues, eigenvectors, 15A18 Eigenvalues, singular values, and eigenvectors     

Conformational analysis of [Met5]-enkephalin: Solvation and ionization considerations

[Met5]-Enkephalin has the sequence Tyr-Gly-Gly-Phe-Met. Only the extended conformation of the peptide has been observed by X-ray crystallography. Nuclear magnetic resonance spectroscopy supports the presence of a turn at Gly 3 and Phe 4 in dimethyl sulfoxide. In this study, the peptide conformational states and thermodynamic properties are understood in terms of ionization state and solvent environment. In the calculation, final conformations obtained from multiple independent Monte Carlo simulated annealing conformational searches are starting points for molecular dynamics simulations. In an aqueous environment given by the use of solvation free energy and the zwitterionic state, dominant structural motifs computed are G-P Type II bend, G-G Type II bend, and G-G Type I bend motifs, in order of increasing free energy. In the calculation of the peptide with neutral N- and C-termini and solvation free energy, the extended conformer dominates (by at least a factor of 2.5), and the conformation of another low free energy conformer superimposes well on the pharmacophoric groups of morphine. Neutralization of charge and solvation induce and stabilize the extended conformation, respectively. A mechanism of inter-conversion between the extended conformer and three bent conformers is supported by /-scatter plots, and by the conformer relative free energies. An estimate of the entropy change of receptor unbinding is 8.3 cal K-1 mol-1, which gives rise to a -2.5 kcal/mol entropy contribution to the free energy of unbinding at 25 °C. The conformational analysis methodology described here should be useful in studies on short peptides and flexible protein surface loops that have important biological implications.     

The Relevance of the Collaborative Effect in Determining the Performances of Photorefractive Polymer Materials

A derivative of 2‐methylindole, 3‐[2‐(4‐nitrophenyl)ethenyl]‐1‐allyl‐2‐methylindole, NPEMI‐A, is studied for its photoconductivity and photorefractivity behaviour. Its blends with the organic polymer poly‐(2,3‐dimethyl‐N‐vinylindole), PVDMI, are also investigated. Due to the expected and devised mutual solubility of the two components of the blends, it is possible to carry out measurements with the weight percent of the chromophore NPEMI‐A changing from zero to 100. Films were produced by a squeezing process between two ITO‐covered glass sheets. No opacity phenomena, that are so common for many other organic blends due to the segregation of the dissolved chromophore, are observed. The photorefractive optical gain Γ₂ is obtained as a function of the chromophore content. Differential scanning calorimetry measurements (DSC) are also carried out to obtain the whole change of the glass transition temperature T_g as a function of the amount of chromophore contained in the blends. From the experimental trend of T_g a meaningful quantitative estimate of the value of the electrostatic interactions acting in the studied blends, is obtained. The importance of the value of T_g, and of the electrostatic interactions, in determining the extent of the photorefractivity is clearly evident. The results are compared for NPEMI‐A (Γ₂=210 cm^?1) and for NPEMI‐E (Γ₂ ≈ 2000 cm^?1) that has a N‐2‐ethylhexyl group instead of a N‐allyl group. The Pockels and Kerr contributions and—for the first time—a “collaborative effect” of the photorefractivity of NPEMI‐A are distinguished and quantitatively evaluated.     

Identification of methylated regions with peak search based on Poisson model from massively parallel methylated DNA immunoprecipitation-sequencing data

DNA methylation is one of the most important epigenetic modification types, which plays a critical role in gene expression. High efficient surveying of whole genome DNA methylation has been aims of many researchers for long. Recently, the rapidly developed massively parallel DNA‐sequencing technologies open the floodgates to vast volumes of sequence data, enabling a paradigm shift in profiling the whole genome methylation. Here, we describe a strategy, combining methylated DNA immunoprecipitation sequencing with peak search to identify methylated regions on a whole‐genome scale. Massively parallel methylated DNA immunoprecipitation sequencing combined with methylation DNA immunoprecipitation was adopted to obtain methylated DNA sequence data from human leukemia cell line K562, and the methylated regions were identified by peak search based on Poison model. From our result, 140 958 non‐overlapping methylated regions have been identified in the whole genome. Also, the credibility of result has been proved by its strong correlation with bisulfite‐sequencing data (Pearson R²=0.92). It suggests that this method provides a reliable and high‐throughput strategy for whole genome methylation identification.