Stochastic dynamics modeling of the protein sequence length distribution in genomes: implications for microbial evolution

In this paper, we report an analysis of the protein sequence length distribution for 13 bacteria, four archaea and one eukaryote whose genomes have been completely sequenced, The frequency distribution of protein sequence length for all the 18 organisms are remarkably similar, independent of genome size and can be described in terms of a lognormal probability distribution function. A simple stochastic model based on multiplicative processes has been proposed to explain the sequence length distribution. The stochastic model supports the random-origin hypothesis of protein sequences in genomes. Distributions of large proteins deviate from the overall lognormal behavior. Their cumulative distribution follows a power-law analogous to Pareto's law used to describe the income distribution of the wealthy. The protein sequence length distribution in genomes of organisms has important implications for microbial evolution and applications. (C) 1999 Elsevier Science B.V. All rights reserved.

Structural characterisation of a uracil containing hairpin DNA by NMR and molecular dynamics

Three-dimensional (3D) structure of a hairpin DNA d-CTAGAGGATCCTTTUGGATCCT (22mer; abbreviated as U4-hairpin), which has a uracil nucleotide unit at the fourth position from the 5' end of the tetra-loop has been solved by NMR spectroscopy. The H-1 resonances of this hairpin have been assigned almost completely. NMR restrained molecular dynamics and energy minimisation procedures have been used to describe the 3D structure of the U4 hairpin. This study establishes that the stem of the hairpin adopts a right handed B-DNA conformation while the T-12 and U-15 nucleotide stack upon 3' and 5' ends of the stem, respectively. Further, T-14 stacks upon both T-12 and U-15 while T-13 partially stacks upon T-14. Very weak stacking interaction is observed between T-13 and T-12. All the individual nucleotide bases adopt 'anti' conformation with respect to their sugar moiety. The turning phosphate in the loop is located between T-13 and T-14. The stereochemistry of U-15 mimics the situation where uracil would stack in a B-DNA conformation. This could be the reason as to why the U4-hairpin is found to be the best substrate for its interaction with uracil DNA glycosylase (UDG) compared to the other substrates in which the uracil is at the first, second and third positions of the tetra-loop from its 5' end, as reported previously.

Subquadratic quantum number dependence and other anomalies of vibrational dephasing in liquid nitrogen: Molecular dynamics simulation study from the triple point to the critical point and beyond

Vibrational phase relaxation near gas-liquid and liquid-solid phase coexistence has been studied by molecular dynamics simulations of N-N stretch in N-2. Experimentally observed pronounced insensitivity of phase relaxation from the triple point to beyond the boiling point is found to originate from a competition between density relaxation and resonant-energy transfer terms. The sharp rise in relaxation rate near the critical point (CP) can be attributed at least partly to the sharp, rise in vibration-rotation coupling contribution. Substantial subquadratic quantum number dependence of overtone dephasing rate is found near the CP and in supercritical fluids. [S0031-9007 (99)09318-7].

Algal Photosynthetic Dynamics in Urban Lakes under Stress Conditions

Urban lakes form vital ecosystems supporting livelihood with social, economic and aesthetic benefits that are essential for quality life. This depends on the biotic and abiotic components in an ecosystem. The structure of an ecosystem forms a decisive factor in sustaining its functional abilities which include nutrient cycling, oxygen production, etc. A community assemblage of primary producers (algae) plays a crucial role in maintaining the balance as they form the base of energy pyramid in the ecosystem. Algae assimilate carbon in the environment via photosynthetic activities and releases oxygen for the next level of biotic elements in an ecosystem. Besides these, algal cells rich in protein serve as food and feed, used as manure and for production of biofuels. Understanding algal photosynthetic dynamics helps in assessing the level of dissolved oxygen (DO), food (fish, etc.), waste assimilation, etc. Algal chlorophyll content, algal biomass, primary productivity and algal photosynthetic quotient are some of the parameters that help in assessing the status of urban lakes. Chlorophyll content gives a measure of the growth, spread and quantity of algae. Unplanned rapid urbanization in Bangalore in recent times has resulted in either disappearance of lake ecosystems or deteriorated the lake water quality impairing the ecological processes. This paper computes algal growth, community structure, primary productivity and composition for three major lakes (T G Halli, Bellandur and Varthur lakes) under contrast levels of anthropogenic influences.

