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.

Rhodium assisted C-H activation of N-(2 '-hydroxyphenyl)benzaldimines. Synthesis, structure and electrochemical properties of a group of organorhodium complexes

Reaction of a group of N-(2'-hydroxyphenyl)benzaldimines, derived from 2-aminophenol and five para-substituted benzaldehydes (the para substituents are OCH3, CH3, H, Cl and NO2), with [Rh(PPh3)(3)Cl] in refluxing toluene in the presence of a base (NEW afforded a family of organometallic complexes of rhodium(III). The crystal structure of one complex has been determined by X-ray crystallography. In these complexes the benzaldimine ligands are coordinated to the metal center, via dissociation of the phenolic proton and the phenyl proton at the ortho position of the phenyl ring in the imine fragment, as dianionic tridentate C,N,O-donors, and the two PPh3 ligands are trans. The complexes are diamagnetic (low-spin d(6), S = 0) and show intense MLCT transitions in the visible region. Cyclic voltammetry shows a Rh(III)-Rh(IV) oxidation within 0.63-0.93 V vs SCE followed by an oxidation of the coordinated benzaldimine ligand. A reduction of the coordinated benzaldimine is also observed within -0.96 to -1.04 V vs SCE. Potential of the Rh(Ill)-Rh(IV) oxidation is found to be sensitive to the nature of the para-substituent. (c) 2006 Elsevier B.V. All rights reserved.

An investigation of the germylene addition reaction, GeH2+C2H2: Time-resolved gas-phase kinetic studies and quantum chemical calculations of the reaction energy surface

Time resolved studies of germylene, GeH2, generated by laser flash photolysis of 3,4-dimethylgermacyclopentene-3, have been carried out to obtain rate constants for its bimolecular reaction with acetylene, C2H2. The reaction was studied in the gas-phase over the pressure range 1-100 Tort, with SF6 as bath gas, at 5 temperatures in the range 297-553 K. The reaction showed a very slight pressure dependence at higher temperatures. The high pressure rate constants (obtained by extrapolation at the three higher temperatures) gave the Arrhenius equation: log(k(infinity)/cm(3) molecule(-1) s(-1)) (-10.94 +/- 0.05) + (6.10 +/- 0.36 kJ mol(-1))/RTln10. These Arrhenius parameters are consistent with a fast reaction occurring at approximately 30% of the collision rate at 298 K. Quantum chemical calculations (both DFT and ab initio G2//B3LYP and G2//QCISD) of the GeC2H4 potential energy surface (PES), show that GeH2 + C2H2 react initially to form germirene which can isomerise to vinylgermylene with a relatively low barrier. RRKM modelling, based on a loose association transition state, but assuming vinylgermylene is the end product (used in combination with a weak collisional deactivation model) predicts a strong pressure dependence using the calculated energies, in conflict with the experimental evidence. The detailed GeC2H4 PES shows considerable complexity with ten other accessible stable minima (B3LYP level), the three most stable of which are all germylenes. Routes through this complex surface were examined in detail. The only product combination which appears capable of satisfying the (P-3) + C2H4.C2H4 was confirmed as a product by GC observed lack of a strong pressure dependence is Ge(P-3) + C2H4. C2H4 was confirmed as a product by GC analysis. Although the formation of these products are shown to be possible by singlet-triplet curve crossing during dissociation of 1-germiranylidene (1-germacyclopropylidene), it seems more likely (on thermochernical grounds) that the triplet biradical, (GeCH2CH2.)-Ge-., is the immediate product precursor. Comparisons are made with the reaction of SiH2 with C2H2.

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.

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

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 NO. 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 299-592 K. The second-order rate constants at 10 Torr fitted the Arrhenius equation log(k/cm(3) molecule(-1) s(-1)) = (- 11.66 +/- 0.01) + (6.20 +/- 0.10 kJ mol(-1))IRT In 10 The rate constants showed a variation with pressure of a factor of ca. 2 over the available range, almost independent of temperature. The data could not be fitted by RRKM calculations to a simple third body assisted association reaction alone. However, a mechanistic model with an additional (pressure independent) side channel gave a reasonable fit to the data. Ab initio calculations at the G3 level supported a mechanism in which the initial adduct, bent H2SiNO, can ring close to form cyclo-H2SiNO, which is partially collisionally stabilized. In addition, bent H2SiNO can undergo a low barrier isomerization reaction leading, via a sequence of steps, ultimately to dissociation products of which the lowest energy pair are NH2 + SiO. The rate controlling barrier for this latter pathway is only 16 kJ mol(-1) below the energy of SiH2 + NO. This is consistent with the kinetic findings. A particular outcome of this work is that, despite the pressure dependence and the effects of the secondary barrier (in the side reaction), the initial encounter of SiH2 with NO occurs at the collision rate. Thus, silylene can be as reactive with odd electron molecules as with many even electron species. Some comparisons are drawn with the reactions of CH2 + NO and SiCl2 + NO.

Reactions of silylene with unreactive molecules. 2. nitrogen: gas-phase kinetic and theoretical studies

Time-resolved studies of the reaction of silylene, SiH2, with N-2 have been attempted at 296, 417, and 484 K, using laser flash photolysis to generate and monitor SiH2. No conclusive evidence for reaction could be found even with pressures of N-2 of 500 Torr. This enables us to set upper limits of ca. 3 x 10(-15) cm(3) molecule(-1) s(-1) for the second-order rate constants. A lower limit for the activation energy, E-a, of ca. 47 kJ mol(-1) is also derived. Ab initio calculations at the G3 level indicate that the only SiH2N2 species of lower energy than the separated reactants is the H2Si...N-2 donor-acceptor (ylid) species with a relative enthalpy of -26 kJ mol(-1), insufficient for observation of reaction under the experimental conditions. Ten bound species on the SiH2N2 surface were found and their energies calculated as well as those of the potential dissociation products: HSiN + NH((3)Sigma(-)) and HNSi + NH((3)Sigma(-)). Additionally two of the transition states involving cyclic-SiH2N2 (siladiazirine) were explored. It appears that siladiazirine is neither thermodynamically nor kinetically stable. The findings indicate that Si-N-d bonds (where N-d is double-bonded nitrogen) are not particularly strong. An unexpected cyclic intermediate was found in the isomerization of silaisocyanamide to silacyanamide.

Direct detection of dimethylstannylene and tetramethyldistannene in solution and the gas phase by laser flash photolysis of 1,1-dimethylstannacyclopent-3-enes

The photochemistry of 1,1-dimethyl- and 1,1,3,4-tetramethylstannacyclopent-3-ene (4a and 4b,respectively) has been studied in the gas phase and in hexane solution by steady-state and 193-nm laser flash photolysis methods. Photolysis of the two compounds results in the formation of 1,3-butadiene (from 4a) and 2,3-dimethyl-1,3-butadiene (from 4b) as the major products, suggesting that cycloreversion to yield dimethylstannylene (SnMe2) is the main photodecomposition pathway of these molecules. Indeed, the stannylene has been trapped as the Sn-H insertion product upon photolysis of 4a in hexane containing trimethylstannane. Flash photolysis of 4a in the gas phase affords a transient absorbing in the 450-520nm range that is assigned to SnMe2 by comparison of its spectrum and reactivity to those previously reported from other precursors. Flash photolysis of 4b in hexane solution affords results consistent with the initial formation of SnMe2 (lambda(max) approximate to 500 nm), which decays over similar to 10 mu s to form tetramethyldistannene (5b; lambda(max) approximate to 470 nm). The distannene decays over the next ca. 50 mu s to form at least two other longer-lived species, which are assigned to higher SnMe2 oligomers. Time-dependent DFT calculations support the spectral assignments for SnMe2 and Sn2Me4, and calculations examining the variation in bond dissociation energy with substituent (H, Me, and Ph) in disilenes, digermenes, and distannenes rule out the possibility that dimerization of SnMe2 proceeds reversibly. Addition of methanol leads to reversible reaction with SnMe2 to form a transient absorbing at lambda(max) approximate to 360 nm, which is assigned to the Lewis acid-base complex between SnMe2 and the alcohol.

Statistical deconvolution of enthalpic energetic contributions to MHC-peptide binding affinity

Background: MHC Class I molecules present antigenic peptides to cytotoxic T cells, which forms an integral part of the adaptive immune response. Peptides are bound within a groove formed by the MHC heavy chain. Previous approaches to MHC Class I-peptide binding prediction have largely concentrated on the peptide anchor residues located at the P2 and C-terminus positions. Results: A large dataset comprising MHC-peptide structural complexes was created by remodelling pre-determined x-ray crystallographic structures. Static energetic analysis, following energy minimisation, was performed on the dataset in order to characterise interactions between bound peptides and the MHC Class I molecule, partitioning the interactions within the groove into van der Waals, electrostatic and total non-bonded energy contributions. Conclusion: The QSAR techniques of Genetic Function Approximation (GFA) and Genetic Partial Least Squares (G/PLS) algorithms were used to identify key interactions between the two molecules by comparing the calculated energy values with experimentally-determined BL50 data. Although the peptide termini binding interactions help ensure the stability of the MHC Class I-peptide complex, the central region of the peptide is also important in defining the specificity of the interaction. As thermodynamic studies indicate that peptide association and dissociation may be driven entropically, it may be necessary to incorporate entropic contributions into future calculations.

Chemical composition and reactivity of water on clean and oxygen-covered Pd{111}

The chemical composition and dissociation behaviour of water adsorbed on clean and oxygen pre-covered Pd{111} was studied using high-resolution time-resolved and temperature-programmed X-ray photoelectron spectroscopy. We find that water remains intact at all temperatures up to desorption on the clean surface and at high oxygen coverage(0.69 ML) when a surface oxide is formed. The highest desorption peaks occur at 163 K from the clean surface and at 172 K from the surface oxide. At the intermediate coverage of 0.20 ML oxygen reacts with coadsorbed water at 155 K, to generate a mixed H2O/OH layer exhibiting a (root 3- x root 3)R30 degrees diffraction pattern, which is stable up to 177 K. The measured ratio between intact water and the hydroxyl species in this layer varies between 1.5 and 2 depending on temperature. (C) 2008 Elsevier B.V. All rights reserved.

Modifying the adsorption characteristics of water on Ru{0001} with preadsorbed oxygen

The coadsorption of water and preadsorbed oxygen on Ru{0001) was studied by synchrotron-based high-resolution x-ray photoelectron spectroscopy. A dramatic change was observed in the interaction of water with oxygen between low and high oxygen precoverages. Low oxygen coverages below 0.18 ML induce partial dissociation, which leads to an adsorbed layer of H2O and OH. Around half the oxygen atoms take part in this reaction. All OH recombines upon heating to 200 K and desorbs together with H2O. Oxygen coverages between 0.20 and 0.50 ML inhibit dissociation, instead a highly stable intact water species is observed, which desorbs at 220 K. This species is significantly more stable than intact water on the clean surface. The stabilization is most likely due to the formation of hydrogen bonds with neighboring oxygen atoms. For intermediate oxygen coverages around 0.18 ML, the dissociation behavior depends on the preparation conditions, which points toward possible mechanisms and pathways for partial dissociation of water on Ru{0001}.

Synergetic effects of the Cu/Pt{110} surface alloy: enhanced reactivity of water and carbon monoxide

We have used synchrotron-based high-resolution X-ray photoelectron spectroscopy in combination with ab initio density functional theory calculations to investigate the characteristics of water and CO adsorption on the bimetallic Cu/Pt{110}-(2 x 1) surface at a Cu coverage near 0.5 ML. Cu fills the troughs of the reconstructed clean surface forming nanowires, which are stable up to 830 K. Their presence dramatically influences the adsorption of water and CO. Water adsorption changes from intact to partially dissociated while the desorption temperature of CO on this surface increases by up to 27 K with respect to the clean Pt{110} surface. Ab initio calculations and experimental valence band spectra reveal that the Cu 3d-band is narrowed and shifted upward with respect to bulk Cu surfaces. This and electron donation to surface Pt atoms cause the increase in the bond strength between CO and the Pt surface atoms. The pathway for water dissociation occurs via Cu surface atoms. The heat of adsorption of water bonding to Cu surface atoms was calculated to be 0.82 eV, which is significantly higher than on the clean Pt{110} surface; the activation energy for partial dissociation is 0.53 eV (not corrected for zero point energy).

Palladium and platinum complexes of 2-(2 '-carboxyphenylazo)-4-methylphenol: synthesis, structure and spectral properties

Reaction of 2-(2'-carboxyphenylazo)-4-methylphenol (H2L) with [M(PPh3)(2)Cl-2] (M = Pd, Pt) affords mixed-ligand complexes of type [M(PPh3)(L)]. Structures of both the complexes have been determined by X-ray crystallography. Both the complexes are square planar, where the 2-(2'-carboxyphenylazo)-4-methylphenol is coordinated to the metal center, via dissociation of the two acidic protons, as a dianionic tridentate O,N,O-donor, and the fourth position is occupied by the triphenylphosphine. These complexes show intense MLCT transitions in the visible region.

Soluble redox-active polymetallic chains [{Ru-0(CO)(L)(bpy)}(m)](n) (bpy=2,2 '-bipyridine, L = PrCN, Cl-; m=0,-1): electrosynthesis and characterization

Electrochemical and spectroelectrochemical techniques were employed to study in detail the formation and so far unreported spectroscopic properties of soluble electroactive molecular chains with nonbridged metal-metal backbones, namely, [{Ru-0(CO)(PrCN)(bpy)}(m)](n) (m = 0, -1) and [{Ru-0(CO)(bpy)Cl}(m)](n) (m = -1, -2; bpy = 2,2'-bipyridine). The precursors cis-(Cl)-[Ru-II(CO)(MeCN)(bpy)Cl-2] (in PrCN) and mer-[Ru-II(CO)(bpy)Cl-3](-) (in tetrahydrofuran (THF) and PrCN) undergo one-electron reductions to reactive radicals [Ru-II(CO)(MeCN)(bpy(center dot-))Cl-2](-) and [Ru-II(CO)(bpy(center dot-))Cl-3](2-), respectively. Both [bpy(center dot-)]-containing species readily electropolymerize on concomitant dissociation of two chloride ligands and consumption of a second electron. Along this path, mer-to-fac isomerization of the bpy-reduced trichlorido complex (supported by density functional theory calculations) and a concentration-dependent oligomerization process contribute to the complex reactivity pattern. In situ spectroelectrochemistry (IR, UV/vis a has revealed that the charged polymer [{Ru-0(CO)(bpy)Cl}(-)](n) is stable in THF, but in PrCN it converts readily to [Ru-0(CO)(PrCN)(bpy)](n). An excess of chloride ions retards this substitution at low temperatures. Both polymetallic chains are completely soluble in the electrolyte solution and can be reduced reversibly to the corresponding [bpy(center dot-)]-containing species.

Interactions of bovine serum albumin with ethylene oxide/butylene oxide copolymers in aqueous solution

The interactions of bovine serum albumin (BSA) with three ethylene oxide/butylene oxide (E/B) copolymers having different block lengths and varying molecular architectures is examined in this study in aqueous solutions. Dynamic light scattering (DLS) indicates the absence of BSA-polymer binding in micellar systems of copolymers with lengthy hydrophilic blocks. On the contrary, stable protein-polyrner aggregates were observed in the case of E18B10 block copolymer. Results from DLS and SAXS suggest the dissociation of E/B copolymer micelles in the presence of protein and the absorption of polymer chains to BSA surface. At high protein loadings, bound BSA adopts a more compact conformation in solution. The secondary structure of the protein remains essentially unaffected even at high polymer concentrations. Raman spectroscopy was used to give insight to the configurations of the bound molecules in concentrated solutions. In the vicinity of the critical gel concentration of E18B10 introduction of BSA can dramatically modify the phase diagram, inducing a gel-sol-gel transition. The overall picture of the interaction diagram of the E18B10-BSA reflects the shrinkage of the suspended particles due to destabilization of micelles induced by BSA and the gelator nature of the globular protein. SAXS and rheology were used to further characterize the structure and flow behavior of the polymer-protein hybrid gels and sols.

Full-dimensional vibrational calculations for H5O2+ using an ab initio potential energy surface

We report quantum diffusion Monte Carlo (DMC) and variational calculations in full dimensionality for selected vibrational states of H5O2+ using a new ab initio potential energy surface [X. Huang, B. Braams, and J. M. Bowman, J. Chem. Phys. 122, 044308 (2005)]. The energy and properties of the zero-point state are focused on in the rigorous DMC calculations. OH-stretch fundamentals are also calculated using "fixed-node" DMC calculations and variationally using two versions of the code MULTIMODE. These results are compared with infrared multiphoton dissociation measurements of Yeh [L. I. Yeh, M. Okumura, J. D. Myers, J. M. Price, and Y. T. Lee, J. Chem. Phys. 91, 7319 (1989)]. Some preliminary results for the energies of several modes of the shared hydrogen are also reported.