A novel Bayesian approach to drug design with simple and complex enzymes: Use with lens aldehyde dehydrogenase and the glyoxalase pathway

Purpose: Acquiring details of kinetic parameters of enzymes is crucial to biochemical understanding, drug development, and clinical diagnosis in ocular diseases. The correct design of an experiment is critical to collecting data suitable for analysis, modelling and deriving the correct information. As classical design methods are not targeted to the more complex kinetics being frequently studied, attention is needed to estimate parameters of such models with low variance. Methods: We have developed Bayesian utility functions to minimise kinetic parameter variance involving differentiation of model expressions and matrix inversion. These have been applied to the simple kinetics of the enzymes in the glyoxalase pathway (of importance in posttranslational modification of proteins in cataract), and the complex kinetics of lens aldehyde dehydrogenase (also of relevance to cataract). Results: Our successful application of Bayesian statistics has allowed us to identify a set of rules for designing optimum kinetic experiments iteratively. Most importantly, the distribution of points in the range is critical; it is not simply a matter of even or multiple increases. At least 60 % must be below the KM (or plural if more than one dissociation constant) and 40% above. This choice halves the variance found using a simple even spread across the range.With both the glyoxalase system and lens aldehyde dehydrogenase we have significantly improved the variance of kinetic parameter estimation while reducing the number and costs of experiments. Conclusions: We have developed an optimal and iterative method for selecting features of design such as substrate range, number of measurements and choice of intermediate points. Our novel approach minimises parameter error and costs, and maximises experimental efficiency. It is applicable to many areas of ocular drug design, including receptor-ligand binding and immunoglobulin binding, and should be an important tool in ocular drug discovery.

Agonist binding, agonist affinity and agonist efficacy at G protein-coupled receptors

Measurements of affinity and efficacy are fundamental for work on agonists both in drug discovery and in basic studies on receptors. In this review I wish to consider methods for measuring affinity and efficacy at G protein coupled receptors (GPCRs). Agonist affinity may be estimated in terms of the dissociation constant for agonist binding to a receptor using ligand binding or functional assays. It has, however, been suggested that measurements of affinity are always contaminated by efficacy so that it is impossible to separate the two parameters. Here I show that for many GPCRs, if receptor/G protein coupling is suppressed, experimental measurements of agonist affinity using ligand binding (K-obs) provide quite accurate measures of the agonist microscopic dissociation constant (K-A). Also in pharmacological functional studies, good estimates of agonist dissociation constants are possible. Efficacy can be quantitated in several ways based on functional data ( maximal effect of the agonist (E-max), ratio of agonist dissociation constant to concentration of agonist giving half maximal effect in functional assay ( K-obs/ EC50), a combined parameter EmaxKobs/EC50). Here I show that EmaxKobs/EC50 provides the best assessment of efficacy for a range of agonists across the full range of efficacy for full to partial agonists. Considerable evidence now suggests that ligand efficacy may be dependent on the pathway used to assess it. The efficacy of a ligand may, therefore, be multidimensional. It is still, however, necessary to have accurate measures of efficacy in different pathways.

Effect of an ancillary ligand on single helix-double helix interconversion in copper complexes. Copper(I)-water bond

Reaction of Cu(ClO(4))(2)center dot 6H(2)O with the 1:2 condensate of benzildihydrazone and 2-acetylpyridine, in methanol in equimolar ratio yields a green compound which upon recrystallisation from 1:1 CH(2)Cl(2)-C(6)H(6) mixture affords [CuL(H(2)O)](ClO(4))(2)center dot 1/2C(6)H(6). The complex crystallises in the space group P-1 with a = 8.028(11) angstrom, b = 12.316(17) angstrom, c = 18.14(3) angstrom, alpha = 97.191(10)degrees, beta = 94.657(10)degrees and gamma = 108.039(10)degrees. It is single helical with the metal having a distorted trigonal bipyramidal N(4)O coordination sphere. The acid dissociation constant of the Cu(I) complex in CH(3)CN is 3.34 +/- 0.19. The X band EPR spectrum of the compound is rhombic with g(1) = 2.43, g(2) = 2.10 g(3) = 2.02 and A(1) = 79.3 x 10(-4) cm(-1). The Cu(II/I) potential of the complex in CH(2)Cl(2) at a glassy carbon electrode is 0.43 V vs SCE. It is argued that the copper-water bond persists in the corresponding copper(I) species. Its implications on the single helix-double helix interconversion in copper helicates are discussed. DFT calculations at the B3LYP/6-311G** level shows that the binding energy of water in the single helicol live-coordinate copper(I) species [CuL(H(2)O)](+) is similar to 40 kJ mol(-1).

