Chain Length Effects on Helix-Hairpin Distribution in Short Peptides with Aib-(D)Ala and Aib-Aib Segments

The Aib-(D)Ala dipeptide segment has a tendency to form both type-I'/III' and type-I/III beta-turns. The occurrence of prime turns facilitates the formation of beta-hairpin conformations, while type-I/III turns can nucleate helix formation. The octapeptide Boc-Leu-Phe-Val-Aib-(D)Ala-Leu-Phe-Val-OMe (1) has been previously shown to form a beta-hairpin in the crystalline state and in solution. The effects of sequence truncation have been examined using the model peptides Boc-Phe-Val-Aib-Xxx-Leu-Phe-NHMe (2, 6), Boc-Val-Aib-Xxx-Leu-NHMe (3, 7), and Boc-Aib-Xxx-NHMe (4, 8), where Xxx = (D)Ala, Aib. For peptides with central Aib-Aib segments, Boc-Phe-Val-Aib-Aib-Leu-Phe-NHMe (6), Boc-Val-Aib-Aib-Leu-NHMe (7), and Boc-Aib-Aib-NHMe (8) helical conformations have been established by NMR studies in both hydrogen bonding (CD(3)OH) and non-hydrogen bonding (CDCl(3)) solvents. In contrast, the corresponding hexapeptide Boc-Phe-Val-Aib-(D)Ala-Leu-Phe-Val-NHMe (2) favors helical conformations in CDCl(3) and beta-hairpin conformations in CD(3)OH. The beta-turn conformations (type-I'/III) stabilized by intramolecular 4 -> 1 hydrogen bonds are observed for the peptide Boc-Aib-(D)Ala-NHMe (4) and Boc-Aib-Aib-NIiMe (8) in crystals. The tetrapeptide Boc-Val-Aib-Aib-Leu-NHMe (7) adopts an incipient 3(10)-helical conformation stabilized by three 4 -> 1 hydrogen bonds. The peptide Boc-Val-Aib-(D)Ala-Leu-NHMe (3) adopts a novel et-turn conformation, stabilized by three intramolecular hydrogen bonds (two 4 -> 1 and one 5 -> 1). The Aib-L(D)Ala segment adopts a type-I' beta-turn conformation. The observation of an NOE between Val (1) NH <-> HNCH(3) (5) in CD(3)OH suggests, that the solid state conformation is maintained in methanol solutions. (C) 2011 Wiley Periodicals, Inc. Biopolymers (Pept Sci) 96: 744-756, 2011.

Asparagine and glutamine differ in their propensities to form specific side chain-backbone hydrogen bonded motifs in proteins

Short range side chain-backbone hydrogen bonded motifs involving Asn and Gln residues have been identified from a data set of 1370 protein crystal structures (resolution = 1.5 angstrom). Hydrogen bonds involving residues i - 5 to i + 5 have been considered. Out of 12,901 Asn residues, 3403 residues (26.4%) participate in such interactions, while out of 10,934 Gln residues, 1780 Gln residues (16.3%) are involved in these motifs. Hydrogen bonded ring sizes (Cn, where n is the number of atoms involved), directionality and internal torsion angles are used to classify motifs. The occurrence of the various motifs in the contexts of protein structure is illustrated. Distinct differences are established between the nature of motifs formed by Asn and Gln residues. For Asn, the most highly populated motifs are the C10 (COdi .NHi + 2), C13 (COdi .NHi + 3) and C17 (NdHi .COi - 4) structures. In contrast, Gln predominantly forms C16 (COei .NHi - 3), C12 (NeHi .COi - 2), C15 (NeHi .COi - 3) and C18 (NeHi .COi - 4) motifs, with only the C18motif being analogous to the Asn C17structure. Specific conformational types are established for the Asn containing motifs, which mimic backbone beta-turns and a-turns. Histidine residues are shown to serve as a mimic for Asn residues in side chain-backbone hydrogen bonded ring motifs. Illustrative examples from protein structures are considered. Proteins 2012; (c) 2011 Wiley Periodicals, Inc.

A novel variable resolution global spectral method on the sphere

A variable resolution global spectral method is created on the sphere using High resolution Tropical Belt Transformation (HTBT). HTBT belongs to a class of map called reparametrisation maps. HTBT parametrisation of the sphere generates a clustering of points in the entire tropical belt; the density of the grid point distribution decreases smoothly in the domain outside the tropics. This variable resolution method creates finer resolution in the tropics and coarser resolution at the poles. The use of FFT procedure and Gaussian quadrature for the spectral computations retains the numerical efficiency available with the standard global spectral method. Accuracy of the method for meteorological computations are demonstrated by solving Helmholtz equation and non-divergent barotropic vorticity equation on the sphere. (C) 2011 Elsevier Inc. All rights reserved.

Hypoxia-Mediated Impairment of the Mitochondrial Respiratory Chain Inhibits the Bactericidal Activity of Macrophages

In infected tissues oxygen tensions are low. As innate immune cells have to operate under these conditions, we analyzed the ability of macrophages (M phi) to kill Escherichia coli or Staphylococcus aureus in a hypoxic microenvironment. Oxygen restriction did not promote intracellular bacterial growth but did impair the bactericidal activity of the host cells against both pathogens. This correlated with a decreased production of reactive oxygen intermediates (ROI) and reactive nitrogen intermediates. Experiments with phagocyte NADPH oxidase (PHOX) and inducible NO synthase (NOS2) double-deficient M phi revealed that in E. coli- or S. aureus-infected cells the reduced antibacterial activity during hypoxia was either entirely or partially independent of the diminished PHOX and NOS2 activity. Hypoxia impaired the mitochondrial activity of infected M phi. Inhibition of the mitochondrial respiratory chain activity during normoxia (using rotenone or antimycin A) completely or partially mimicked the defective antibacterial activity observed in hypoxic E. coli-or S. aureus-infected wild-type M phi, respectively. Accordingly, inhibition of the respiratory chain of S. aureus-infected, normoxic PHOX-/- NOS2(-/-) M phi further raised the bacterial burden of the cells, which reached the level measured in hypoxic PHOX-/- NOS2(-/-) M phi cultures. Our data demonstrate that the reduced killing of S. aureus or E. coli during hypoxia is not simply due to a lack of PHOX and NOS2 activity but partially or completely results from an impaired mitochondrial antibacterial effector function. Since pharmacological inhibition of the respiratory chain raised the generation of ROI but nevertheless phenocopied the effect of hypoxia, ROI can be excluded as the mechanism underlying the antimicrobial activity of mitochondria.

Study of Low Temperature Magnetic Properties of a Single Chain Magnet with Alternate Isotropic and Non-collinear Anisotropic Units

Here we study thermodynamic properties of an important class of single-chain magnets (SCMs), where alternate units are isotropic and anisotropic with anisotropy axes being non-collinear. This class of SCMs shows slow relaxation at low temperatures which results from the interplay of two different relaxation mechanisms, namely dynamical and thermal. Here anisotropy is assumed to be large and negative, as a result, anisotropic units behave like canted spins at low temperatures; but even then simple Ising-type model does not capture the essential physics of the system due to quantum mechanical nature of the isotropic units. We here show how statistical behavior of this class of SCMs can be studied using a transfer matrix (TM) method. We also, for the first time, discuss in detail how weak inter-chain interactions can be treated by a TM method. The finite size effect is also discussed which becomes important for low temperature dynamics. At the end of this paper, we apply this technique to study a real helical chain magnet.

Network properties of protein-decoy structures

Convergence of the vast sequence space of proteins into a highly restricted fold/conformational space suggests a simple yet unique underlying mechanism of protein folding that has been the subject of much debate in the last several decades. One of the major challenges related to the understanding of protein folding or in silico protein structure prediction is the discrimination of non-native structures/decoys from the native structure. Applications of knowledge-based potentials to attain this goal have been extensively reported in the literature. Also, scoring functions based on accessible surface area and amino acid neighbourhood considerations were used in discriminating the decoys from native structures. In this article, we have explored the potential of protein structure network (PSN) parameters to validate the native proteins against a large number of decoy structures generated by diverse methods. We are guided by two principles: (a) the PSNs capture the local properties from a global perspective and (b) inclusion of non-covalent interactions, at all-atom level, including the side-chain atoms, in the network construction accommodates the sequence dependent features. Several network parameters such as the size of the largest cluster, community size, clustering coefficient are evaluated and scored on the basis of the rank of the native structures and the Z-scores. The network analysis of decoy structures highlights the importance of the global properties contributing to the uniqueness of native structures. The analysis also exhibits that the network parameters can be used as metrics to identify the native structures and filter out non-native structures/decoys in a large number of data-sets; thus also has a potential to be used in the protein `structure prediction' problem.

Distribution of transverse chain fluctuations in harmonically confined semiflexible polymers

Two different experimental studies of polymer dynamics based on single-molecule fluorescence imaging have recently found evidence of heterogeneities in the widths of the putative tubes that surround filaments of F-actin during their motion in concentrated solution. In one J. Glaser, D. Chakraborty, K. Kroy, I. Lauter, M. Degawa, N. Kirchesner, B. Hoffmann, R. Merkel, and M. Giesen, Phys. Rev. Lett. 105, 037801 (2010)], the observations were explained in terms of the statistics of a worm-like chain confined to a potential determined self-consistently by a binary collision approximation, and in the other B. Wang, J. Guan, S. M. Anthony, S. C. Bae, K. S. Schweizer, and S. Granick, Phys. Rev. Lett. 104, 118301 (2010)], they were explained in terms of the scaling properties of a random fluid of thin rods. In this paper, we show, using an exact path integral calculation, that the distribution of the length-averaged transverse fluctuations of a harmonically confined weakly bendable rod (one possible realization of a semiflexible chain in a tube), is in good qualitative agreement with the experimental data, although it is qualitatively different in analytic structure from the earlier theoretical predictions. We also show that similar path integral techniques can be used to obtain an exact expression for the time correlation function of fluctuations in the tube cross section. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4712306]

Dynamics of end-to-end loop formation: A flexible chain in the presence of hydrodynamic interaction

Based on the Wilemski-Fixman approach G. Wilemski, M. Fixman, J. Chem. Phys. 60 (1974) 866], we show that, for a flexible chain in theta solvent, hydrodynamic interaction treated with a pre-averaging approximation makes ring closing faster if the chain is not very short. We also show that the ring closing time for a long chain with hydrodynamic interaction in theta solvent scales with the chain length (N) as N-1.5, in agreement with the previous renormalization group calculation based prediction by Freidman and O'Shaughnessy B. Friedman, B. O'Shaughnessy, Phys. Rev. A 40 (1989) 5950]. (C) 2012 Elsevier B.V. All rights reserved.

The Contribution of Rare Species to Community Phylogenetic Diversity across a Global Network of Forest Plots

Niche differentiation has been proposed as an explanation for rarity in species assemblages. To test this hypothesis requires quantifying the ecological similarity of species. This similarity can potentially be estimated by using phylogenetic relatedness. In this study, we predicted that if niche differentiation does explain the co-occurrence of rare and common species, then rare species should contribute greatly to the overall community phylogenetic diversity (PD), abundance will have phylogenetic signal, and common and rare species will be phylogenetically dissimilar. We tested these predictions by developing a novel method that integrates species rank abundance distributions with phylogenetic trees and trend analyses, to examine the relative contribution of individual species to the overall community PD. We then supplement this approach with analyses of phylogenetic signal in abundances and measures of phylogenetic similarity within and between rare and common species groups. We applied this analytical approach to 15 long-term temperate and tropical forest dynamics plots from around the world. We show that the niche differentiation hypothesis is supported in six of the nine gap-dominated forests but is rejected in the six disturbance-dominated and three gap-dominated forests. We also show that the three metrics utilized in this study each provide unique but corroborating information regarding the phylogenetic distribution of rarity in communities.

Confinement and viscoelastic effects on chain closure dynamics

Chemical reactions inside cells are typically subject to the effects both of the cell's confining surfaces and of the viscoelastic behavior of its contents. In this paper, we show how the outcome of one particular reaction of relevance to cellular biochemistry - the diffusion-limited cyclization of long chain polymers - is influenced by such confinement and crowding effects. More specifically, starting from the Rouse model of polymer dynamics, and invoking the Wilemski-Fixman approximation, we determine the scaling relationship between the mean closure time t(c) of a flexible chain (no excluded volume or hydrodynamic interactions) and the length N of its contour under the following separate conditions: (a) confinement of the chain to a sphere of radius d and (b) modulation of its dynamics by colored Gaussian noise. Among other results, we find that in case (a) when d is much smaller than the size of the chain, t(c) similar to Nd-2, and that in case (b), t(c) similar to N-2/(2 (2H)), H being a number between 1/2 and 1 that characterizes the decay of the noise correlations. H is not known a priori, but values of about 0.7 have been used in the successful characterization of protein conformational dynamics. At this value of H (selected for purposes of illustration), t(c) similar to N-3.4, the high scaling exponent reflecting the slow relaxation of the chain in a viscoelastic medium. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4729041]

ENTANGLEMENT IN A 3-SPIN HEISENBERG-XY CHAIN WITH NEAREST-NEIGHBOR INTERACTIONS, SIMULATED IN AN NMR QUANTUM SIMULATOR

The evolution of entanglement in a 3-spin chain with nearest-neighbor Heisenberg-XY interactions for different initial states is investigated here. In an NMR experimental implementation, we generate multipartite entangled states starting from initial separable pseudo-pure states by simulating nearest-neighbor XY interactions in a 3-spin linear chain of nuclear spin qubits. For simulating XY interactions, we follow algebraic method of Zhang et al. Phys. Rev. A 72 (2005) 012331]. Bell state between end qubits has been generated by using only the unitary evolution of the XY Hamiltonian. For generating W-state and GHZ-state a single qubit rotation is applied on second and all the three qubits, respectively after the unitary evolution of the XY Hamiltonian.

Retrobiosynthetic Approach Delineates the Biosynthetic Pathway and the Structure of the Acyl Chain of Mycobacterial Glycopeptidolipids

Glycopeptidolipids (GPLs) are dominant cell surface molecules present in several non-tuberculous and opportunistic mycobacterial species. GPLs from Mycobacterium smegmatis are composed of a lipopeptide core unit consisting of a modified C-26-C-34 fatty acyl chain that is linked to a tetrapeptide (Phe-Thr-Ala-alaninol). The hydroxyl groups of threonine and terminal alaninol are further modified by glycosylations. Although chemical structures have been reported for 16 GPLs from diverse mycobacteria, there is still ambiguity in identifying the exact position of the hydroxyl group on the fatty acyl chain. Moreover, the enzymes involved in the biosynthesis of the fatty acyl component are unknown. In this study we show that a bimodular polyketide synthase in conjunction with a fatty acyl-AMP ligase dictates the synthesis of fatty acyl chain of GPL. Based on genetic, biochemical, and structural investigations, we determine that the hydroxyl group is present at the C-5 position of the fatty acyl component. Our retrobiosynthetic approach has provided a means to understand the biosynthesis of GPLs and also resolve the long-standing debate on the accurate structure of mycobacterial GPLs.

Involvement of the Global Regulator H-NS in the Survival of Escherichia coli in Stationary Phase

Long-term batch cultures of Escherichia coli grown in nutrient-rich medium accumulate mutations that provide a growth advantage in the stationary phase (GASP). We have examined the survivors of prolonged stationary phase to identify loci involved in conferring a growth advantage and show that a mutation in the hns gene causing reduced activity of the global regulator H-NS confers a GASP phenotype under specific conditions. The hns-66 allele bears a point mutation within the termination codon of the H-NS open reading frame, resulting in a longer protein that is partially functional. Although isolated from a long-term stationary-phase culture of the parent carrying the rpoS819 allele that results in reduced RpoS activity, the hns-66 survivor showed a growth disadvantage in the early stationary phase (24 to 48 h) when competed against the parent. The hns-66 mutant is also unstable and reverts at a high frequency in the early stationary phase by accumulating second-site suppressor mutations within the ssrA gene involved in targeting aberrant proteins for proteolysis. The mutant was more stable and showed a moderate growth advantage in combination with the rpoS819 allele when competed against a 21-day-old parent. These studies show that H-NS is a target for mutations conferring fitness gain that depends on the genetic background as well as on the stage of the stationary phase.

Albedo enhancement over land to counteract global warming: impacts on hydrological cycle

A recent modelling study has shown that precipitation and runoff over land would increase when the reflectivity of marine clouds is increased to counter global warming. This implies that large scale albedo enhancement over land could lead to a decrease in runoff over land. In this study, we perform simulations using NCAR CAM3.1 that have implications for Solar Radiation Management geoengineering schemes that increase the albedo over land. We find that an increase in reflectivity over land that mitigates the global mean warming from a doubling of CO2 leads to a large residual warming in the southern hemisphere and cooling in the northern hemisphere since most of the land is located in northern hemisphere. Precipitation and runoff over land decrease by 13.4 and 22.3%, respectively, because of a large residual sinking motion over land triggered by albedo enhancement over land. Soil water content also declines when albedo over land is enhanced. The simulated magnitude of hydrological changes over land are much larger when compared to changes over oceans in the recent marine cloud albedo enhancement study since the radiative forcing over land needed (-8.2 W m(-2)) to counter global mean radiative forcing from a doubling of CO2 (3.3 W m(-2)) is approximately twice the forcing needed over the oceans (-4.2 W m(-2)). Our results imply that albedo enhancement over oceans produce climates closer to the unperturbed climate state than do albedo changes on land when the consequences on land hydrology are considered. Our study also has important implications for any intentional or unintentional large scale changes in land surface albedo such as deforestation/afforestation/reforestation, air pollution, and desert and urban albedo modification.

Markov chain modeling of evolution of strains in reinforced concrete flexural beams

From the analysis of experimentally observed variations in surface strains with loading in reinforced concrete beams, it is noted that there is a need to consider the evolution of strains (with loading) as a stochastic process. Use of Markov Chains for modeling stochastic evolution of strains with loading in reinforced concrete flexural beams is studied in this paper. A simple, yet practically useful, bi-level homogeneous Gaussian Markov Chain (BLHGMC) model is proposed for determining the state of strain in reinforced concrete beams. The BLHGMC model will be useful for predicting behavior/response of reinforced concrete beams leading to more rational design.