Scaling law characterizing the dynamics of the transition of HIV-1 to error catastrophe

Increasing the mutation rate, mu, of viruses above a threshold, mu(c), has been predicted to trigger a catastrophic loss of viral genetic information and is being explored as a novel intervention strategy. Here, we examine the dynamics of this transition using stochastic simulations mimicking within-host HIV-1 evolution. We find a scaling law governing the characteristic time of the transition: tau approximate to 0.6/(mu - mu(c)). The law is robust to variations in underlying evolutionary forces and presents guidelines for treatment of HIV-1 infection with mutagens. We estimate that many years of treatment would be required before HIV-1 can suffer an error catastrophe.

Importance of carbon dioxide physiological forcing to future climate change

An increase in atmospheric carbon dioxide (CO2) concentration influences climate both directly through its radiative effect (i.e., trapping longwave radiation) and indirectly through its physiological effect (i.e., reducing transpiration of land plants). Here we compare the climate response to radiative and physiological effects of increased CO2 using the National Center for Atmospheric Research (NCAR) coupled Community Land and Community Atmosphere Model. In response to a doubling of CO2, the radiative effect of CO2 causes mean surface air temperature over land to increase by 2.86 ± 0.02 K (± 1 standard error), whereas the physiological effects of CO2 on land plants alone causes air temperature over land to increase by 0.42 ± 0.02 K. Combined, these two effects cause a land surface warming of 3.33 ± 0.03 K. The radiative effect of doubling CO2 increases global runoff by 5.2 ± 0.6%, primarily by increasing precipitation over the continents. The physiological effect increases runoff by 8.4 ± 0.6%, primarily by diminishing evapotranspiration from the continents. Combined, these two effects cause a 14.9 ± 0.7% increase in runoff. Relative humidity remains roughly constant in response to CO2-radiative forcing, whereas relative humidity over land decreases in response to CO2-physiological forcing as a result of reduced plant transpiration. Our study points to an emerging consensus that the physiological effects of increasing atmospheric CO2 on land plants will increase global warming beyond that caused by the radiative effects of CO2.

Rank dependence of orientational relaxation in dipolar systems

