Bayesian design and conduct of phase II single-arm clinical trials with binary outcomes: a tutorial

The aim of phase II single-arm clinical trials of a new drug is to determine whether it has sufficient promising activity to warrant its further development. For the last several years Bayesian statistical methods have been proposed and used. Bayesian approaches are ideal for earlier phase trials as they take into account information that accrues during a trial. Predictive probabilities are then updated and so become more accurate as the trial progresses. Suitable priors can act as pseudo samples, which make small sample clinical trials more informative. Thus patients have better chances to receive better treatments. The goal of this paper is to provide a tutorial for statisticians who use Bayesian methods for the first time or investigators who have some statistical background. In addition, real data from three clinical trials are presented as examples to illustrate how to conduct a Bayesian approach for phase II single-arm clinical trials with binary outcomes.

A Bayesian approach for dose-escalation in a phase I clinical trial incorporating pharmacodynamic endpoints

Bayesian decision procedures have already been proposed for and implemented in Phase I dose-escalation studies in healthy volunteers. The procedures have been based on pharmacokinetic responses reflecting the concentration of the drug in blood plasma and are conducted to learn about the dose-response relationship while avoiding excessive concentrations. However, in many dose-escalation studies, pharmacodynamic endpoints such as heart rate or blood pressure are observed, and it is these that should be used to control dose-escalation. These endpoints introduce additional complexity into the modeling of the problem relative to pharmacokinetic responses. Firstly, there are responses available following placebo administrations. Secondly, the pharmacodynamic responses are related directly to measurable plasma concentrations, which in turn are related to dose. Motivated by experience of data from a real study conducted in a conventional manner, this paper presents and evaluates a Bayesian procedure devised for the simultaneous monitoring of pharmacodynamic and pharmacokinetic responses. Account is also taken of the incidence of adverse events. Following logarithmic transformations, a linear model is used to relate dose to the pharmacokinetic endpoint and a quadratic model to relate the latter to the pharmacodynamic endpoint. A logistic model is used to relate the pharmacokinetic endpoint to the risk of an adverse event.

Direct detection of chlorosilylene and time resolved study of some of its reactions in the gas-phase using a new photochemical precursor

Chlorosilylene, ClSiH, was prepared by 193 nm laser flash photolysis of 1-chloro-1-silacyclopent-3-ene in the gas phase. ClSiH was monitored in real time at 457.9 nm using a CW argon ion laser. The kinetics of reactions of ClSiH with C2H4, CH2 = CHCMe3, C2H2, Me2O, SO2, HCl, MeSiH3, Me2SiH2, Me3SiH, MeGeH3, MeGeH3 and precursor have been studied at ambient temperatures for the first time. Addition reactions of ClSiH and reactions with lone pair donors are faster than insertion reactions. Surprisingly ClSiH inserts faster into Si-H than Ge-H bonds. ClSiH is intermediate in reactivity between SiH2 and SiCI2. Relative reactivities of CISiH and SiH2 vary considerably. (c) 2005 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.

Gas phase kinetic and quantum chemical studies of the reactions of silylene with the methylsilanes. Absolute rate constants, temperature dependences, RRKM modelling and potential energy surfaces

Time resolved studies of silylene, SiH2, generated by the 193 nm laser. ash photolysis of phenylsilane, have been carried out to obtain rate coefficients for its bimolecular reactions with methyl-, dimethyl- and trimethyl-silanes in the gas phase. The reactions were studied over the pressure range 3 - 100 Torr with SF6 as bath gas and at five temperatures in the range 300 - 625 K. Only slight pressure dependences were found for SiH2 + MeSiH3 ( 485 and 602 K) and for SiH2 + Me2SiH2 ( 600 K). The high pressure rate constants gave the following Arrhenius parameters: [GRAPHICS] These are consistent with fast, near to collision-controlled, association processes. RRKM modelling calculations are consistent with the observed pressure dependences ( and also the lack of them for SiH2 + Me3SiH). Ab initio calculations at both second order perturbation theory (MP2) and coupled cluster (CCSD(T)) levels, showed the presence of weakly-bound complexes along the reaction pathways. In the case of SiH2 + MeSiH3 two complexes, with different geometries, were obtained consistent with earlier studies of SiH2 + SiH4. These complexes were stabilised by methyl substitution in the substrate silane, but all had exceedingly low barriers to rearrangement to product disilanes. Although methyl groups in the substrate silane enhance the intrinsic SiH2 insertion rates, it is doubtful whether the intermediate complexes have a significant effect on the kinetics. A further calculation on the reaction MeSiH + SiH4 shows that the methyl substitution in the silylene should have a much more significant kinetic effect ( as observed in other studies).

