Estimating successive pK(a) values of polyprotic acids from ab initio molecular dynamics using metadynamics: the dissociation of phthalic acid and its isomers

Estimation of the dissociation constant, or pK(a), of weak acids continues to be a central goal in theoretical chemistry. Here we show that ab initio Car-Parrinello molecular dynamics simulations in conjunction with metadynamics calculations of the free energy profile of the dissociation reaction can provide reasonable estimates of the successive pK(a) values of polyprotic acids. We use the distance-dependent coordination number of the protons bound to the hydroxyl oxygen of the carboxylic group as the collective variable to explore the free energy profile of the dissociation process. Water molecules, sufficient to complete three hydration shells surrounding the acid molecule, were included explicitly in the computation procedure. Two distinct minima corresponding to the dissociated and un-dissociated states of the acid are observed and the difference in their free energy values provides the estimate for pK(a), the acid dissociation constant. We show that the method predicts the pK(a) value of benzoic acid in good agreement with experiment and then show using phthalic acid (benzene dicarboxylic acid) as a test system that both the first and second pK(a) values as well, as the subtle difference in their values for different isomers can be predicted in reasonable agreement with experimental data.

First principles study of structural stability and site preference in Co-3 (W,X)

Since the discovery 1] of gamma' precipitate (L1(2) - Co-3 (Al, W)) in the Co-Al-W ternary system, there has been an increased interest in Co-based superalloys. Since these alloys have two phase microstructures (gamma + gamma') similar to Ni-based superalloys 2], they are viable candidates in high temperature applications, particularly in land-based turbines. The role of alloying on stability of the gamma' phase has been an active area of research. In this study, electronic structure calculations were done to probe the effect of alloying in Co3W with L1(2) structure. Compositions of type Co-3(W, X), (where X/Y = Mn, Fe, Ni, Pt, Cr, Al, Si, V, W, Ta, Ti, Nb, Hf, Zr and Mo) were studied. Effect of alloying on equilibrium lattice parameters and ground state energies was used to calculate Vegard's coefficients and site preference related data. The effect of alloying on the stability of the L1(2) structure vis a vis other geometrically close packed ordered structures was also studied for a range of Co3X compounds. Results suggest that the penchant of element for the W sublattice can be predicted by comparing heats of formation of Co3X in different structures.

Design of a Portable Compliant Device for Estimating the Failure-Load of Mesoscale Cemented Sand Specimens

In this paper, we present the design and development of a portable, hand-operated composite compliant mechanism for estimating the failure-load of cm-sized stiff objects whose stiffness is of the order of 10 s of kN/m. The motivation for the design comes from the need to estimate the failure-load of mesoscale cemented sand specimens in situ, which is not possible with traditional devices used for large specimens or very small specimens. The composite compliant device, developed in this work, consists of two compliant mechanisms: a force-amplifying compliant mechanism (FaCM) to amplify sufficiently the force exerted by hand in order to break the specimen and a displacement-amplifying compliant mechanism (DaCM) to enable measurement of the force using a proximity sensor. The two mechanisms are designed using the selection-maps technique to amplify the force up to 100N by about a factor of 3 and measure the force with a resolution of 15 mN. The composite device, made using a FaCM, a DaCM, and a Hall effect-based proximity sensor, was tested on mesoscale cemented sand specimens that were 10mm in diameter and 20mm in length. The results are compared with those of a large commercial instrument. Through the experiments, it was observed that the failure-load of the cemented sand specimens varied from 0.95N to 24.33 N, depending on the percentage of cementation and curing period. The estimation of the failure-load using the compliant device was found to be within 1.7% of the measurements obtained using the commercial instrument and thus validating the design. The details of the design, prototyping, specimen preparation, testing, and the results comprise the paper.

Analytical approach to quantile estimation in regional frequency analysis based on fuzzy framework

