Correlation between local structural dynamics of proteins inferred from NMR ensembles and evolutionary dynamics of homologues of known structure

Conformational changes in proteins are extremely important for their biochemical functions. Correlation between inherent conformational variations in a protein and conformational differences in its homologues of known structure is still unclear. In this study, we have used a structural alphabet called Protein Blocks (PBs). PBs are used to perform abstraction of protein 3-D structures into a 1-D strings of 16 alphabets (a-p) based on dihedral angles of overlapping pentapeptides. We have analyzed the variations in local conformations in terms of PBs represented in the ensembles of 801 protein structures determined using NMR spectroscopy. In the analysis of concatenated data over all the residues in all the NMR ensembles, we observe that the overall nature of inherent local structural variations in NMR ensembles is similar to the nature of local structural differences in homologous proteins with a high correlation coefficient of .94. High correlation at the alignment positions corresponding to helical and beta-sheet regions is only expected. However, the correlation coefficient by considering only the loop regions is also quite high (.91). Surprisingly, segregated position-wise analysis shows that this high correlation does not hold true to loop regions at the structurally equivalent positions in NMR ensembles and their homologues of known structure. This suggests that the general nature of local structural changes is unique; however most of the local structural variations in loop regions of NMR ensembles do not correlate to their local structural differences at structurally equivalent positions in homologues.

Correlation between thermodynamic anomalies and pathways of ice nucleation in supercooled water

The well-known classical nucleation theory (CNT) for the free energy barrier towards formation of a nucleus of critical size of the new stable phase within the parent metastable phase fails to take into account the influence of other metastable phases having density/order intermediate between the parent metastable phase and the final stable phase. This lacuna can be more serious than capillary approximation or spherical shape assumption made in CNT. This issue is particularly significant in ice nucleation because liquid water shows rich phase diagram consisting of two (high and low density) liquid phases in supercooled state. The explanations of thermodynamic and dynamic anomalies of supercooled water often invoke the possible influence of a liquid-liquid transition between two metastable liquid phases. To investigate both the role of thermodynamic anomalies and presence of distinct metastable liquid phases in supercooled water on ice nucleation, we employ density functional theoretical approach to find nucleation free energy barrier in different regions of phase diagram. The theory makes a number of striking predictions, such as a dramatic lowering of nucleation barrier due to presence of a metastable intermediate phase and crossover in the dependence of free energy barrier on temperature near liquid-liquid critical point. These predictions can be tested by computer simulations as well as by controlled experiments. (C) 2014 AIP Publishing LLC.

Energy funneling and macromolecular conformational dynamics: a 2D Raman correlation study of PEG melting

Cascading energy landscapes through funneling has been postulated as a mechanistic route for achieving the lowest energy configuration of a macromolecular system (such as proteins and polymers). In particular, understanding the molecular mechanism for the melting and crystallization of polymers is a challenging fundamental question. The structural modifications that lead to the melting of poly(ethylene glycol) (PEG) are investigated here. Specific Raman bands corresponding to different configurations of the PEG chain have been identified, and the molecular structural dynamics of PEG melting have been addressed using a combination of Raman spectroscopy, 2D Raman correlation and density functional theory (DFT) calculations. The melting dynamics of PEG have been unambiguously explained along the C-O bond rotation coordinate.

Solubility-Hardness Correlation in Molecular Crystals: Curcumin and Sulfathiazole Polymorphs

Curcumin and sulfathiazole exist as three and five polymorphs, respectively. We correlate solubility and mechanical properties in these polymorphic systems. It is seen that hardness (H) is inversely proportional to the solubility of a polymorph. H of the polymorphs is explained on the basis of slip planes in the crystal structure, the Schmid factor (m), and the relative orientation of molecules with respect to the nanoindenter direction. Effectively, H is a useful parameter (compared to melting point, T-m, and density, rho) that correlates well with the solubility of a polymorph. Such a correlation is helpful in systems like curcumin and sulfathiazole in which the Gibbs free energy of the polymorphs are close to one another. To summarize, a softer polymorph is more soluble.

Correlation between Optical Properties and Nanomorphology of Fluoranthene-Based Conjugated Copolymer

Nanoparticles of conjugated polymers are receiving attention due to their interesting optical properties. Here we report nanoparticles of fluoranthene-based conjugated copolymer prepared by the Suzuki coupling reaction. The copolymer forms nanoparticles by the spontaneous self-assembly after evaporation of organic solvent. The mean diameter of the nanoparticles can be manipulated by varying solvent composition. We investigated the parameters that govern the nanostructured morphology of polymer by systematic variation of good and poor solvent. The UV vis and time-resolved fluorescence spectroscopy measurement reveal the use of poor solvent in the organization of nanostructures. Furthermore, transmission electron microscopy highlights the importance of rigidity of the polymer backbone in morphological development.

