How optimization of potential functions affects protein folding.

The relationship between the optimization of the potential function and the foldability of theoretical protein models is studied based on investigations of a 27-mer cubic-lattice protein model and a more realistic lattice model for the protein crambin. In both the simple and the more complicated systems, optimization of the energy parameters achieves significant improvements in the statistical-mechanical characteristics of the systems and leads to foldable protein models in simulation experiments. The foldability of the protein models is characterized by their statistical-mechanical properties--e.g., by the density of states and by Monte Carlo folding simulations of the models. With optimized energy parameters, a high level of consistency exists among different interactions in the native structures of the protein models, as revealed by a correlation function between the optimized energy parameters and the native structure of the model proteins. The results of this work are relevant to the design of a general potential function for folding proteins by theoretical simulations.

Potential energy functions for ground state surfaces of HCO and HNO

Analytical potential functions are reported for the ground state surfaces of HCO and HNO, the functions being derived from spectroscopic and ab initio data. Harmonized force fields have been deduced for the stable configurations of both molecules and vibration frequencies predicted for the metastable species COH and NOH.

Fluorescence Microscopy Imaging Denoising with Log-Euclidean Priors and Photobleaching Compensation

Fluorescent protein microscopy imaging is nowadays one of the most important tools in biomedical research. However, the resulting images present a low signal to noise ratio and a time intensity decay due to the photobleaching effect. This phenomenon is a consequence of the decreasing on the radiation emission efficiency of the tagging protein. This occurs because the fluorophore permanently loses its ability to fluoresce, due to photochemical reactions induced by the incident light. The Poisson multiplicative noise that corrupts these images, in addition with its quality degradation due to photobleaching, make long time biological observation processes very difficult. In this paper a denoising algorithm for Poisson data, where the photobleaching effect is explicitly taken into account, is described. The algorithm is designed in a Bayesian framework where the data fidelity term models the Poisson noise generation process as well as the exponential intensity decay caused by the photobleaching. The prior term is conceived with Gibbs priors and log-Euclidean potential functions, suitable to cope with the positivity constrained nature of the parameters to be estimated. Monte Carlo tests with synthetic data are presented to characterize the performance of the algorithm. One example with real data is included to illustrate its application.

Culture of mouse peritoneal macrophages with mouse serum induces lipid bodies that associate with the parasitophorous vacuole and decrease their microbicidal capacity against Toxoplasma gondii

Lipid bodies [lipid droplets (LBs)] are lipid-rich organelles involved in lipid metabolism, signalling and inflammation. Recent findings suggest a role for LBs in host response to infection; however, the potential functions of this organelle in Toxoplasma gondii infection and how it alters macrophage microbicidal capacity during infection are not well understood. Here, we investigated the role of host LBs in T. gondii infection in mouse peritoneal macrophages in vitro. Macrophages cultured with mouse serum (MS) had higher numbers of LBs than those cultured in foetal bovine serum and can function as a model to study the role of LBs during intracellular pathogen infection. LBs were found in association with the parasitophorous vacuole, suggesting that T. gondii may benefit from this lipid source. Moreover, increased numbers of macrophage LBs correlated with high prostaglandin E2 (PGE2) production and decreased nitric oxide (NO) synthesis. Accordingly, LB-enriched macrophages cultured with MS were less efficient at controlling T. gondii growth. Treatment of macrophages cultured with MS with indomethacin, an inhibitor of PGE2 production, increased the microbicidal capacity against T. gondii. Collectively, these results suggest that culture with MS caused a decrease in microbicidal activity of macrophages against T. gondii by increasing PGE2 while lowering NO production.

Characterization of the expression, regulation and function of Sorting Nexin 9 and Sorting Nexin 18