Regional frequency analysis is widely used for estimating quantiles of hydrological extreme events at sparsely gauged/ungauged target sites in river basins. It involves identification of a region (group of watersheds) resembling watershed of the target site, and use of information pooled from the region to estimate quantile for the target site. In the analysis, watershed of the target site is assumed to completely resemble watersheds in the identified region in terms of mechanism underlying generation of extreme event. In reality, it is rare to find watersheds that completely resemble each other. Fuzzy clustering approach can account for partial resemblance of watersheds and yield region(s) for the target site. Formation of regions and quantile estimation requires discerning information from fuzzy-membership matrix obtained based on the approach. Practitioners often defuzzify the matrix to form disjoint clusters (regions) and use them as the basis for quantile estimation. The defuzzification approach (DFA) results in loss of information discerned on partial resemblance of watersheds. The lost information cannot be utilized in quantile estimation, owing to which the estimates could have significant error. To avert the loss of information, a threshold strategy (TS) was considered in some prior studies. In this study, it is analytically shown that the strategy results in under-prediction of quantiles. To address this, a mathematical approach is proposed in this study and its effectiveness in estimating flood quantiles relative to DFA and TS is demonstrated through Monte-Carlo simulation experiments and case study on Mid-Atlantic water resources region, USA. (C) 2015 Elsevier B.V. All rights reserved.

Nest site characterization of sympatric hornbills in a tropical dry forest

Hornbills, among the largest and most threatened tropical frugivores, provide important seed dispersal services. Hornbill nest site characteristics are known primarily from wet tropical forests. Nests of the Indian grey hornbill Ocyceros birostris and Oriental pied hornbill Anthracoceros albirostris were characterized in a tropical dry forest. Despite A. albirostris being twice the size of O. birostris, few of the nest cavity attributes were different. A. albirostris nests were surrounded by higher proportion of mixed forest and lower sal forest compared to O. birostris. In this landscape, the larger A. albirostris may prefer to nest in sites with more food plants compared to the smaller O. birostris.

Use of permeable reactive barrier to mitigate groundwater nitrate contamination from on-site sanitation

Nearly 50% of India's population depends on variants of pit-toilet systems for human waste disposal. Nitrate contamination of groundwater by pit-toilet leachate is a major environmental concern in the country as it sources a major proportion (50-80%) of potable water from aquifers. Therefore, minimizing nitrate contamination of groundwater due to leachate infiltration from pit-toilet systems is essential. Batch and column experiments demonstrated the capability of bentonite-enhanced sand (BES) specimens to reduce nitrate concentrations in synthetic solutions (initial NO3-N concentration = 22.7 mg/L, C/N = 3) by about 85-90% in 10 to 24 hour by a heterotrophic denitrification process. Based on the laboratory results, it is recommended that use of a BES-permeable reactive barrier layer at the base of pit-toilets will facilitate heterotrophic denitrification and mitigate nitrate contamination of the underlying aquifer.

Efficient Site-Specific Incorporation of Thioamides into Peptides on a Solid Support

Designing bioactive peptides containing thioamide functionality to modulate their pharmacological properties has been thwarted so far because of various synthetic challenges. The fast, efficient, and inexpensive synthesis and incorporation of a wide range of thionated amino acids into a growing peptide chain on a solid support is reported using standard Fmoc-based chemistry. The commonly employed methodology is comprehensively investigated and optimized with significant improvements regarding the quantity of reagents and reaction conditions. The utility of the protocol is further demonstrated in the synthesis of dithionated linear and monothionated cyclic peptides, which has been a daunting task.

Energy Harvesting WSNs for Accurately Estimating the Maximum Sensor Reading: Trade-Offs and Optimal Design

Computing the maximum of sensor readings arises in several environmental, health, and industrial monitoring applications of wireless sensor networks (WSNs). We characterize the several novel design trade-offs that arise when green energy harvesting (EH) WSNs, which promise perpetual lifetimes, are deployed for this purpose. The nodes harvest renewable energy from the environment for communicating their readings to a fusion node, which then periodically estimates the maximum. For a randomized transmission schedule in which a pre-specified number of randomly selected nodes transmit in a sensor data collection round, we analyze the mean absolute error (MAE), which is defined as the mean of the absolute difference between the maximum and that estimated by the fusion node in each round. We optimize the transmit power and the number of scheduled nodes to minimize the MAE, both when the nodes have channel state information (CSI) and when they do not. Our results highlight how the optimal system operation depends on the EH rate, availability and cost of acquiring CSI, quantization, and size of the scheduled subset. Our analysis applies to a general class of sensor reading and EH random processes.

Effects of site stiffness and source to receiver distance on surface wave tests' results

