Two-Phase Thermodynamic Model for Efficient and Accurate Absolute Entropy of Water from Molecular Dynamics Simulations

Presented here is the two-phase thermodynamic (2PT) model for the calculation of energy and entropy of molecular fluids from the trajectory of molecular dynamics (MD) simulations. In this method, the density of state (DoS) functions (including the normal modes of translation, rotation, and intramolecular vibration motions) are determined from the Fourier transform of the corresponding velocity autocorrelation functions. A fluidicity parameter (f), extracted from the thermodynamic state of the system derived from the same MD, is used to partition the translation and rotation modes into a diffusive, gas-like component (with 3Nf degrees of freedom) and a nondiffusive, solid-like component. The thermodynamic properties, including the absolute value of entropy, are then obtained by applying quantum statistics to the solid component and applying hard sphere/rigid rotor thermodynamics to the gas component. The 2PT method produces exact thermodynamic properties of the system in two limiting states: the nondiffusive solid state (where the fluidicity is zero) and the ideal gas state (where the fluidicity becomes unity). We examine the 2PT entropy for various water models (F3C, SPC, SPC/E, TIP3P, and TIP4P-Ew) at ambient conditions and find good agreement with literature results obtained based on other simulation techniques. We also validate the entropy of water in the liquid and vapor phases along the vapor-liquid equilibrium curve from the triple point to the critical point. We show that this method produces converged liquid phase entropy in tens of picoseconds, making it an efficient means for extracting thermodynamic properties from MD simulations.

Diversity of DNA methyltransferases that recognize asymmetric target sequences

DNA methyltransferases (MTases) are a group of enzymes that catalyze the methyl group transfer from S-adenosyl-L-methionine in a sequence-specific manner. Orthodox Type II DNA MTases usually recognize palindromic DNA sequences and add a methyl group to the target base (either adenine or cytosine) on both strands. However, there are a number of MTases that recognize asymmetric target sequences and differ in their subunit organization. In a bacterial cell, after each round of replication, the substrate for any MTase is hemimethylated DNA, and it therefore needs only a single methylation event to restore the fully methylated state. This is in consistent with the fact that most of the DNA MTases studied exist as monomers in solution. Multiple lines of evidence suggest that some DNA MTases function as dimers. Further, functional analysis of many restriction-modification systems showed the presence of more than one or fused MTase genes. It was proposed that presence of two MTases responsible for the recognition and methylation of asymmetric sequences would protect the nascent strands generated during DNA replication from cognate restriction endonuclease. In this review, MTases recognizing asymmetric sequences have been grouped into different subgroups based on their unique properties. Detailed characterization of these unusual MTases would help in better understanding of their specific biological roles and mechanisms of action. The rapid progress made by the genome sequencing of bacteria and archaea may accelerate the identification and study of species- and strain-specific MTases of host-adapted bacteria and their roles in pathogenic mechanisms.

Study of magnetic field effect on the specific heat and entropy in DyBa2Cu3O7-x

The change in the specific heat by the application of magnetic field up to 161 for high temperature superconductor system for DyBa2Cu3O7-x by Revaz et al. [23] is examined through the phenomenological Ginzburg-Landau(G-L) theory of anisotropic Type-II superconductors. The observed specific heat anomaly near T-c with magnetic field is explained qualitatively through the expression <Delta C > = (B-a/T-c) t/(1 - t)(alpha Theta(gamma)lambda(2)(m)(0)), which is the anisotropic formulation of the G-L theory in the London limit developed by Kogan and coworkers; relating to the change in specific heat Delta C for the variation of applied magnetic field for different orientations with c-axis. The analysis of this equation explains satisfactorily the specific heat anomaly near T-c and determines the anisotropic ratio gamma as 5.608, which is close to the experimental value 5.3 +/- 0.5given in the paper of Revaz et al. for this system. (C) 2010 Elsevier B.V. All rights reserved.

A Low ML-Decoding Complexity, High Coding Gain, Full-Rate, Full-Diversity STBC for 4 x 2 MIMO System

This paper proposes a full-rate, full-diversity space-time block code(STBC) with low maximum likelihood (ML) decoding complexity and high coding gain for the 4 transmit antenna, 2 receive antenna (4 x 2) multiple-input multiple-output (MIMO) system that employs 4/16-QAM. For such a system, the best code known is the DjABBA code and recently, Biglieri, Hong and Viterbo have proposed another STBC (BHV code) for 4-QAM which has lower ML-decoding complexity than the DjABBA code but does not have full-diversity like the DjABBA code. The code proposed in this paper has the same ML-decoding complexity as the BHV code for any square M-QAM but has full-diversity for 4- and 16-QAM. Compared with the DjABBA code, the proposed code has lower ML-decoding complexity for square M-QAM constellation, higher coding gain for 4- and 16-QAM, and hence a better codeword error rate (CER) performance. Simulation results confirming this are presented.