Summary : Sorting nexin (SNX) family members play important roles in intracellular protein and membrane trafficking, The membrane-tubulating SNX9 protein has been shown to interact with multiple components of the endocytic machinery and to participate in clathrin-mediated endocytosis of cell surface receptors. It has not been investigated if SNX9 may also participate in other protein sorting pathways that involve vesicular transport, specifically the biogenesis of lysosome-related organelles (LROs). Closely related to SNX9 is SNXl8, whose function is largely unknown. In this work, we have characterized the expression of SNX9 and SNXl8 in LRO-containing cells and investigated their role in protein trafficking during the formation of LROs. Our results indicate that SNX9 and SNXl8 are not essential for the formation of LROs, nor for the sorting of melanosomal proteins. We investigated how the level of intracellular SNX9 protein is regulated and found that it is a substrate of the ubiquitin ligase Itch, a member of the NEDD4 family of E3 ubiquitin ligases. Itch ubiquitylates SNX9 and regulates SNX9 levels by enhancing its degradation. Using ? truncated proteins we found that the interaction with SNX9 is mediated by the proline-rich domain of Itch, a domain distinct from the conventional WW recognition domain, and the SH3 domain of SNX9. Interaction with the PRD of Itch is essential for SNX9 ubiquitylation and degradation. We further showed that Itch binding is not affected by tyrosine phosphorylation of SNX9. Using lentivector-mediated siRNA techniques, we found that Itch regulates the level of melanosomal proteins, while knock-down of SNX9 does not alter their level. Interestingly, we revealed that silencing of SNXIS affects the amount of the melanosomal protein Melan-A, but also of SNX9, and that SNXl8 can interact with SNX9. Taken together, our results highlight that the pool of substrates of NEDD4 family E3 ligases extends to proteins containing SH3 domains and provide insight into the potential functions of SNXI8. Résumé : Les membres de la famille des Sorting Nexins (SNX) jouent des rôles importants dans le trafic intracellulaire de protéines et membranes. Il a été démontré que la protéine SNX9, qui génère les tubules membranaires, interagit avec plusieurs composants de la machinerie d'endocytose et participe à l'endocytose des récepteurs de surface mediée par la clathrine. Aucune étude n'a investigué si SNX9 pourrait aussi participer à d'autres voies de trafic de protéines tel que le transport vésiculaire, et plus particulièrement la biogenèse des organites lysosomaux ("lysosome-related organelles", LR©s). SNXl8 est similaire à SNX9, mais sa fonction est largement inconnue. Dans ce travail, nous avons caractérisé l'expression de SNX9 et SNX18 dans des cellules contenants des LROs et investigué leur rôle dans le trafic de protéines pendant la formation des LROS. Nos résultats indiquent que SNX9 et SNXI8 ne sont essentiels ni pour la formation des LR©s, ni pour le trafic de protéines mélanosomales. Nous avons examiné la régulation du niveau intracellulaire de la protéine SNX9 et avons trouvé qu'elle est un substrat de l'ubiquitine ligase Itch, un membre de la famille NEDD4 des ubiquitine ligases E3. Itch ubiquitine SNX9 et régule les niveaux de SNX9 en augmentant sa dégradation. En utilisant des protéines mutées nous avons découvert que l'interaction avec SNX9 est médiée par le domaine riche en proline de Itch, qui est différent du domaine conventionnel de reconnaissance WW, et par le domaine SH3 de SNX9. L'interaction avec le domaine riche en proline de Itch est essentielle pour l'ubiquitination et la dégradation de SNX9. De plus, nous avons montré que cette liaison n'est pas affectée par la phosphorylation des résidus tyrosine de SNX9. En utilisant des vecteurs lentiviraux exprimant des siARN, nous avons trouvé que Itch régule les niveaux de protéines mélanosomales, alors que l'extinction de l'expression de SNX9 ne change pas leurs niveaux. En autre, nous avons révélé que la diminution de SNXl8 affecte le niveau de la protéine mélanosomale Melan-A et de SNX9, et aussi que SNXl8 peut interagir avec SNX9. En résumé, nos résultats démontrent que l'ensemble des substrats de la famille NEDD4 des ubiquitine ligases E3 s'élargit aux protéines contenant des domaines SH3 et ouvrent des perspectives sur les fonctions potentielles de SNXl8.

Computer simulation study of liquid lithium at 470 and 843 K

Both structural and dynamical properties of 7Li at 470 and 843 K are studied by molecular dynamics simulation and the results are comapred with the available experimental data. Two effective interatomic potentials are used, i.e., a potential derived from the Ashcroft pseudopotential [Phys. Lett. 23, 48 (1966)] and a recently proposed potential deduced from the neutral pseudoatom method [J. Phys.: Condens. Matter 5, 4283 (1993)]. Although the shape of the two potential functions is very different, the majority of the properties calculated from them are very similar. The differences among the results using the two interaction models are carefully discussed.

The evolution of lncRNA repertoires and expression patterns in tetrapods.

Only a very small fraction of long noncoding RNAs (lncRNAs) are well characterized. The evolutionary history of lncRNAs can provide insights into their functionality, but the absence of lncRNA annotations in non-model organisms has precluded comparative analyses. Here we present a large-scale evolutionary study of lncRNA repertoires and expression patterns, in 11 tetrapod species. We identify approximately 11,000 primate-specific lncRNAs and 2,500 highly conserved lncRNAs, including approximately 400 genes that are likely to have originated more than 300 million years ago. We find that lncRNAs, in particular ancient ones, are in general actively regulated and may function predominantly in embryonic development. Most lncRNAs evolve rapidly in terms of sequence and expression levels, but tissue specificities are often conserved. We compared expression patterns of homologous lncRNA and protein-coding families across tetrapods to reconstruct an evolutionarily conserved co-expression network. This network suggests potential functions for lncRNAs in fundamental processes such as spermatogenesis and synaptic transmission, but also in more specific mechanisms such as placenta development through microRNA production.

