Function inferences from a molecular structural model of bacterial ParE toxin

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

The massless two-loop two-point function and zeta functions in counterterms of Feynman diagrams

The main part of this thesis describes a method of calculating the massless two-loop two-point function which allows expanding the integral up to an arbitrary order in the dimensional regularization parameter epsilon by rewriting it as a double Mellin-Barnes integral. Closing the contour and collecting the residues then transforms this integral into a form that enables us to utilize S. Weinzierl's computer library nestedsums. We could show that multiple zeta values and rational numbers are sufficient for expanding the massless two-loop two-point function to all orders in epsilon. We then use the Hopf algebra of Feynman diagrams and its antipode, to investigate the appearance of Riemann's zeta function in counterterms of Feynman diagrams in massless Yukawa theory and massless QED. The class of Feynman diagrams we consider consists of graphs built from primitive one-loop diagrams and the non-planar vertex correction, where the vertex corrections only depend on one external momentum. We showed the absence of powers of pi in the counterterms of the non-planar vertex correction and diagrams built by shuffling it with the one-loop vertex correction. We also found the invariance of some coefficients of zeta functions under a change of momentum flow through these vertex corrections.

An integrative transcranial magnetic stimulation mapping technique using non-linear curve fitting

The aim of this study is to develop a new simple method for analyzing one-dimensional transcranial magnetic stimulation (TMS) mapping studies in humans. Motor evoked potentials (MEP) were recorded from the abductor pollicis brevis (APB) muscle during stimulation at nine different positions on the scalp along a line passing through the APB hot spot and the vertex. Non-linear curve fitting according to the Levenberg-Marquardt algorithm was performed on the averaged amplitude values obtained at all points to find the best-fitting symmetrical and asymmetrical peak functions. Several peak functions could be fitted to the experimental data. Across all subjects, a symmetric, bell-shaped curve, the complementary error function (erfc) gave the best results. This function is characterized by three parameters giving its amplitude, position, and width. None of the mathematical functions tested with less or more than three parameters fitted better. The amplitude and position parameters of the erfc were highly correlated with the amplitude at the hot spot and with the location of the center of gravity of the TMS curve. In conclusion, non-linear curve fitting is an accurate method for the mathematical characterization of one-dimensional TMS curves. This is the first method that provides information on amplitude, position and width simultaneously.

Fitting 3D morphable models using implicit representations

We consider the problem of approximating the 3D scan of a real object through an affine combination of examples. Common approaches depend either on the explicit estimation of point-to-point correspondences or on 2-dimensional projections of the target mesh; both present drawbacks. We follow an approach similar to [IF03] by representing the target via an implicit function, whose values at the vertices of the approximation are used to define a robust cost function. The problem is approached in two steps, by approximating first a coarse implicit representation of the whole target, and then finer, local ones; the local approximations are then merged together with a Poisson-based method. We report the results of applying our method on a subset of 3D scans from the Face Recognition Grand Challenge v.1.0.

Effects of coronary sinus occlusion on myocardial ischaemia in humans: role of coronary collateral function

OBJECTIVE This study tested the hypotheses that intermittent coronary sinus occlusion (iCSO) reduces myocardial ischaemia, and that the amount of ischaemia reduction is related to coronary collateral function. DESIGN Prospective case-control study with intraindividual comparison of myocardial ischaemia during two 2-min coronary artery balloon occlusions with and without simultaneous iCSO by a balloon-tipped catheter. SETTING University Hospital. PATIENTS 35 patients with chronic stable coronary artery disease. INTERVENTION 2-min iCSO. MAIN OUTCOME MEASURES Myocardial ischaemia as assessed by intracoronary (i.c.) ECG ST shift at 2 min of coronary artery balloon occlusion. Collateral flow index (CFI) without iCSO, that is, the ratio between mean distal coronary occlusive (Poccl) and mean aortic pressure (Pao) both minus central venous pressure. RESULTS I.c. ECG ST segment shift (elevation in all) at the end of the procedure with iCSO versus without iCSO was 1.33±1.25 mV versus 1.85±1.45 mV, p<0.0001. Regression analysis showed that the degree of i.c. ECG ST shift reduction during iCSO was related to CFI, best fitting a Lorentzian function (r(2)=0.61). Ischaemia reduction with iCSO was greatest at a CFI of 0.05-0.20, whereas in the low and high CFI range the effect of iCSO was absent. CONCLUSIONS ICSO reduces myocardial ischaemia in patients with chronic coronary artery disease. Ischaemia reduction by iCSO depends on coronary collateral function. A minimal degree of collateral function is necessary to render iCSO effective. ICSO cannot manifest an effect when collateral function prevents ischaemia in the first place.

A fitting algorithm for random modeling the PLC channel