The gas-phase reactions Of SiCl4 and Si2Cl6 with CH3OH and C2H5OH: an investigation by mass spectrometry and matrix-isolation infrared spectroscopy

The gas phase reactions Of SiCl4 and Si2Cl6 With CH3OH and C2H5OH have been investigated using both mass spectrometry and matrix isolation techniques. SiCl4 reacts with both CH3OH and C2H5OH upon mixing of the vapours for times in excess of 3 h to generate the HCl-elimination products SiCl3OR (R = CH3 or C2H5). The identity of these products is confirmed by deuteration experiments and by ab initio calculations at the HF/6-31G(d) level. Further products are generated when the mixture is passed through a tube heated to 750degreesC. Si2Cl6 reacts with CH3OH and C2H5OH via a different mechanism in which the Si-Si bond is cleaved to yield SiCl3OR and HCl. Other products of the type SiCl4-n(OCH3)(n) are tentatively identified by a combination of mass spectrometric and matrix isolation measurements. These latter products indicate further replacement of Cl atoms by OR groups as a result of reaction of CH3OH or C2H5OH with the initial product.

The kinetics of the gas-phase decomposition of the 2-methyl-2-butoxyl and 2-methyl-2-pentoxyl radicals

The kinetics of the title reactions have been studied by relative-rate methods as a function of temperature. Relative-rate coefficients for the two decomposition channels of 2-methyl-2-butoxyl have been measured at five different temperatures between 283 and 345 K and the observed temperature dependence is consistent with the results of some previous experimental studies. The kinetics of the two decomposition channels of 2-methyl-2-pentoxyl have also been investigated, as a function of temperature, relative to the estimated rate of isomerisation of this radical. Room-temperature rate coefficient data for the two decomposition channels of both 2-methyl-2-pentoxyl and 2-methyl-2-butxoyl (after combining the relative rate coefficient for this latter with a value for the rate coefficient of the major channel, extrapolated from the data presented by Batt et al., Int. J. Chem. Kinet., 1978, 10, 931) are shown to be consistent with a non-linear kinetic correlation, for alkoxyl radical decomposition rate data, previously presented by this laboratory (Johnson et al., Atmos. Environ., 2004, 38, 1755-1765).

Organic acid formation in the gas-phase ozonolysis of alpha-pinene

The mechanism of formation of pinonic and norpinonic acids from alpha-pinene ozonolysis has been investigated by studying the products of the ozonolysis of an enone derived from alpha-pinene using gas chromatography coupled to mass spectrometry.

Multifunctional acid formation from the gas-phase ozonolysis of beta-pinene

The gas-phase ozonolysis of beta-pinene was studied in static chamber experiments, using gas chromatography coupled to mass spectrometric and flame ionisation detection to separate and detect products. A range of multifunctional organic acids-including pinic acid, norpinic acid, pinalic-3- acid, pinalic-4-acid, norpinalic acid and OH-pinalic acid-were identified in the condensed phase after derivatisation. Formation yields for these products under systematically varying reaction conditions (by adding different OH radical scavengers and Criegee intermediate scavengers) were investigated and compared with those observed from alpha-pinene ozonolysis, allowing detailed information on product formation mechanisms to be elucidated. In addition, branching ratios for the initial steps of the reaction were inferred from quantitative measurements of primary carbonyl formation. Atmospheric implications of this work are discussed.

Formation of organic acids from the gas-phase ozonolysis of terpinolene

Gas-phase ozonolysis of terpinolene was studied in static chamber experiments using gas chromatography coupled to mass spectrometric and flame ionisation detection to separate and detect products. Two isomers of C-7-diacids and three isomers of C-7-aldehydic acids were identified in the condensed phase after derivatisation. Possible mechanisms of formation of these acids were investigated using different OH radical scavengers and relative humidities, and were compared to those reported earlier for the ozonolysis of beta-pinene. In addition, branching ratios for some of the individual reaction steps, e. g. the branching ratio between the two hydroperoxide channels of the C-7-CI, were deduced from the quantitative product yield data. Branching ratios for POZ decomposition and the stabilisation/decomposition of the C-7-CI were also obtained from measurements of the C-7 primary carbonyl product.

Mechanisms for the formation of organic acids in the gas-phase ozonolysis of 3-carene