Diffusing-wave spectroscopy in an inhomogeneous object: Local viscoelastic spectra from ultrasound-assisted measurement of correlation decay arising from the ultrasound focal volume

We demonstrate diffusing-wave spectroscopy (DWS) in a localized region of a viscoelastically inhomogeneous object by measurement of the intensity autocorrelation g(2)(tau)] that captures only the decay introduced by the temperature-induced Brownian motion in the region. The region is roughly specified by the focal volume of an ultrasound transducer which introduces region specific mechanical vibration owing to insonification. Essential characteristics of the localized non-Markovian dynamics are contained in the decay of the modulation depth M(tau)], introduced by the ultrasound forcing in the focal volume selected, on g(2)(tau). The modulation depth M(tau(i)) at any delay time tau(i) can be measured by short-time Fourier transform of g(2)(tau) and measurement of the magnitude of the spectrum at the ultrasound drive frequency. By following the established theoretical framework of DWS, we are able to connect the decay in M(tau) to the mean-squared displacement (MSD) of scattering centers and the MSD to G*(omega), the complex viscoelastic spectrum. A two-region composite polyvinyl alcohol phantom with different viscoelastic properties is selected for demonstrating local DWS-based recovery of G*(omega) corresponding to these regions from the measured region specific M(tau(i))vs tau(i). The ultrasound-assisted measurement of MSD is verified by simulating, using a generalized Langevin equation (GLE), the dynamics of the particles in the region selected as well as by the usual DWS experiment without the ultrasound. It is shown that whereas the MSD obtained by solving the GLE without the ultrasound forcing agreed with its experimental counterpart covering small and large values of tau, the match was good only in the initial transients in regard to experimental measurements with ultrasound.

Nonlinear optical absorption in a graphene infrared photodetector

The photoresponse of the graphene photodetector elucidated strong dependence on several optical parameters, such as the angle of incidence and the incident power of infrared exposure at room temperature. The sinusoidal dependence of the photoresponse on incidence angle, which had not been realized before, has now been revealed. The combined effect of the photo excited charge carrier and the photon drag effect explain this nonlinear optical absorption in graphene at lower incident power. The nonlinear dependence of the charge carrier generation on the incident power revealed that this process contributed to the nonlinear photoresponse. However, a deviation is observed at a higher incident power due to the induction of thermal effects in the graphene lattice. This work demonstrates the tunability of the graphene photodetector under a systematic variation that involves both parameters.

Processing-microstructure-yield strength correlation in a near beta Ti alloy, Ti-5Al-5Mo-5V-3Cr

A combined set of thermo-mechanical steps recommended for high strength beta Ti alloy are homogenization, deformation, recrystallization, annealing and ageing steps in sequence. Recrystallization carried out above or below beta transus temperature generates either beta annealed (lath type morphology of alpha) or bimodal (lath+globular morphology of alpha) microstructure. Through variations in heat treatment parameters at these processing steps, wide ranges of length scales of features have been generated in both types of microstructures in a near beta Ti alloy, Ti-5Al-5Mo-5V-3Cr (Ti-5553). 0.2% Yield strength (YS) has been correlated to various microstructural features and associated heat treatment parameters. Relative importance of microstructural features in influencing YS has been identified. Process parameters at different steps have been identified and recommended for attaining different levels of YS for this near beta Ti alloy. (C) 2014 Elsevier B.V. All rights reserved.

Multimodal therapy for complete regression of malignant melanoma using constrained nonlinear optimal dynamic inversion

Using a realistic nonlinear mathematical model for melanoma dynamics and the technique of optimal dynamic inversion (exact feedback linearization with static optimization), a multimodal automatic drug dosage strategy is proposed in this paper for complete regression of melanoma cancer in humans. The proposed strategy computes different drug dosages and gives a nonlinear state feedback solution for driving the number of cancer cells to zero. However, it is observed that when tumor is regressed to certain value, then there is no need of external drug dosages as immune system and other therapeutic states are able to regress tumor at a sufficiently fast rate which is more than exponential rate. As model has three different drug dosages, after applying dynamic inversion philosophy, drug dosages can be selected in optimized manner without crossing their toxicity limits. The combination of drug dosages is decided by appropriately selecting the control design parameter values based on physical constraints. The process is automated for all possible combinations of the chemotherapy and immunotherapy drug dosages with preferential emphasis of having maximum possible variety of drug inputs at any given point of time. Simulation study with a standard patient model shows that tumor cells are regressed from 2 x 107 to order of 105 cells because of external drug dosages in 36.93 days. After this no external drug dosages are required as immune system and other therapeutic states are able to regress tumor at greater than exponential rate and hence, tumor goes to zero (less than 0.01) in 48.77 days and healthy immune system of the patient is restored. Study with different chemotherapy drug resistance value is also carried out. (C) 2014 Elsevier Ltd. All rights reserved.