The molecular dynamics study of lithium ion conduction in phosphate glasses and the role of non-bridging oxygen

Molecular dynamics (MD) simulation of lithium phosphate (Li2O-P2O5) glasses with varying Li2O content has been carried out. Two different P-O distances corresponding to phosphorus coordination with bridging oxygen (BO) and non-bridging oxygen (NBO) were identified in the simulated glasses. NBO-BO interconversion or bond switching was noted, which results in a dynamic equilibration of the tetrahedral phosphate units (P-n, n = 1,3 indicates the number of bridging oxygen atoms in the coordination of phosphorus). The NBO-BO bond switching is mildly activated with an effective activation barrier of 0.03-0.05 eV. Lithium ion jumps do not appear to be strongly coupled to bond switching. But the number of Li+ ions coordinated to an optimum number of NBOs and the number of Li+ ions jumping out of their sites appear to be correlated. Detailed analysis was made of the dynamics of P-n species and new insights have been obtained regarding ion migration in network-modified phosphate glasses.

Biological water: Femtosecond dynamics of macromolecular hydration

The unique features of a macromolecule and water as a solvent make the issue of solvation unconventional, with questions about the static versus dynamic nature of hydration and the, physics of orientational and translational diffusion at the boundary. For proteins, the hydration shell that covers the surface is critical to the stability of its structure and function. Dynamically speaking, the residence time of water at the surface is a signature of its mobility and binding. With femtosecond time resolution it is possible to unravel the shortest residence times which are key for the description of the hydration layer, static or dynamic. In this article we review these issues guided by experimental studies, from this laboratory, of polar hydration dynamics at the surfaces of two proteins (Subtilisin Carlsberg (SC) and Monellin). The natural probe tryptophan amino acid was used for the interrogation of the dynamics, and for direct comparison we also studied the behavior in bulk water - a complete hydration in 1 ps. We develop a theoretical description of solvation and relate the theory to the experimental observations. In this - theoretical approach, we consider the dynamical equilibrium in the hydration shell, defining the rate processes for breaking and making the transient hydrogen bonds, and the effective friction in the layer which is defined by the translational and orientational motions of water molecules. The relationship between the residence time of water molecules and the observed slow component in solvation dynamics is a direct one. For the two proteins studied, we observed a "bimodal decay" for the hydration correlation function, with two primary relaxation times: ultrafast, typically 1 ps or less, and longer, typically 15-40 ps, and both are related to the residence time at the protein surface, depending on the binding energies. We end by making extensions to studies of the denatured state of the protein, random coils, and the biomimetic micelles, and conclude with our thoughts on the relevance of the dynamics of native structures to their functions.

Comparison of the ultrafast to slow time scale dynamics of three liquid crystals in the isotropic phase

The dynamics of three liquid crystals, 4'(pentyloxy)-4-biphenylcarbonitrile (5-OCB), 4'-pentyl-4-biphenylcarbonitrile (5-CB), and 1-isothiocyanato-(4-propylcyclohexyl)benzene (3-CHBT), are investigated from very short time (similar to1 ps) to very long time (>100 ns) as a function of temperature using optical heterodyne detected optical Kerr effect experiments. For all three liquid crystals, the data decay exponentially only on the longest time scale (> several ns). The temperature dependence of the long time scale exponential decays is described well by the Landau-de Gennes theory of the randomization of pseudonematic domains that exist in the isotropic phase of liquid crystals near the isotropic to nematic phase transition. At short time, all three liquid crystals display power law decays. Over the full range of times, the data for all three liquid crystals are fit with a model function that contains a short time power law. The power law exponents for the three liquid crystals range between 0.63 and 0.76, and the power law exponents are temperature independent over a wide range of temperatures. Integration of the fitting function gives the empirical polarizability-polarizability (orientational) correlation function. A preliminary theoretical treatment of collective motions yields a correlation function that indicates that the data can decay as a power law at short times. The power law component of the decay reflects intradomain dynamics. (C) 2002 American Institute of Physics.