Entropy and dynamics of water in hydration layers of a bilayer

We compute the entropy and transport properties of water in the hydration layer of dipalmitoylphosphatidylcholine bilayer by using a recently developed theoretical scheme two-phase thermodynamic model, termed as 2PT method; S.-T. Lin et al., J. Chem. Phys. 119, 11792 (2003)] based on the translational and rotational velocity autocorrelation functions and their power spectra. The weights of translational and rotational power spectra shift from higher to lower frequency as one goes from the bilayer interface to the bulk. Water molecules near the bilayer head groups have substantially lower entropy (48.36 J/mol/K) than water molecules in the intermediate region (51.36 J/mol/K), which have again lower entropy than the molecules (60.52 J/mol/K) in bulk. Thus, the entropic contribution to the free energy change (T Delta S) of transferring an interface water molecule to the bulk is 3.65 kJ/mol and of transferring intermediate water to the bulk is 2.75 kJ/mol at 300 K, which is to be compared with 6.03 kJ/mol for melting of ice at 273 K. The translational diffusion of water in the vicinity of the head groups is found to be in a subdiffusive regime and the rotational diffusion constant increases going away from the interface. This behavior is supported by the slower reorientational relaxation of the dipole vector and OH bond vector of interfacial water. The ratio of reorientational relaxation time for Legendre polynomials of order 1 and 2 is approximately 2 for interface, intermediate, and bulk water, indicating the presence of jump dynamics in these water molecules. (C) 2010 American Institute of Physics. doi:10.1063/1.3494115]

Maximum-entropy calculation of the electron density at 4 A resolution of Pf1 filamentous bacteriophage

A 4 A electron-density map of Pf1 filamentous bacterial virus has been calculated from x-ray fiber diffraction data by using the maximum-entropy method. This method produces a map that is free of features due to noise in the data and enables incomplete isomorphous-derivative phase information to be supplemented by information about the nature of the solution. The map shows gently curved (banana-shaped) rods of density about 70 A long, oriented roughly parallel to the virion axis but slewing by about 1/6th turn while running from a radius of 28 A to one of 13 A. Within these rods, there is a helical periodicity with a pitch of 5 to 6 A. We interpret these rods to be the helical subunits of the virion. The position of strongly diffracted intensity on the x-ray fiber pattern shows that the basic helix of the virion is right handed and that neighboring nearly parallel protein helices cross one another in an unusual negative sense.

DNA Binding, Coprotease, and Strand Exchange Activities of Mycobacterial RecA Proteins: Implications for Functional Diversity among RecA Nucleoprotein Filaments

One of the fundamental questions concerning homologous recombination is how RecA or its homologues recognize several DNA sequences with high affinity and catalyze all the diverse biological activities. In this study, we show that the extent of single-stranded DNA binding and strand exchange (SE) promoted by mycobacterial RecA proteins with DNA substrates having various degrees of GC content was comparable with that observed for Escherichia coli RecA. However, the rate and extent of SE promoted by these recombinases showed a strong negative correlation with increasing amounts of sequence divergence embedded at random across the length of the donor strand. Conversely, a positive correlation was seen between SE efficiency and the degree of sequence divergence in the recipient duplex DNA. The extent of heteroduplex formation was not significantly affected when both the pairing partners contained various degrees of sequence divergence, although there was a moderate decrease in the case of mycobacterial RecA proteins with substrates containing larger amounts of sequence divergence. Whereas a high GC content had no discernible effect on E. coli RecA coprotease activity, a negative correlation was apparent between mycobacterial RecA proteins and GC content. We further show clear differences in the extent of SE promoted by E. coli and mycobacterial RecA proteins in the presence of a wide range of ATP:ADP ratios. Taken together, our findings disclose the existence of functional diversity among E. coli and mycobacterial RecA nucleoprotein filaments, and the milieu of sequence divergence (i.e., in the donor or recipient) exerts differential effects on heteroduplex formation, which has implications for the emergence of new genetic variants.

Diversity in Functional Organization of Class I and Class II Biotin Protein Ligase