The prolyl-aminodipeptidases and their inhibitors as therapeutic targets for fibrogenic disorders

Many biologically active peptides are protected from general proteolytic degradation by evolutionary conserved prolines (Pro), due to conformational constraints imposed by the Pro residue. Thus the biological importance of prolyl-specific peptidases points to a high potential for drug discovery for this family of enzymes. Panels of inhibitors have been synthesized and their effects, determined in biological models, suggest the inhibition of families of enzymes with similar activities. Prolyl-specific aminodipeptidases include dipeptidyl-aminodipeptidase IV (DPP IV)/CD26, DPP8, DPP9 and fibroblast activation protease-alpha (FAP-alpha)/seprase, able to release X-Pro dipeptides from the N-terminus of peptides. DPP IV inhibitors are in clinical use for type 2 diabetes. In this review, the expression and the potential functions of prolyl-aminodipeptidases are reviewed in diseases, and the inhibitors developed for these enzymes are discussed, with a specific focus on inhibitors able to discriminate between DPP IV and fibroblast activation protease-alpha (FAPalpha)/seprase as potential leads for the treatment of fibrogenic diseases.

Termodinâmica estatística de líquidos com o método de Monte Carlo. I. Metodologia

Statistical mechanics Monte Carlo simulation is reviewed as a formalism to study thermodynamic properties of liquids. Considering the importance of free energy changes in chemical processes, the thermodynamic perturbation theory implemented in the Monte Carlo method is discussed. The representation of molecular interaction by the Lennard-Jones and Coulomb potential functions is also discussed. Charges derived from quantum molecular electrostatic potential are also discussed as an useful methodology to generate an adequate set of partial charges to be used in liquid simulation.

Termodinâmica estatística de líquidos com o método de Monte Carlo. II. O líquido clorofórmio

Thermodynamic properties and radial distribution functions for liquid chloroform were calculated using the Monte Carlo method implemented with Metropolis algorithm in the NpT ensemble at 298 K and 1 atm. A five site model was developed to represent the chloroform molecules. A force field composed by Lennard-Jones and Coulomb potential functions was used to calculate the intermolecular energy. The partial charges needed to represent the Coulombic interactions were obtained from quantum chemical ab initio calculations. The Lennard-Jones parameters were adjusted to reproduce experimental values for density and enthalpy of vaporization for pure liquid. All thermodynamic results are in excelent agreement with experimental data. The correlation functions calculated are in good accordance with theoretical results avaliable in the literature. The free energy for solvating one chloroform molecule into its own liquid at 298 K and 1 atm was computed as an additional test of the potential model. The result obtained compares well with the experimental value. The medium effects on cis/trans convertion of a hypotetical solute in water TIP4P and chloroform solvents were also accomplished. The results obtained from this investigation are in agreement with estimates of the continuous theory of solvation.

On the equation of state and atomic mean square displacement of crystals

We have presented a Green's function method for the calculation of the atomic mean square displacement (MSD) for an anharmonic Hamil toni an . This method effectively sums a whole class of anharmonic contributions to MSD in the perturbation expansion in the high temperature limit. Using this formalism we have calculated the MSD for a nearest neighbour fcc Lennard Jones solid. The results show an improvement over the lowest order perturbation theory results, the difference with Monte Carlo calculations at temperatures close to melting is reduced from 11% to 3%. We also calculated the MSD for the Alkali metals Nat K/ Cs where a sixth neighbour interaction potential derived from the pseudopotential theory was employed in the calculations. The MSD by this method increases by 2.5% to 3.5% over the respective perturbation theory results. The MSD was calculated for Aluminum where different pseudopotential functions and a phenomenological Morse potential were used. The results show that the pseudopotentials provide better agreement with experimental data than the Morse potential. An excellent agreement with experiment over the whole temperature range is achieved with the Harrison modified point-ion pseudopotential with Hubbard-Sham screening function. We have calculated the thermodynamic properties of solid Kr by minimizing the total energy consisting of static and vibrational components, employing different schemes: The quasiharmonic theory (QH), ).2 and).4 perturbation theory, all terms up to 0 ().4) of the improved self consistent phonon theory (ISC), the ring diagrams up to o ().4) (RING), the iteration scheme (ITER) derived from the Greens's function method and a scheme consisting of ITER plus the remaining contributions of 0 ().4) which are not included in ITER which we call E(FULL). We have calculated the lattice constant, the volume expansion, the isothermal and adiabatic bulk modulus, the specific heat at constant volume and at constant pressure, and the Gruneisen parameter from two different potential functions: Lennard-Jones and Aziz. The Aziz potential gives generally a better agreement with experimental data than the LJ potential for the QH, ).2, ).4 and E(FULL) schemes. When only a partial sum of the).4 diagrams is used in the calculations (e.g. RING and ISC) the LJ results are in better agreement with experiment. The iteration scheme brings a definitive improvement over the).2 PT for both potentials.

