On the interaction of bovine serum albumin with ionic surfactants: Temperature induced EPR changes of a maleimide nitroxide reflect local protein dynamics and probe solvent accessibility

The interaction of bovine serum albumin (BSA) with the ionic surfactants sodium dodecylsulfate (SDS, anionic), cetyltrimethylammonium chloride (CTAC, cationic) and N-hexadecyl-N,N-dimethyl-3-ammonio-1-propanesulfonate (HPS, zwitterionic) was studied by electron paramagnetic resonance (EPR) spectroscopy of spin label covalently bound to the single free thiol group of the protein. EPR spectra simulation allows to monitor the protein dynamics at the labeling site and to estimate the changes in standard Gibbs free energy, enthalpy and entropy for transferring the nitroxide side chain from the more motionally restricted to the less restricted component. Whereas SDS and CTAC showed similar increases in the dynamics of the protein backbone for all measured concentrations. HPS presented a smaller effect at concentrations above 1.5 mM. At 10 mM of surfactants and 0.15 mM BSA, the standard Gibbs free energy change was consistent with protein backbone conformations more expanded and exposed to the solvent as compared to the native protein, but with a less pronounced effect for HPS. In the presence of the surfactants, the enthalpy change, related to the energy required to dissociate the nitroxide side chain from the protein, was greater, suggesting a lower water activity. The nitroxide side chain also detected a higher viscosity environment in the vicinity of the paramagnetic probe induced by the addition of the surfactants. The results suggest that the surfactant-BSA interaction, at higher surfactant concentration, is affected by the affinities of the surfactant to its own micelles and micelle-like aggregates. Complementary DLS data suggests that the temperature induced changes monitored by the nitroxide probe reflects local changes in the vicinity of the single thiol group of Cys-34 BSA residue. (C) 2011 Elsevier B.V. All rights reserved.

Structure, Dynamics, and RNA Interaction Analysis of the Human SBDS Protein

Shwachman-Bodian-Diamond syndrome is an autosomal recessive genetic syndrome with pleiotropic phenotypes, including pancreatic deficiencies, bone marrow dysfunctions with increased risk of myelodysplasia or leukemia, and skeletal abnormalities. This syndrome has been associated with mutations in the SBDS gene, which encodes a conserved protein showing orthologs in Archaea and eukaryotes. The Shwachman-Bodian-Diamond syndrome pleiotropic phenotypes may be an indication of different cell type requirements for a fully functional SBDS protein. RNA-binding activity has been predicted for archaeal and yeast SBDS orthologs, with the latter also being implicated in ribosome biogenesis. However, full-length SBDS orthologs function in a species-specific manner, indicating that the knowledge obtained from model systems may be of limited use in understanding major unresolved issues regarding SBDS function, namely, the effect of mutations in human SBDS on its biochemical function and the specificity of RNA interaction. We determined the solution structure and backbone dynamics of the human SBDS protein and describe its RNA binding site using NMR spectroscopy. Similarly to the crystal structures of Archaea, the overall structure of human SBDS comprises three well-folded domains. However, significant conformational exchange was observed in NMR dynamics experiments for the flexible linker between the N-terminal domain and the central domain, and these experiments also reflect the relative motions of the domains. RNA titrations monitored by heteronuclear correlation experiments and chemical shift mapping analysis identified a classic RNA binding site at the N-terminal FYSH (fungal, Yhr087wp, Shwachman) domain that concentrates most of the mutations described for the human SBDS. (C) 2010 Elsevier Ltd. All rights reserved.

Combining structure and dynamics: non-denaturing high-pressure effect on lysozyme in solution

Small-angle X-ray scattering (SAXS) and elastic and quasi-elastic neutron scattering techniques were used to investigate the high-pressure-induced changes on interactions, the low-resolution structure and the dynamics of lysozyme in solution. SAXS data, analysed using a global-fit procedure based on a new approach for hydrated protein form factor description, indicate that lysozyme completely maintains its globular structure up to 1500 bar, but significant modi. cations in the protein-protein interaction potential occur at approximately 600-1000 bar. Moreover, the mass density of the protein hydration water shows a clear discontinuity within this pressure range. Neutron scattering experiments indicate that the global and the local lysozyme dynamics change at a similar threshold pressure. A clear evolution of the internal protein dynamics from diffusing to more localized motions has also been probed. Protein structure and dynamics results have then been discussed in the context of protein-water interface and hydration water dynamics. According to SAXS results, the new configuration of water in the first hydration layer induced by pressure is suggested to be at the origin of the observed local mobility changes.

