Glycogen phosphorylase inhibitors: a free energy perturbation analysis of glucopyranose spirohydantoin analogues

GP catalyzes the phosphorylation of glycogen to Glc-1-P. Because of its fundamental role in the metabolism of glycogen, GP has been the target for a systematic structure-assisted design of inhibitory compounds, which could be of value in the therapeutic treatment of type 2 diabetes mellitus. The most potent catalytic-site inhibitor of GP identified to date is spirohydantoin of glucopyranose (hydan). In this work, we employ MD free energy simulations to calculate the relative binding affinities for GP of hydan and two spirohydantoin analogues, methyl-hydan and n-hydan, in which a hydrogen atom is replaced by a methyl- or amino group, respectively. The results are compared with the experimental relative affinities of these ligands, estimated by kinetic measurements of the ligand inhibition constants. The calculated binding affinity for methyl-hydan (relative to hydan) is 3.75 +/- 1.4 kcal/mol, in excellent agreement with the experimental value (3.6 +/- 0.2 kcal/mol). For n-hydan, the calculated value is 1.0 +/- 1.1 kcal/mol, somewhat smaller than the experimental result (2.3 +/- 0.1 kcal/mol). A free energy decomposition analysis shows that hydan makes optimum interactions with protein residues and specific water molecules in the catalytic site. In the other two ligands, structural perturbations of the active site by the additional methyl- or amino group reduce the corresponding binding affinities. The computed binding free energies are sensitive to the preference of a specific water molecule for two well-defined positions in the catalytic site. The behavior of this water is analyzed in detail, and the free energy profile for the translocation of the water between the two positions is evaluated. The results provide insights into the role of water molecules in modulating ligand binding affinities. A comparison of the interactions between a set of ligands and their surrounding groups in X-ray structures is often used in the interpretation of binding free energy differences and in guiding the design of new ligands. For the systems in this work, such an approach fails to estimate the order of relative binding strengths, in contrast to the rigorous free energy treatment.

Multiple-time higher-order perturbation analysis of the regularized long-wavelength equation

By considering the long-wavelength limit of the regularized long wave (RLW) equation, we study its multiple-time higher-order evolution equations. As a first result, the equations of the Korteweg-de Vries hierarchy are shown to play a crucial role in providing a secularity-free perturbation theory in the specific case of a solitary-wave solution. Then, as a consequence, we show that the related perturbative series can be summed and gives exactly the solitary-wave solution of the RLW equation. Finally, some comments and considerations are made on the N-soliton solution, as well as on the limitations of applicability of the multiple-scale method in obtaining uniform perturbative series.

Kinematic analysis and modelling of biped locomotion systems

This paper presents the kinematic study of robotic biped locomotion systems. The main purpose is to determine the kinematic characteristics and the system performance during walking. For that objective, the prescribed motion of the biped is completely characterised in terms of five locomotion variables: step length, hip height, maximum hip ripple, maximum foot clearance and link lengths. In this work, we propose four methods to quantitatively measure the performance of the walking robot: energy analysis, perturbation analysis, lowpass frequency response and locomobility measure. These performance measures are discussed and compared in determining the robustness and effectiveness of the resulting locomotion.

A new approach to evaluate imperfection sensitivity in asymmetric bifurcation buckling analysis

A direct procedure for the evaluation of imperfection sensitivity in bifurcation problems is presented. The problems arise in the context of the general theory of elastic stability for discrete structural systems, in which the energy criterion of stability of structures and the total potential energy formulation are employed. In cases of bifurcation buckling the sensitivity of the critical load with respect to an imperfection parameter e is singular at the state given by epsilon =0, so that, a regular perturbation expansion of the solution is not possible. In this work we describe a direct procedure to obtain the relations between the critical loads, the generalized coordinates at the critical state, the eigenvector, and the amplitude of the imperfection, using singular perturbation analysis. The expansions are assumed in terms of arbitrary powers of the imperfection parameter, so that both exponents and coefficients of the expansion are unknown. The solution of the series exponents is obtained by searching the least degenerate solution. The formulation is here applied to asymmetric bifurcations, for which explicit expressions of the coefficients are obtained. The use of the method is illustrated by a simple example, which allows consideration of the main features of the formulation.