An analysis of the methyl rotation dynamics in the So (x' A) and T (a A) states of thioacetaldehyde from Pyrolysis jet spectra /

Jet-cooled, laser-induced phosphorescence excitation spectra (LIP) of thioacetaldehyde CH3CHS, CH3CDS, CD3CHS and CD3CDS have been observed over the region 15800 - 17300 cm"^ in a continuous pyrolysis jet. The vibronic band structure of the singlet-triplet n -* n* transition were attributed to the strong coupling of the methyl torsion and aldehydic hydrogen wagging modes . The vibronic peaks have been assigned in terms of two upper electronic state (T^) vibrations; the methyl torsion mode v^g, and the aldehydic hydrogen wagging mode v^^. The electronic origin O^a^ is unequivocally assigned as follows: CH3CHS (16294.9 cm"'' ), CH3CDS (16360.9 cm"'' ), CD3CHS (16299.7 cm"^ ), and CD3CDS (16367.2 cm"'' ). To obtain structural and dynamical information about the two electronic states, potential surfaces V(e,a) for the 6 (methyl torsion) and a (hydrogen wagging) motions were generated by ab initio quantum mechanical calculations with a 6-3 IG* basis in which the structural parameters were fully relaxed. The kinetic energy coefficients BQ(a,e) , B^(a,G) , and the cross coupling term B^(a,e) , were accurately represented as functions of the two active coordinates, a and 9. The calculations reveal that the molecule adopts an eclipsed conformation for the lower Sq electronic state (a=0°,e=0"') with a barrier height to internal rotation of 541.5 cm"^ which is to be compared to 549.8 cm"^ obtained from the microwave experiment. The conformation of the upper T^ electronic state was found to be staggered (a=24 . 68° ,e=-45. 66° ) . The saddle point in the path traced out by the aldehyde wagging motion was calculated to be 175 cm"^ above the equilibrium configuration. The corresponding maxima in the path taken by methyl torsion was found to be 322 cm'\ The small amplitude normal vibrational modes were also calculated to aid in the assignment of the spectra. Torsional-wagging energy manifolds for the two states were derived from the Hamiltonian H(a,e) which was solved variationally using an extended two dimensional Fourier expansion as a basis set. A torsionalinversion band spectrum was derived from the calculated energy levels and Franck-Condon factors, and was compared with the experimental supersonic-jet spectra. Most of the anomalies which were associated with the interpretation of the observed spectrum could be accounted for by the band profiles derived from ab initio SCF calculations. A model describing the jet spectra was derived by scaling the ab initio potential functions. The global least squares fitting generates a triplet state potential which has a minimum at (a=22.38° ,e=-41.08°) . The flatter potential in the scaled model yielded excellent agreement between the observed and calculated frequency intervals.

Analyse de la corrélation conditionnelle dérivée de la coévolution d’un système de trois gènes par un modèle du maximum de vraisemblance

Les gènes codant pour des protéines peuvent souvent être regroupés et intégrés en modules fonctionnels par rapport à un organelle. Ces modules peuvent avoir des composantes qui suivent une évolution corrélée pouvant être conditionnelle à un phénotype donné. Les gènes liés à la motilité possèdent cette caractéristique, car ils se suivent en cascade en réponse à des stimuli extérieurs. L’hyperthermophilie, d’autre part, est interreliée à la reverse gyrase, cependant aucun autre élément qui pourrait y être associé avec certitude n’est connu. Ceci peut être dû à un déplacement de gènes non orthologues encore non résolu. En utilisant une approche bio-informatique, une modélisation mathématique d’évolution conditionnelle corrélée pour trois gènes a été développée et appliquée sur des profils phylétiques d’archaea. Ceci a permis d’établir des théories quant à la fonction potentielle du gène du flagelle FlaD/E ainsi que l’histoire évolutive des gènes lui étant liés et ayant contribué à sa formation. De plus, une histoire évolutive théorique a été établie pour une ligase liée à l’hyperthermophilie.