Evolutionary dynamics of the immunodominant repeats of the Plasmodium vivax malaria-vaccine candidate circumsporozoite protein (CSP)

The circumsporozoite protein (CSP) of Plasmodium vivax, a major target for malaria vaccine development, has immunodominant B-cell epitopes mapped to central nonapeptide repeat arrays. To determine whether rearrangements of repeat motifs during mitotic DNA replication of parasites create significant CSP diversity under conditions of low effective meiotic recombination rates, we examined csp alleles from sympatric P. vivax isolates systematically sampled from an area of low malaria endemicity in Brazil over a period of 14 months. Nine unique csp types, comprising six different nona peptide repeats, were observed in 45 isolates analyzed. Identical or nearly identical repeats predominated in most arrays, consistent with their recent expansion. We found strong linkage disequilibrium at sites across the chromosome 8 segment flanking the csp locus, consistent with rare meiotic recombination in this region. We conclude that CSP repeat diversity may not be severely constrained by rare meiotic recombination in areas of low malaria endemicity. New repeat variants may be readily created by nonhomologous recombination even when meiotic recombination is rare, with potential implications for CSP-based vaccine development. (C) 2010 Elsevier B.V. All rights reserved.

Sequence diversity and evolutionary dynamics of the dimorphic antigen merozoite surface protein-6 and other Msp genes of Plasmodium falciparum

Immune evasion by Plasmodium falciparum is favored by extensive allelic diversity of surface antigens. Some of them, most notably the vaccine-candidate merozoite surface protein (MSP)-1, exhibit a poorly understood pattern of allelic dimorphism, in which all observed alleles group into two highly diverged allelic families with few or no inter-family recombinants. Here we describe contrasting levels and patterns of sequence diversity in genes encoding three MSP-1-associated surface antigens of P. falciparum, ranging from an ancient allelic dimorphism in the Msp-6 gene to a near lack of allelic divergence in Msp-9 to a more classical multi-allele polymorphism in Msp-7 Other members of the Msp-7 gene family exhibit very little polymorphism in non-repetitive regions. A comparison of P. falciparum Msp-6 sequences to an orthologous sequence from P. reichenowi provided evidence for distinct evolutionary histories of the 5` and 3` segments of the dimorphic region in PfMsp-6, consistent with one dimorphic lineage having arisen from recombination between now-extinct ancestral alleles. In addition. we uncovered two surprising patterns of evolution in repetitive sequence. Firsts in Msp-6, large deletions are associated with (nearly) identical sequence motifs at their borders. Second, a comparison of PfMsp-9 with the P. reichenowi ortholog indicated retention of a significant inter-unit diversity within an 18-base pair repeat within the coding region of P. falciparum, but homogenization in P. reichenowi. (C) 2009 Elsevier B.V. All rights reserved.

Only subtle protein conformational adaptations are required for ligand binding to thyroid hormone receptors: Simulations using a novel multipoint steered molecular dynamics approach