Heteronuclear proton double quantum-carbon single quantum scalar correlation in solids

A new NMR experiment that exploits the advantages of proton double quantum (DQ) NMR through a proton DQ-carbon single quantum (SQ) correlation experiment in the solid state is proposed. Analogous to the previously proposed 2D H-1 (DQ)-C-13 refocused INEPT experiment (Webber et al., 2010), the correlation between H-1 and C-13 is achieved through scalar coupling evolution, while the double quantum coherence among protons is generated through dipolar couplings. However, the new experiment relies on C-13 transverse coherence for scalar transfer. The new experiment dubbed MAS-J-H-1 (DQ)-C-13-HMQC, is particularly suited for unlabeled molecules and can provide higher sensitivity than its INEPT counterpart. The experiment is applied to four different samples. (C) 2014 Elsevier Inc. All rights reserved.

Nonlinear preconditioning for efficient and accurate interface capturing in simulation of multicomponent compressible flows

Single fluid schemes that rely on an interface function for phase identification in multicomponent compressible flows are widely used to study hydrodynamic flow phenomena in several diverse applications. Simulations based on standard numerical implementation of these schemes suffer from an artificial increase in the width of the interface function owing to the numerical dissipation introduced by an upwind discretization of the governing equations. In addition, monotonicity requirements which ensure that the sharp interface function remains bounded at all times necessitate use of low-order accurate discretization strategies. This results in a significant reduction in accuracy along with a loss of intricate flow features. In this paper we develop a nonlinear transformation based interface capturing method which achieves superior accuracy without compromising the simplicity, computational efficiency and robustness of the original flow solver. A nonlinear map from the signed distance function to the sigmoid type interface function is used to effectively couple a standard single fluid shock and interface capturing scheme with a high-order accurate constrained level set reinitialization method in a way that allows for oscillation-free transport of the sharp material interface. Imposition of a maximum principle, which ensures that the multidimensional preconditioned interface capturing method does not produce new maxima or minima even in the extreme events of interface merger or breakup, allows for an explicit determination of the interface thickness in terms of the grid spacing. A narrow band method is formulated in order to localize computations pertinent to the preconditioned interface capturing method. Numerical tests in one dimension reveal a significant improvement in accuracy and convergence; in stark contrast to the conventional scheme, the proposed method retains its accuracy and convergence characteristics in a shifted reference frame. Results from the test cases in two dimensions show that the nonlinear transformation based interface capturing method outperforms both the conventional method and an interface capturing method without nonlinear transformation in resolving intricate flow features such as sheet jetting in the shock-induced cavity collapse. The ability of the proposed method in accounting for the gravitational and surface tension forces besides compressibility is demonstrated through a model fully three-dimensional problem concerning droplet splash and formation of a crownlike feature. (C) 2014 Elsevier Inc. All rights reserved.

Dynamic Site Characterization and Correlation of Shear Wave Velocity with Standard Penetration Test `N' Values for the City of Agartala, Tripura State, India

Seismic site characterization is the basic requirement for seismic microzonation and site response studies of an area. Site characterization helps to gauge the average dynamic properties of soil deposits and thus helps to evaluate the surface level response. This paper presents a seismic site characterization of Agartala city, the capital of Tripura state, in the northeast of India. Seismically, Agartala city is situated in the Bengal Basin zone which is classified as a highly active seismic zone, assigned by Indian seismic code BIS-1893, Indian Standard Criteria for Earthquake Resistant Design of Structures, Part-1 General Provisions and Buildings. According to the Bureau of Indian Standards, New Delhi (2002), it is the highest seismic level (zone-V) in the country. The city is very close to the Sylhet fault (Bangladesh) where two major earthquakes (M (w) > 7) have occurred in the past and affected severely this city and the whole of northeast India. In order to perform site response evaluation, a series of geophysical tests at 27 locations were conducted using the multichannel analysis of surface waves (MASW) technique, which is an advanced method for obtaining shear wave velocity (V (s)) profiles from in situ measurements. Similarly, standard penetration test (SPT-N) bore log data sets have been obtained from the Urban Development Department, Govt. of Tripura. In the collected data sets, out of 50 bore logs, 27 were selected which are close to the MASW test locations and used for further study. Both the data sets (V (s) profiles with depth and SPT-N bore log profiles) have been used to calculate the average shear wave velocity (V (s)30) and average SPT-N values for the upper 30 m depth of the subsurface soil profiles. These were used for site classification of the study area recommended by the National Earthquake Hazard Reduction Program (NEHRP) manual. The average V (s)30 and SPT-N classified the study area as seismic site class D and E categories, indicating that the city is susceptible to site effects and liquefaction. Further, the different data set combinations between V (s) and SPT-N (corrected and uncorrected) values have been used to develop site-specific correlation equations by statistical regression, as `V (s)' is a function of SPT-N value (corrected and uncorrected), considered with or without depth. However, after considering the data set pairs, a probabilistic approach has also been presented to develop a correlation using a quantile-quantile (Q-Q) plot. A comparison has also been made with the well known published correlations (for all soils) available in the literature. The present correlations closely agree with the other equations, but, comparatively, the correlation of shear wave velocity with the variation of depth and uncorrected SPT-N values provides a more suitable predicting model. Also the Q-Q plot agrees with all the other equations. In the absence of in situ measurements, the present correlations could be used to measure V (s) profiles of the study area for site response studies.

