Molecular Dynamics Study of the Structure, Flexibility, and Hydrophilicity of PETIM Dendrimers: A Comparison with PAMAM Dendrimers

A new class of dendrimers, the poly(propyl ether imine) (PETIM) dendrimer, has been shown to be a novel hyperbranched polymer having potential applications as a drug delivery vehicle. Structure and dynamics of the amine terminated PETIM dendrimer and their changes with respect to the dendrimer generation are poorly understood. Since most drugs are hydrophobic in nature, the extent of hydrophobicity of the dendrimer core is related to its drug encapsulation and retention efficacy. In this study, we carry out fully atomistic molecular dynamics (MD) simulations to characterize the structure of PETIM (G2-G6) dendrimers in salt solution as a function of dendrimer generation at different protonation levels. Structural properties such as radius of gyration (R-g), radial density distribution, aspect ratio, and asphericity are calculated. In order to assess the hydrophilicity of the dendrimer, we compute the number of bound water molecules in the interior of dendrirner as well as the number of dendrimer-water hydrogen bonds. We conclude that PETIM dendrimers have relatively greater hydrophobicity and flexibility when compared with their extensively investigated PAMAM counterparts. Hence PETIM dendrimers are expected to have stronger interactions with lipid membranes as well as improved drug encapsulation and retention properties when compared with PAMAM dendrimers. We compute the root-mean-square fluctuation of dendrimers as well as their entropy to quantify the flexibility of the dendrimer. Finally we note that structural and solvation properties computed using force field parameters derived based on the CHARMM general purpose force field were in good quantitative agreement with those obtained using the generalized Amber force field (GAFF).

Seismic intensity map of South India for estimated future earthquakes

In this study, an attempt has been made to prepare the seismic intensity map for south India considering the probable earthquakes in the region. Anbazhagan et al. (Nat Hazards 60:1325-1345, 2012) have identified eight probable future earthquake zones in south India based on rupture-based seismic hazard analysis. Anbazhagan et al. (Eng Geol 171:81-95, 2014) has estimated the maximum future earthquake magnitude at these eight zones using regional rupture character. In this study, the whole south India is divided into several grids of size 1(o) x 1(o) and the intensity at each grid point is calculated using the regional intensity model for the maximum earthquake magnitude at each of the eight zones. The intensity due to earthquakes at these zones is mapped and thus eight seismic intensity maps are prepared. The final seismic intensity map of south India is obtained by considering the maximum intensity at each grid point due to the estimated earthquakes. By looking at the seismic intensity map, one can expect slight to heavy damage due to the probable earthquake magnitudes. Heavy damage may happen close to the probable earthquake zones.

Interlaboratory comparison of magnesium isotopic compositions of 12 felsic to ultramafic igneous rock standards analyzed by MC-ICPMS

To evaluate the interlaboratory mass bias for high-precision stable Mg isotopic analysis of natural materials, a suite of silicate standards ranging in composition from felsic to ultramafic were analyzed in five laboratories by using three types of multicollector inductively coupled plasma mass spectrometer (MC-ICPMS). Magnesium isotopic compositions from all labs are in agreement for most rocks within quoted uncertainties but are significantly (up to 0.3 parts per thousand in Mg-26/Mg-24, > 4 times of uncertainties) different for some mafic samples. The interlaboratory mass bias does not correlate with matrix element/Mg ratios, and the mechanism for producing it is uncertain but very likely arises from column chemistry. Our results suggest that standards with different matrices are needed to calibrate the efficiency of column chemistry and caution should be taken when dealing with samples with complicated matrices. Well-calibrated standards with matrix elements matching samples should be used to reduce the interlaboratory mass bias.

H2O-CH4 and H2S-CH4 complexes: a direct comparison through molecular beam experiments and ab initio calculations