Thyroid hormone receptors (TR) are hormone-dependent transcription regulators that play a major role in human health, development, and metabolic functions. The thyroid hormone resistance syndrome, diabetes, obesity, and some types of cancer are just a few examples of important diseases that are related to TR malfunctioning, particularly impaired hormone binding. Ligand binding to and dissociation from the receptor ultimately control gene transcription and, thus, detailed knowledge of binding and release mechanisms are fundamental for the comprehension of the receptor`s biological function and development of pharmaceuticals. In this work, we present the first computational study of ligand entry into the ligand binding domain (LBD) of a nuclear receptor. We report molecular dynamics simulations of ligand binding to TRs using a generalization of the steered molecular dynamics technique designed to perform single-molecule pulling simulations along arbitrarily nonlinear driving pathways. We show that only gentle protein movements and conformational adaptations are required for ligand entry into the LBDs and that the magnitude of the forces applied to assist ligand binding are of the order of the forces involved in ligand dissociation. Our simulations suggest an alternative view for the mechanisms ligand binding and dissociation of ligands from nuclear receptors in which ligands can simply diffuse through the protein surface to reach proper positioning within the binding pocket. The proposed picture indicates that the large-amplitude protein motions suggested by the apo- and holo-RXR alpha crystallographic structures are not required, reconciling conformational changes of LBDs required for ligand entry with other nuclear receptors apo-structures that resemble the ligand-bound LBDs.

Dynamics of biochemical and morphophysiological changes during zygotic embryogenesis in Acca sellowiana (Berg.) Burr.

Acca sellowiana (Berg.) Burr. is a native Myrtaceae from southern Brazil and Uruguay, now the subject of a domestication and breeding program. Biotechnological tools have been used to assist in this program. The establishment of a reliable protocol of somatic embryogenesis has been pursued, with a view to capturing and fixing genetic gains. The rationale behind this work relies on the fact that deepening comprehension of the general metabolism of zygotic embryogenesis may certainly improve the protocol for somatic embryogenesis. Thus, in the present work we studied the accumulation of protein, total sugars, starch, amino acids, polyamines (PAs), IAA and ABA, in different stages of A. sellowiana zygotic embryogenesis. Starch is the predominant storage compound during zygotic embryo development. Increased synthesis of amino acids in the cotyledonary stage, mainly of asparagine, was observed throughout development. Total free PAs showed increased synthesis, whereas total conjugated PAs were mainly observed in the early developmental stages. IAA decreased and ABA increased with the progression from early to late embryogenesis. Besides providing basic information on the morphophysiological and biochemical changes of zygotic embryogenesis, the results here obtained may provide adequate strategies towards the modulation of somatic embryogenesis in this species as well as in other woody angiosperms.

Thiol and Sulfenic Acid Oxidation of AhpE, the One-Cysteine Peroxiredoxin from Mycobacterium tuberculosis: Kinetics, Acidity Constants, and Conformational Dynamics

Drug resistance and virulence of Mycobacterium tuberculosis are partially related to the pathogen`s antioxidant systems. Peroxide detoxification in this bacterium is achieved by the heme-containing catalase peroxidase and different two-cysteine peroxiredoxins. M. tuberculosis genome also codifies for a putative one-cysteine peroxiredoxin, alkyl hydroperoxide reductase E (MtAhpE). Its expression was previously demonstrated at a transcriptional level, and the crystallographic structure of the recombinant protein was resolved under reduced and oxidized states. Herein, we report that the conformation of MtAhpE changed depending on its single cysteine redox state, as reflected by different tryptophan fluorescence properties and changes in quaternary structure. Dynamics of fluorescence changes, complemented by competition kinetic assays, were used to perform protein functional studies. MtAhE reduced peroxynitrite 2 orders of magnitude faster than hydrogen peroxide (1.9 x 10(7) M(-1) s(-1) vs 8.2 x 10(4) M(-1) s(-1) at pH 7.4 and 25 degrees C, respectively). The latter also caused cysteine overoxidation to sulfinic acid, but at much slower rate constant (40 M(-1) s(-1)). The pK(a) of the thiol in the reduced enzyme was 5.2, more than one unit lower than that of the sulfenic acid in the oxidized enzyme. The pH profile of hydrogen peroxide-mediated thiol and sulfenic acid oxidations indicated thiolate and sulfenate as the reacting species. The formation of sulfenic acid as well as the catalytic peroxidase activity of MtAhpE was demonstrated using the artificial reducing substrate thionitrobenzoate. Taken together, our results indicate that MtAhpE is a relevant component in the antioxidant repertoire of M. tuberculosis probably involved in peroxide and specially peroxynitrite detoxification.

Population dynamics of genetically diverse Plasmodium falciparum lineages: community-based prospective study in rural Amazonia