The cell envelope of Mycobacterium tuberculosis (M. tuberculosis) is composed of a variety of lipids including mycolic acids, sulpholipids, lipoarabinomannans, etc., which impart rigidity crucial for its survival and pathogenesis. Acyl CoA carboxylase (ACC) provides malonyl-CoA and methylmalonyl-CoA, committed precursors for fatty acid and essential for mycolic acid synthesis respectively. Biotin Protein Ligase (BPL/BirA) activates apo-biotin carboxyl carrier protein (BCCP) by biotinylating it to an active holo-BCCP. A minimal peptide (Schatz), an efficient substrate for Escherichia coli BirA, failed to serve as substrate for M. tuberculosis Biotin Protein Ligase (MtBPL). MtBPL specifically biotinylates homologous BCCP domain, MtBCCP87, but not EcBCCP87. This is a unique feature of MtBPL as EcBirA lacks such a stringent substrate specificity. This feature is also reflected in the lack of self/promiscuous biotinylation by MtBPL. The N-terminus/HTH domain of EcBirA has the selfbiotinable lysine residue that is inhibited in the presence of Schatz peptide, a peptide designed to act as a universal acceptor for EcBirA. This suggests that when biotin is limiting, EcBirA preferentially catalyzes, biotinylation of BCCP over selfbiotinylation. R118G mutant of EcBirA showed enhanced self and promiscuous biotinylation but its homologue, R69A MtBPL did not exhibit these properties. The catalytic domain of MtBPL was characterized further by limited proteolysis. Holo-MtBPL is protected from proteolysis by biotinyl-59 AMP, an intermediate of MtBPL catalyzed reaction. In contrast, apo-MtBPL is completely digested by trypsin within 20 min of co-incubation. Substrate selectivity and inability to promote self biotinylation are exquisite features of MtBPL and are a consequence of the unique molecular mechanism of an enzyme adapted for the high turnover of fatty acid biosynthesis.

A Low-Complexity, Full-Rate, Full-Diversity 2 x 2 STBC with Golden Code's Coding Gain

This paper presents a low-ML-decoding-complexity, full-rate, full-diversity space-time block code (STBC) for a 2 transmit antenna, 2 receive antenna multiple-input multiple-output (MIMO) system, with coding gain equal to that of the best and well known Golden code for any QAM constellation. Recently, two codes have been proposed (by Paredes, Gershman and Alkhansari and by Sezginer and Sari), which enjoy a lower decoding complexity relative to the Golden code, but have lesser coding gain. The 2 x 2 STBC presented in this paper has lesser decoding complexity for non-square QAM constellations, compared with that of the Golden code, while having the same decoding complexity for square QAM constellations. Compared with the Paredes-Gershman-Alkhansari and Sezginer-Sari codes, the proposed code has the same decoding complexity for non-rectangular QAM constellations. Simulation results, which compare the codeword error rate (CER) performance, are presented.

On the improvement of diversity-multiplexing gain tradeoff in a training based TDD-simo system

This paper investigates the diversity-multiplexing gain tradeoff (DMT) of a time-division duplex (TDD) single-input multiple-output (SIMO) system with perfect channel state information (CSI) at the receiver (CSIR) and partial CSI at the transmitter (CSIT). The partial CSIT is acquired through a training sequence from the receiver to the transmitter. The training sequence is chosen in an intelligent manner based on the CSIR, to reduce the training length by a factor of r, the number of receive antennas. We show that, for the proposed training scheme and a given channel coherence time, the diversity order increases linearly with r for nonzero multiplexing gain. This is a significant improvement over conventional orthogonal training schemes.

Molecular interpretation of the linear relationship between the entropy and the enthalpy of activation of charge transfer reactions in polar liquids

The well-known linear relationship (T?S# =??H# +?, where 1 >? > 0,? > 0) between the entropy (?S#) and the enthalpy (?H#) of activation for reactions in polar liquids is investigated by using a molecular theory. An explicit derivation of this linear relation from first principles is presented for an outersphere charge transfer reaction. The derivation offers microscopic interpretation for the quantities? and?. It has also been possible to make connection with and justify the arguments of Bell put forward many years ago.

Demographic history of India and mtDNA-sequence diversity

The demographic history of India was examined by comparing mtDNA sequences obtained from members of three culturally divergent Indian subpopulations (endogamous caste groups). While an inferred tree revealed some clustering according to caste affiliation, there was no clear separation into three genetically distinct groups along caste lines. Comparison of pairwise nucleotide difference distributions, however, did indicate a difference in growth patterns between two of the castes. The Brahmin population appears to have undergone either a rapid expansion or steady growth. The low-ranking Mukri caste, however, may have either maintained a roughly constant population size or undergone multiple bottlenecks during that period. Comparison of the Indian sequences to those obtained from other populations, using a tree, revealed that the Indian sequences, along with ah other non-African samples, form a starlike cluster. This cluster may represent a major expansion, possibly originating in southern Asia, taking place at some point after modern humans initially left Africa.

Parasite abundance and diversity in mammals: correlates with host ecology

Fecally dispersed parasites of 12 wild mammal species in Mudumalai Sanctuary, southern India, were studied, Fecal propagule densities and parasite diversity measures were correlated with host ecological variables. Host species with higher predatory pressure had lower parasite loads and parasite diversity. Host body weight, home range, population density, gregariousness, and diet did not show predicted effects on parasite loads. Measures of a! diversity were positively correlated with parasite abundance and were negatively correlated with beta diversity, Based on these data, hypotheses regarding determinants of parasite community are discussed.