By using six 4.5 Hz geophones, surface wave tests were performed on four different sites by dropping freely a 65 kg mass from a height of 5 m. The receivers were kept far away from the source to eliminate the arrival of body waves. Three different sources to nearest receiver distances (S), namely, 46 m, 56 m and 66 m, were chosen. Dispersion curves were drawn for all the sites. The maximum wavelength (lambda(max)), the maximum depth (d(max)) up to which exploration can be made and the frequency content of the signals depends on the site stiffness and the value of S. A stiffer site yields greater values of lambda(max) and d(max). For stiffer sites, an increase in S leads to an increase in lambda(max). The predominant time durations of the signals increase from stiffer to softer sites. An inverse analysis was also performed based on the stiffness matrix approach in conjunction with the maximum vertical flexibility coefficient of ground surface to establish the governing mode of excitation. For the Site 2, the results from the surface wave tests were found to compare reasonably well with that determined on the basis of cross boreholes seismic tests. (C) 2015 Elsevier Ltd. All rights reserved.

Ab Initio Molecular Dynamics Simulations of Amino Acids in Aqueous Solutions: Estimating pK(a) Values from Metadynamics Sampling

Changes in the protonation and deprotonation of amino acid residues in proteins play a key role in many biological processes and pathways. Here, we report calculations of the free-energy profile for the protonation deprotonation reaction of the 20 canonical alpha amino acids in aqueous solutions using ab initio Car-Parrinello molecular dynamics simulations coupled with metad-ynamics sampling. We show here that the calculated change in free energy of the dissociation reaction provides estimates of the multiple pK(a) values of the amino acids that are in good agreement with experiment. We use the bond-length-dependent number of the protons coordinated to the hydroxyl oxygen of the carboxylic and the amine groups as the collective variables to explore the free-energy profiles of the Bronsted acid-base chemistry of amino acids in aqueous solutions. We ensure that the amino acid undergoing dissociation is solvated by at least three hydrations shells with all water molecules included in the simulations. The method works equally well for amino acids with neutral, acidic and basic side chains and provides estimates of the multiple pK(a) values with a mean relative error, with respect to experimental results, of 0.2 pK(a) units.

Electronic Structure Evolution across the Peierls Metal-Insulator Transition in a Correlated Ferromagnet

Transition metal compounds often undergo spin-charge-orbital ordering due to strong electron-electron correlations. In contrast, low-dimensional materials can exhibit a Peierls transition arising from low-energy electron-phonon-coupling-induced structural instabilities. We study the electronic structure of the tunnel framework compound K2Cr8O16, which exhibits a temperature-dependent (T-dependent) paramagnetic-to-ferromagnetic- metal transition at T-C = 180 K and transforms into a ferromagnetic insulator below T-MI = 95 K. We observe clear T-dependent dynamic valence (charge) fluctuations from above T-C to T-MI, which effectively get pinned to an average nominal valence of Cr+3.75 (Cr4+:Cr3+ states in a 3:1 ratio) in the ferromagnetic-insulating phase. High-resolution laser photoemission shows a T-dependent BCS-type energy gap, with 2G(0) similar to 3.5(k(B)T(MI)) similar to 35 meV. First-principles band-structure calculations, using the experimentally estimated on-site Coulomb energy of U similar to 4 eV, establish the necessity of strong correlations and finite structural distortions for driving the metal-insulator transition. In spite of the strong correlations, the nonintegral occupancy (2.25 d-electrons/Cr) and the half-metallic ferromagnetism in the t(2g) up-spin band favor a low-energy Peierls metal-insulator transition.

Site-Selective Addition of Maleimide to Indole at the C-2 Position: Ru(II)-Catalyzed C-H Activation

Synthesis of 3-(indol-2-yl)succinimide derivatives is presented using a directing group strategy. Selective functionalization of C-2 in the presence of highly reactive C-3 in indole derivatives has been achieved. A conjugate addition product instead of Heck-type product has been brought about by careful selection of the alkene partner (maleimides and maleate esters) such that a beta-hydride elimination is avoided.

Novel pppGpp binding site at the C-terminal region of the Rel enzyme from Mycobacterium smegmatis

