Simulating maize phenology as a function of air temperature with a linear and a nonlinear model

The objective of this study was to adapt a nonlinear model (Wang and Engel - WE) for simulating the phenology of maize (Zea mays L.), and to evaluate this model and a linear one (thermal time), in order to predict developmental stages of a field-grown maize variety. A field experiment, during 2005/2006 and 2006/2007 was conducted in Santa Maria, RS, Brazil, in two growing seasons, with seven sowing dates each. Dates of emergence, silking, and physiological maturity of the maize variety BRS Missões were recorded in six replications in each sowing date. Data collected in 2005/2006 growing season were used to estimate the coefficients of the two models, and data collected in the 2006/2007 growing season were used as independent data set for model evaluations. The nonlinear WE model accurately predicted the date of silking and physiological maturity, and had a lower root mean square error (RMSE) than the linear (thermal time) model. The overall RMSE for silking and physiological maturity was 2.7 and 4.8 days with WE model, and 5.6 and 8.3 days with thermal time model, respectively.

Improving node number simulation in soybean

The objective of this study was to improve the simulation of node number in soybean cultivars with determinate stem habits. A nonlinear model considering two approaches to input daily air temperature data (daily mean temperature and daily minimum/maximum air temperatures) was used. The node number on the main stem data of ten soybean cultivars was collected in a three-year field experiment (from 2004/2005 to 2006/2007) at Santa Maria, RS, Brazil. Node number was simulated using the Soydev model, which has a nonlinear temperature response function [f(T)]. The f(T) was calculated using two methods: using daily mean air temperature calculated as the arithmetic average among daily minimum and maximum air temperatures (Soydev tmean); and calculating an f(T) using minimum air temperature and other using maximum air temperature and then averaging the two f(T)s (Soydev tmm). Root mean square error (RMSE) and deviations (simulated minus observed) were used as statistics to evaluate the performance of the two versions of Soydev. Simulations of node number in soybean were better with the Soydev tmm version, with a 0.5 to 1.4 node RMSE. Node number can be simulated for several soybean cultivars using only one set of model coefficients, with a 0.8 to 2.4 node RMSE.

EADock : design of a new molecular docking algorithm and some of its applications