The characteristics of the power-line communication (PLC) channel are difficult to model due to the heterogeneity of the networks and the lack of common wiring practices. To obtain the full variability of the PLC channel, random channel generators are of great importance for the design and testing of communication algorithms. In this respect, we propose a random channel generator that is based on the top-down approach. Basically, we describe the multipath propagation and the coupling effects with an analytical model. We introduce the variability into a restricted set of parameters and, finally, we fit the model to a set of measured channels. The proposed model enables a closed-form description of both the mean path-loss profile and the statistical correlation function of the channel frequency response. As an example of application, we apply the procedure to a set of in-home measured channels in the band 2-100 MHz whose statistics are available in the literature. The measured channels are divided into nine classes according to their channel capacity. We provide the parameters for the random generation of channels for all nine classes, and we show that the results are consistent with the experimental ones. Finally, we merge the classes to capture the entire heterogeneity of in-home PLC channels. In detail, we introduce the class occurrence probability, and we present a random channel generator that targets the ensemble of all nine classes. The statistics of the composite set of channels are also studied, and they are compared to the results of experimental measurement campaigns in the literature.

Uncertainty in ellipse fitting using a flatbed scanner: development and experimental verification

In the field of dimensional metrology, the use of optical measuring machines requires the handling of a large number of measurement points, or scanning points, taken from the image of the measurand. The presence of correlation between these measurement points has a significant influence on the uncertainty of the result. The aim of this work is the development of an estimation procedure for the uncertainty of measurement in a geometrically elliptical shape, taking into account the correlation between the scanning points. These points are obtained from an image produced using a commercial flat bed scanner. The characteristic parameters of the ellipse (coordinates of the center, semi-axes and the angle of the semi-major axis with regard to the horizontal) are determined using a least squares fit and orthogonal distance regression. The uncertainty is estimated using the information from the auto-correlation function of the residuals and is propagated through the fitting algorithm according to the rules described in Evaluation of Measurement Data—Supplement 2 to the ‘Guide to the Expression of Uncertainty in Measurement’—Extension to any number of output quantities. By introducing the concept of cut-off length, it can be observed how it is possible to take into account the presence of the correlation in the estimation of uncertainty in a very simple way while avoiding underestimation.

Effect of curve fitting parameters on quantitative analysis of Fe0.94O and Fe2O3 using XPS

The XPS peaks of Fe 3p for Fe2+ and Fe3+ in FeO and Fe2O3, respectively, have been measured and the effects of curve fitting parameters on interpretation of the data have been analysed. Firstly, the peak fit parameters, i.e. (1) the number of peaks to be deconvoluted, (2) the range of the peak for back ground subtraction, (3) straight line (Li) or the Shirley (Sh) background subtraction method, (4) GL ratio (the ratio of Gaussian and Lorentzian contribution to the peak shape) and (5) asymmetry factor (AS), are manually selected. Secondly, the standard peak fit parameters were systematically investigated. The peak shape was fitted to a Voigt function by changing the peak position, the peak height and the full width at half maximum (FWHM) to minimize the chi(2). The recommended peak positions and peak parameters for Fe2+ and Fe3+ in iron oxides have been determined. (c) 2006 Elsevier B.V. All rights reserved.

Extracellular loops 1 and 3 and their associated transmembrane regions of the calcitonin receptor-like receptor are needed for CGRP receptor function

The first and third extracellular loops (ECL) of G protein-coupled receptors (GPCRs) have been implicated in ligand binding and receptor function. This study describes the results of an alanine/leucine scan of ECLs 1 and 3 and loop-associated transmembrane (TM) domains of the secretin-like GPCR calcitonin receptor-like receptor which associates with receptor activity modifying protein 1 to form the CGRP receptor. Leu195Ala, Val198Ala and Ala199Leu at the top of TM2 all reduced aCGRP-mediated cAMP production and internalization; Leu195Ala and Ala199Leu also reduced aCGRP binding. These residues form a hydrophobic cluster within an area defined as the "minor groove" of rhodopsin-like GPCRs. Within ECL1, Ala203Leu and Ala206Leu influenced the ability of aCGRP to stimulate adenylate cyclase. In TM3, His219Ala, Leu220Ala and Leu222Ala have influences on aCGRP binding and cAMP production; they are likely to indirectly influence the binding site for aCGRP as well as having an involvement in signal transduction. On the exofacial surfaces of TMs 6 and 7, a number of residues were identified that reduced cell surface receptor expression, most noticeably Leu351Ala and Glu357Ala in TM6. The residues may contribute to the RAMP1 binding interface. Ile360Ala impaired aCGRP-mediated cAMP production. Ile360 is predicted to be located close to ECL2 and may facilitate receptor activation. Identification of several crucial functional loci gives further insight into the activation mechanism of this complex receptor system and may aid rational drug design.

The slope of the psychometric function and non-stationarity of thresholds in spatiotemporal contrast vision

