Sequential designs for phase III clinical trials incorporating treatment selection

Most statistical methodology for phase III clinical trials focuses on the comparison of a single experimental treatment with a control. An increasing desire to reduce the time before regulatory approval of a new drug is sought has led to development of two-stage or sequential designs for trials that combine the definitive analysis associated with phase III with the treatment selection element of a phase II study. In this paper we consider a trial in which the most promising of a number of experimental treatments is selected at the first interim analysis. This considerably reduces the computational load associated with the construction of stopping boundaries compared to the approach proposed by Follman, Proschan and Geller (Biometrics 1994; 50: 325-336). The computational requirement does not exceed that for the sequential comparison of a single experimental treatment with a control. Existing methods are extended in two ways. First, the use of the efficient score as a test statistic makes the analysis of binary, normal or failure-time data, as well as adjustment for covariates or stratification straightforward. Second, the question of trial power is also considered, enabling the determination of sample size required to give specified power. Copyright © 2003 John Wiley & Sons, Ltd.

Origin of avian genome size and structure in non-avian dinosaurs

Avian genomes are small and streamlined compared with those of other amniotes by virtue of having fewer repetitive elements and less non-coding DNA(1,2). This condition has been suggested to represent a key adaptation for flight in birds, by reducing the metabolic costs associated with having large genome and cell sizes(3,4). However, the evolution of genome architecture in birds, or any other lineage, is difficult to study because genomic information is often absent for long-extinct relatives. Here we use a novel bayesian comparative method to show that bone-cell size correlates well with genome size in extant vertebrates, and hence use this relationship to estimate the genome sizes of 31 species of extinct dinosaur, including several species of extinct birds. Our results indicate that the small genomes typically associated with avian flight evolved in the saurischian dinosaur lineage between 230 and 250 million years ago, long before this lineage gave rise to the first birds. By comparison, ornithischian dinosaurs are inferred to have had much larger genomes, which were probably typical for ancestral Dinosauria. Using comparative genomic data, we estimate that genome-wide interspersed mobile elements, a class of repetitive DNA, comprised 5 - 12% of the total genome size in the saurischian dinosaur lineage, but was 7 - 19% of total genome size in ornithischian dinosaurs, suggesting that repetitive elements became less active in the saurischian lineage. These genomic characteristics should be added to the list of attributes previously considered avian but now thought to have arisen in non-avian dinosaurs, such as feathers(5), pulmonary innovations 6, and parental care and nesting

Likelihood-based estimation of microsatellite mutation rates

Microsatellites are widely used in genetic analyses, many of which require reliable estimates of microsatellite mutation rates, yet the factors determining mutation rates are uncertain. The most straightforward and conclusive method by which to study mutation is direct observation of allele transmissions in parent-child pairs, and studies of this type suggest a positive, possibly exponential, relationship between mutation rate and allele size, together with a bias toward length increase. Except for microsatellites on the Y chromosome, however, previous analyses have not made full use of available data and may have introduced bias: mutations have been identified only where child genotypes could not be generated by transmission from parents' genotypes, so that the probability that a mutation is detected depends on the distribution of allele lengths and varies with allele length. We introduce a likelihood-based approach that has two key advantages over existing methods. First, we can make formal comparisons between competing models of microsatellite evolution; second, we obtain asymptotically unbiased and efficient parameter estimates. Application to data composed of 118,866 parent-offspring transmissions of AC microsatellites supports the hypothesis that mutation rate increases exponentially with microsatellite length, with a suggestion that contractions become more likely than expansions as length increases. This would lead to a stationary distribution for allele length maintained by mutational balance. There is no evidence that contractions and expansions differ in their step size distributions.

A SAXS study of flow alignment of thermotropic liquid crystal mixtures

We report on the capillary flow behaviour of thermotropic liquid crystal mixtures containing 4-n-octyl-4'-cyanobiphenyl (8CB) and 4-n-pentyl-4'-cyanobiphenyl (5CB). The liquid crystal mixtures are studied in the Nematic (N) and Smectic (SA) phases at room temperature. Polarised optical microscopy (POM), rheology and simultaneous X-ray diffraction (XRD)/capillary flow experiments are performed to characterise the system. Polarised optical microscopy reveals a dramatic change in optical texture when the 5CB content is increased from 20 to 30% in the mixtures. X-ray diffraction results show that the system goes through a SA-N phase transition, such that the mixtures are smectic for 10-20% 5CB and nematic for 30-90% 5CB. Smectic mixtures flow with the layers aligned along the flow direction (mesogens perpendicular to flow) while nematic mixtures flow with the mesogens aligned in the flow direction. Simultaneous XRD/shear flow experiments show that the SA-N transition is independent of the flow rate in the range 1-6 ml min-1. The correlation length of the liquid crystal order decreases with increasing 5CB content. Rheology is used to prove that the correlation length behaviour is related to a reduction in the viscosity of the mixture.