3 Summary 3. 1 English The pharmaceutical industry has been facing several challenges during the last years, and the optimization of their drug discovery pipeline is believed to be the only viable solution. High-throughput techniques do participate actively to this optimization, especially when complemented by computational approaches aiming at rationalizing the enormous amount of information that they can produce. In siiico techniques, such as virtual screening or rational drug design, are now routinely used to guide drug discovery. Both heavily rely on the prediction of the molecular interaction (docking) occurring between drug-like molecules and a therapeutically relevant target. Several softwares are available to this end, but despite the very promising picture drawn in most benchmarks, they still hold several hidden weaknesses. As pointed out in several recent reviews, the docking problem is far from being solved, and there is now a need for methods able to identify binding modes with a high accuracy, which is essential to reliably compute the binding free energy of the ligand. This quantity is directly linked to its affinity and can be related to its biological activity. Accurate docking algorithms are thus critical for both the discovery and the rational optimization of new drugs. In this thesis, a new docking software aiming at this goal is presented, EADock. It uses a hybrid evolutionary algorithm with two fitness functions, in combination with a sophisticated management of the diversity. EADock is interfaced with .the CHARMM package for energy calculations and coordinate handling. A validation was carried out on 37 crystallized protein-ligand complexes featuring 11 different proteins. The search space was defined as a sphere of 15 R around the center of mass of the ligand position in the crystal structure, and conversely to other benchmarks, our algorithms was fed with optimized ligand positions up to 10 A root mean square deviation 2MSD) from the crystal structure. This validation illustrates the efficiency of our sampling heuristic, as correct binding modes, defined by a RMSD to the crystal structure lower than 2 A, were identified and ranked first for 68% of the complexes. The success rate increases to 78% when considering the five best-ranked clusters, and 92% when all clusters present in the last generation are taken into account. Most failures in this benchmark could be explained by the presence of crystal contacts in the experimental structure. EADock has been used to understand molecular interactions involved in the regulation of the Na,K ATPase, and in the activation of the nuclear hormone peroxisome proliferatoractivated receptors a (PPARa). It also helped to understand the action of common pollutants (phthalates) on PPARy, and the impact of biotransformations of the anticancer drug Imatinib (Gleevec®) on its binding mode to the Bcr-Abl tyrosine kinase. Finally, a fragment-based rational drug design approach using EADock was developed, and led to the successful design of new peptidic ligands for the a5ß1 integrin, and for the human PPARa. In both cases, the designed peptides presented activities comparable to that of well-established ligands such as the anticancer drug Cilengitide and Wy14,643, respectively. 3.2 French Les récentes difficultés de l'industrie pharmaceutique ne semblent pouvoir se résoudre que par l'optimisation de leur processus de développement de médicaments. Cette dernière implique de plus en plus. de techniques dites "haut-débit", particulièrement efficaces lorsqu'elles sont couplées aux outils informatiques permettant de gérer la masse de données produite. Désormais, les approches in silico telles que le criblage virtuel ou la conception rationnelle de nouvelles molécules sont utilisées couramment. Toutes deux reposent sur la capacité à prédire les détails de l'interaction moléculaire entre une molécule ressemblant à un principe actif (PA) et une protéine cible ayant un intérêt thérapeutique. Les comparatifs de logiciels s'attaquant à cette prédiction sont flatteurs, mais plusieurs problèmes subsistent. La littérature récente tend à remettre en cause leur fiabilité, affirmant l'émergence .d'un besoin pour des approches plus précises du mode d'interaction. Cette précision est essentielle au calcul de l'énergie libre de liaison, qui est directement liée à l'affinité du PA potentiel pour la protéine cible, et indirectement liée à son activité biologique. Une prédiction précise est d'une importance toute particulière pour la découverte et l'optimisation de nouvelles molécules actives. Cette thèse présente un nouveau logiciel, EADock, mettant en avant une telle précision. Cet algorithme évolutionnaire hybride utilise deux pressions de sélections, combinées à une gestion de la diversité sophistiquée. EADock repose sur CHARMM pour les calculs d'énergie et la gestion des coordonnées atomiques. Sa validation a été effectuée sur 37 complexes protéine-ligand cristallisés, incluant 11 protéines différentes. L'espace de recherche a été étendu à une sphère de 151 de rayon autour du centre de masse du ligand cristallisé, et contrairement aux comparatifs habituels, l'algorithme est parti de solutions optimisées présentant un RMSD jusqu'à 10 R par rapport à la structure cristalline. Cette validation a permis de mettre en évidence l'efficacité de notre heuristique de recherche car des modes d'interactions présentant un RMSD inférieur à 2 R par rapport à la structure cristalline ont été classés premier pour 68% des complexes. Lorsque les cinq meilleures solutions sont prises en compte, le taux de succès grimpe à 78%, et 92% lorsque la totalité de la dernière génération est prise en compte. La plupart des erreurs de prédiction sont imputables à la présence de contacts cristallins. Depuis, EADock a été utilisé pour comprendre les mécanismes moléculaires impliqués dans la régulation de la Na,K ATPase et dans l'activation du peroxisome proliferatoractivated receptor a (PPARa). Il a également permis de décrire l'interaction de polluants couramment rencontrés sur PPARy, ainsi que l'influence de la métabolisation de l'Imatinib (PA anticancéreux) sur la fixation à la kinase Bcr-Abl. Une approche basée sur la prédiction des interactions de fragments moléculaires avec protéine cible est également proposée. Elle a permis la découverte de nouveaux ligands peptidiques de PPARa et de l'intégrine a5ß1. Dans les deux cas, l'activité de ces nouveaux peptides est comparable à celles de ligands bien établis, comme le Wy14,643 pour le premier, et le Cilengitide (PA anticancéreux) pour la seconde.