The slope of the two-interval, forced-choice psychometric function (e.g. the Weibull parameter, ß) provides valuable information about the relationship between contrast sensitivity and signal strength. However, little is known about how or whether ß varies with stimulus parameters such as spatiotemporal frequency and stimulus size and shape. A second unresolved issue concerns the best way to estimate the slope of the psychometric function. For example, if an observer is non-stationary (e.g. their threshold drifts between experimental sessions), ß will be underestimated if curve fitting is performed after collapsing the data across experimental sessions. We measured psychometric functions for 2 experienced observers for 14 different spatiotemporal configurations of pulsed or flickering grating patches and bars on each of 8 days. We found ß ˜ 3 to be fairly constant across almost all conditions, consistent with a fixed nonlinear contrast transducer and/or a constant level of intrinsic stimulus uncertainty (e.g. a square law transducer and a low level of intrinsic uncertainty). Our analysis showed that estimating a single ß from results averaged over several experimental sessions was slightly more accurate than averaging multiple estimates from several experimental sessions. However, the small levels of non-stationarity (SD ˜ 0.8 dB) meant that the difference between the estimates was, in practice, negligible.

Determination of excitation profile and dielectric function spatial nonuniformity in porous silicon by using WKB approach

We develop an analytical model based on the WKB approach to evaluate the experimental results of the femtosecond pump-probe measurements of the transmittance and reflectance obtained on thin membranes of porous silicon. The model allows us to retrieve a pump-induced nonuniform complex dielectric function change along the membrane depth. We show that the model fitting to the experimental data requires a minimal number of fitting parameters while still complying with the restriction imposed by the Kramers-Kronig relation. The developed model has a broad range of applications for experimental data analysis and practical implementation in the design of devices involving a spatially nonuniform dielectric function, such as in biosensing, wave-guiding, solar energy harvesting, photonics and electro-optical devices.

Gut hormone GPCRs: structure, function, drug discovery

Crystallization and determination of the high resolution three-dimensional structure of the β2-adrenergic receptor in 2007 was followed by structure elucidation of a number of other receptors, including those for neurotensin and glucagon. These major advances foster the understanding of structure-activity relationship of these receptors and structure-based rational design of new ligands having more predictable activity. At present, structure determination of gut hormone receptors in complex with their ligands (natural, synthetic) and interacting signalling proteins, for example, G-proteins, arrestins, represents a challenge which promises to revolutionize gut hormone endocrinonology.

A signaling visualization toolkit to support rational design of combination therapies and biomarker discovery: SiViT

Targeted cancer therapy aims to disrupt aberrant cellular signalling pathways. Biomarkers are surrogates of pathway state, but there is limited success in translating candidate biomarkers to clinical practice due to the intrinsic complexity of pathway networks. Systems biology approaches afford better understanding of complex, dynamical interactions in signalling pathways targeted by anticancer drugs. However, adoption of dynamical modelling by clinicians and biologists is impeded by model inaccessibility. Drawing on computer games technology, we present a novel visualisation toolkit, SiViT, that converts systems biology models of cancer cell signalling into interactive simulations that can be used without specialist computational expertise. SiViT allows clinicians and biologists to directly introduce for example loss of function mutations and specific inhibitors. SiViT animates the effects of these introductions on pathway dynamics, suggesting further experiments and assessing candidate biomarker effectiveness. In a systems biology model of Her2 signalling we experimentally validated predictions using SiViT, revealing the dynamics of biomarkers of drug resistance and highlighting the role of pathway crosstalk. No model is ever complete: the iteration of real data and simulation facilitates continued evolution of more accurate, useful models. SiViT will make accessible libraries of models to support preclinical research, combinatorial strategy design and biomarker discovery.

Sensitivity analysis in optimized parametric curve fitting

Purpose – Curve fitting from unordered noisy point samples is needed for surface reconstruction in many applications -- In the literature, several approaches have been proposed to solve this problem -- However, previous works lack formal characterization of the curve fitting problem and assessment on the effect of several parameters (i.e. scalars that remain constant in the optimization problem), such as control points number (m), curve degree (b), knot vector composition (U), norm degree (k), and point sample size (r) on the optimized curve reconstruction measured by a penalty function (f) -- The paper aims to discuss these issues -- Design/methodology/approach - A numerical sensitivity analysis of the effect of m, b, k and r on f and a characterization of the fitting procedure from the mathematical viewpoint are performed -- Also, the spectral (frequency) analysis of the derivative of the angle of the fitted curve with respect to u as a means to detect spurious curls and peaks is explored -- Findings - It is more effective to find optimum values for m than k or b in order to obtain good results because the topological faithfulness of the resulting curve strongly depends on m -- Furthermore, when an exaggerate number of control points is used the resulting curve presents spurious curls and peaks -- The authors were able to detect the presence of such spurious features with spectral analysis -- Also, the authors found that the method for curve fitting is robust to significant decimation of the point sample -- Research limitations/implications - The authors have addressed important voids of previous works in this field -- The authors determined, among the curve fitting parameters m, b and k, which of them influenced the most the results and how -- Also, the authors performed a characterization of the curve fitting problem from the optimization perspective -- And finally, the authors devised a method to detect spurious features in the fitting curve -- Practical implications – This paper provides a methodology to select the important tuning parameters in a formal manner -- Originality/value - Up to the best of the knowledge, no previous work has been conducted in the formal mathematical evaluation of the sensitivity of the goodness of the curve fit with respect to different possible tuning parameters (curve degree, number of control points, norm degree, etc.)