Design and analysis of distributed utility maximization algorithm for multihop wireless network with inaccurate feedback

Distributed network utility maximization (NUM) is receiving increasing interests for cross-layer optimization problems in multihop wireless networks. Traditional distributed NUM algorithms rely heavily on feedback information between different network elements, such as traffic sources and routers. Because of the distinct features of multihop wireless networks such as time-varying channels and dynamic network topology, the feedback information is usually inaccurate, which represents as a major obstacle for distributed NUM application to wireless networks. The questions to be answered include if distributed NUM algorithm can converge with inaccurate feedback and how to design effective distributed NUM algorithm for wireless networks. In this paper, we first use the infinitesimal perturbation analysis technique to provide an unbiased gradient estimation on the aggregate rate of traffic sources at the routers based on locally available information. On the basis of that, we propose a stochastic approximation algorithm to solve the distributed NUM problem with inaccurate feedback. We then prove that the proposed algorithm can converge to the optimum solution of distributed NUM with perfect feedback under certain conditions. The proposed algorithm is applied to the joint rate and media access control problem for wireless networks. Numerical results demonstrate the convergence of the proposed algorithm. © 2013 John Wiley & Sons, Ltd.

Design and analysis of distributed utility maximization algorithm for multihop wireless network with inaccurate feedback

Distributed network utility maximization (NUM) is receiving increasing interests for cross-layer optimization problems in multihop wireless networks. Traditional distributed NUM algorithms rely heavily on feedback information between different network elements, such as traffic sources and routers. Because of the distinct features of multihop wireless networks such as time-varying channels and dynamic network topology, the feedback information is usually inaccurate, which represents as a major obstacle for distributed NUM application to wireless networks. The questions to be answered include if distributed NUM algorithm can converge with inaccurate feedback and how to design effective distributed NUM algorithm for wireless networks. In this paper, we first use the infinitesimal perturbation analysis technique to provide an unbiased gradient estimation on the aggregate rate of traffic sources at the routers based on locally available information. On the basis of that, we propose a stochastic approximation algorithm to solve the distributed NUM problem with inaccurate feedback. We then prove that the proposed algorithm can converge to the optimum solution of distributed NUM with perfect feedback under certain conditions. The proposed algorithm is applied to the joint rate and media access control problem for wireless networks. Numerical results demonstrate the convergence of the proposed algorithm. © 2013 John Wiley & Sons, Ltd.

Reactivity of chars and carbons: New insights through molecular modeling

A new model proposed for the gasification of chars and carbons incorporates features of the turbostratic nanoscale structure that exists in such materials. The model also considers the effect of initial surface chemistry and different reactivities perpendicular to the edges and to the faces of the underlying crystallite planes comprising the turbostratic structure. It may be more realistic than earlier models based on pore or grain structure idealizations when the carbon contains large amounts of crystallite matter. Shrinkage of the carbon particles in the chemically controlled regime is also possible due to the random complete gasification of crystallitic planes. This mechanism can explain observations in the literature of particle size reduction. Based on the model predictions, both initial surface chemistry and the number of stacked planes in the crystallites strongly influence the reactivity and particle shrinkage. Its test results agree well with literature data on the air-oxidation of Spherocarb and show that it accurately predicts the variation of particle size with conversion. Model parameters are determined entirely from rate measurements.

Performance Issues in Biped Walking Robots

