Photochemical cross-linking of plastically compressed collagen gel produces an optimal scaffold for corneal tissue engineering

The experiments were designed to use photochemically cross-linked plastically compressed collagen (PCPCC) gel to support corneal epithelial cells. A plastically compressed collagen (PCC) scaffold was photo cross-linked by UVA in the presence of riboflavin to form a biomaterial with optimal mechanical properties. The breaking force, rheology, surgical suture strength, transparency, ultrastructure, and cell-based biocompatibility were compared between PCPCC and PCC gels. The breaking force increased proportionally with an increased concentration of riboflavin. The stress required to reach breaking point of the PCPCC scaffolds was over two times higher compared to the stress necessary to break PCC scaffolds in the presence of 0.1% riboflavin. Rheology results indicated that the structural properties of PCC remain unaltered after UVA cross-linking. The PCC gels were more easily broken than PCPCC gels when sutured on to bovine corneas. The optical density values of PCPCC and PCC showed no significant differences (p > 0.05). SEM analyses showed that the collagen fibres within the PCPCC gels were similar in morphology to PCC gels. No difference in cell-based biocompatibility was seen between the PCPCC and PCC scaffolds in terms of their ability to support the ex vivo expansion of corneal epithelial cells or their subsequent differentiation evidenced by similar levels of cytokeratin 14. In conclusion, PCPCC scaffold is an optimal biomaterial for use in therapeutic tissue engineering of the cornea.

A tutorial on pseudospectral methods for computational optimal control

[English] This paper is a tutorial introduction to pseudospectral optimal control. With pseudospectral methods, a function is approximated as a linear combination of smooth basis functions, which are often chosen to be Legendre or Chebyshev polynomials. Collocation of the differential-algebraic equations is performed at orthogonal collocation points, which are selected to yield interpolation of high accuracy. Pseudospectral methods directly discretize the original optimal control problem to recast it into a nonlinear programming format. A numerical optimizer is then employed to find approximate local optimal solutions. The paper also briefly describes the functionality and implementation of PSOPT, an open source software package written in C++ that employs pseudospectral discretization methods to solve multi-phase optimal control problems. The software implements the Legendre and Chebyshev pseudospectral methods, and it has useful features such as automatic differentiation, sparsity detection, and automatic scaling. The use of pseudospectral methods is illustrated in two problems taken from the literature on computational optimal control. [Portuguese] Este artigo e um tutorial introdutorio sobre controle otimo pseudo-espectral. Em metodos pseudo-espectrais, uma funcao e aproximada como uma combinacao linear de funcoes de base suaves, tipicamente escolhidas como polinomios de Legendre ou Chebyshev. A colocacao de equacoes algebrico-diferenciais e realizada em pontos de colocacao ortogonal, que sao selecionados de modo a minimizar o erro de interpolacao. Metodos pseudoespectrais discretizam o problema de controle otimo original de modo a converte-lo em um problema de programa cao nao-linear. Um otimizador numerico e entao empregado para obter solucoes localmente otimas. Este artigo tambem descreve sucintamente a funcionalidade e a implementacao de um pacote computacional de codigo aberto escrito em C++ chamado PSOPT. Tal pacote emprega metodos de discretizacao pseudo-spectrais para resolver problemas de controle otimo com multiplas fase. O PSOPT permite a utilizacao de metodos de Legendre ou Chebyshev, e possui caractersticas uteis tais como diferenciacao automatica, deteccao de esparsidade e escalonamento automatico. O uso de metodos pseudo-espectrais e ilustrado em dois problemas retirados da literatura de controle otimo computacional.

Retrieval of mesospheric electron densities using an optimal estimation inverse method

We present a new method to determine mesospheric electron densities from partially reflected medium frequency radar pulses. The technique uses an optimal estimation inverse method and retrieves both an electron density profile and a gradient electron density profile. As well as accounting for the absorption of the two magnetoionic modes formed by ionospheric birefringence of each radar pulse, the forward model of the retrieval parameterises possible Fresnel scatter of each mode by fine electronic structure, phase changes of each mode due to Faraday rotation and the dependence of the amplitudes of the backscattered modes upon pulse width. Validation results indicate that known profiles can be retrieved and that χ2 tests upon retrieval parameters satisfy validity criteria. Application to measurements shows that retrieved electron density profiles are consistent with accepted ideas about seasonal variability of electron densities and their dependence upon nitric oxide production and transport.