Self-organisation in the assembly of gels from mixtures of different dendritic peptide building blocks

This paper investigates dendritic peptides capable of assembling into nanostructured gels, and explores the effect on self-assembly of mixing different molecular building blocks. Thermal measurements, small angle Xray scattering (SAXS) and circular dichroism (CD) spectroscopy are used to probe these materials on macroscopic, nanoscopic and molecular length scales. The results from these investigations demonstrate that in this case, systems with different "size" and "chirality" factors can self-organise, whilst systems with different "shape" factors cannot. The "size" and "chirality" factors are directly connected with the molecular information programmed into the dendritic peptides, whilst the shape factor depends on the group linking these peptides together-this is consistent with molecular recognition hydrogen bond pathways between the peptidic building blocks controlling the ability of these systems to self-recognise. These results demonstrate that mixtures of relatively complex peptides, with only subtle differences on the molecular scale, can self-organise into nanoscale structures, an important step in the spontaneous assembly of ordered systems from complex mixtures.

Heterogeneously catalyzed asymmetric hydrogenation of C = C bonds directed by surface-tethered chiral modifiers

Asymmetric hydrogenation of C=C bonds is of the highest importance in organic synthesis, and such reactions are currently carried out with organometallic homogeneous catalysts. Achieving heterogeneous metal-catalyzed hydrogenation, a highly desirable goal, necessitates forcing the crucial enantiodifferentiating step to take place at the metal surface. By synthesis and application of six chiral sulfide ligands that anchor robustly to Pd nanoparticles and resist displacement, we have for the first time accomplished heterogeneous enantioselective catalytic hydrogenation of isophorone. High resolution XPS data established that ligand adsorption from solution occurred exclusively on the Pd nanoparticles and not on the carbon support. All ligands contained a pyrrolidine nitrogen to enable their interaction with the isophorone substrate while the sulfide functionality provided the required interaction with the Pd surface. Enantioselective turnover numbers of up to similar to 100 product molecules per ligand molecule were found with a very large variation in asymmetric induction between ligands: observed enantiomeric excesses increased with increasing size of the alkyl group in the sulfide. This likely reflects varying degrees of ligand dispersion on the surface: bulky substituent groups hinder close approach of ligand molecules to each other, inhibiting close-packed island formation, favoring dispersion as separate molecules, and leading to effective asymmetric induction. Conversely, small substituents favor island formation leading to very low asymmetric induction. Enantioselective reaction most likely involves initial formation of an enamine or iminium species, confirmed by use of an analogous tertiary amine, which leads to racemic product. Ligand rigidity and resistance to self-assembled monolayer formation are important attributes that should be designed into improved chiral modifiers.

Cannabinoids stimulate fibroblastic colony formation by bone marrow cells indirectly via CB2 receptors

Recently, the cannabinoid receptors CB1 and CB2 were shown to modulate bone formation and resorption in vivo, although little is known of the mechanisms underlying this. The effects of cannabinoids on mesenchymal stem cell (MSC) recruitment in whole bone marrow were investigated using either the fibroblastic colony-forming unit (CFU-f) assay or high-density cultures of whole bone marrow. Levels of the CB1 and CB2 receptors were assessed by flow cytometry. Treatment of CFU-f cultures with the endocannabinoid 2-arachidonylglycerol (2-AG) dose-dependently increased fibroblastic and differentiated colony formation along with colony size. The nonspecific agonists CP 55,940 and WIN 55,212 both increased colony numbers, as did the CB2 agonists BML190 and JWH015. The CB1-specific agonist ACEA had no effect, whereas the CB2 antagonist AM630 blocked the effect of the natural cannabinoid tetrahydrocannabivarin, confirming mediation via the CB2 receptor. Treatment of primary bone marrow cultures with 2-AG stimulated proliferation and collagen accumulation, whereas treatment of subcultures of MSC had no effect, suggesting that the target cell is not the MSC but an accessory cell present in bone marrow. Subcultures of MSCs were negative for CB1 and CB2 receptors as shown by flow cytometry, whereas whole bone marrow contained a small population of cells positive for both receptors. These data suggest that cannabinoids may stimulate the recruitment of MSCs from the bone marrow indirectly via an accessory cell and mediated via the CB2 receptor. This recruitment may be one mechanism responsible for the increased bone formation seen after cannabinoid treatment in vivo.