Nonlinear analysis of a thin pre-twisted and delaminated anisotropic strip

The present work is aimed at the development of an efficient mathematical model to assess the degradation in the stiffness properties of an anisotropic strip due to delamination. In particular, the motive is to capture those nonlinear effects in a strip that arise due to the geometry of the structure, in the presence of delamination. The variational asymptotic method (VAM) is used as a mathematical tool to simplify the original 3D problem to a 1D problem. Further simplification is achieved by modeling the delaminated structure by a sublaminate approach. By VAM, a 2D nonlinear sectional analysis is carried out to determine compact expression for the stiffness terms. The stiffness terms, both linear and nonlinear, are derived as functions of delamination length and location in closed form. In general, the results from the analysis include fully coupled nonlinear 1D stiffness coefficients, 3D strain field, 3D stress field, and in-plane and warping fields. In this work, the utility of the model is demonstrated for a static case, and its capability to capture the trapeze effect in the presence of delamination is investigated and compared with results available in the literature.

Regimes of nonlinear depletion and regularity in the 3D Navier-Stokes equations

The periodic 3D Navier-Stokes equations are analyzed in terms of dimensionless, scaled, L-2m-norms of vorticity D-m (1 <= m <= infinity). The first in this hierarchy, D-1, is the global enstrophy. Three regimes naturally occur in the D-1-D-m plane. Solutions in the first regime, which lie between two concave curves, are shown to be regular, owing to strong nonlinear depletion. Moreover, numerical experiments have suggested, so far, that all dynamics lie in this heavily depleted regime 1]; new numerical evidence for this is presented. Estimates for the dimension of a global attractor and a corresponding inertial range are given for this regime. However, two more regimes can theoretically exist. In the second, which lies between the upper concave curve and a line, the depletion is insufficient to regularize solutions, so no more than Leray's weak solutions exist. In the third, which lies above this line, solutions are regular, but correspond to extreme initial conditions. The paper ends with a discussion on the possibility of transition between these regimes.

Microstructure-hardness-fretting wear resistance correlation in ultrafine grained Cu-TiB2-Pb composites

The retention of the desired combination of mechanical/tribological properties in ultrafine grained materials presents important challenges in the field of bulk metallic composites. In order to address this aspect, the present work demonstrates how one can achieve a good combination of hardness and wear resistance in Cu-Pb-TiB2 composites, consolidated by spark plasma sintering at low temperatures ( < 500 degrees C). Transmission electron microscope (TEM) studies reveal ultrafine grains of Cu (100-400 nm) with coarser TiB2 particles (1-2 mu m) along with fine scale Pb dispersoid at triple junctions or at the grain boundaries of Cu. Importantly, a high hardness of around 2.2 GPa and relative density of close to 90% relative density (rho(theo)) have been achieved for Cu-15 wt% TiB2-10 wt% Pb composite. Such property theo, combination has never been reported for any Cu-based nanocomposite, by conventional processing route. In reference to the tribological performance, fretting wear tests were conducted on the sintered nanocomposites and a good combination of steady state COF (0.6-0.7) and wear rate (10-4 mm(3)/N m) were measured. An inverse relationship between wear rate and hardness was recorded and this commensurates well with Archard's relationship of abrasive wear. The formation of a wear-resistant delaminated tribolayer consisting of TiB2 particles and ultrafine oxide debris, (Cu, Fe, Ti)(x)O-y as confirmed from subsurface imaging using focused ion beam microscopy has been identified as the key factors for the low wear rate of these composites. (C) 2014 Elsevier B.V. All rights reserved.