An important yet unsolved problem in the field of orientational relaxation in dipolar liquids is the dependence of the correlation functions C(l)(t), C(l)(t) = [4pi/(2l + 1)SIGMA(m = -l)l [Y(lm)(OMEGA(0)Y(lm)(OMEGA(t))] on the rank l (where Y(lm)(OMEGA) are the usual spherical harmonics). The existing theories on this effect differ in their predictions. To investigate this, we have carried out extensive computer simulations of a Brownian dipolar lattice. The dielectric friction was found to decrease rapidly with increasing l, in qualitative agreement with the predictions of Hubbard-Wolynes. However, the observed effect is much stronger than the predictions of the existing theories.

Non-Gaussianity and direction-dependent systematics in HST Key Project data

Two new statistics, namely Delta(chi 2) and Delta(chi), based on the extreme value theory, were derived by Gupta et al. We use these statistics to study the direction dependence in the HST Key Project data, which provides one of the most precise measurements of the Hubble constant. We also study the non-Gaussianity in this data set using these statistics. Our results for Delta(chi 2) show that the significance of direction-dependent systematics is restricted to well below the 1 sigma confidence limit; however, the presence of non-Gaussian features is subtle. On the other hand, the Delta(chi). statistic, which is more sensitive to direction dependence, shows direction dependence systematics to be at a slightly higher confidence level, and the presence of non-Gaussian features at a level similar to the Delta(chi 2) statistic.

Performance analysis of adaptive physical layer network coding for wireless two-way Relaying

The analysis of modulation schemes for the physical layer network-coded two way relaying scenario is presented which employs two phases: Multiple access (MA) phase and Broadcast (BC) phase. Depending on the signal set used at the end nodes, the minimum distance of the effective constellation seen at the relay becomes zero for a finite number of channel fade states referred as the singular fade states. The singular fade states fall into the following two classes: (i) the ones which are caused due to channel outage and whose harmful effect cannot be mitigated by adaptive network coding called the non-removable singular fade states and (ii) the ones which occur due to the choice of the signal set and whose harmful effects can be removed called the removable singular fade states. In this paper, we derive an upper bound on the average end-to-end Symbol Error Rate (SER), with and without adaptive network coding at the relay, for a Rician fading scenario. It is shown that without adaptive network coding, at high Signal to Noise Ratio (SNR), the contribution to the end-to-end SER comes from the following error events which fall as SNR-1: the error events associated with the removable and nonremovable singular fade states and the error event during the BC phase. In contrast, for the adaptive network coding scheme, the error events associated with the removable singular fade states fall as SNR-2, thereby providing a coding gain over the case when adaptive network coding is not used. Also, it is shown that for a Rician fading channel, the error during the MA phase dominates over the error during the BC phase. Hence, adaptive network coding, which improves the performance during the MA phase provides more gain in a Rician fading scenario than in a Rayleigh fading scenario. Furthermore, it is shown that for large Rician factors, among those removable singular fade states which have the same magnitude, those which have the least absolute value of the phase - ngle alone contribute dominantly to the end-to-end SER and it is sufficient to remove the effect of only such singular fade states.

A Rotor Flux Estimation During Zero and Active Vector Periods Using Current Error Space Vector From a Hysteresis Controller for a Sensorless Vector Control of IM Drive

This paper proposes a sensorless vector control scheme for general-purpose induction motor drives using the current error space phasor-based hysteresis controller. In this paper, a new technique for sensorless operation is developed to estimate rotor voltage and hence rotor flux position using the stator current error during zero-voltage space vectors. It gives a comparable performance with the vector control drive using sensors especially at a very low speed of operation (less than 1 Hz). Since no voltage sensing is made, the dead-time effect and loss of accuracy in voltage sensing at low speed are avoided here, with the inherent advantages of the current error space phasor-based hysteresis controller. However, appropriate device on-state drops are compensated to achieve a steady-state operation up to less than 1 Hz. Moreover, using a parabolic boundary for current error, the switching frequency of the inverter can be maintained constant for the entire operating speed range. Simple sigma L-s estimation is proposed, and the parameter sensitivity of the control scheme to changes in stator resistance, R-s is also investigated in this paper. Extensive experimental results are shown at speeds less than 1 Hz to verify the proposed concept. The same control scheme is further extended from less than 1 Hz to rated 50 Hz six-step operation of the inverter. Here, the magnetic saturation is ignored in the control scheme.

Stochastic Simulations Suggest that HIV-1 Survives Close to Its Error Threshold

The use of mutagenic drugs to drive HIV-1 past its error threshold presents a novel intervention strategy, as suggested by the quasispecies theory, that may be less susceptible to failure via viral mutation-induced emergence of drug resistance than current strategies. The error threshold of HIV-1, mu(c), however, is not known. Application of the quasispecies theory to determine mu(c) poses significant challenges: Whereas the quasispecies theory considers the asexual reproduction of an infinitely large population of haploid individuals, HIV-1 is diploid, undergoes recombination, and is estimated to have a small effective population size in vivo. We performed population genetics-based stochastic simulations of the within-host evolution of HIV-1 and estimated the structure of the HIV-1 quasispecies and mu(c). We found that with small mutation rates, the quasispecies was dominated by genomes with few mutations. Upon increasing the mutation rate, a sharp error catastrophe occurred where the quasispecies became delocalized in sequence space. Using parameter values that quantitatively captured data of viral diversification in HIV-1 patients, we estimated mu(c) to be 7 x 10(-5) -1 x 10(-4) substitutions/site/replication, similar to 2-6 fold higher than the natural mutation rate of HIV-1, suggesting that HIV-1 survives close to its error threshold and may be readily susceptible to mutagenic drugs. The latter estimate was weakly dependent on the within-host effective population size of HIV-1. With large population sizes and in the absence of recombination, our simulations converged to the quasispecies theory, bridging the gap between quasispecies theory and population genetics-based approaches to describing HIV-1 evolution. Further, mu(c) increased with the recombination rate, rendering HIV-1 less susceptible to error catastrophe, thus elucidating an added benefit of recombination to HIV-1. Our estimate of mu(c) may serve as a quantitative guideline for the use of mutagenic drugs against HIV-1.

Co-evolution of specific amino acid in sigma 1.2 region and nucleotide base in the discriminator to act as sensors of small molecule effectors of transcription initiation in mycobacteria

The transcription from rrn and a number of other promoters is regulated by initiating ribonucleotides (iNTPs) and guanosine tetra/penta phosphate (p)ppGpp], either by strengthening or by weakening of the RNA polymerase (RNAP)-promoter interactions during initiation. Studies in Escherichia coli revealed the importance of a sequence termed discriminator, located between -10 and the transcription start site of the responsive promoters in this mode of regulation. Instability of the open complex at these promoters is attributed to the lack of stabilizing interactions between the suboptimal discriminator and the 1.2 region of sigma 70 (Sig70) in RNAP holoenzyme. We demonstrate a different pattern of interaction between the promoters and sigma A (SigA) of Mycobacterium tuberculosis to execute similar regulation. Instead of cytosine and methionine, thymine at three nucleotides downstream to -10 element and leucine 232 in SigA are found to be essential for iNTPs and pppGpp mediated response at the rrn and gyr promoters of the organism. The specificity of the interaction is substantiated by mutational replacements, either in the discriminator or in SigA, which abolish the nucleotide mediated regulation in vitro or in vivo. Specific yet distinct bases and the amino acids appear to have co-evolved' to retain the discriminator-sigma 1.2 region regulatory switch operated by iNTPs/pppGpp during the transcription initiation in different bacteria.

X-ray Structure of Gelonin at 1.8 Å Resolution

Gelonin is a single chain ribosome inactivating protein (RIP) with potential application in the treatment of cancer and AIDS. Diffraction quality crystals grown using PEG3350, belong to the space group P2(1), with it a = 49.4 Angstrom b = 44.9 Angstrom, c = 137.4 Angstrom and beta = 98.4 degrees, and contain two molecules in the asymmetric unit. Diffraction data collected to 1.8 Angstrom resolution has a R(m) value of 7.3%. Structure of gelonin has been solved by the molecular replacement method, using ricin A chain as the search model. Crystallographic refinement using X-PLOR resulted in a model for which the r.m.s deviations from ideal bond lengths and bond angles are 0.012 Angstrom and 2.7 degrees, respectively The final R-factor is 18.4% for 39,806 reflections for which I > 1.0 sigma(I).The C-alpha atoms of the two molecules in the asymmetric unit superpose to within 0.38 Angstrom for 247 atom pairs. The overall fold of gelonin is similar to that of other RIPs such as ricin A chain and alpha-momorcharin, the r.m.s.d. for C-alpha superpositions being 1.3 and 1.4 Angstrom, respectively The-catalytic residues (Glu166, Arg169 and Tyr113) in the active site form a hydrogen bond scheme similar to that observed in other RIPs. The conformation of Tyr74 in the active site, however, is significantly different from that in alpha-momorcharin. Three well defined water molecules are located in the active site cavity and one of them, X319, superposes to within 0.2 Angstrom of a corresponding water molecule in the structure of alpha-momorcharin. Any of the three could be the substrate water molecule in the hydrolysis reaction catalysed by gelonin.Difference electron density for a N-linked sugar moiety has been observed near only one of the two potential glycosylation sites in the sequence. The amino acid at position 239 has been established as Lys by calculation of omit electron density maps.The two cysteine residues in the sequence, Cys44 and Cys50, form a disulphide bond, and are therefore not available for disulphide conjugation with antibodies. Based on the structure, the region of the molecule that is involved in intradimer interactions is suggested to be suitable for introducing a Cys residue for purposes of conjugation with an antibody to produce useful immunotoxins.

Axial binding of 2-methylimidazole to the tetraarylcarboxylatodiruthenium (II, III) core: X-ray crystal structure of [Ru-2(O2CC6H4-p-OMe)(4)(2-mimH)(2)](ClO4). 1.75CH(2)Cl(2).H2O

Diruthenium (II. III) complexes of the type [Ru-2(O2CAr)(4) (2-mimH)(2)](ClO4) (Ar = C6H4-p-X : X=OMe,1, X=Me, 2, 2-mimH=2-methylimidazole) have been isolated from the reaction of Ru2Cl(O2CAr)(4) with 2-mimH in CH2Cl2 followed by the addition of NaClO4. The crystal structure of 1.1.75CH(2)Cl(2).H2O has been determined. The crystal belongs to the monoclinic space group p2(1)/c with the following unit cell dimensions for the C40H40N4O16ClRu2.1.75CH(2)Cl(2).H2O (M = 1237.0) : a = 12.347(3)Angstrom, b = 17.615(5)Angstrom, c = 26.148(2)Angstrom,beta = 92.88(1)degrees. v = 5679(2)Angstrom(3). Z=4, D-c = 1.45 g cm(-3). lambda(Mo-K-alpha) = 0.7107 Angstrom, mu(Mo-K-alpha) = 8.1 cm(-1), T = 293 K, R = 0.0815 (wR(2) = 0.2118) for 5834 reflections with 1 > 2 sigma(I). The complex has a tetracarboxylatodiruthenium (II, III) core and two axially bound 2-methylimidazole ligands. The Ru-Ru bond length is 2.290(1)Angstrom. The Ru-Ru bond order is 2.5 and the complex is three-electron paramagnetic. The complex shows an irreversible Ru-2(II,III)-->Ru-2(Il,II) reduction near -0.2 V vs SCE in CH2Cl2-0. 1 MTBAP. The complexes exemplify the first adduct of the tetracarboxylatodiruthenium (II,III) core having N-donor ligands

Mesh-free approximations via the error reproducing kernel method and applications to nonlinear systems developing shocks

Error estimates for the error reproducing kernel method (ERKM) are provided. The ERKM is a mesh-free functional approximation scheme [A. Shaw, D. Roy, A NURBS-based error reproducing kernel method with applications in solid mechanics, Computational Mechanics (2006), to appear (available online)], wherein a targeted function and its derivatives are first approximated via non-uniform rational B-splines (NURBS) basis function. Errors in the NURBS approximation are then reproduced via a family of non-NURBS basis functions, constructed using a polynomial reproduction condition, and added to the NURBS approximation of the function obtained in the first step. In addition to the derivation of error estimates, convergence studies are undertaken for a couple of test boundary value problems with known exact solutions. The ERKM is next applied to a one-dimensional Burgers equation where, time evolution leads to a breakdown of the continuous solution and the appearance of a shock. Many available mesh-free schemes appear to be unable to capture this shock without numerical instability. However, given that any desired order of continuity is achievable through NURBS approximations, the ERKM can even accurately approximate functions with discontinuous derivatives. Moreover, due to the variation diminishing property of NURBS, it has advantages in representing sharp changes in gradients. This paper is focused on demonstrating this ability of ERKM via some numerical examples. Comparisons of some of the results with those via the standard form of the reproducing kernel particle method (RKPM) demonstrate the relative numerical advantages and accuracy of the ERKM.

R-parameter: A local truncation error based adaptive framework for finite volume compressible flow solvers

A residual-based strategy to estimate the local truncation error in a finite volume framework for steady compressible flows is proposed. This estimator, referred to as the -parameter, is derived from the imbalance arising from the use of an exact operator on the numerical solution for conservation laws. The behaviour of the residual estimator for linear and non-linear hyperbolic problems is systematically analysed. The relationship of the residual to the global error is also studied. The -parameter is used to derive a target length scale and consequently devise a suitable criterion for refinement/derefinement. This strategy, devoid of any user-defined parameters, is validated using two standard test cases involving smooth flows. A hybrid adaptive strategy based on both the error indicators and the -parameter, for flows involving shocks is also developed. Numerical studies on several compressible flow cases show that the adaptive algorithm performs excellently well in both two and three dimensions.

Error correction in nonbinary words

A simple error detecting and correcting procedure is described for nonbinary symbol words; here, the error position is located using the Hamming method and the correct symbol is substituted using a modulo-check procedure.

Construction of Stability Multipliers with Prescribed Phase Characteristics: an Improved Value for $\sigma

For a feedback system consisting of a transfer function $G(s)$ in the forward path and a time-varying gain $n(t)(0 \leqq n(t) \leqq k)$ in the feedback loop, a stability multiplier $Z(s)$ has been constructed (and used to prove stability) by Freedman [2] such that $Z(s)(G(s) + {1 / K})$ and $Z(s - \sigma )(0 < \sigma < \sigma _ * )$ are strictly positive real, where $\sigma _ * $ can be computed from a knowledge of the phase-angle characteristic of $G(i\omega ) + {1 / k}$ and the time-varying gain $n(t)$ is restricted by $\sigma _ * $ by means of an integral inequality. In this note it is shown that an improved value for $\sigma _ * $ is possible by making some modifications in his derivation. ©1973 Society for Industrial and Applied Mathematics.

Multicarrier analysis of magnetotransport data at low and high electric fields

We present some results on multicarrier analysis of magnetotransport data, Both synthetic as well as data from narrow gap Hg0.8Cd0.2Te samples are used to demonstrate applicability of various algorithms vs. nonlinear least square fitting, Quantitative Mobility Spectrum Analysis (QMSA) and Maximum Entropy Mobility Spectrum Analysis (MEMSA). Comments are made from our experience oil these algorithms, and, on the inversion procedure from experimental R/sigma-B to S-mu specifically with least square fitting as an example. Amongst the conclusions drawn are: (i) Experimentally measured resistivity (R-xx, R-xy) should also be used instead of just the inverted conductivity (sigma(xx), sigma(xy)) to fit data to semiclassical expressions for better fits especially at higher B. (ii) High magnetic field is necessary to extract low mobility carrier parameters. (iii) Provided the error in data is not large, better estimates to carrier parameters of remaining carrier species can be obtained at any stage by subtracting highest mobility carrier contribution to sigma from the experimental data and fitting with the remaining carriers. (iv)Even in presence of high electric field, an approximate multicarrier expression can be used to guess the carrier mobilities and their variations before solving the full Boltzmann equation.