Temporal changes in the prevalence of antigenic variants in Plasmodium falciparum populations have been interpreted as evidence of immune-mediated frequency-dependent selection, but evolutively neutral processes may generate similar patterns of serotype replacement. Over 4 years, we investigated the population dynamics of P. falciparum polymorphisms the community level by using 11 putatively neutral microsatellite markers. Plasmodium falciparum Populations were less diverse than sympatric P. vivax isolates, with less multiple-clone infections, lower number of alleles per locus and lower Virtual heterozygosity, but both species showed significant multilocus linkage disequilibrium. Evolutively neutral P. falciparum polymorphisms showed a high turnover rate, with few lineages persisting for several months in the population. Similar results had previously been obtained, in the same community, for sympatric P. vivax isolates. In contrast, the prevalence of the 2 dimorphic types of a major antigen, MSP-2, remained remarkably stable throughout the Study period. We Suggest that the relatively fast turnover of parasite lineages represents the typical population dynamics of neutral polymorphisms in small populations, with clear implications for the detection of frequency-dependent selection of polymorphisms.

Protein lethality investigated in terms of long range dynamical interactions

The relationship between network structure/dynamics and biological function constitutes a fundamental issue in systems biology. However, despite many related investigations, the correspondence between structure and biological functions is not yet fully understood. A related subject that has deserved particular attention recently concerns how essentiality is related to the structure and dynamics of protein interactions. In the current work, protein essentiality is investigated in terms of long range influences in protein-protein interaction networks by considering simulated dynamical aspects. This analysis is performed with respect to outward activations, an approach which models the propagation of interactions between proteins by considering self-avoiding random walks. The obtained results are compared to protein local connectivity. Both the connectivity and the outward activations were found to be strongly related to protein essentiality.

Molecular Basis of the Thermostability and Thermophilicity of Laminarinases: X-ray Structure of the Hyperthermostable Laminarinase from Rhodothermus marinus and Molecular Dynamics Simulations

Glycosyl hydrolases are enzymes capable of breaking the glycosidic linkage of polysaccharides and have considerable industrial and biotechnological applications. Driven by the later applications, it is frequently desirable that glycosyl hydrolases display stability and activity under extreme environment conditions, such as high temperatures and extreme pHs. Here, we present X-ray structure of the hyperthermophilic laminarinase from Rhodothermus marinus (RmLamR) determined at 1.95 angstrom resolution and molecular dynamics simulation studies aimed to comprehend the molecular basis, for the thermal stability of this class of enzymes. As most thermostable proteins, RmLamR contains a relatively large number of salt bridges, which are not randomly distributed on the structure. On the contrary, they form clusters interconnecting beta-sheets of the catalytic domain. Not all salt bridges, however, are beneficial for the protein thermostability: the existence of charge-charge interactions permeating the hydrophobic core of the enzymes actually contributes to destabilize the structure by facilitating water penetration into hydrophobic cavities, as can be seen in the case of mesophilic enzymes. Furthermore, we demonstrate that the mobility of the side-chains is perturbed differently in each class of enzymes. The side-chains of loop residues surrounding the catalytic cleft in the mesophilic laminarinase gain mobility and obstruct the active site at high temperature. By contrast, thermophilic laminarinases preserve their active site flexibility, and the active-site cleft remains accessible for recognition of polysaccharide substrates even at high temperatures. The present results provide structural insights into the role played by salt-bridges and active site flexibility on protein thermal stability and may be relevant for other classes of proteins, particularly glycosyl hydrolases.

The Catalytic Acid in the Dephosphorylation of the Cdk2-pTpY/CycA Protein Complex by Cdc25B Phosphatase

The development of anticancer therapeutics that target Cdc25 phosphatases is now an active area of research. A complete understanding of the Cdc25 catalytic mechanism would certainly allow a more rational inhibitor design. However, the identity of the catalytic acid used by Cdc25 has been debated and not established unambiguously. Results of molecular dynamics simulations with a calibrated hybrid potential for the first reaction step catalyzed by Cdc25B in complex with its natural substrate, the Cdk2-pTpY/CycA protein complex, are presented here. The calculated reaction free-energy profiles are in very good agreement with experimental measurements and are used to discern between different proposals for the general acid. In addition, the simulations give useful insight on interactions that can be explored for the design of inhibitors specific to Cdc25.