This paper presents the dynamic analysis of robotic biped systems. The main goal is to gain insight into the phenomena of walking and to evaluate its performance. In this study, we propose three methods to quantitatively measure the dynamic efficiency of walking: energy analysis, perturbation analysis and lowpass frequency analysis. In order to accomplish this goal, the prescribed motion of the biped is completely characterised in terms of a set of locomotion variables, namely: step lenght, hip height, hip ripple, hip offset, foot clearance and link lenghts. based on these variables and their influence, the performance measures aer discussed and the results compared with those observed in nature.

Speed of reaction-diffusion fronts in spatially heterogeneous media

The front speed problem for nonuniform reaction rate and diffusion coefficient is studied by using singular perturbation analysis, the geometric approach of Hamilton-Jacobi dynamics, and the local speed approach. Exact and perturbed expressions for the front speed are obtained in the limit of large times. For linear and fractal heterogeneities, the analytic results have been compared with numerical results exhibiting a good agreement. Finally we reach a general expression for the speed of the front in the case of smooth and weak heterogeneities

Modulating the electronic structure of amino acids: interaction of model lewis acids with anthranilic acid

On the basis of theoretical B3LYP calculations, Yáñez and co-workers (J. Chem. Theory Comput. 2012, 8, 2293) illustrated that beryllium ions are capable of significantly modulating (changing) the electronic structures of imidazole. In this computational organic chemistry study, the interaction of this β-amino acid and five model Lewis acids (BeF1+, Be2+, AlF2(1+), AlF2+, and Al3+) were investigated. Several aspects were addressed: natural bond orbitals, including second order perturbation analysis of intra-molecular charge delocalization and the natural population analysis atomic charges; molecular geometries; selected infrared stretching frequencies (C-N, C-O, and N-H), and selected ¹H-NMR chemical shifts. The data illustrate that this interaction can weaken the H-O bond and goes beyond strengthening the intra-molecular hydrogen bond (N...H-O) to cause a spontaneous transfer of the proton to the nitrogen atom in five cases generating zwitterion structures. Many new features are observed. Most importantly, the zwitterion structures include a stabilizing hydrogen bond (N-H...O) that varies in relative strength according to the Lewis acid. These findings explain the experimental observations of α-amino acids (for example: J. Am. Chem. Soc. 2001, 123, 3577) and are the first reported fundamental electronic structure characterization of β-amino acids in zwitterion form.

Lesser Scaup Population Dynamics: What Can Be Learned from Available Data?

Populations of Lesser Scaup (Aythya affinis) have declined markedly in North America since the early 1980s. When considering alternatives for achieving population recovery, it would be useful to understand how the rate of population growth is functionally related to the underlying vital rates and which vital rates affect population growth rate the most if changed (which need not be those that influenced historical population declines). To establish a more quantitative basis for learning about life history and population dynamics of Lesser Scaup, we summarized published and unpublished estimates of vital rates recorded between 1934 and 2005, and developed matrix life-cycle models with these data for females breeding in the boreal forest, prairie-parklands, and both regions combined. We then used perturbation analysis to evaluate the effect of changes in a variety of vital-rate statistics on finite population growth rate and abundance. Similar to Greater Scaup (Aythya marila), our modeled population growth rate for Lesser Scaup was most sensitive to unit and proportional change in adult female survival during the breeding and non-breeding seasons, but much less so to changes in fecundity parameters. Interestingly, population growth rate was also highly sensitive to unit and proportional changes in the mean of nesting success, duckling survival, and juvenile survival. Given the small samples of data for key aspects of the Lesser Scaup life cycle, we recommend additional research on vital rates that demonstrate a strong effect on population growth and size (e.g., adult survival probabilities). Our life-cycle models should be tested and regularly updated in the future to simultaneously guide science and management of Lesser Scaup populations in an adaptive context.

Evaluating the Small Population Paradigm for Rare Large-Bodied Woodpeckers, with Implications for the Ivory-billed Woodpecker