Oxidation of low density lipoprotein by iron at lysosomal pH: implications for atherosclerosis

Low density lipoprotein (LDL) has recently been shown to be oxidised by iron within the lysosomes of macrophages and this is a novel potential mechanism for LDL oxidation in atherosclerosis. Our aim was to characterise the chemical and physical changes induced in LDL by iron at lysosomal pH and to investigate the effects of iron chelators and α-tocopherol on this process. LDL was oxidised by iron at pH 4.5 and 37°C and its oxidation monitored by spectrophotometry and HPLC. LDL was oxidised effectively by FeSO4 (5-50 µM) and became highly aggregated at pH 4.5, but not at pH 7.4. Cholesteryl esters decreased and after a pronounced lag 7-ketocholesterol increased greatly. Total hydroperoxides (measured by tri-iodide assay) increased up to 24 h and then decreased only slowly. The lipid composition after 12 h at pH 4.5 and 37°C was similar to that of LDL oxidised by copper at pH 7.4 and 4°C, i.e. rich in hydroperoxides but low in oxysterols. Previously oxidised LDL aggregated rapidly and spontaneously at pH 4.5, but not at pH 7.4. Ferrous was much more effective than ferric iron at oxidising LDL when added after the oxidation was already underway. The iron chelators diethylenetriaminepentaacetic acid and, to a lesser extent, desferrioxamine inhibited LDL oxidation when added during its initial stages, but were unable to prevent LDL aggregating after it had been partially oxidised. Surprisingly, desferrioxamine increased the rate of LDL modification when added late in the oxidation process. α-Tocopherol enrichment of LDL initially increased the oxidation of LDL, but inhibited it later. The presence of oxidised and highly aggregated lipid within lysosomes has the potential to perturb the function of these organelles and to promote atherosclerosis.

Effects of oxidised low density lipoprotein on dendritic cells: a possible immunoregulatory component of the atherogenic micro-environment?

Objective: The objective of this study was to explore the relationship between low density lipoprotein (LDL) and dendritic cell (DC) activation, based upon the hypothesis that reactive oxygen species (ROS)-mediated modification of proteins that may be present in local DC microenvironments could be important as mediators of this activation. Although LDL are known to be oxidised in vivo, and taken up by macrophages during atherogenesis; their effect on DC has not been explored previously. Methods: Human DCs were prepared from peripheral blood monocytes using GM-CSF and IL-4. Plasma LDLs were isolated by sequential gradient centrifugation, oxidised in CuSO4, and oxidation arrested to yield mild, moderate and highly oxidised LDL forms. DCs exposed to these LDLs were investigated using combined phenotypic, functional (autologous T cell activation), morphological and viability assays. Results: Highly-oxidised LDL increased DC HLA-DR, CD40 and CD86 expression, corroborated by increased DC-induced T cell proliferation. Both native and oxidised LDL induced prominent DC clustering. However, high concentrations of highly-oxidised LDL inhibited DC function, due to increased DC apoptosis. Conclusions: This study supports the hypothesis that oxidised LDL are capable of triggering the transition from sentinel to messenger DC. Furthermore, the DC clustering–activation–apoptosis sequence in the presence of different LDL forms is consistent with a regulatory DC role in immunopathogenesis of atheroma. A sequence of initial accumulation of DC, increasing LDL oxidation, and DC-induced T cell activation, may explain why local breach of tolerance can occur. Above a threshold level, however, supervening DC apoptosis limits this, contributing instead to the central plaque core.

Ethnic differences in beta-cell function, dietary intake and expression of the metabolic syndrome among UK adults of south Asian, black African-Caribbean and white-European origin at high risk of metabolic syndrome