Topological phases, Majorana modes and quench dynamics in a spin ladder system

We explore the salient features of the `Kitaev ladder', a two-legged ladder version of the spin-1/2 Kitaev model on a honeycomb lattice, by mapping it to a one-dimensional fermionic p-wave superconducting system. We examine the connections between spin phases and topologically non-trivial phases of non-interacting fermionic systems, demonstrating the equivalence between the spontaneous breaking of global Z(2) symmetry in spin systems and the existence of isolated Majorana modes. In the Kitaev ladder, we investigate topological properties of the system in different sectors characterized by the presence or absence of a vortex in each plaquette of the ladder. We show that vortex patterns can yield a rich parameter space for tuning into topologically non-trivial phases. We introduce and employ a new topological invariant for explicitly determining the presence of zero energy Majorana modes at the boundaries of such phases. Finally, we discuss dynamic quenching between topologically non-trivial phases in the Kitaev ladder and, in particular, the post-quench dynamics governed by tuning through a quantum critical point.

DNA binding properties of novel dansylated distamycin analogues in which the fluorophore is directly conjugated to the N-methyl-pyrrole carboxamide backbone

Polyamides that are structural analogues of the naturally occurring DNA minor groove binding antibiotic distamycin (Dst) are promising candidates as gene modulators. Developing strategies for the large scale screening and monitoring of the cellular distribution of such ligands would aid the faster discovery of molecules, which would have eventual utility in molecular biology and medicine. Attachment of fluorescent tags would be a useful step towards this end. A fundamental question in this connection is whether the tag modifies the DNA binding affinity of the parent compounds. Towards answering this question, we have developed two oligopeptides that bear the dansyl (N, N-dimethylaminonaphthalene sulfonamido fluorophore) coupled directly to the N-terminus of the conjugated N-methylpyrrole carboxamide network, and possess three or four N-methyl pyrrole carboxamide units (abbreviated as Dn3 and Dn4 respectively). DNA binding abilities of these molecules were assessed from fluorescence titration experiments, duplex-DNA T-m analysis (employing both UV and fluorescence spectroscopy), induced circular dichroism measurements (ICD), salt dependence of ICD and apparent binding constant measurements (K-app) employing ethidium bromide (EtBr) displacement assay. Both these molecules 'reported' DNA binding in the form of an enhanced fluorescence emission. As judged from the ICD measurements, salt dependence of ICD, T-m analysis and K-app measurements, the binding affinities of the molecules that possessed dansyl group at their N-termini were lower than the ones with equivalent number of amide units, but possessed N-methylpyrrole carboxamide unit at their N-termini. These results would have implications in the future design of fluorescent polyamides.

Dynamics of crystal size distributions with size-dependent rates

The growth and dissolution dynamics of nonequilibrium crystal size distributions (CSDs) can be determined by solving the governing population balance equations (PBEs) representing reversible addition or dissociation. New PBEs are considered that intrinsically incorporate growth dispersion and yield complete CSDs. We present two approaches to solving the PBEs, a moment method and a numerical scheme. The results of the numerical scheme agree with the moment technique, which can be solved exactly when powers on mass-dependent growth and dissolution rate coefficients are either zero or one. The numerical scheme is more general and can be applied when the powers of the rate coefficients are non-integers or greater than unity. The influence of the size dependent rates on the time variation of the CSDs indicates that as equilibrium is approached, the CSDs become narrow when the exponent on the growth rate is less than the exponent on the dissolution rate. If the exponent on the growth rate is greater than the exponent on the dissolution rate, then the polydispersity continues to broaden. The computation method applies for crystals large enough that interfacial stability issues, such as ripening, can be neglected. (C) 2002 Elsevier Science B.V. All rights reserved.

Essential dynamics and sidechain hydrogen bond cluster studies on eosinophil cationic protein

Eosinophil Cationic Protein (ECP) is a member of RNase A superfamily which carries out the obligatory catalytic role of cleaving RNA. It is involved in a variety of biological functions. Molecular dynamics simulations followed by essential dynamics analysis on this protein are carried out with the goal of gaining insights into the dynamical properties at atomic level. The top essential modes contribute to subspaces and to the transition phase. Further, the sidechain-sidechain/sidechain-mainchain hydrogen bond clusters are analyzed in the top modes, and compared with those of crystal structure. The role of residues identified by these methods is discussed in the context of concerted motion, structure and stability of the protein.

Excited state structure and dynamics of p-benzoquinone and bromanil from time-resolved resonance Raman spectra and simulation

p-Benzoquinone and its halogen substituted derivatives are known to have differing reactivities in the triplet excited state. While bromanil catalyzes the reduction of octaethylporphyrin most efficiently among the halogenated p-benzoquinones, the reaction does not take place in presence of the unsubstituted p-benzoquinone (T. Nakano and Y. Mori, Bull. Chem. Soc. Jpn., 67, 2627 (1994)). Understanding of such differences requires a detailed knowledge of the triplet state structures, normal mode compositions and excited state dynamics. In this paper, we apply a recently presented scheme (M. Puranik, S. Umapathy, J. G. Snijders, and J. Chandrasekhar, J. Chem, Phys., 115, 6106 (2001)) that combines parameters from experiment and computation in a wave packet dynamics simulation to the triplet states of p-benzoquinone and bromanil. The absorption and resonance Raman spectra of both the molecules have been simulated. The normal mode compositions and mode specific excited state displacements have been presented and compared. Time-dependent evolution of the absorption and Raman overlaps for all the observed modes has been discussed in detail. In p-benzoquinone, the initial dynamics is along the C=C stretching and C-H bending modes whereas in bromanil nearly equal displacements are observed along all the stretching coordinates.

Diperoxovanadate can substitute for H(2)O(2) at much lower concentration in inducing features of premature cellular senescence in mouse fibroblasts (NIH3T3)

Stress induced premature senescence (SIPS) in mammalian cells is an accelerated ageing response and experimentally obtained on treatment of cells with high concentrations of H(2)O(2), albeit at sub-lethal doses, because H(2)O(2) gets depleted by abundant cellular catalase. In the present study diperoxovanadate (DPV) was used as it is known to be stable at physiological pH, to be catalase-resistant and to substitute for H(2)O(2) in its activities at concentrations order of magnitudes lower. On treating NIH3T3 cells with DPV, SIPS-like morphology was observed along with an immediate response of rounding of the cells by disruption of actin cytoskeleton and transient G2/M arrest. DPV could bring about growth arrest and senescence associated features at 25 mu M dose, which were not seen with similar doses of either H(2)O(2) or vanadate. A minimal dose of 150 mu M of H(2)O(2) was required to induce similar affects as 25 mu M DPV. Increase in senescent associated markers such as p21, HMGA2 and PAI-1 was more prominent in DPV treated cells compared to similar dose of H(2)O(2). DPV-treated cells showed marked relocalization of Cyclin D1 from nucleus to cytoplasm. These results indicate that DPV, stable inorganic peroxide, is more efficient in inducing SIPS at lower concentrations compared to H(2)O(2). (C) 2011 Elsevier Ireland Ltd. All rights reserved.

Hydrogen-bond dynamics near a micellar surface: Origin of the universal slow relaxation at complex aqueous interfaces

The dynamics of hydrogen bonds among water molecules themselves and with the polar head groups (PHG) at a micellar surface have been investigated by long molecular dynamics simulations. The lifetime of the hydrogen bond between a PHG and a water molecule is found to be much longer than that between any two water molecules, and is likely to be a general feature of hydrophilic surfaces of organized assemblies. Analyses of individual water trajectories suggest that water molecules can remain bound to the micellar surface for more than 100 ps. The activation energy for such a transition from the bound to a free state for the water molecules is estimated to be about 3.5 kcal/mol.