Simulation of soybean growth and yield under northeastern Amazon climatic conditions

The objective of this work was to parameterize, calibrate, and validate a new version of the soybean growth and yield model developed by Sinclair, under natural field conditions in northeastern Amazon. The meteorological data and the values of soybean growth and leaf area were obtained from an agrometeorological experiment carried out in Paragominas, PA, Brazil, from 2006 to 2009. The climatic conditions during the experiment were very distinct, with a slight reduction in rainfall in 2007, due to the El Niño phenomenon. There was a reduction in the leaf area index (LAI) and in biomass production during this year, which was reproduced by the model. The simulation of the LAI had root mean square error (RMSE) of 0.55 to 0.82 m² m-2, from 2006 to 2009. The simulation of soybean yield for independent data showed a RMSE of 198 kg ha-1, i.e., an overestimation of 3%. The model was calibrated and validated for Amazonian climatic conditions, and can contribute positively to the improvement of the simulations of the impacts of land use change in the Amazon region. The modified version of the Sinclair model is able to adequately simulate leaf area formation, total biomass, and soybean yield, under northeastern Amazon climatic conditions.

CO2-response function of radiation use efficiency in rice for climate change scenarios

The objective of this work was to evaluate a generalized response function to the atmospheric CO2 concentration [f(CO2)] by the radiation use efficiency (RUE) in rice. Experimental data on RUE at different CO2 concentrations were collected from rice trials performed in several locations around the world. RUE data were then normalized, so that all RUE at current CO2 concentration were equal to 1. The response function was obtained by fitting normalized RUE versus CO2 concentration to a Morgan-Mercer-Flodin (MMF) function, and by using Marquardt's method to estimate the model coefficients. Goodness of fit was measured by the standard deviation of the estimated coefficients, the coefficient of determination (R²), and the root mean square error (RMSE). The f(CO2) describes a nonlinear sigmoidal response of RUE in rice, in function of the atmospheric CO2 concentration, which has an ecophysiological background, and, therefore, renders a robust function that can be easily coupled to rice simulation models, besides covering the range of CO2 emissions for the next generation of climate scenarios for the 21st century.

JPEG standard uniform quantization error modeling with applications to sequential and progressive operation modes

In this paper we propose a method for computing JPEG quantization matrices for a given mean square error or PSNR. Then, we employ our method to compute JPEG standard progressive operation mode definition scripts using a quantization approach. Therefore, it is no longer necessary to use a trial and error procedure to obtain a desired PSNR and/or definition script, reducing cost. Firstly, we establish a relationship between a Laplacian source and its uniform quantization error. We apply this model to the coefficients obtained in the discrete cosine transform stage of the JPEG standard. Then, an image may be compressed using the JPEG standard under a global MSE (or PSNR) constraint and a set of local constraints determined by the JPEG standard and visual criteria. Secondly, we study the JPEG standard progressive operation mode from a quantization based approach. A relationship between the measured image quality at a given stage of the coding process and a quantization matrix is found. Thus, the definition script construction problem can be reduced to a quantization problem. Simulations show that our method generates better quantization matrices than the classical method based on scaling the JPEG default quantization matrix. The estimation of PSNR has usually an error smaller than 1 dB. This figure decreases for high PSNR values. Definition scripts may be generated avoiding an excessive number of stages and removing small stages that do not contribute during the decoding process with a noticeable image quality improvement.

Electromyographic, cerebral, and muscle hemodynamic responses during intermittent, isometric contractions of the biceps brachii at three submaximal intensities.

