Dynamics of ribonuclease A and ribonuclease S: Computational and experimental studies

RNase S is a complex consisting of two proteolytic fragments of RNase A: the S peptide (residues 1-20) and S protein (residues 21-124). RNase S and RNase A have very similar X-ray structures and enzymatic activities. previous experiments have shown increased rates of hydrogen exchange and greater sensitivity to tryptic cleavage for RNase S relative to RNase A. It has therefore been asserted that the RNase S complex is considerably more dynamically flexible than RNase A. In the present study we examine the differences in the dynamics of RNaseS and RNase A computationally, by MD simulations, and experimentally, using trypsin cleavage as a probe of dynamics. The fluctuations around the average solution structure during the simulation were analyzed by measuring the RMS deviation in coordinates. No significant differences between RNase S and RNase A dynamics were observed in the simulations. We were able to account for the apparent discrepancy between simulation and experiment by a simple model, According to this model, the experimentally observed differences in dynamics can be quantitatively explained by the small amounts of free S peptide and S protein that are present in equilibrium with the RNase S complex. Thus, folded RNase A and the RNase S complex have identical dynamic behavior, despite the presence of a break in polypeptide chain between residues 20 and 21 in the latter molecule. This is in contrast to what has been widely believed for over 30 years about this important fragment complementation system.

Crossover from Ising to mean-field critical behavior in an aqueous electrolyte solution

The near-critical behavior of the susceptibility deduced from light-scattering measurements in a ternary liquid mixture of 3-methylpyridine, water, and sodium bromide has been determined. The measurements have been performed in the one-phase region near the lower consolute points of samples with different concentrations of sodium bromide. A crossover from Ising asymptotic behavior to mean-field behavior has been observed. As the concentration of sodium bromide increases, the crossover becomes more pronounced, and the crossover temperature shifts closer to the critical temperature. The data are well described by a model that contains two independent crossover parameters. The crossover of the susceptibility critical exponent γ from its Ising value γ=1.24 to the mean-field value γ=1 is sharp and nonmonotonic. We conclude that there exists an additional length scale in the system due to the presence of the electrolyte which competes with the correlation length of the concentration fluctuations. An analogy with crossover phenomena in polymer solutions and a possible connection with multicritical phenomena is discussed.

Intermittent Dynamics, Stochastic Resonance and Dynamical Heterogeneity in Supercooled Liquid Water

Here we find through computer simulations and theoretical analysis that the low temperature thermodynamic anomalies of liquid water arises from the intermittent fluctuation between its high density and low density forms, consisting largely of 5-coordinated and 4-coordinated water molecules, respectively. The fluctuations exhibit strong dynamic heterogeneity (defined by the four point time correlation function), accompanied by a divergence like growth of the dynamic correlation length, of the type encountered in fragile supercooled liquids. The intermittency has been explained by invoking a two state model often employed to understand stochastic resonance, with the relevant periodic perturbation provided here by the fluctuation of the total volume of the system.

Precipitation In Small Systems--II. Mean Field Equations More Effective Than Population Balance

Part I (Manjunath et al., 1994, Chem. Engng Sci. 49, 1451-1463) of this paper showed that the random particle numbers and size distributions in precipitation processes in very small drops obtained by stochastic simulation techniques deviate substantially from the predictions of conventional population balance. The foregoing problem is considered in this paper in terms of a mean field approximation obtained by applying a first-order closure to an unclosed set of mean field equations presented in Part I. The mean field approximation consists of two mutually coupled partial differential equations featuring (i) the probability distribution for residual supersaturation and (ii) the mean number density of particles for each size and supersaturation from which all average properties and fluctuations can be calculated. The mean field equations have been solved by finite difference methods for (i) crystallization and (ii) precipitation of a metal hydroxide both occurring in a single drop of specified initial supersaturation. The results for the average number of particles, average residual supersaturation, the average size distribution, and fluctuations about the average values have been compared with those obtained by stochastic simulation techniques and by population balance. This comparison shows that the mean field predictions are substantially superior to those of population balance as judged by the close proximity of results from the former to those from stochastic simulations. The agreement is excellent for broad initial supersaturations at short times but deteriorates progressively at larger times. For steep initial supersaturation distributions, predictions of the mean field theory are not satisfactory thus calling for higher-order approximations. The merit of the mean field approximation over stochastic simulation lies in its potential to reduce expensive computation times involved in simulation. More effective computational techniques could not only enhance this advantage of the mean field approximation but also make it possible to use higher-order approximations eliminating the constraints under which the stochastic dynamics of the process can be predicted accurately.