Anharmonic force constant calculations

The relationship of the anharmonic force constants in curvilinear internal coordinates to the observed vibration-rotation spectrum of a molecule is reviewed. A simplified method of setting up the required non-linear coordinate transformations is described: this makes use of an / tensor, which is a straightforward generalization of the / matrix used in the customary description of harmonic force constant calculations. General formulae for the / tensor elements, in terms of the familiar L matrix elements, are presented. The use of non-linear symmetry coordinates and redundancies are described. Sample calculations on the water and ammonia molecules are reported.

Analytical potentials for triatomic molecules IX. The prediction of anharmonic force constants from potential energy surfaces based on harmonic force fields and dissociation energies for SO2 and O3

Analytical potential energy functions which are valid at all dissociation limits have been derived for the ground states of SO2 and O3. The procedure involves minimizing the errors between the observed vibrational spectra and spectra calculated by a variational procedure. Good agreement is obtained between the observed and calculated spectra for both molecules. Comparisons are made between anharmonic force fields, previously determined from the spectral data, and the force fields obtained by differentiating the derived analytical functions at the equilibrium configurations.

ASYM20: a program for force constant and normal coordinate calculations, with a critical review of the theory involved

The theory of harmonic force constant refinement calculations is reviewed, and a general-purpose program for force constant and normal coordinate calculations is described. The program, called ASYM20. is available through Quantum Chemistry Program Exchange. It will work on molecules of any symmetry containing up to 20 atoms and will produce results on a series of isotopomers as desired. The vibrational secular equations are solved in either nonredundant valence internal coordinates or symmetry coordinates. As well as calculating the (harmonic) vibrational wavenumbers and normal coordinates, the program will calculate centrifugal distortion constants, Coriolis zeta constants, harmonic contributions to the α′s. root-mean-square amplitudes of vibration, and other quantities related to gas electron-diffraction studies and thermodynamic properties. The program will work in either a predict mode, in which it calculates results from an input force field, or in a refine mode, in which it refines an input force field by least squares to fit observed data on the quantities mentioned above. Predicate values of the force constants may be included in the data set for a least-squares refinement. The program is written in FORTRAN for use on a PC or a mainframe computer. Operation is mainly controlled by steering indices in the input data file, but some interactive control is also implemented.

Change in the effect of constant photoperiods on the rate of sexual maturation in modern genotypes of domestic pullet

1. Data for modern egg-type hybrids reared on constant daylengths show that, as expected, they mature more quickly than earlier genotypes. However, the constant photoperiod which gives earliest sexual maturity has not changed as a result of selection and is 10 h for both early and modern genotypes. 2. Further analysis showed that the rate of delay in sexual maturity for constant photoperiods above 10 h is similar for modern and for early hybrids ( +0.29 d for each incremental one hour of photoperiod), the response of modern hybrids below 10 h ( +4.22 d for each one-hour reduction in photoperiod) is more than double that of early hybrids ( +1.71 d/h).

Time-resolved gas-phase kinetic and quantum chemical studies of the reaction of silylene with oxygen

Time-resolved kinetic studies of the reaction of silylene, SiH2, generated by laser flash photolysis of phenylsilane, have been carried out to obtain rate constants for its bimolecular reaction with O-2. The reaction was studied in the gas phase over the pressure range 1-100 Torr in SF6 bath gas, at five temperatures in the range 297-600 K. The second order rate constants at 10 Torr were fitted to the Arrhenius equation: log(k/cm(3) molecule(-1) s(-1)) = (-11.08 +/- 0.04) + (1.57 +/- 0.32 kJ mol(-1))/RT ln10 The decrease in rate constant values with increasing temperature, although systematic is very small. The rate constants showed slight increases in value with pressure at each temperature, but this was scarcely beyond experimental uncertainty. From estimates of Lennard-Jones collision rates, this reaction is occurring at ca. 1 in 20 collisions, almost independent of pressure and temperature. Ab initio calculations at the G3 level backed further by multi-configurational (MC) SCF calculations, augmented by second order perturbation theory (MRMP2), support a mechanism in which the initial adduct, H2SiOO, formed in the triplet state (T), undergoes intersystem crossing to the more stable singlet state (S) prior to further low energy isomerisation processes leading, via a sequence of steps, ultimately to dissociation products of which the lowest energy pair are H2O + SiO. The decomposition of the intermediate cyclo-siladioxirane, via O-O bond fission, plays an important role in the overall process. The bottleneck for the overall process appears to be the T -> S process in H2SiOO. This process has a small spin orbit coupling matrix element, consistent with an estimate of its rate constant of 1 x 10(9) s(-1) obtained with the aid of RRKM theory. This interpretation preserves the idea that, as in its reactions in general, SiH2 initially reacts at the encounter rate with O-2. The low values for the secondary reaction barriers on the potential energy surface account for the lack of an observed pressure dependence. Some comparisons are drawn with the reactions of CH2 + O-2 and SiCl2 + O-2.

