Comparative analysis of thermophilic and mesophilic proteins using Protein Energy Networks

Background: Thermophilic proteins sustain themselves and function at higher temperatures. Despite their structural and functional similarities with their mesophilic homologues, they show enhanced stability. Various comparative studies at genomic, protein sequence and structure levels, and experimental works highlight the different factors and dominant interacting forces contributing to this increased stability. Methods: In this comparative structure based study, we have used interaction energies between amino acids, to generate structure networks called as Protein Energy Networks (PENs). These PENs are used to compute network, sub-graph, and node specific parameters. These parameters are then compared between the thermophile-mesophile homologues. Results: The results show an increased number of clusters and low energy cliques in thermophiles as the main contributing factors for their enhanced stability. Further more, we see an increase in the number of hubs in thermophiles. We also observe no community of electrostatic cliques forming in PENs. Conclusion: In this study we were able to take an energy based network approach, to identify the factors responsible for enhanced stability of thermophiles, by comparative analysis. We were able to point out that the sub-graph parameters are the prominent contributing factors. The thermophiles have a better-packed hydrophobic core. We have also discussed how thermophiles, although increasing stability through higher connectivity retains conformational flexibility, from a cliques and communities perspective.

Comparative analysis of thermophilic and mesophilic proteins using Protein Energy Networks

Background: Thermophilic proteins sustain themselves and function at higher temperatures. Despite their structural and functional similarities with their mesophilic homologues, they show enhanced stability. Various comparative studies at genomic, protein sequence and structure levels, and experimental works highlight the different factors and dominant interacting forces contributing to this increased stability. Methods: In this comparative structure based study, we have used interaction energies between amino acids, to generate structure networks called as Protein Energy Networks (PENs). These PENs are used to compute network, sub-graph, and node specific parameters. These parameters are then compared between the thermophile-mesophile homologues. Results: The results show an increased number of clusters and low energy cliques in thermophiles as the main contributing factors for their enhanced stability. Further more, we see an increase in the number of hubs in thermophiles. We also observe no community of electrostatic cliques forming in PENs. Conclusion: In this study we were able to take an energy based network approach, to identify the factors responsible for enhanced stability of thermophiles, by comparative analysis. We were able to point out that the sub-graph parameters are the prominent contributing factors. The thermophiles have a better-packed hydrophobic core. We have also discussed how thermophiles, although increasing stability through higher connectivity retains conformational flexibility, from a cliques and communities perspective.

Studies on the conformation of amino acids XII. Energy calculations on prolyl residue

The favoured conformations of the prolyl residue have been obtained by calculating their potential energies arising from bond-angle strain, torsion-angle strain, non-bonded and electrostatic interatomic energies. In addition to the five membered ring, the peptide unit at the amino end (with ω = 180°) and the C′ atom at the carboxyl end have been taken into account. It is found that there are two local minima in the configurational space of the parameters defining the conformation, as is actually observed-one (denoted by B) with Cγ displaced on the same side as C′, which is lower in energy than the other (denoted by A) with Cγ displaced on the opposite side of C′. The other four atoms Cδ, N, Cα, Cβ are nearly in a plane. The conformations of minimum energy (for both A and B) have bond angles very close to the mean observed values while the torsion angles are well within the range observed in various structures for each type. Taking into account the fact that the influence of neighbouring molecules in a crystal structure may make the conformation of a molecule different from the minimal one, the ranges of the conformational parameters for which the energy is within 0.6 kcal/mole above the minimum value (called the "most probable range") and within 1.2 kcal/mole (called the "probable range") have been determined. The ranges thus obtained, agree well with observation, and most of the observed data lie within the most probable ranges, although differing appreciably from the conformation of minimum energy. The study has been extended, in a limited way, to the conformation of the ring in the amino acid proline. Since the nitrogen is tetrahedral in this (as contrasted with being planar in the prolyl residue), it is found that any one of the five atoms can be out of plane (either way), with the other four lying nearly in a plane. These correspond to low energy conformations (up to 1.2 kcal/mole above the minimum). One such example, in which the Cα atom is out of plane is known for dl-proline · HCl. It is also shown that in these calculations energies due to bond length distortions can be neglected to a good degree of approximation, provided the 'best' values of the bond lengths for the particular compound are used in the theoretical calculations.

Top Quark Mass Measurement in the lepton+jets Channel Using a Matrix Element Method and in situ Jet Energy Calibration

A precision measurement of the top quark mass m_t is obtained using a sample of ttbar events from ppbar collisions at the Fermilab Tevatron with the CDF II detector. Selected events require an electron or muon, large missing transverse energy, and exactly four high-energy jets, at least one of which is tagged as coming from a b quark. A likelihood is calculated using a matrix element method with quasi-Monte Carlo integration taking into account finite detector resolution and jet mass effects. The event likelihood is a function of m_t and a parameter DJES to calibrate the jet energy scale /in situ/. Using a total of 1087 events, a value of m_t = 173.0 +/- 1.2 GeV/c^2 is measured.

Search for Anomalous Production of Events with a Photon, Jet, b-quark Jet, and Missing Transverse Energy

We present a signature-based search for anomalous production of events containing a photon, two jets, of which at least one is identified as originating from a b quark, and missing transverse energy. The search uses data corresponding to 2.0/fb of integrated luminosity from p-pbar collisions at a center-of-mass energy of sqrt(s)=1.96 TeV, collected with the CDF II detector at the Fermilab Tevatron. From 6,697,466 events with a photon candidate with transverse energy ET> 25 GeV, we find 617 events with missing transverse energy > 25 GeV and two or more jets with ET> 15 GeV, at least one identified as originating from a b quark, versus an expectation of 607+- 113 events. Increasing the requirement on missing transverse energy to 50 GeV, we find 28 events versus an expectation of 30+-11 events. We find no indications of non-standard-model phenomena.