A cross-sectional analysis of ethnic differences in dietary intake, insulin sensitivity and beta-cell function, using the intravenous glucose tolerance test (IVGTT), was conducted on 497 healthy adult participants of the ‘Reading, Imperial, Surrey, Cambridge, and Kings’ (RISCK) study. Insulin sensitivity (Si) was significantly lower in African-Caribbean (AC) and South Asian (SA) participants [IVGTT-Si; AC: 2.13 vs SA: 2.25 vs white-European (WE): 2.84 (×10−4 mL µU min)2, p < 0.001]. AC participants had a higher prevalence of anti-hypertensive therapy (AC: 19.7% vs SA: 7.5%), the most cardioprotective lipid profile [total:high-density lipoprotein (HDL); AC: 3.52 vs SA: 4.08 vs WE: 3.83, p = 0.03] and more pronounced hyperinsulinaemia [IVGTT–acute insulin response (AIR)] [AC: 575 vs SA: 428 vs WE: 344 mL/µU/min)2, p = 0.002], specifically in female participants. Intake of saturated fat and carbohydrate was lower and higher in AC (10.9% and 50.4%) and SA (11.1% and 52.3%), respectively, compared to WE (13.6% and 43.8%, p < 0.001). Insulin resistance in ACs is characterised by ‘normal’ lipid profiles but high rates of hypertension and pronounced hyperinsulinaemia.

Estimation of optimal gravity wave parameters for climate models using data assimilation

There is a current need to constrain the parameters of gravity wave drag (GWD) schemes in climate models using observational information instead of tuning them subjectively. In this work, an inverse technique is developed using data assimilation principles to estimate gravity wave parameters. Because mostGWDschemes assume instantaneous vertical propagation of gravity waves within a column, observations in a single column can be used to formulate a one-dimensional assimilation problem to estimate the unknown parameters. We define a cost function that measures the differences between the unresolved drag inferred from observations (referred to here as the ‘observed’ GWD) and the GWD calculated with a parametrisation scheme. The geometry of the cost function presents some difficulties, including multiple minima and ill-conditioning because of the non-independence of the gravity wave parameters. To overcome these difficulties we propose a genetic algorithm to minimize the cost function, which provides a robust parameter estimation over a broad range of prescribed ‘true’ parameters. When real experiments using an independent estimate of the ‘observed’ GWD are performed, physically unrealistic values of the parameters can result due to the non-independence of the parameters. However, by constraining one of the parameters to lie within a physically realistic range, this degeneracy is broken and the other parameters are also found to lie within physically realistic ranges. This argues for the essential physical self-consistency of the gravity wave scheme. A much better fit to the observed GWD at high latitudes is obtained when the parameters are allowed to vary with latitude. However, a close fit can be obtained either in the upper or the lower part of the profiles, but not in both at the same time. This result is a consequence of assuming an isotropic launch spectrum. The changes of sign in theGWDfound in the tropical lower stratosphere, which are associated with part of the quasi-biennial oscillation forcing, cannot be captured by the parametrisation with optimal parameters.

Optimal frame length for keeping normalized goodput with lowest requirement on BER

This paper presents an adaptive frame length mechanism based on a cross-layer analysis of intrinsic relations between the MAC frame length, bit error rate (BER) of the wireless link and normalized goodput. The proposed mechanism selects the optimal frame length that keeps the service normalized goodput at required levels while satisfying the lowest requirement on the BER, thus increasing the transmission reliability. Numerical results are provided and show that an optimal frame length satisfying the lowest BER requirement does indeed exist. The performance of BER requirement as a function of the MAC frame length is evaluated and compared for transmission scenarios with and without automatic repeat request (ARQ). Furthermore, issues related to the MAC overhead length are also discussed to illuminate the functionality and performance of the proposed mechanism.

Surface water temperature observations of large lakes by optimal estimation