Acid dissociation of 2,2 '-dipyridylamine in non-aqueous medium when chelated to some Ru(II)N-4 cores

[Ru(2,2'-bipyridine)(2)(Hdpa)](BF4)(2) center dot 2H(2)O (1), [Ru(1,10-phenanthroline)(2)(Hdpa)] (PF6)(2) center dot CH2Cl2 (2) and [Ru(4,4,4',4'-tetramethyl-2,2'- bisoxazoline)(2)(Hdpa)] (PF6)(2) (3) are synthesized where Hdpa is 2,2'-dipyridylamine. The X-ray crystal structures of 1 and 2 have been determined. Hdpa in 1 and 2 is found to bind the metal via the two pyridyl N ends. Comparing the NMR spectra in DMSO-d(6), it is concluded that 3 has a similar structure. The pK(a) values (for the dissociation of the NH proton in Hdpa) of free Hdpa and its complexes are determined in acetonitrile by exploiting molar conductance. These correlate linearly with the chemical shift of the NH proton in the respective entities. (C) 2007 Elsevier B.V. All rights reserved.

A monoaqua zinc complex. Unique acid dissociation behaviour

Using the I : 2 condensate of benzil dihydrazone and 2-acetylpyridine as the ligand L, two complexes of zinc, [ZnL(CH3COO)]PF6 (1) and [ZnL(H2O)CIO4]CIO4 H2O (2), are synthesised from Zn(CH3COO)(2).2H(2)O and Zn(CIO4)(2).6H(2)O, respectively. From X-ray crystallography, both the complexes are found to be single helical with the metal in distorted octahedral N4O2 environment. In 1, the two oxygen atoms come from the bidentate acetate while 2 is a monoaqua complex with a perchlorate anion bound to the metal in monodentate fashion. The perchlorate in 2 is not at all weakly bound [Zn-O(perchlorate) 2.256(4) A]. Still in acetonitrile solution, the coordinated perchlorate ion dissociates upon deprotonation [reaction (i)].

Dissociation of water on oxygen-covered Rh{111}

The adsorption of water and coadsorption with oxygen on Rh{111} under ultrahigh vacuum conditions was studied using synchrotron-based photoemission and photoabsorption spectroscopy. Water adsorbs intact on the clean surface at temperatures below 154 K. Irradiation with x-rays, however, induces fast dissociation and the formation of a mixed OH+H(2)O layer indicating that the partially dissociated layer is thermodynamically more stable. Coadsorption of water and oxygen at a coverage below 0.3 monolayers has a similar effect, leading to the formation of a hydrogen-bonded network of water and hydroxyl molecules at a ratio of 3:2. The partially dissociated layers are more stable than chemisorbed intact water with the maximum desorption temperatures up to 30 K higher. For higher oxygen coverage, up to 0.5 monolayers, water does not dissociate and an intact water species is observed above 160 K, which is characterized by an O 1s binding energy 0.6 eV higher than that of chemisorbed water and a high desorption temperature similar to the partially dissociated layer. The extra stabilization is most likely due to hydrogen bonds with atomic oxygen.

Positive schizotypy and trait dissociation as vulnerability factors for post-traumatic distress

Objective. This study investigated whether trait positive schizotypy or trait dissociation was associated with increased levels of data-driven processing and symptoms of post-traumatic distress following a road traffic accident. Methods. Forty-five survivors of road traffic accidents were recruited from a London Accident and Emergency service. Each completed measures of trait positive schizotypy, trait dissociation, data-driven processing, and post-traumatic stress. Results. Trait positive schizotypy was associated with increased levels of data-driven processing and post-traumatic symptoms during a road traffic accident, whereas trait dissociation was not. Conclusions. Previous results which report a significant relationship between trait dissociation and post-traumatic symptoms may be an artefact of the relationship between trait positive schizotypy and trait dissociation.