Variable step LMS algorithm using the accumulated instantaneous error concept

The convergence speed of the standard Least Mean Square adaptive array may be degraded in mobile communication environments. Different conventional variable step size LMS algorithms were proposed to enhance the convergence speed while maintaining low steady state error. In this paper, a new variable step LMS algorithm, using the accumulated instantaneous error concept is proposed. In the proposed algorithm, the accumulated instantaneous error is used to update the step size parameter of standard LMS is varied. Simulation results show that the proposed algorithm is simpler and yields better performance than conventional variable step LMS.

Blind equalisation of multilevel PAM data for nonminimum phase channels via second- and fourth-order cumulants

This paper introduces a new blind equalisation algorithm for the pulse amplitude modulation (PAM) data transmitted through nonminimum phase (NMP) channels. The algorithm itself is based on a noncausal AR model of communication channels and the second- and fourth-order cumulants of the received data series, where only the diagonal slices of cumulants are used. The AR parameters are adjusted at each sample by using a successive over-relaxation (SOR) scheme, a variety of the ordinary LMS scheme, but with a faster convergence rate and a greater robustness to the selection of the ‘step-size’ in iterations. Computer simulations are implemented for both linear time-invariant (LTI) and linear time-variant (LTV) NMP channels, and the results show that the algorithm proposed in this paper has a fast convergence rate and a potential capability to track the LTV NMP channels.

Use of water-miscible retinyl palmitate as markers of chylomicrons gives earlier peak response of plasma retinyl esters compared with oil-soluble retinyl palmitate

Delayed peak response of plasma retinyl esters (RE) relative to plasma triacylglycerols (TAG) and apolipoprotein (Apo) B-48 responses following a fat load supplemented with vitamin A raised doubts about the use of vitamin A to label dietary-derived lipids and lipoproteins. The present study compared the use of water-miscible and oil-soluble retinyl palmitate (RP) as markers of dietary-derived lipoproteins in healthy subjects along with the measurements of postprandial plasma TAG and ApoB-48 responses to investigate whether the delayed peak response observed was due to delayed intestinal output of RE from oil-based solutions. Nine healthy female subjects were given a standard test meal containing a dose (112 mg) of RP in either water-miscible or oil-soluble form in random order, on two separate occasions after a 12 h overnight fast. The results showed that the mean plasma RE concentrations reached a peak significantly later than mean plasma TAG and ApoB-48 concentrations when oil-soluble RP was consumed, whereas plasma RE peaked earlier relative to plasma TAG and ApoB-48 responses when water-miscible RP was used. The results suggested a more rapid absorption with a significantly higher and earlier peak response of plasma RE when water-miscible RP was consumed. This was in contrast to the delayed initial appearance and later sustained higher concentrations of plasma RE during the late postprandial period when oil-soluble RP was consumed. The RE response to the water-miscible RP showed better concordance with plasma TAG response than that of oil-soluble RP.

The relationship between size, diversification and risk

Property portfolio diversification takes many forms, most of which can be associated with asset size. In other words larger property portfolios are assumed to have greater diversification potential than small portfolios. In addition, since greater diversification is generally associated with lower risk it is assumed that larger property portfolios will also have reduced return variability compared with smaller portfolios. If large property portfolios can simply be regarded as scaled-up, better-diversified versions of small property portfolios, then the greater a portfolio’s asset size, the lower its risk. This suggests a negative relationship between asset size and risk. However, if large property portfolios are not simply scaled-up versions of small portfolios, the relationship between asset size and risk may be unclear. For instance, if large portfolios hold riskier assets or pursue more volatile investment strategies, it may be that a positive relationship between asset size and risk would be observed, even if large property portfolios are more diversified. This paper tests the empirical relationship between property portfolio size, diversification and risk, in Institutional portfolios in the UK, during the period from 1989 to 1999 to determine which of these two characterisations is more appropriate.

On the momentum fluxes associated with mountain waves in directionally sheared flows