New molecular beam scattering experiments have been performed to measure the total ( elastic plus inelastic) cross sections as a function of the velocity in collisions between water and hydrogen sulfide projectile molecules and the methane target. Measured data have been exploited to characterize the range and strength of the intermolecular interaction in such systems, which are of relevance as they drive the gas phase molecular dynamics and the clathrate formation. Complementary information has been obtained by rotational spectra, recorded for the hydrogen sulfide-methane complex, with a pulsed nozzle Fourier transform microwave spectrometer. Extensive ab initio calculations have been performed to rationalize all the experimental findings. The combination of experimental and theoretical information has established the ground for the understanding of the nature of the interaction and allows for its basic components to be modelled, including charge transfer, in these weakly bound systems. The intermolecular potential for H2S-CH4 is significantly less anisotropic than for H2O-CH4, although both of them have potential minima that can be characterized as `hydrogen bonded'.

Use of polydispersity index as control parameter to study melting/freezing of Lennard-Jones system: Comparison among predictions of bifurcation theory with Lindemann criterion, inherent structure analysis and Hansen-Verlet rule

Using polydispersity index as an additional order parameter we investigate freezing/melting transition of Lennard-Jones polydisperse systems (with Gaussian polydispersity in size), especially to gain insight into the origin of the terminal polydispersity. The average inherent structure (IS) energy and root mean square displacement (RMSD) of the solid before melting both exhibit quite similar polydispersity dependence including a discontinuity at solid-liquid transition point. Lindemann ratio, obtained from RMSD, is found to be dependent on temperature. At a given number density, there exists a value of polydispersity index (delta (P)) above which no crystalline solid is stable. This transition value of polydispersity(termed as transition polydispersity, delta (P) ) is found to depend strongly on temperature, a feature missed in hard sphere model systems. Additionally, for a particular temperature when number density is increased, delta (P) shifts to higher values. This temperature and number density dependent value of delta (P) saturates surprisingly to a value which is found to be nearly the same for all temperatures, known as terminal polydispersity (delta (TP)). This value (delta (TP) similar to 0.11) is in excellent agreement with the experimental value of 0.12, but differs from hard sphere transition where this limiting value is only 0.048. Terminal polydispersity (delta (TP)) thus has a quasiuniversal character. Interestingly, the bifurcation diagram obtained from non-linear integral equation theories of freezing seems to provide an explanation of the existence of unique terminal polydispersity in polydisperse systems. Global bond orientational order parameter is calculated to obtain further insights into mechanism for melting.

Comparative sequence analyses of genome and transcriptome reveal novel transcripts and variants in the Asian elephant Elephas maximus

The Asian elephant Elephas maximus and the African elephant Loxodonta africana that diverged 5-7 million years ago exhibit differences in their physiology, behaviour and morphology. A comparative genomics approach would be useful and necessary for evolutionary and functional genetic studies of elephants. We performed sequencing of E. maximus and map to L. africana at similar to 15X coverage. Through comparative sequence analyses, we have identified Asian elephant specific homozygous, non-synonymous single nucleotide variants (SNVs) that map to 1514 protein coding genes, many of which are involved in olfaction. We also present the first report of a high-coverage transcriptome sequence in E. maximus from peripheral blood lymphocytes. We have identified 103 novel protein coding transcripts and 66-long non-coding (lnc)RNAs. We also report the presence of 181 protein domains unique to elephants when compared to other Afrotheria species. Each of these findings can be further investigated to gain a better understanding of functional differences unique to elephant species, as well as those unique to elephantids in comparison with other mammals. This work therefore provides a valuable resource to explore the immense research potential of comparative analyses of transcriptome and genome sequences in the Asian elephant.

A comparison of the fine-scale structure of the diurnal cycle of tropical rain and lightning

In this study, the fine-scale structure of the diurnal variability of ground-based lightning is systematically compared with satellite-based rain. At the outset, it is shown that tropical variability of lightning exhibits a prominent diurnal mode, much like rain. A comparison of the geographical distribution of the timing of the diurnal maximum shows that there is very good agreement between the two observables over continental and coastal regions throughout the tropics. Following this global tropical comparison, we focus on two regions, Borneo and equatorial South America, both of which show the interplay between oceanward and landward propagations of the phase of the diurnal maximum. Over Borneo, both rain and lightning clearly show a climatological cycle of ``breathing in'' (afternoon to early morning) and ``breathing out'' (morning to early afternoon). Over the equatorial east coast of South America, landward propagation is noticed in rain and lightning from early afternoon to early morning. Along the Pacific coast of South America, both rain and lightning show oceanward propagation. Though qualitatively consistent, over both regions the propagation is seen to extend further in rainfall. Additionally, given that lightning highlights vigorous convection, the timing of its diurnal maximum often precedes that of rainfall in the convective life cycle. (C) 2015 Elsevier B.V. All rights reserved.

Inter-comparison and performance evaluation of chemistry transport models over Indian region

Aerosol loading over the South Asian region has the potential to affect the monsoon rainfall, Himalayan glaciers and regional air-quality, with implications for the billions in this region. While field campaigns and network observations provide primary data, they tend to be location/season specific. Numerical models are useful to regionalize such location-specific data. Studies have shown that numerical models underestimate the aerosol scenario over the Indian region, mainly due to shortcomings related to meteorology and the emission inventories used. In this context, we have evaluated the performance of two such chemistry-transport models: WRF-Chem and SPRINTARS over an India-centric domain. The models differ in many aspects including physical domain, horizontal resolution, meteorological forcing and so on etc. Despite these differences, both the models simulated similar spatial patterns of Black Carbon (BC) mass concentration, (with a spatial correlation of 0.9 with each other), and a reasonable estimates of its concentration, though both of them under-estimated vis-a-vis the observations. While the emissions are lower (higher) in SPRINTARS (WRF-Chem), overestimation of wind parameters in WRF-Chem caused the concentration to be similar in both models. Additionally, we quantified the under-estimations of anthropogenic BC emissions in the inventories used these two models and three other widely used emission inventories. Our analysis indicates that all these emission inventories underestimate the emissions of BC over India by a factor that ranges from 1.5 to 2.9. We have also studied the model simulations of aerosol optical depth over the Indian region. The models differ significantly in simulations of AOD, with WRF-Chem having a better agreement with satellite observations of AOD as far as the spatial pattern is concerned. It is important to note that in addition to BC, dust can also contribute significantly to AOD. The models differ in simulations of the spatial pattern of mineral dust over the Indian region. We find that both meteorological forcing and emission formulation contribute to these differences. Since AOD is column integrated parameter, description of vertical profiles in both models, especially since elevated aerosol layers are often observed over Indian region, could be also a contributing factor. Additionally, differences in the prescription of the optical properties of BC between the models appear to affect the AOD simulations. We also compared simulation of sea-salt concentration in the two models and found that WRF-Chem underestimated its concentration vis-a-vis SPRINTARS. The differences in near-surface oceanic wind speeds appear to be the main source of this difference. In-spite of these differences, we note that there are similarities in their simulation of spatial patterns of various aerosol species (with each other and with observations) and hence models could be valuable tools for aerosol-related studies over the Indian region. Better estimation of emission inventories could improve aerosol-related simulations. (C) 2015 Elsevier Ltd. All rights reserved.

Photometric metallicity map of the Large Magellanic Cloud

We have estimated a metallicity map of the Large Magellanic Cloud (LMC) using the Magellanic Cloud Photometric Survey (MCPS) and Optical Gravitational Lensing Experiment (OGLE III) photometric data. This is a first of its kind map of metallicity up to a radius of 4 degrees-5 degrees, derived using photometric data and calibrated using spectroscopic data of Red Giant Branch (RGB) stars. We identify the RGB in the V, (V - I) colour-magnitude diagrams of small subregions of varying sizes in both data sets. We use the slope of the RGB as an indicator of the average metallicity of a subregion, and calibrate the RGB slope to metallicity using spectroscopic data for field and cluster red giants in selected subregions. The average metallicity of the LMC is found to be Fe/H] = -0.37 dex (sigmaFe/H] = 0.12) from MCPS data, and Fe/H] = -0.39 dex (sigmaFe/H] = 0.10) from OGLE III data. The bar is found to be the most metal-rich region of the LMC. Both the data sets suggest a shallow radial metallicity gradient up to a radius of 4 kpc (-0.049 +/- 0.002 dex kpc(-1) to -0.066 +/- 0.006 dex kpc(-1)). Subregions in which the mean metallicity differs from the surrounding areas do not appear to correlate with previously known features; spectroscopic studies are required in order to assess their physical significance.

Comparison of Convective and Radiative Heating Modes on the Thermophysical Changes of a Cerium Nitrate Droplet

In this paper, we try to establish the equivalence or similarity in the thermal and physiochemical changes in precursor droplets (cerium nitrate) in convective and radiative fields. The radiative field is created through careful heating of the droplet using a monochromatic light source (CO2 laser). The equivalence is also established for different modes of convection like droplet injected into a high-speed flow and droplet experiencing a convective flow due to acoustic streaming (levitated) only. The thermophysical changes are studied in an aqueous cerium nitrate droplet, and the dissociation of cerium nitrate to ceria is modeled using modified Kramers' reaction rate formulation. It is observed that vaporization, species accumulation, and chemical characteristics obtained in a convectively heated droplet are retained in a radiatively heated droplet by careful adjustment of the laser intensity. The timescales and ceria yield match reasonably well for both the cases. It is also noted that similar conclusions are drawn in both levitated droplet and a nonlevitated droplet.

Spatio-temporal Map Formation Based on a Potential Function

We revisit the problem of temporal self organization using activity diffusion based on the neural gas (NGAS) algorithm. Using a potential function formulation motivated by a spatio-temporal metric, we derive an adaptation rule for dynamic vector quantization of data. Simulations results show that our algorithm learns the input distribution and time correlation much faster compared to the static neural gas method over the same data sequence under similar training conditions.

Comparison of Matching Pursuit Algorithm with Other Signal Processing Techniques for Computation of the Time-Frequency Power Spectrum of Brain Signals

Signals recorded from the brain often show rhythmic patterns at different frequencies, which are tightly coupled to the external stimuli as well as the internal state of the subject. In addition, these signals have very transient structures related to spiking or sudden onset of a stimulus, which have durations not exceeding tens of milliseconds. Further, brain signals are highly nonstationary because both behavioral state and external stimuli can change on a short time scale. It is therefore essential to study brain signals using techniques that can represent both rhythmic and transient components of the signal, something not always possible using standard signal processing techniques such as short time fourier transform, multitaper method, wavelet transform, or Hilbert transform. In this review, we describe a multiscale decomposition technique based on an over-complete dictionary called matching pursuit (MP), and show that it is able to capture both a sharp stimulus-onset transient and a sustained gamma rhythm in local field potential recorded from the primary visual cortex. We compare the performance of MP with other techniques and discuss its advantages and limitations. Data and codes for generating all time-frequency power spectra are provided.

Three-dimensional plotted hydroxyapatite scaffolds with predefined architecture: comparison of stabilization by alginate cross-linking versus sintering

Scaffolds for bone tissue engineering are essentially characterized by porous three-dimensional structures with interconnected pores to facilitate the exchange of nutrients and removal of waste products from cells, thereby promoting cell proliferation in such engineered scaffolds. Although hydroxyapatite is widely being considered for bone tissue engineering applications due to its occurrence in the natural extracellular matrix of this tissue, limited reports are available on additive manufacturing of hydroxyapatite-based materials. In this perspective, hydroxyapatite-based three-dimensional porous scaffolds with two different binders (maltodextrin and sodium alginate) were fabricated using the extrusion method of three-dimensional plotting and the results were compared in reference to the structural properties of scaffolds processed via chemical stabilization and sintering routes, respectively. With the optimal processing conditions regarding to pH and viscosity of binder-loaded hydroxyapatite pastes, scaffolds with parallelepiped porous architecture having up to 74% porosity were fabricated. Interestingly, sintering of the as-plotted hydroxyapatite-sodium alginate (cross-linked with CaCl2 solution) scaffolds led to the formation of chlorapatite (Ca9.54P5.98O23.8Cl1.60(OH)(2.74)). Both the sintered scaffolds displayed progressive deformation and delayed fracture under compressive loading, with hydroxyapatite-alginate scaffolds exhibiting a higher compressive strength (9.5 +/- 0.5MPa) than hydroxyapatite-maltodextrin scaffolds (7.0 +/- 0.6MPa). The difference in properties is explained in terms of the phase assemblage and microstructure.

Visualisation of charge-transfer excitations in donor-acceptor molecules using the particle-hole map: a case study

Charge-transfer (CT) excitations are essential for photovoltaic phenomena in organic solar cells. Owing to the complexity of molecular geometries and orbital coupling, a detailed analysis and spatial visualisation of CT processes can be challenging. In this paper, a new detail-oriented visualisation scheme, the particle-hole map (PHM), is applied and explained for the purpose of spatial analysis of excitations in organic molecules. The PHM can be obtained from the output of a time-dependent density-functional theory calculation with negligible additional computational cost, and provides a useful physical picture for understanding the origins and destinations of electrons and holes during an excitation process. As an example, we consider intramolecular CT excitations in Diketopyrrolopyrrole-based molecules, and relate our findings to experimental results.