This study examined the electromyographic, cerebral and muscle hemodynamic responses during intermittent isometric contractions of biceps brachii at 20, 40, and 60% of maximal voluntary contraction (MVC). Eleven volunteers completed 2 min of intermittent isometric contractions (12/min) at an elbow angle of 90° interspersed with 3 min rest between intensities in systematic order. Surface electromyography (EMG) was recorded from the right biceps brachii and near infrared spectroscopy (NIRS) was used to simultaneously measure left prefrontal and right biceps brachii oxyhemoglobin (HbO2), deoxyhemoglobin (HHb), and total hemoglobin (Hbtot). Transcranial Doppler ultrasound was used to measure middle cerebral artery velocity (MCAv) bilaterally. Finger photoplethysmography was used to record beat-to-beat blood pressure and heart rate. EMG increased with force output from 20 to 60% MVC (P < 0.05). Cerebral HbO2 and Hbtot increased while HHb decreased during contractions with differences observed between 60% vs. 40% and 20% MVC (P < 0.05). Muscle HbO2 decreased while HHb increased during contractions with differences being observed among intensities (P < 0.05). Muscle Hbtot increased from rest at 20% MVC (P < 0.05), while no further change was observed at 40 and 60% MVC (P > 0.05). MCAv increased from rest to exercise but was not different among intensities (P > 0.05). Force output correlated with the root mean square EMG and changes in muscle HbO2 (P < 0.05), but not changes in cerebral HbO2 (P > 0.05) at all three intensities. Force output declined by 8% from the 1st to the 24th contraction only at 60% MVC and was accompanied by systematic increases in RMS, cerebral HbO2 and Hbtot with a leveling off in muscle HbO2 and Hbtot. These changes were independent of alterations in mean arterial pressure. Since cerebral blood flow and oxygenation were elevated at 60% MVC, we attribute the development of fatigue to reduced muscle oxygen availability rather than impaired central neuronal activation.

Robust power allocation algorithms for MIMO OFDM systems with imperfect CSI

This paper presents a Bayesian approach to the design of transmit prefiltering matrices in closed-loop schemes robust to channel estimation errors. The algorithms are derived for a multiple-input multiple-output (MIMO) orthogonal frequency division multiplexing (OFDM) system. Two different optimizationcriteria are analyzed: the minimization of the mean square error and the minimization of the bit error rate. In both cases, the transmitter design is based on the singular value decomposition (SVD) of the conditional mean of the channel response, given the channel estimate. The performance of the proposed algorithms is analyzed,and their relationship with existing algorithms is indicated. As withother previously proposed solutions, the minimum bit error rate algorithmconverges to the open-loop transmission scheme for very poor CSI estimates.

Motion estimation using higher-order statistics

The objective of this paper is to introduce a fourth-order cost function of the displaced frame difference (DFD) capable of estimatingmotion even for small regions or blocks. Using higher than second-orderstatistics is appropriate in case the image sequence is severely corruptedby additive Gaussian noise. Some results are presented and compared to those obtained from the mean kurtosis and the mean square error of the DFD.

Maximum likelihood estimation of position in GNSS

In this letter, we obtain the Maximum LikelihoodEstimator of position in the framework of Global NavigationSatellite Systems. This theoretical result is the basis of a completelydifferent approach to the positioning problem, in contrastto the conventional two-steps position estimation, consistingof estimating the synchronization parameters of the in-viewsatellites and then performing a position estimation with thatinformation. To the authors’ knowledge, this is a novel approachwhich copes with signal fading and it mitigates multipath andjamming interferences. Besides, the concept of Position–basedSynchronization is introduced, which states that synchronizationparameters can be recovered from a user position estimation. Weprovide computer simulation results showing the robustness ofthe proposed approach in fading multipath channels. The RootMean Square Error performance of the proposed algorithm iscompared to those achieved with state-of-the-art synchronizationtechniques. A Sequential Monte–Carlo based method is used todeal with the multivariate optimization problem resulting fromthe ML solution in an iterative way.

Cardiorespiratory and Cardiac Autonomic Responses to 30-15 Intermittent Fitness Test in Team Sport Players