The direct impact of mountain waves on the atmospheric circulation is due to the deposition of wave momentum at critical levels, or levels where the waves break. The first process is treated analytically in this study within the framework of linear theory. The variation of the momentum flux with height is investigated for relatively large shears, extending the authors’ previous calculations of the surface gravity wave drag to the whole atmosphere. A Wentzel–Kramers–Brillouin (WKB) approximation is used to treat inviscid, steady, nonrotating, hydrostatic flow with directional shear over a circular mesoscale mountain, for generic wind profiles. This approximation must be extended to third order to obtain momentum flux expressions that are accurate to second order. Since the momentum flux only varies because of wave filtering by critical levels, the application of contour integration techniques enables it to be expressed in terms of simple 1D integrals. On the other hand, the momentum flux divergence (which corresponds to the force on the atmosphere that must be represented in gravity wave drag parameterizations) is given in closed analytical form. The momentum flux expressions are tested for idealized wind profiles, where they become a function of the Richardson number (Ri). These expressions tend, for high Ri, to results by previous authors, where wind profile effects on the surface drag were neglected and critical levels acted as perfect absorbers. The linear results are compared with linear and nonlinear numerical simulations, showing a considerable improvement upon corresponding results derived for higher Ri.

Asymptotic gravity wave drag expressions for non-hydrostatic rotating flow over a ridge

Asymptotic expressions are derived for the mountain wave drag in flow with constant wind and static stability over a ridge when both rotation and non-hydrostatic effects are important. These expressions, which are much more manageable than the corresponding exact drag expressions (when these do exist) are found to provide accurate approximations to the drag, even when non-hydrostatic and rotation effects are strong, despite having been developed for cases where these effects are weak. The derived expressions are compared with approximations to the drag found previously, and their asymptotic behaviour in various limits is studied.

An analytical model of mountain wave drag for wind profiles with shear and curvature

An analytical model is developed to predict the surface drag exerted by internal gravity waves on an isolated axisymmetric mountain over which there is a stratified flow with a velocity profile that varies relatively slowly with height. The model is linear with respect to the perturbations induced by the mountain, and solves the Taylor–Goldstein equation with variable coefficients using a Wentzel–Kramers–Brillouin (WKB) approximation, formally valid for high Richardson numbers, Ri. The WKB solution is extended to a higher order than in previous studies, enabling a rigorous treatment of the effects of shear and curvature of the wind profile on the surface drag. In the hydrostatic approximation, closed formulas for the drag are derived for generic wind profiles, where the relative magnitude of the corrections to the leading-order drag (valid for a constant wind profile) does not depend on the detailed shape of the orography. The drag is found to vary proportionally to Ri21, decreasing as Ri decreases for a wind that varies linearly with height, and increasing as Ri decreases for a wind that rotates with height maintaining its magnitude. In these two cases the surface drag is predicted to be aligned with the surface wind. When one of the wind components varies linearly with height and the other is constant, the surface drag is misaligned with the surface wind, especially for relatively small Ri. All these results are shown to be in fairly good agreement with numerical simulations of mesoscale nonhydrostatic models, for high and even moderate values of Ri.

Comparative biology of mammalian telomeres: hypotheses on ancestral states and the roles of telomeres in longevity determination

Progressive telomere shortening from cell division (replicative aging) provides a barrier for human tumor progression. This program is not conserved in laboratory mice, which have longer telomeres and constitutive telomerase. Wild species that do ⁄ do not use replicative aging have been reported, but the evolution of different phenotypes and a conceptual framework for understanding their uses of telomeres is lacking. We examined telomeres ⁄ telomerase in cultured cells from > 60 mammalian species to place different uses of telomeres in a broad mammalian context. Phylogeny-based statistical analysis reconstructed ancestral states. Our analysis suggested that the ancestral mammalian phenotype included short telomeres (< 20 kb, as we now see in humans) and repressed telomerase. We argue that the repressed telomerase was a response to a higher mutation load brought on by the evolution of homeothermy. With telomerase repressed, we then see the evolution of replicative aging. Telomere length inversely correlated with lifespan, while telomerase expression co-evolved with body size. Multiple independent times smaller, shorter-lived species changed to having longer telomeres and expressing telomerase. Trade-offs involving reducing the energetic ⁄ cellular costs of specific oxidative protection mechanisms (needed to protect < 20 kb telomeres in the absence oftelomerase) could explain this abandonment of replicative aging. These observations provide a conceptual framework for understanding different uses of telomeres in mammals, support a role for human-like telomeres in allowing longer lifespans to evolve, demonstrate the need to include telomere length in the analysis of comparative studies of oxidative protection in the biology of aging, and identify which mammals can be used as appropriate model organisms for the study of the role of telomeres in human cancer and aging. Key words: evolution of telomeres; immortalization; telomerase; replicative aging; senescence.