Search for standard model Higgs boson production in association with a W boson using a neural network discriminant at CDF

We present a search for standard model Higgs boson production in association with a W boson in proton-antiproton collisions at a center of mass energy of 1.96 TeV. The search employs data collected with the CDF II detector that correspond to an integrated luminosity of approximately 1.9 inverse fb. We select events consistent with a signature of a single charged lepton, missing transverse energy, and two jets. Jets corresponding to bottom quarks are identified with a secondary vertex tagging method, a jet probability tagging method, and a neural network filter. We use kinematic information in an artificial neural network to improve discrimination between signal and background compared to previous analyses. The observed number of events and the neural network output distributions are consistent with the standard model background expectations, and we set 95% confidence level upper limits on the production cross section times branching fraction ranging from 1.2 to 1.1 pb or 7.5 to 102 times the standard model expectation for Higgs boson masses from 110 to $150 GeV/c^2, respectively.

Search for Anomalous Production of Events with Two Photons and Additional Energetic Objects at CDF

We present results of a search for anomalous production of two photons together with an electron, muon, $\tau$ lepton, missing transverse energy, or jets using $p\bar{p}$ collision data from 1.1-2.0 fb$^{-1}$ of integrated luminosity collected by the Collider Detector at Fermilab (CDF). The event yields and kinematic distributions are examined for signs of new physics without favoring a specific model of new physics. The results are consistent with the standard model expectations. The search employs several new analysis techniques that significantly reduce instrumental backgrounds in channels with an electron and missing transverse energy.

Free energy of mixing of divalent basic oxides with silica

An expression derived for the free energy of mixing of a divalent basic oxide (MO) with SiO2 based on a model of silicate structure, takes into account the distribution of O2- (from MO) into the silica network, the mixing of silicate ions with O2- and the enthalpy of mixing. The resulting expression is ΔGmix=RT{N11n (2N1-N)2/4N1(1-N)+N21n N 2-N/1-N}, where N={(β+N1)-√(β+N 1)2-8βN1N2}/2β β=characteristic constant for the system N1=mol fraction of silica N2=mol fraction of MO. For the proper choice of β, calculated values of the activity of MO for the system PbO-SiO2, MnO-SiO2, FeO-SiO2 and CaO-SiO2 are in good agreement with experiment. The model predicts that the activity of the basic oxide decreases with increase in temperature.

Potential functions for hydrogen bond interactions - IV. Minimum Energy Conformation of the α-Helical Structure of PoIy-L-Alanine

Making use of the empirical potential functions for peptide NH .. O bonds, developed in this laboratory, the relative stabilities of the rightand left-handed α-helical structures of poly-L-alanine have been investigated, by calculating their conformational energies (V). The value of Vmin of the right-handed helix (αP) is about - 10.4 kcal/mole, and that of the left-handed helix (αM) is about - 9.6 kcal/mole, showing that the former is lower in energy by 0.8 kcal/mole. The helical parameters of the stable conformation of αP are n ∼ 3.6 and h ∼ 1.5 Å. The hydrogen bond of length 2.85 Å and nonlinearity of about 10° adds about 4.0 kcal/ mole to the stabilising energy of the helix in the minimum enregy region. The energy minimum is not sharply defined, but occurs over a long valley, suggesting that a distribution of conformations (φ{symbol}, ψ) of nearly the same energy may occur for the individual residues in a helix. The experimental data of a-helical fibres of poly-L-alanine are in good agreement with the theoretical results for αP. In the case of proteins, the mean values of (φ{symbol}, ψ) for different helices are distributed, but they invariably occur within the contour for V = Vmin + 2 kcal/mole for αP.

Electron energy equation for an atomic radiating plasma

The electron-energy equation for an atomic radiating plasma is considered in this work. Using the atomic model of Bates, Kingston and McWhirter, the radiation loss-term valid for all optical thicknesses is obtained. A study of the energy gained by electrons in inelastic collisions shows that the radiation loss term can be neglected only for rapidly-decaying or fast-growing plasmas. Emission from optically thin plasmas is considered next and an exact expression is given for the total radiation loss in a recombination continuum. A derivation of the Kramers-Unsöld approximation is presented and the error involved in estimating the total emitted recombination radiation by this approximation is shown to be small.

Energy landscape view of phase transitions and slow dynamics in thermotropic liquid crystals

Thermotropic liquid crystals are known to display rich phase behavior on temperature variation. Although the nematic phase is orientationally ordered but translationally disordered, a smectic phase is characterized by the appearance of a partial translational order in addition to a further increase in orientational order. In an attempt to understand the interplay between orientational and translational order in the mesophases that thermotropic liquid crystals typically exhibit upon cooling from the high-temperature isotropic phase, we investigate the potential energy landscapes of a family of model liquid crystalline systems. The configurations of the system corresponding to the local potential energy minima, known as the inherent structures, are determined from computer simulations across the mesophases. We find that the depth of the potential energy minima explored by the system along an isochor grows through the nematic phase as temperature drops in contrast to its insensitivity to temperature in the isotropic and smectic phases. The onset of the growth of the orientational order in the parent phase is found to induce a translational order, resulting in a smectic-like layer in the underlying inherent structures; the inherent structures, surprisingly, never seem to sustain orientational order alone if the parent nematic phase is sandwiched between the high-temperature isotropic phase and the low-temperature smectic phase. The Arrhenius temperature dependence of the orientational relaxation time breaks down near the isotropic-nematic transition. We find that this breakdown occurs at a temperature below which the system explores increasingly deeper potential energy minima.