Buchheit, M, Al Haddad, H, Millet GP, Lepretre, PM, Newton, M, and Ahmaidi, S. Cardiorespiratory and cardiac autonomic responses to 30-15 Intermittent Fitness Test in team sport players. J Strength Cond Res 23(1): xxx-xxx, 2009-The 30-15 Intermittent Fitness Test (30-15IFT) is an attractive alternative to classic continuous incremental field tests for defining a reference velocity for interval training prescription in team sport athletes. The aim of the present study was to compare cardiorespiratory and autonomic responses to 30-15IFT with those observed during a standard continuous test (CT). In 20 team sport players (20.9 +/- 2.2 years), cardiopulmonary parameters were measured during exercise and for 10 minutes after both tests. Final running velocity, peak lactate ([La]peak), and rating of perceived exertion (RPE) were also measured. Parasympathetic function was assessed during the postexercise recovery phase via heart rate (HR) recovery time constant (HRRtau) and HR variability (HRV) vagal-related indices. At exhaustion, no difference was observed in peak oxygen uptake (&OV0312;o2peak), respiratory exchange ratio, HR, or RPE between 30-15IFT and CT. In contrast, 30-15IFT led to significantly higher minute ventilation, [La]peak, and final velocity than CT (p < 0.05 for all parameters). All maximal cardiorespiratory variables observed during both tests were moderately to well correlated (e.g., r = 0.76, p = 0.001 for &OV0312;o2peak). Regarding ventilatory thresholds (VThs), all cardiorespiratory measurements were similar and well correlated between the 2 tests. Parasympathetic function was lower after 30-15IFT than after CT, as indicated by significantly longer HHRtau (81.9 +/- 18.2 vs. 60.5 +/- 19.5 for 30-15IFT and CT, respectively, p < 0.001) and lower HRV vagal-related indices (i.e., the root mean square of successive R-R intervals differences [rMSSD]: 4.1 +/- 2.4 and 7.0 +/- 4.9 milliseconds, p < 0.05). In conclusion, the 30-15IFT is accurate for assessing VThs and &OV0312;o2peak, but it alters postexercise parasympathetic function more than a continuous incremental protocol.

Critical flux and fouling in ultrafiltration of proteins

In this thesis different parameters influencing critical flux in protein ultrafiltration and membrane foul-ing were studied. Short reviews of proteins, cross-flow ultrafiltration, flux decline and criticalflux and the basic theory of Partial Least Square analysis (PLS) are given at the beginning. The experiments were mainly performed using dilute solutions of globular proteins, commercial polymeric membranes and laboratory scale apparatuses. Fouling was studied by flux, streaming potential and FTIR-ATR measurements. Critical flux was evaluated by different kinds of stepwise procedures and by both con-stant pressure and constant flux methods. The critical flux was affected by transmembrane pressure, flow velocity, protein concentration, mem-brane hydrophobicity and protein and membrane charges. Generally, the lowest critical fluxes were obtained at the isoelectric points of the protein and the highest in the presence of electrostatic repulsion between the membrane surface and the protein molecules. In the laminar flow regime the critical flux increased with flow velocity, but not any more above this region. An increase in concentration de-creased the critical flux. Hydrophobic membranes showed fouling in all charge conditionsand, furthermore, especially at the beginning of the experiment even at very low transmembrane pressures. Fouling of these membranes was thought to be due to protein adsorption by hydrophobic interactions. The hydrophilic membranes used suffered more from reversible fouling and concentration polarisation than from irreversible foul-ing. They became fouled at higher transmembrane pressures becauseof pore blocking. In this thesis some new aspects on critical flux are presented that are important for ultrafiltration and fractionation of proteins.

Neuro-mechanical and metabolic adjustments to the repeated anaerobic sprint test in professional football players.