Mycobacterium tuberculosis elicits the stringent response under unfavorable growth conditions, such as those encountered by the pathogen inside the host. The hallmark of this response is production of guanosine tetra-and pentaphosphates, collectively termed (p)ppGpp, which have pleiotropic effects on the bacterial physiology. As the stringent response is connected to survival under stress, it is now being targeted for developing inhibitors against bacterial persistence. The Rel enzyme in mycobacteria has two catalytic domains at its N-terminus that are involved in the synthesis and hydrolysis of (p)ppGpp, respectively. However, the function of the C-terminal region of the protein remained unknown. Here, we have identified a binding site for pppGpp in the C-terminal region of Rel. The binding affinity of pppGpp was quantified by isothermal titration calorimetry. The binding site was determined by crosslinking using the nucleotide analog azido-pppGpp, and examining the crosslink product by mass spectrometry. Additionally, mutations in the Rel protein were created to confirm the site of pppGpp binding by isothermal titration calorimetry. These mutants showed increased pppGpp synthesis and reduced hydrolytic activity. We believe that binding of pppGpp to Rel provides a feedback mechanism that allows the protein to detect and adjust the (p)ppGpp level in the cell. Our work suggests that such sites should also be considered while designing inhibitors to target the stringent response.

Probing the role of highly conserved residues in triosephosphate isomerase-analysis of site specific mutants at positions 64 and 75 in the Plasmodial enzyme

Highly conserved residues in enzymes are often found to be clustered close to active sites, suggesting that functional constraints dictate the nature of amino acid residues accommodated at these sites. Using the Plasmodiumfalciparum triosephosphate isomerase (PfTIM) enzyme () as a template, we have examined the effects of mutations at positions 64 and 75, which are not directly involved in the proton transfer cycle. Thr (T) occurring at position 75 is completely conserved, whereas only Gln (Q) and Glu (E) are accommodated at position 64. Biophysical and kinetic data are reported for four T75 (T75S/V/C/N) and two Q64 (Q64N/E) mutants. The dimeric structure is weakened in the Q64E and Q64N mutants, whereas dimer integrity is unimpaired in all four T75 mutants. Measurement of the concentration dependence of enzyme activity permits an estimate of K-d values for dimer dissociation (Q64N=73.79.2nm and Q64E=44.6 +/- 8.4nm). The T75S/V/C mutants have activities comparable to the wild-type enzyme, whereas a fourfold drop is observed for T75N. All four T75 mutants show a dramatic fall in activity between 35 degrees C and 45 degrees C. Crystal structure determination of the T75S/V/N mutants provides insights into the variations in local interactions, with the T75N mutant showing the largest changes. Hydrogen-bond interactions determine dimer stability restricting the choice of residues at position 64 to Gln (Q) and Glu (E). At position 75, the overwhelming preference for Thr (T) may be dictated by the imperative of maintaining temperature stability of enzyme activity.

Plant reproductive traits mediate tritrophic feedback effects within an obligate brood-site pollination mutualism

Plants, herbivores and parasitoids affect each other directly and indirectly; however, feedback effects mediated by host plant traits have rarely been demonstrated in these tritrophic interactions. Brood-site pollination mutualisms (e.g. those involving figs and fig wasps) represent specialised tritrophic communities where the progeny of mutualistic pollinators and of non-mutualistic gallers (both herbivores) together with that of their parasitoids develop within enclosed inflorescences called syconia (hence termed brood-sites or microcosms). Plant reproductive phenology (which affects temporal brood-site availability) and inflorescence size (representing brood-site size) are plant traits that could affect reproductive resources, and hence relationships between trees, pollinators and non-pollinating wasps. Analysing wasp and seed contents of syconia, we examined direct, indirect, trophic and non-trophic relationships within the interaction web of the fig-fig wasp community of Ficus racemosa in the context of brood site size and availability. We demonstrate that in addition to direct resource competition and predator-prey (host-parasitoid) interactions, these communities display exploitative or apparent competition and trait-mediated indirect interactions. Inflorescence size and plant reproductive phenology impacted plant-herbivore and plant-parasitoid associations. These plant traits also influenced herbivore-herbivore and herbivore-parasitoid relationships via indirect effects. Most importantly, we found a reciprocal effect between within-tree reproductive asynchrony and fig wasp progeny abundances per syconium that drives a positive feedback cycle within the system. The impact of a multitrophic feedback cycle within a community built around a mutualistic core highlights the need for a holistic view of plant-herbivore-parasitoid interactions in the community ecology of mutualisms.