Axially homogeneous, zero mean flow buoyancy-driven turbulence in a vertical pipe

We report an experimental study of a new type of turbulent flow that is driven purely by buoyancy. The flow is due to an unstable density difference, created using brine and water, across the ends of a long (length/diameter = 9) vertical pipe. The Schmidt number Sc is 670, and the Rayleigh number (Ra) based on the density gradient and diameter is about 10(8). Under these conditions the convection is turbulent, and the time-averaged velocity at any point is `zero'. The Reynolds number based on the Taylor microscale, Re-lambda, is about 65. The pipe is long enough for there to be an axially homogeneous region, with a linear density gradient, about 6-7 diameters long in the midlength of the pipe. In the absence of a mean flow and, therefore, mean shear, turbulence is sustained just by buoyancy. The flow can be thus considered to be an axially homogeneous turbulent natural convection driven by a constant (unstable) density gradient. We characterize the flow using flow visualization and particle image velocimetry (PIV). Measurements show that the mean velocities and the Reynolds shear stresses are zero across the cross-section; the root mean squared (r.m.s.) of the vertical velocity is larger than those of the lateral velocities (by about one and half times at the pipe axis). We identify some features of the turbulent flow using velocity correlation maps and the probability density functions of velocities and velocity differences. The flow away from the wall, affected mainly by buoyancy, consists of vertically moving fluid masses continually colliding and interacting, while the flow near the wall appears similar to that in wall-bound shear-free turbulence. The turbulence is anisotropic, with the anisotropy increasing to large values as the wall is approached. A mixing length model with the diameter of the pipe as the length scale predicts well the scalings for velocity fluctuations and the flux. This model implies that the Nusselt number would scale as (RaSc1/2)-Sc-1/2, and the Reynolds number would scale as (RaSc-1/2)-Sc-1/2. The velocity and the flux measurements appear to be consistent with the Ra-1/2 scaling, although it must be pointed out that the Rayleigh number range was less than 10. The Schmidt number was not varied to check the Sc scaling. The fluxes and the Reynolds numbers obtained in the present configuration are Much higher compared to what would be obtained in Rayleigh-Benard (R-B) convection for similar density differences.

Response And Stability Of A Stochastic Beam-Column Using Stochastic Fem

A beam-column resting on continuous Winkler foundation and discrete elastic supports is considered. The beam-column is of variable cross-section and the variation of sectional properties along the axis of the beam-column is deterministic. Young's modulus, mass per unit length and distributed axial loadings of the beam-column have a stochastic distribution. The foundation stiffness coefficient of the Winkler model, the stiffnesses of discrete elastic supports, stiffnesses of end springs and the end thrust, are all considered as random parameters. The material property fluctuations and distributed axial loadings are considered to constitute independent, one-dimension uni-variate homogeneous real stochastic fields in space. The foundation stiffness coefficient, stiffnesses of the discrete elastic supports, stiffnesses of end springs and the end thrust are considered to constitute independent random variables. Static response, free vibration and stability behaviour of the beam-column are studied. Hamilton's principle is used to formulate the problem using stochastic FEM. Sensitivity vectors of the response and stability parameters are evaluated. Using these statistics of free vibration frequencies, mode shapes, buckling parameters, etc., are evaluated. A numerical example is given.

Ion conductance in electrolyte solutions

We develop a new theoretical formulation to study ion conductance in electrolyte solutions, based on a mode coupling theory treatment of the electrolyte friction. The new theory provides expressions for both the ion atmosphere relaxation and electrophoretic contributions to the total electrolyte friction that acts on a moving ion. While the ion atmosphere relaxation term arises from the time-dependent microscopic interaction of the moving ion with the surrounding ions in the solution, the electrophoretic term originates from the coupling of the ion's velocity to the collective current mode of the ion atmosphere. Mode coupling theory, combined with time-dependent density functional theory of ion atmosphere fluctuations, leads to self-consistent expressions for these two terms which also include the effects of self-motion of the ion under consideration. These expressions have been solved for the concentration dependence of electrolyte friction and ion conductance. It is shown that in the limit of very low ion concentration, the present theory correctly reduces to the well-known Debye-Huckel-Onsager limiting law which predicts a linear dependence of conductance on the square root of ion concentration (c). At moderate and high concentrations, the present theory predicts a significant nonlinear and weaker dependence on root c which is in very good agreement with experimental results. The present theory is self-contained and does not involve any adjustable parameter.

Resistance Noise in Electrically Biased Bilayer Graphene

We demonstrate that the low-frequency resistance fluctuations, or noise, in bilayer graphene are strongly connected to its band structure and display a minimum when the gap between the conduction and valence band is zero. Using double-gated bilayer graphene devices we have tuned the zero gap and charge neutrality points independently, which offers a versatile mechanism to investigate the low-energy band structure, charge localization, and screening properties of bilayer graphene.

LISA as a dark energy probe

Recently, it has been shown that the inclusion of higher signal harmonics in the inspiral signals of binary supermassive black holes (SMBH) leads to dramatic improvements in the parameter estimation with Laser Interferometer Space Antenna (LISA). In particular, the angular resolution becomes good enough to identify the host galaxy or galaxy cluster, in which case the redshift can be determined by electromagnetic means. The gravitational wave signal also provides the luminosity distance with high accuracy, and the relationship between this and the redshift depends sensitively on the cosmological parameters, such as the equation-of-state parameter w = p(DE)/rho(DE) of dark energy. Using binary SMBH events at z < 1 with appropriate masses and orientations, one would be able to constrain w to within a few per cent. We show that, if the measured sky location is folded into the error analysis, the uncertainty on w goes down by an additional factor of 2-3, leaving weak lensing as the only limiting factor in using LISA as a dark energy probe.

Structural mobility and transformations in globular proteins

Although globular proteins are endowed with well defined three-dimensional structures, they exhibit substantial mobility within the framework of the given threedimensional structure. The different types of mobility found in proteins by and large correspond to the different levels of organisational hierarchy in protein architecture. They are of considerable structural and functional significance, and can be broadly classified into(a) thermal and conformational fluctuations, (b) segmental mobility, (c) interdomain mobility and (d) intersubunit mobility. Protein crystallographic studies has provided a wealth of information on all of them. The temperature factors derived from X-ray diffraction studies provide a measure of atomic displacements caused by thermal and conformational fluctuations. The variation of displacement along the polypeptide chain have provided functionally significant information on the flexibility of different regions of the molecule in proteins such as myoglobin, lysozyme and prealbumin. Segmental mobility often involves the movement of a region or a segment of a molecule with respect to the rest, as in the transition between the apo and the holo structures of lactate dehydrogenase. It may also involve rigidification of a disordered region of the molecule as in the activation of the zymogens of serine proteases. Transitions between the apo and the holo structures of alcohol dehydrogenase,and between the free and the sugar bound forms of hexokinase, are good examples of interdomain mobility caused by hinge-bending. The capability of different domains to move semi-independently contributes greatly to the versatility of immunoglobulin molecules. Interdomain mobility in citrate synthase appears to be more complex and its study has led to an alternative description of domain closure. The classical and the most thoroughly studied case of intersubunit mobility is that in haemoglobin. The stereochemical mechanism of the action of this allosteric protein clearly brings out the functional subtilities that could be achieved through intersubunit movements. In addition to ligand binding and activation,environmental changes also often cause structural transformations. The reversible transformation between 2 Zn insulin and 4 Zn insulin is caused by changes in the ionic strength of the medium. Adenylate Kinase provides a good example for functionally significant reversible conformational transitions induced by variation in pH. Available evidences indicate that reversible structural transformations in proteins could also be caused by changes in the aqueous environment, including those in the amount of water surrounding protein molecules.

Transient state work fluctuation theorem for a classical harmonic oscillator linearly coupled to a harmonic bath

Based on a Hamiltonian description we present a rigorous derivation of the transient state work fluctuation theorem and the Jarzynski equality for a classical harmonic oscillator linearly coupled to a harmonic heat bath, which is dragged by an external agent. Coupling with the bath makes the dynamics dissipative. Since we do not assume anything about the spectral nature of the harmonic bath the derivation is not restricted only to the Ohmic bath, rather it is more general, for a non-Ohmic bath. We also derive expressions of the average work done and the variance of the work done in terms of the two-time correlation function of the fluctuations of the position of the harmonic oscillator. In the case of an Ohmic bath, we use these relations to evaluate the average work done and the variance of the work done analytically and verify the transient state work fluctuation theorem quantitatively. Actually these relations have far-reaching consequences. They can be used to numerically evaluate the average work done and the variance of the work done in the case of a non-Ohmic bath when analytical evaluation is not possible.

The Milky Way's Dark Matter Halo Appears To Be Lopsided

The atomic hydrogen gas (H I) disk in the outer region (beyond similar to 10 kpc from the center) of Milky Way can provide valuable information about the structure of the dark matter halo. The recent three-dimensional thickness map of the outer H I disk from the all sky 21 cm line Leiden/Argentine/Bonn survey, gives us a unique opportunity to investigate the structure of the dark matter halo of Milky Way in great detail. A striking feature of this new survey is the north-south (N-S) asymmetry in the thickness map of the atomic hydrogen gas. Assuming vertical hydrostatic equilibrium under the total potential of the Galaxy, we derive the model thickness map of the H I gas. We show that simple axisymmetric halo models, such as softened isothermal halo (producing a flat rotation curve with V-c similar to 220 km s(-1)) or any halo with density falling faster than the isothermal one, are not able to explain the observed radial variation of the gas thickness. We also show that such axisymmetric halos along with different H I velocity dispersion in the two halves, cannot explain the observed asymmetry in the thickness map. Amongst the nonaxisymmetric models, it is shown that a purely lopsided (m = 1, first harmonic) dark matter halo with reasonable H I velocity dispersion fails to explain the N-S asymmetry satisfactorily. However, we show that by superposing a second harmonic (m = 2) out of phase onto a purely lopsided halo, e. g., our best fit and more acceptable model A (with parameters epsilon(1)(h) = 0.2, epsilon(2)(h) = 0.18, and sigma(H I) = 8.5 km s(-1)) can provide an excellent fit to the observation and reproduce the N-S asymmetry naturally. The emerging picture of the asymmetric dark matter halo is supported by the. cold dark matter halos formed in the cosmological N-body simulation.

Contribution from normal and starburst galaxies to the extragalactic gamma-ray background (EGRB)

The extragalactic diffuse emission at gamma-ray energies has interesting cosmological implications since these photons suffer little or no attenuation during their propagation from the site of origin. The emission could originate from either truly diffuse processes or from unresolved point sources such as AGNs, normal galaxies and starburst galaxies. Here, we examine the unresolved point source origin of the extragalactic gamma-ray background emission from normal galaxies and starburst galaxies. gamma-ray emission from normal galaxies is primarily coming from cosmic-ray interactions with interstellar matter and radiation (similar to 90%) along with a small contribution from discrete point sources (similar to 10%). Starburst galaxies are expected to have enhanced supernovae activity which leads to higher cosmic-ray densities, making starburst galaxies sufficiently luminous at gamma-ray energies to be detected by the current gamma-ray mission(Fermi Gamma-ray Space Telescope).

New temperature fluctuation method for direct determination of thermal activation energy of deep levels in semiconductors

A new method is suggested where the thermal activation energy is measured directly and not as a slope of an Arrhenius plot. The sample temperature T is allowed to fluctuate about a temperature T0. The reverse-biased sample diode is repeatedly pulsed towards zero bias and the transient capacitance C1 at time t1 is measured The activation energy is obtained by monitoring the fluctuations in C1 and T. The method has been used to measure the activation energy of the gold acceptor level in silicon.

Alamethicin and related membrane channel forming polypeptides

Alamethicin and several related microbial polypeptides, which contain a high proportion of agr-aminoisobutyric acid (Aib) residues, possess the ability to modify the permeability properties of phospholipid bilayer membranes. Alamethicin induces excitability phenomena in model membranes and has served as an excellent model for the study of voltage sensitive transmembrane channels. This review summarizes various aspects of the structural chemistry and membrane modifying properties of alamethicin and related Alb containing peptides. The presence of Aib residues in these sequences, constrains the polypeptides to 310 or agr-helical conformations. Functional membrane channels are formed by aggregation of cylindrical peptide helices, which span the lipid bilayer, forming a scaffolding for an aqueous column across the membrane. After consideration of the available data on the conductance characteristics of alamethicin channels, a working, hypothesis for a channel model is outlined. Channel aggregates in the lipid phase may be stabilized by intermolecular hydrogen bonding, involving a central glutamine residue and also by interactions between the macro-dipoles of proximate peptide helices. Fluctuations between different conductance states are rationalized by transitions between states of different aggregation and hence altered dimensions of the aqueous core or by changes in net dipole moment of the aggregate. Ion fluxes through the channel may also be affected by the electric field within the aqueous core.