Six large-bodied, ≥ 120 g, woodpecker species are listed as near-threatened to critically endangered by the International Union for Conservation of Nature (IUCN). The small population paradigm assumes that these populations are likely to become extinct without an increase in numbers, but the combined influences of initial population size and demographic rates, i.e., annual adult survival and fecundity, may drive population persistence for these species. We applied a stochastic, stage-based single-population model to available demographic rates for Dryocopus and Campephilus woodpeckers. In particular, we determined the change in predicted extinction rate, i.e., proportion of simulated populations that went extinct within 100 yr, to concomitant changes in six input parameters. To our knowledge, this is the first study to evaluate the combined importance of initial population size and demographic rates for the persistence of large-bodied woodpeckers. Under a worse-case scenario, the median time to extinction was 7 yr (range: 1–32). Across the combinations of other input values, increasing initial population size by one female induced, on average, 0.4%–3.2% (range: 0%–28%) reduction in extinction rate. Increasing initial population size from 5–30 resulted in extinction rates < 0.05 under limited conditions: (1) all input values were intermediate, or (2) Allee effect present and annual adult survival ≥ 0.8. Based on our model, these species can persist as rare, as few as five females, and thus difficult-to-detect, populations provided they maintain ≥ 1.1 recruited females annually per adult female and an annual adult survival rate ≥ 0.8. Athough a demographic-based population viability analysis (PVA) is useful to predict how extinction rate changes across scenarios for life-history attributes, the next step for modeling these populations should incorporate more easily acquired data on changes in patch occupancy to make predictions about patch colonization and extinction rates.

Evaluation of peak-picking algorithms for protein mass spectrometry

Peak picking is an early key step in MS data analysis. We compare three commonly used approaches to peak picking and discuss their merits by means of statistical analysis. Methods investigated encompass signal-to-noise ratio, continuous wavelet transform, and a correlation-based approach using a Gaussian template. Functionality of the three methods is illustrated and discussed in a practical context using a mass spectral data set created with MALDI-TOF technology. Sensitivity and specificity are investigated using a manually defined reference set of peaks. As an additional criterion, the robustness of the three methods is assessed by a perturbation analysis and illustrated using ROC curves.

Evaluation of two "integrated" polarimetric Quantitative Precipitation Estimation (QPE) algorithms at C-band

Two so-called “integrated” polarimetric rate estimation techniques, ZPHI (Testud et al., 2000) and ZZDR (Illingworth and Thompson, 2005), are evaluated using 12 episodes of the year 2005 observed by the French C-band operational Trappes radar, located near Paris. The term “integrated” means that the concentration parameter of the drop size distribution is assumed to be constant over some area and the algorithms retrieve it using the polarimetric variables in that area. The evaluation is carried out in ideal conditions (no partial beam blocking, no ground-clutter contamination, no bright band contamination, a posteriori calibration of the radar variables ZH and ZDR) using hourly rain gauges located at distances less than 60 km from the radar. Also included in the comparison, for the sake of benchmarking, is a conventional Z = 282R1.66 estimator, with and without attenuation correction and with and without adjustment by rain gauges as currently done operationally at Météo France. Under those ideal conditions, the two polarimetric algorithms, which rely solely on radar data, appear to perform as well if not better, pending on the measurements conditions (attenuation, rain rates, …), than the conventional algorithms, even when the latter take into account rain gauges through the adjustment scheme. ZZDR with attenuation correction is the best estimator for hourly rain gauge accumulations lower than 5 mm h−1 and ZPHI is the best one above that threshold. A perturbation analysis has been conducted to assess the sensitivity of the various estimators with respect to biases on ZH and ZDR, taking into account the typical accuracy and stability that can be reasonably achieved with modern operational radars these days (1 dB on ZH and 0.2 dB on ZDR). A +1 dB positive bias on ZH (radar too hot) results in a +14% overestimation of the rain rate with the conventional estimator used in this study (Z = 282R^1.66), a -19% underestimation with ZPHI and a +23% overestimation with ZZDR. Additionally, a +0.2 dB positive bias on ZDR results in a typical rain rate under- estimation of 15% by ZZDR.