PURPOSE: This study aimed to determine the neuro-mechanical and metabolic adjustments in the lower limbs induced by the running anaerobic sprint test (the so-called RAST). METHODS: Eight professional football players performed 6 × 35 m sprints interspersed with 10 s of active recovery on artificial turf with their football shoes. Sprinting mechanics (plantar pressure insoles), root mean square activity of the vastus lateralis (VL), rectus femoris (RF), and biceps femoris (BF) muscles (surface electromyography, EMG) and VL muscle oxygenation (near-infrared spectroscopy) were monitored continuously. RESULTS: Sprint time, contact time and total stride duration increased from the first to the last repetition (+17.4, +20.0 and +16.6 %; all P < 0.05), while flight time and stride length remained constant. Stride frequency (-13.9 %; P < 0.001) and vertical stiffness decreased (-27.2 %; P < 0.001) across trials. Root mean square EMG activities of RF and BF (-18.7 and -18.1 %; P < 0.01 and 0.001, respectively), but not VL (-1.2 %; P > 0.05), decreased over sprint repetitions and were correlated with the increase in running time (r = -0.82 and -0.90; both P < 0.05). Together with a better maintenance of RF and BF muscles activation levels over sprint repetitions, players with a better repeated-sprint performance (lower cumulated times) also displayed faster muscle de- (during sprints) and re-oxygenation (during recovery) rates (r = -0.74 and -0.84; P < 0.05 and 0.01, respectively). CONCLUSION: The repeated anaerobic sprint test leads to substantial alterations in stride mechanics and leg-spring behaviour. Our results also strengthen the link between repeated-sprint ability and the change in neuromuscular activation as well as in muscle de- and re-oxygenation rates.

Plantar flexor neuromuscular adjustments following match-play football in hot and cool conditions.

We assessed neuromuscular fatigue and recovery of the plantar flexors after playing football with or without severe heat stress. Neuromuscular characteristics of the plantar flexors were assessed in 17 male players at baseline and ∼30 min, 24, and 48 h after two 90-min football matches in temperate (∼20 °C and 55% rH) and hot (∼43 °C and 20% rH) environments. Measurements included maximal voluntary strength, muscle activation, twitch contractile properties, and rate of torque development and soleus EMG (i.e., root mean square activity) rise from 0 to 30, -50, -100, and -200 ms during maximal isometric contractions for plantar flexors. Voluntary activation and peak twitch torque were equally reduced (-1.5% and -16.5%, respectively; P < 0.05) post-matches relative to baseline in both conditions, the latter persisting for at least 48 h, whereas strength losses (∼5%) were not significant. Absolute explosive force production declined (P < 0.05) 30 ms after contraction onset independently of condition, with no change at any other epochs. Globally, normalized rate of force development and soleus EMG activity rise values remained unchanged. In football, match-induced alterations in maximal and rapid torque production capacities of the plantar flexors are moderate and do not differ after competing in temperate and hot environments.

Kohinan poisto spektrikuvista pehmeillä morfologisilla suodattimilla

Tämän diplomityön tavoitteena oli tutkia kohinan poistoa spektrikuvista käyttäen pehmeitä morfologisia suodattimia. Työssä painotettiin impulssimaisen kohinan suodattamista. Suodattimien toimintaa arvioitiin numeerisesti keskimääräisen itseisarvovirheen, neliövirheen sekä signaali-kohinasuhteen avulla ja visuaalisesti tarkastelemalla suodatettuja kuvia sekä niiden yksittäisiä spektritasoja. Käytettyjä suodatusmenetelmiä olivat suodatus kuvapisteittäin spektrin suunnassa, suodatus koko spektrissä sekä kuutiomenetelmä ja komponenteittainen suodatus. Suodatettavat kuvat sisälsivät joko suola ja pippuri- tai bittivirhekohinaa. Parhaimmat suodatustulokset sekä numeeristen virhekriteerien että visuaalisen tarkastelun perusteella saatiin komponenteittaisella sekä kuvapisteittäisellä menetelmällä. Työssä käytetyt menetelmät on esitetty algoritmimuodossa. Suodatinalgoritmien toteutukset ja suodatuskokeet tehtiin Matlab-ohjelmistolla.