Cryptic post-transition state bifurcations that reduce the efficiency of lactone-forming Rh-carbenoid C-H insertions.

Byproducts of chemical reactions are generally thought to result from the competition between two reaction pathways, each with its own rate-determining transition state structure. We show here, however, that pathways with a single transition state structure followed by a post-transition state bifurcation may also be a source of undesired products, especially those whose appearance is unexpected. The viability of this scenario for intramolecular C-H insertion reactions affording β-lactones via Rh-carbenoid intermediates is assessed through quantum chemical calculations on potential energy surfaces and quasi-classical molecular dynamics simulations. It appears that, in these cases, the rhodium catalyst is to blame for the accessibility of a second, unintended, pathway following the transition state structure for β-lactone formation that leads to fragmentation to a ketene and carbonyl compound. If an unexpected product is formed via a post-transition state bifurcation, conventional strategies for suppressing its formation are unlikely to succeed. Guidelines for recognizing the presence of a post-transition state bifurcation are described here, along with hints at means for controlling product distributions

Interactive molecular dynamics in virtual reality for accurate flexible protein-ligand docking

Simulating drug binding and unbinding is a challenge, as the rugged energy landscapes that separate bound and unbound states require extensive sampling that consumes significant computational resources. Here, we describe the use of interactive molecular dynamics in virtual reality (iMD-VR) as an accurate low-cost strategy for flexible protein-ligand docking. We outline an experimental protocol which enables expert iMD-VR users to guide ligands into and out of the binding pockets of trypsin, neuraminidase, and HIV-1 protease, and recreate their respective crystallographic protein-ligand binding poses within 5 - 10 minutes. Following a brief training phase, our studies shown that iMD-VR novices were able to generate unbinding and rebinding pathways on similar timescales as iMD-VR experts, with the majority able to recover binding poses within 2.15 Angstrom RMSD of the crystallographic binding pose. These results indicate that iMD-VR affords sufficient control for users to carry out the detailed atomic manipulations required to dock flexible ligands into dynamic enzyme active sites and recover crystallographic poses, offering an interesting new approach for simulating drug docking and generating binding hypotheses.Comment: PLOS ON

Sur un théorème de Poncelet et sa généralisation par M. Horvarth

Most chemical transformations (reactions or conformational changes) that are of interest to researchers have many degrees of freedom, usually too many to visualize without reducing the dimensionality of the system to include only the most important atomic motions. In this article, we describe a method of using Principal Component Analysis (PCA) for analyzing a series of molecular geometries (e.g., a reaction pathway or molecular dynamics trajectory) and determining the reduced dimensional space that captures the most structural variance in the fewest dimensions. The software written to carry out this method is called PathReducer, which permits (1) visualizing the geometries in a reduced dimensional space, (2) determining the axes that make up the reduced dimensional space, and (3) projecting the series of geometries into the low-dimensional space for visualization. We investigated two options to represent molecular structures within PathReducer: aligned Cartesian coordinates and matrices of interatomic distances. We found that interatomic distance matrices better captured non-linear motions in a smaller number of dimensions. To demonstrate the utility of PathReducer, we have carried out a number of applications where we have projected molecular dynamics trajectories into a reduced dimensional space defined by an intrinsic reaction coordinate. The visualizations provided by this analysis show that dynamic paths can differ greatly from the minimum energy pathway on a potential energy surface. Viewing intrinsic reaction coordinates and trajectories in this way provides a quick way to gather qualitative information about the pathways trajectories take relative to a minimum energy path. Given that the outputs from PCA are linear combinations of the input molecular structure coordinates (i.e., Cartesian coordinates or interatomic distances), they can be easily transferred to other types of calculations that require the definition of a reduced dimensional space (e.g., biased molecular dynamics simulations)

Measurement of the Lifetime Difference Between B_s Mass Eigenstates

We present measurements of the lifetimes and polarization amplitudes for B_s --> J/psi phi and B_d --> J/psi K*0 decays. Lifetimes of the heavy (H) and light (L) mass eigenstates in the B_s system are separately measured for the first time by determining the relative contributions of amplitudes with definite CP as a function of the decay time. Using 203 +/- 15 B_s decays, we obtain tau_L = (1.05 +{0.16}/-{0.13} +/- 0.02) ps and tau_H = (2.07 +{0.58}/-{0.46} +/- 0.03) ps. Expressed in terms of the difference DeltaGamma_s and average Gamma_s, of the decay rates of the two eigenstates, the results are DeltaGamma_s/Gamma_s = (65 +{25}/-{33} +/- 1)%, and DeltaGamma_s = (0.47 +{0.19}/-{0.24} +/- 0.01) inverse ps.Comment: 8 pages, 3 figures, 2 tables; as published in Physical Review Letters on 16 March 2005; revisions are for length and typesetting only, no changes in results or conclusion

HIVToolbox, an Integrated Web Application for Investigating HIV

Many bioinformatic databases and applications focus on a limited domain of knowledge federating links to information in other databases. This segregated data structure likely limits our ability to investigate and understand complex biological systems. To facilitate research, therefore, we have built HIVToolbox, which integrates much of the knowledge about HIV proteins and allows virologists and structural biologists to access sequence, structure, and functional relationships in an intuitive web application. HIV-1 integrase protein was used as a case study to show the utility of this application. We show how data integration facilitates identification of new questions and hypotheses much more rapid and convenient than current approaches using isolated repositories. Several new hypotheses for integrase were created as an example, and we experimentally confirmed a predicted CK2 phosphorylation site. Weblink: [http://hivtoolbox.bio-toolkit.com

Tracking Cats: Problems with Placing Feline Carnivores on δ18O, δD Isoscapes

Several felids are endangered and threatened by the illegal wildlife trade. Establishing geographic origin of tissues of endangered species is thus crucial for wildlife crime investigations and effective conservation strategies. As shown in other species, stable isotope analysis of hydrogen and oxygen in hair (δD(h), δ(18)O(h)) can be used as a tool for provenance determination. However, reliably predicting the spatial distribution of δD(h) and δ(18)O(h) requires confirmation from animal tissues of known origin and a detailed understanding of the isotopic routing of dietary nutrients into felid hair.We used coupled δD(h) and δ(18)O(h) measurements from the North American bobcat (Lynx rufus) and puma (Puma concolor) with precipitation-based assignment isoscapes to test the feasibility of isotopic geo-location of felidae. Hairs of felid and rabbit museum specimens from 75 sites across the United States and Canada were analyzed. Bobcat and puma lacked a significant correlation between H/O isotopes in hair and local waters, and also exhibited an isotopic decoupling of δ(18)O(h) and δD(h). Conversely, strong δD and δ(18)O coupling was found for key prey, eastern cottontail rabbit (Sylvilagus floridanus; hair) and white-tailed deer (Odocoileus virginianus; collagen, bone phosphate).Puma and bobcat hairs do not adhere to expected pattern of H and O isotopic variation predicted by precipitation isoscapes for North America. Thus, using bulk hair, felids cannot be placed on δ(18)O and δD isoscapes for use in forensic investigations. The effective application of isotopes to trace the provenance of feline carnivores is likely compromised by major controls of their diet, physiology and metabolism on hair δ(18)O and δD related to body water budgets. Controlled feeding experiments, combined with single amino acid isotope analysis of diets and hair, are needed to reveal mechanisms and physiological traits explaining why felid hair does not follow isotopic patterns demonstrated in many other taxa