IDSite: An accurate approach to predict P450-mediated drug metabolism

Accurate prediction of drug metabolism is crucial for drug design. Since a large majority of drugs metabolism involves P450 enzymes, we herein describe a computational approach, IDSite, to predict P450-mediated drug metabolism. To model induced-fit effects, IDSite samples the conformational space with flexible docking in Glide followed by two refinement stages using the Protein Local Optimization Program (PLOP). Sites of metabolism (SOMs) are predicted according to a physical-based score that evaluates the potential of atoms to react with the catalytic iron center. As a preliminary test, we present in this paper the prediction of hydroxylation and O-dealkylation sites mediated by CYP2D6 using two different models: a physical-based simulation model, and a modification of this model in which a small number of parameters are fit to a training set. Without fitting any parameters to experimental data, the Physical IDSite scoring recovers 83% of the experimental observations for 56 compounds with a very low false positive rate. With only 4 fitted parameters, the Fitted IDSite was trained with the subset of 36 compounds and successfully applied to the other 20 compounds, recovering 94% of the experimental observations with high sensitivity and specificity for both sets.     

Single-molecule conductance through multiple π-π-stacked benzene rings determined with direct electrode-to-benzene ring connections

Understanding electron transport across π-π-stacked systems will help to answer fundamental questions about biochemical redox processes and benefit the design of new materials and molecular devices. Herein we employed the STM break-junction technique to measure the single-molecule conductance of multiple π-π-stacked aromatic rings. We studied electron transport through up to four stacked benzene rings held together in an eclipsed fashion via a paracyclophane scaffold. We found that the strained hydrocarbons studied herein couple directly to gold electrodes during the measurements; hence, we did not require any heteroatom binding groups as electrical contacts. Density functional theory-based calculations suggest that the gold atoms of the electrodes bind to two neighboring carbon atoms of the outermost cyclophane benzene rings in η(2) fashion. Our measurements show an exponential decay of the conductance with an increasing number of stacked benzene rings, indicating a nonresonant tunneling mechanism. Furthermore, STM tip-substrate displacement data provide additional evidence that the electrodes bind to the outermost benzene rings of the π-π-stacked molecular wires.     

Insights into the different dioxygen activation pathways of methane and toluene monooxygenase hydroxylases

The methane and toluene monooxygenase hydroxylases (MMOH and TMOH, respectively) have almost identical active sites, yet the physical and chemical properties of their oxygenated intermediates, designated P*, H(peroxo), Q, and Q* in MMOH and ToMOH(peroxo) in a subclass of TMOH, ToMOH, are substantially different. We review and compare the structural differences in the vicinity of the active sites of these enzymes and discuss which changes could give rise to the different behavior of H(peroxo) and Q. In particular, analysis of multiple crystal structures reveals that T213 in MMOH and the analogous T201 in TMOH, located in the immediate vicinity of the active site, have different rotatory configurations. We study the rotational energy profiles of these threonine residues with the use of molecular mechanics (MM) and quantum mechanics/molecular mechanics (QM/MM) computational methods and put forward a hypothesis according to which T213 and T201 play an important role in the formation of different types of peroxodiiron(III) species in MMOH and ToMOH. The hypothesis is indirectly supported by the QM/MM calculations of the peroxodiiron(III) models of ToMOH and the theoretically computed Mo?ssbauer spectra. It also helps explain the formation of two distinct peroxodiiron(III) species in the T201S mutant of ToMOH. Additionally, a role for the ToMOD regulatory protein, which is essential for intermediate formation and protein functioning in the ToMO system, is advanced. We find that the low quadrupole splitting parameter in the Mo?ssbauer spectrum observed for a ToMOH(peroxo) intermediate can be explained by protonation of the peroxo moiety, possibly stabilized by the T201 residue. Finally, similarities between the oxygen activation mechanisms of the monooxygenases and cytochrome P450 are discussed.     

DYNGA: a general purpose QM-MM-MD program. I. Application to water

We present a general purpose QM-MM-MD engine (DYNGA) designed to test alternative hybrid Hamiltonians geared towards the treatment of problems of interest in structural biology including the use of experimental data constraints. In this first presentation we use DYNGA to explore the behaviour of a traditional QM-MM approach in the treatment of the water—water interaction. We find the potential energy hypersurface for the water dimer computed with the HF 4–31G*/TIP3P hybrid Hamiltonian tends to be too flat. We also explore the effect of using traditional QM-MM techniques on proton wires and conclude there is a need for improvement, possibly addressed by using polarizable force fields.     

PAST AND FUTURE MANAGEMENT OF A COLLAPSED FISHERY: THE BAY OF BISCAY ANCHOVY

Abstract We analyze the efficiency of the international management of the Bay of Biscay anchovy. While a sharing agreement between France and Spain has been in place since 1992, the fish stock collapsed in 2005 and the fishery closed from 2005 to spring 2010. We consider differences in production technologies between both countries and calibrate our model using data from 1987 to 2009. Our results suggest two sources of rent dissipation under the existing sharing agreement: inefficient quota allocation and production inefficiencies due to inflexible national regulations. We discuss several alternatives to improve management.     

ConFlo III - an interface for high precision delta(13)C and delta(15)N analysis with an extended dynamic range

A newly developed interface coupling a CHN combustion device (elemental analyser 'EA') to an isotope ratio mass spectrometer is described and evaluated. The purpose of the device is to extend the dynamic range of delta(13)C and delta(15)N analysis from less than 2 orders of magnitude to more than 3 orders of magnitude. Carbon isotope ratio measurements of atropine as a model compound have been performed analysing between 1 μg to 5 mg C with acceptable to excellent precision (0.6 to 0.06 per thousand, delta-notation). The correction due to the blank signal is critical for sample amounts smaller than 4 μg C. The maximum sample weight is determined by the combustion capacity of the EA. Larger sample amounts are measured using dilution of a small part of the EA effluent with helium. The dilution mechanism works virtually free of isotope fractionation. Copyright 1999 John Wiley & Sons, Ltd.