This paper describes experimental studies aimed at elucidating mechanisms for the formation of low-volatility organic acids in the gas-phase ozonolysis of 3-carene. Experiments were carried out in a static chamber under 'OH-free' conditions. A range of multifunctional acids-which are analogous to those observed from alpha-pinene ozonolysis-were identified in the condensed phase using gas chromatography coupled to mass spectrometry after derivation. Product yields were determined as a function of different OH radical scavengers and relative humidities to give mechanistic information about their routes of formation. Furthermore, an enone and an enal derived from 3-carene were ozonised in order to probe the early mechanistic steps in the reaction and, in particular, which of the two initially formed Criegee intermediates gives rise to which products. Branching ratios for the formation of the two Criegee Intermediates are determined. Similarities and differences in product formation from 3-carene and alpha-pinene ozonolysis are discussed and possible mechanisms-supported by experimental evidence-are developed for all acids investigated.

Mechanisms for the formation of secondary organic aerosol components from the gas-phase ozonolysis of alpha-pinene

Gas-phase ozonolysis of alpha-pinene was studied in static chamber experiments under 'OH-free' conditions. A range of multifunctional products-in particular low-volatility carboxylic acids-were identified in the condensed phase using gas chromatography coupled to mass spectrometry after derivatisation. The dependence of product yields on reaction conditions (humidity, choice of OH radical scavengers, added Criegee intermediate scavengers, NO2 etc.) was investigated to probe the mechanisms of formation of these products; additional information was obtained by studying the ozonolysis of an enal and an enone derived from alpha-pinene. On the basis of experimental findings, previously suggested mechanisms were evaluated and detailed gas-phase mechanisms were developed to explain the observed product formation. Atmospheric implications of this work are discussed.

Quantitative model for the kinetics of lyotropic phase transitions involving changes in monolayer curvature

We present a kinetic model for transformations between different self-assembled lipid structures. The model shows how data on the rates of phase transitions between mesophases of different geometries can be used to provide information on the mechanisms of the transformations and the transition states involved. This can be used, for example, to gain an insight into intermediate structures in cell membrane fission or fusion. In cases where the monolayer curvature changes on going from the initial to the final mesophase, we consider the phase transition to be driven primarily by the change in the relaxed curvature with pressure or temperature, which alters the relative curvature elastic energies of the two mesophase structures. Using this model, we have analyzed previously published kinetic data on the inter-conversion of inverse bicontinuous cubic phases in the 1-monoolein-30 wt% water system. The data are for a transition between QII(G) and QII(D) phases, and our analysis indicates that the transition state more closely resembles the QII(D) than the QII(G) phase. Using estimated values for the monolayer mean curvatures of the QII(G) and QII(D) phases of -0.123 nm(-1) and -0.133 nm(-1), respectively, gives values for the monolayer mean curvature of the transition state of between -0.131 nm(-1) and -0.132 nm(-1). Furthermore, we estimate that several thousand molecules undergo the phase transition cooperatively within one "cooperative unit", equivalent to 1-2 unit cells of QII(G) or 4-10 unit cells of QII(D).

Detection of neural correlates of self-paced motor activity using empirical mode decomposition phase locking analysis

Transient episodes of synchronisation of neuronal activity in particular frequency ranges are thought to underlie cognition. Empirical mode decomposition phase locking (EMDPL) analysis is a method for determining the frequency and timing of phase synchrony that is adaptive to intrinsic oscillations within data, alleviating the need for arbitrary bandpass filter cut-off selection. It is extended here to address the choice of reference electrode and removal of spurious synchrony resulting from volume conduction. Spline Laplacian transformation and independent component analysis (ICA) are performed as pre-processing steps, and preservation of phase synchrony between synthetic signals. combined using a simple forward model, is demonstrated. The method is contrasted with use of bandpass filtering following the same preprocessing steps, and filter cut-offs are shown to influence synchrony detection markedly. Furthermore, an approach to the assessment of multiple EEG trials using the method is introduced, and the assessment of statistical significance of phase locking episodes is extended to render it adaptive to local phase synchrony levels. EMDPL is validated in the analysis of real EEG data, during finger tapping. The time course of event-related (de)synchronisation (ERD/ERS) is shown to differ from that of longer range phase locking episodes, implying different roles for these different types of synchronisation. It is suggested that the increase in phase locking which occurs just prior to movement, coinciding with a reduction in power (or ERD) may result from selection of the neural assembly relevant to the particular movement. (C) 2009 Elsevier B.V. All rights reserved.