Optimal estimation (OE) and probabilistic cloud screening were developed to provide lake surface water temperature (LSWT) estimates from the series of (advanced) along-track scanning radiometers (ATSRs). Variations in physical properties such as elevation, salinity, and atmospheric conditions are accounted for through the forward modelling of observed radiances. Therefore, the OE retrieval scheme developed is generic (i.e., applicable to all lakes). LSWTs were obtained for 258 of Earth's largest lakes from ATSR-2 and AATSR imagery from 1995 to 2009. Comparison to in situ observations from several lakes yields satellite in situ differences of −0.2 ± 0.7 K for daytime and −0.1 ± 0.5 K for nighttime observations (mean ± standard deviation). This compares with −0.05 ± 0.8 K for daytime and −0.1 ± 0.9 K for nighttime observations for previous methods based on operational sea surface temperature algorithms. The new approach also increases coverage (reducing misclassification of clear sky as cloud) and exhibits greater consistency between retrievals using different channel–view combinations. Empirical orthogonal function (EOF) techniques were applied to the LSWT retrievals (which contain gaps due to cloud cover) to reconstruct spatially and temporally complete time series of LSWT. The new LSWT observations and the EOF-based reconstructions offer benefits to numerical weather prediction, lake model validation, and improve our knowledge of the climatology of lakes globally. Both observations and reconstructions are publically available from http://hdl.handle.net/10283/88.

Density functional theory study of rutile VO2 surfaces

We present the results of a density functional theory (DFT) investigation of the surfaces of rutile-like vanadium dioxide, VO2(R). We calculate the surface energies of low Miller index planes, and find that the most stable surface orientation is the (110). The equilibrium morphology of a VO2(R) particle has an acicular shape, laterally confined by (110) planes and topped by (011) planes. The redox properties of the (110) surface are investigated by calculating the relative surface free energies of the non-stoichiometric compositions as a function of oxygen chemical potential. It is found that the VO2(110) surface is oxidized with respect to the stoichiometric composition, not only at ambient conditions but also at the more reducing conditions under which bulk VO2 is stable in comparison with bulk V2O5. The adsorbed oxygen forms surface vanadyl species much more favorably than surface peroxo species.

On the application of optimal wavelet filter banks for ECG signal classification

This paper discusses ECG signal classification after parametrizing the ECG waveforms in the wavelet domain. Signal decomposition using perfect reconstruction quadrature mirror filter banks can provide a very parsimonious representation of ECG signals. In the current work, the filter parameters are adjusted by a numerical optimization algorithm in order to minimize a cost function associated to the filter cut-off sharpness. The goal consists of achieving a better compromise between frequency selectivity and time resolution at each decomposition level than standard orthogonal filter banks such as those of the Daubechies and Coiflet families. Our aim is to optimally decompose the signals in the wavelet domain so that they can be subsequently used as inputs for training to a neural network classifier.

Solar causes of the long-term increase in geomagnetic activity

We analyze the causes of the century-long increase in geomagnetic activity, quantified by annual means of the aa index, using observations of interplanetary space, galactic cosmic rays, the ionosphere, and the auroral electrojet, made during the last three solar cycles. The effects of changes in ionospheric conductivity, the Earth's dipole tilt, and magnetic moment are shown to be small; only changes in near-Earth interplanetary space make a significant contribution to the long-term increase in activity. We study the effects of the interplanetary medium by applying dimensional analysis to generate the optimum solar wind-magnetosphere energy coupling function, having an unprecedentedly high correlation coefficient of 0.97. Analysis of the terms of the coupling function shows that the largest contributions to the drift in activity over solar cycles 20-22 originate from rises in the average interplanetary magnetic field (IMF) strength, solar wind concentration, and speed; average IMF orientation has grown somewhat less propitious for causing geomagnetic activity. The combination of these factors explains almost all of the 39% rise in aa observed over the last three solar cycles. Whereas the IMF strength varies approximately in phase with sunspot numbers, neither its orientation nor the solar wind density shows any coherent solar cycle variation. The solar wind speed peaks strongly in the declining phase of even-numbered cycles and can be identified as the chief cause of the phase shift between the sunspot numbers and the aa index. The rise in the IMF magnitude, the largest single contributor to the drift in geomagnetic activity, is shown to be caused by a rise in the solar coronal magnetic field, consistent with a rise in the coronal source field, modeled from photospheric observations, and an observed decay in cosmic ray fluxes.

Radar observations of non-thermal plasmas at different aspect angles

Data are presented from the EISCAT CP-3-E experiment which show the presence of non-thermal plasma over a range of latitudes. The O+ ion-velocity distribution function is almost toroidal when the electric field reaches values of 125 mV m−1. The ion temperature derived from such data assuming a Maxwellian distribution function will overestimate the true ion temperature when the observing angle is large with respect to the magnetic field, and underestimate the temperature when the aspect angle is small. When the expressions for the distribution function are extended to include mixed ion composition, an improvement is sometimes found in fitting the observed data, and estimates of the composition can be made. Such an analysis suggests that N2+ can occasionally form a significant part of the total ion density in a narrow height region centred at 275 km.

Synaptic vesicle glycoprotein 2A modulates vesicular release and calcium channel function at peripheral sympathetic synapses

Synaptic vesicle glycoprotein (SV)2A is a transmembrane protein found in secretory vesicles and is critical for Ca2+-dependent exocytosis in central neurons, although its mechanism of action remains uncertain. Previous studies have proposed, variously, a role of SV2 in the maintenance and formation of the readily releasable pool (RRP) or in the regulation of Ca2+ responsiveness of primed vesicles. Such previous studies have typically used genetic approaches to ablate SV2 levels; here, we used a strategy involving small interference RNA (siRNA) injection to knockdown solely presynaptic SV2A levels in rat superior cervical ganglion (SCG) neuron synapses. Moreover, we investigated the effects of SV2A knockdown on voltage-dependent Ca2+ channel (VDCC) function in SCG neurons. Thus, we extended the studies of SV2A mechanisms by investigating the effects on vesicular transmitter release and VDCC function in peripheral sympathetic neurons. We first demonstrated an siRNA-mediated SV2A knockdown. We showed that this SV2A knockdown markedly affected presynaptic function, causing an attenuated RRP size, increased paired-pulse depression and delayed RRP recovery after stimulus-dependent depletion. We further demonstrated that the SV2A–siRNA-mediated effects on vesicular release were accompanied by a reduction in VDCC current density in isolated SCG neurons. Together, our data showed that SV2A is required for correct transmitter release at sympathetic neurons. Mechanistically, we demonstrated that presynaptic SV2A: (i) acted to direct normal synaptic transmission by maintaining RRP size, (ii) had a facilitatory role in recovery from synaptic depression, and that (iii) SV2A deficits were associated with aberrant Ca2+ current density, which may contribute to the secretory phenotype in sympathetic peripheral neurons.

Nonlinear identification using orthogonal forward regression with nested optimal regularization

An efficient data based-modeling algorithm for nonlinear system identification is introduced for radial basis function (RBF) neural networks with the aim of maximizing generalization capability based on the concept of leave-one-out (LOO) cross validation. Each of the RBF kernels has its own kernel width parameter and the basic idea is to optimize the multiple pairs of regularization parameters and kernel widths, each of which is associated with a kernel, one at a time within the orthogonal forward regression (OFR) procedure. Thus, each OFR step consists of one model term selection based on the LOO mean square error (LOOMSE), followed by the optimization of the associated kernel width and regularization parameter, also based on the LOOMSE. Since like our previous state-of-the-art local regularization assisted orthogonal least squares (LROLS) algorithm, the same LOOMSE is adopted for model selection, our proposed new OFR algorithm is also capable of producing a very sparse RBF model with excellent generalization performance. Unlike our previous LROLS algorithm which requires an additional iterative loop to optimize the regularization parameters as well as an additional procedure to optimize the kernel width, the proposed new OFR algorithm optimizes both the kernel widths and regularization parameters within the single OFR procedure, and consequently the required computational complexity is dramatically reduced. Nonlinear system identification examples are included to demonstrate the effectiveness of this new approach in comparison to the well-known approaches of support vector machine and least absolute shrinkage and selection operator as well as the LROLS algorithm.