Interaction Energy Based Protein Structure Networks

The three dimensional structure of a protein is formed and maintained by the noncovalent interactions among the amino acid residues of the polypeptide chain These interactions can be represented collectively in the form of a network So far such networks have been investigated by considering the connections based on distances between the amino acid residues Here we present a method of constructing the structure network based on interaction energies among the amino acid residues in the protein We have investigated the properties of such protein energy based networks (PENs) and have shown correlations to protein structural features such as the clusters of residues involved in stability formation of secondary and super secondary structural units Further we demonstrate that the analysis of PENs in terms of parameters such as hubs and shortest paths can provide a variety of biologically important information such as the residues crucial for stabilizing the folded units and the paths of communication between distal residues in the protein Finally the energy regimes for different levels of stabilization in the protein structure have clearly emerged from the PEN analysis

Insights into the Fold Organization of TIM Barrel from Interaction Energy Based Structure Networks

There are many well-known examples of proteins with low sequence similarity, adopting the same structural fold. This aspect of sequence-structure relationship has been extensively studied both experimentally and theoretically, however with limited success. Most of the studies consider remote homology or ``sequence conservation'' as the basis for their understanding. Recently ``interaction energy'' based network formalism (Protein Energy Networks (PENs)) was developed to understand the determinants of protein structures. In this paper we have used these PENs to investigate the common non-covalent interactions and their collective features which stabilize the TIM barrel fold. We have also developed a method of aligning PENs in order to understand the spatial conservation of interactions in the fold. We have identified key common interactions responsible for the conservation of the TIM fold, despite high sequence dissimilarity. For instance, the central beta barrel of the TIM fold is stabilized by long-range high energy electrostatic interactions and low-energy contiguous vdW interactions in certain families. The other interfaces like the helix-sheet or the helix-helix seem to be devoid of any high energy conserved interactions. Conserved interactions in the loop regions around the catalytic site of the TIM fold have also been identified, pointing out their significance in both structural and functional evolution. Based on these investigations, we have developed a novel network based phylogenetic analysis for remote homologues, which can perform better than sequence based phylogeny. Such an analysis is more meaningful from both structural and functional evolutionary perspective. We believe that the information obtained through the ``interaction conservation'' viewpoint and the subsequently developed method of structure network alignment, can shed new light in the fields of fold organization and de novo computational protein design.

Selectivity in the Interaction of Electron Donor and Acceptor Molecules with Graphene and Single-Walled Carbon Nanotubes

Interaction of electron donor and acceptor molecules with graphene samples prepared by different methods as well as with single-walled carbon nanotubes (SWNTs) has been investigated by isothermal titration calorimetry (ITC). The ITC interaction energies of the graphene samples and SWNTs with electron acceptor molecules are higher than those with electron donor molecules. Thus, tetracyanoethylene (TCNE) shows the highest interaction energy with both graphene and SWNTs. The interaction energy with acceptor molecules varies with the electron affinity as well as with the charge-transfer transition energy for different aromatics. Metallic SWNTs interact reversibly with electron acceptor molecules, resulting in the opening of a gap.

Solvation dynamics, energy distribution and trapping of a light solute ion

In the theoretical treatments of the dynamics of solvation of a newly created ion in a dipolar solvent, the self-motion of the solute is usually ignored. Recently, it has been shown that for a light ion the translational motion of the ion can significantly enhance its own rate of solvation. Therefore, solvation itself may not be the rate determining step in the equilibration. Instead, the rate determining step is the search of the low energy configuration which serves to localize the light ion. In this article a microscopic calculation of the probability distribution of the interaction energy of the nascent charge with the dipolar solvent molecules is presented in order to address this problem of solute trapping. It is found that to a good approximation, this distribution is Gaussian and the second moment of this distribution is exactly equal to the half of its own solvation energy. It is shown that this is in excellent agreement with the simulation results that are available for the model Brownian dipolar lattice and for liquid acetonitrile. If the distortion of the solvent by the ion is negligible then the same relation gives the energy distribution for the solvated ion, with the average centered at the final equilibrium solvation energy. These results are expected to be useful in understanding various chemical processes in dipolar liquids. Another interesting outcome of the present study is a simple dynamic argument that supports Onsager's ''inverse snow-ball'' conjecture of solvation of a light ion. A simple derivation of the semi-phenomenological relation between the solvation time correlation function and the single particle orientation, reported recently by Maroncelli et al. (J. Phys. Chem. 97 (1993) 13), is also presented.

Estimation of surface free energy of pyrites by contact angle measurements

The surface properties of coal-pyrite play a major role in determining its separation from coal in processes such as flotation. The solution pH is an important parameter in determining the surface properties of both coal and coal-pyrite such as surface free energy and zeta-potential. In the present investigation, the effect of pH on the surface free energy of pyrites from different sources was studied. The surface free energy of solids is made up of two components, i.e. the dispersive surface free energy and the acid-base interaction energy. Various methods have been used by previous researchers to evaluate these two components for different solids. In the present study, a new approach was developed and used to study the surface free energy of pyrite surfaces. Results indicate that the dispersion surface free energy of various pyrites is independent of pH while the acid-base interaction energy is strongly dependent on the pH. The acid-base interaction energy is different for each pyrite sample and also the change with pH varies with the type of pyrite. Coal-pyrite was found to be more hydrophobic than ore-pyrite which may be attributed to the presence of carbon in coal-pyrites. The acid-base interaction energy varied little with pH for coal pyrites than ore-pyrite. Comparison of acid-base interaction energy with zeta-potential measurements shows a good correlation between the minimum in acid-base interaction energy and the pHpzc.

A resistometric study of solute–solute interaction in dilute aluminium–silver alloys

The experimentally determined apparent vacancy formation energy values in dilute aluminium—silver alloys showed a divergence from calculated values at higher solute fractions. This is explained in terms of a solute—solute interaction energy of the order of 0.10 ev which exists when the binding energy between a vacancy and a solute atom pair is reduced to zero.

Theoretical studies of lithium bonding in lithium chloride/aliphatic amine complexes

In the systematic study of amine … LiCl [amines = NH3, CH3NH2, (CH3)2NH] complexes the possibility of an ion-pair structure and the effect of methylation on the stabilization energy is investigated. ΔEis evaluated by the SCF/4-31G method and augmented by the approximate dispersion energy calculated perturbationally. The interaction energy decreases with the increasing number of methyl groups in the amine. The dispersion energy plays a negligible role in the stabilization of complexes. None of the systems studied are ion pairs; their Li bonds are of a so-called molecular type. Due to the divergence of the multipole expansion, the attempt to correct the 4-31G stabilization energies via the electrostatic energy fails. The relative order of the ΔE in the series of complexes is verified instead in the extended basis set calculation. The lithium bonds are compared with their H-bonded analogues.

A semi-empirical theory of molten salts

A semi-empirical model is presented for describing the interionic interactions in molten salts using the experimentally available structure data. An extension of Bertaut's method of non-overlapping charges is used to estimate the electrostatic interaction energy in ionic melts. It is shown, in agreement with earlier computer simulation studies, that this energy increases when an ionic salt melts. The repulsion between ions is described using a compressible ion theory which uses structure-independent parameters. The van der Waals interactions and the thermal free energy are also included in the total energy, which is minimised with respect to isostructural volume variations to calculate the equilibrium density. Detailed results are presented for three molten systems, NaCl, CaCl2 and ZnCl2, and are shown to be in satisfactory agreement with experiments. With reliable structural data now being reported for several other molten salts, the present study gains relevance.

Ion channel formation by zervamicin-IIB. A molecular modelling study.

Zervamicin-IIB (Zrv-IIB) is a 16 residue peptaibol which forms voltage-activated, multiple conductance level channels in planar lipid bilayers. A molecular model of Zrv-IIB channels is presented. The structure of monomeric Zrv-IIB is based upon the crystal structure of Zervamicin-Leu. The helical backbone is kinked by a hydroxyproline residue at position 10. Zrv-IIB channels are modelled as helix bundles of from 4 to 8 parallel helices surrounding a central pore. The monomers are packed with their C-terminal helical segments in close contact, and the bundles are stabilized by hydrogen bonds between glutamine 11 and hydroxyproline 10 of adjacent helices. Interaction energy profiles for movement of three different probes species (K+, Cl- and water) through the central pore are analyzed. The conformations of: (a) the sidechain of glutamine 3; (b) the hydroxyl group of hydroxyproline 10; and (c) the C-terminal hydroxyl group are "optimized" in order to maximize favourable interactions between the channel and the probes, resulting in favourable interaction energy profiles for all three. This suggests that conformational flexibility of polar sidechains enables the channel lining to mimic an aqueous environment.

Effect of coordinated ions on structure and flexibility of parallel G-quadruplexes: A molecular dynamics study

ingle tract guanine residues can associate to form stable parallel quadruplex structures in the presence of certain cations. Nanosecond scale molecular dynamics simulations have been performed on fully solvated fibre model of parallel d(G(7)) quadruplex structures with Na+ or K+ ions coordinated in the cavity formed by the O6 atoms of the guanine bases. The AMBER 4.1 force field and Particle Mesh Ewald technique for electrostatic interactions have been used in all simulations. There quadruplex structures are stable during the simulation, with the middle four base tetrads showing root mean square deviation values between 0.5 to 0.8 Angstrom from the initial structure as well the high resolution crystal structure. Even in the absence of any coordinated ion in the initial structure, the G-quadruplex structure remains intact throughout the simulation. During the 1.1 ns MD simulation, one Nai counter ion from the solvent as well as several water molecules enter the central cavity to occupy the empty coordination sites within the parallel quadruplex and help stabilize the structure. Hydrogen bonding pattern depends on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate cations of different sizes. in the absence of any coordinated ion. due to strong mutual repulsion, O6 atoms within G-tetrad are forced farther apart from each other, which leads to a considerably different hydrogen bonding scheme within the G-tetrads and very favourable interaction energy between the guanine bases constituting a G-tetrad. However, a coordinated ion between G-tetrads provides extra stacking energy for the G-tetrads and makes the quadruplex structure more rigid. Na+ ions, within the quadruplex cavity, are more mobile than coordinated K+ ions. A number of hydrogen bonded water molecules are observed within the grooves of all quadruplex structures.

Theoretical studies on alpha-helix--DNA interactions

The interaction energies between (Ala)10 and alpha-helix fragment and different nucleotide sequences in right-handed B-form have been optimized using semi-empirical potential energy functions. The energies are calculated for two different orientations of the alpha-helix, viz., when the alpha-helix axis taken in the N----C direction is (i) parallel and (ii) antiparallel to the 5'-3' ascending strand of DNA, proximal to it. When both the DNA molecule as well as the alpha-helix are treated as rigid molecules it is found that a polyalanine alpha-helix has slightly more favourable contacts when it is in the proximity of a four nucleotide sequence of 5'-(N-A-T-N)-3' type, where N is either a purine or a pyrimidine. However, when the two interacting molecules are allowed to undergo local structural variations then the interaction energy appears to be independent of the base sequence confirming the non-specific nature of these interactions.

Volume effects and associations in liquid alloys

The regular associated solution model for binary systems has been modified by incorporating the size of the complex as an explicit variable. The thermodynamic properties of the liquid alloy and the interactions between theA ?B type of complex and the unassociated atoms in anA-B binary have been evaluated as a function of relative size of the complex using the activity coefficients at infinite dilution and activity data at one other composition in the binary. The computational procedure adopted for determining the concentration of clusters and interaction energies in the associated liquid is similar to that proposed by Lele and Rao. The analysis has been applied to the thermodynamic mixing functions of liquid Al-Ca alloys believed to contain Al2Ca associates. It is found that the size of the cluster significantly affects the interaction energies between the complex and the unassociated atoms, while the equilibrium constant and enthalpy change for the association reaction exhibit only minor variation, when the equations are fitted to experimental data. The interaction energy between unassociated free atoms remains virtually unaltered as the size of the complex is varied between extreme values. Accurate data on free energy, enthalpy, and volume of mixing at the same temperature on alloy systems with compound forming tendency would permit a rigorous test of the proposed model.

Effect of coordinated ions on structure and flexiblity of parallel G-quandruplexes: a molecular dynamics study

Single tract guanine residues can associate to form stable parallel quadruplex structures in the presence of certain cations. Nanosecond scale molecular dynamics simulations have been performed on fully solvated fibre model of parallel d(G7) quadruplex structures with Na+ or K+ ions coordinated in the cavity formed by the 06 atoms of the guanine bases. The AMBER 4.1 force field and Particle Mesh Ewald technique for electrostatic interactions have been used in all simulations. These quadruplex structures are stable during the simulation, with the middle four base tetrads showing root mean square deviation values between 0.5 to 0.8 A from the initial structure as well the high resolution crystal structure. Even in the absence of any coordinated ion in the initial structure, the G-quadruplex structure remains intact throughout the simulation. During the 1.1 ns MD simulation, one Na+ counter ion from the solvent as well as several water molecules enter the central cavity to occupy the empty coordination sites within the parallel quadruplex and help stabilize the structure. Hydrogen bonding pattern depends on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate cations of different sizes. In the absence of any coordinated ion, due to strong mutual repulsion, 06 atoms within G-tetrad are forced farther apart from each other, which leads to a considerably different hydrogen bonding scheme within the G-tetrads and very favourable interaction energy between the guanine bases constituting a G-tetrad. However, a coordinated ion between G-tetrads provides extra stacking energy for the G-tetrads and makes the quadruplex structure more rigid. Na+ ions, within the quadruplex cavity, are more mobile than coordinated K+ ions. A number of hydrogen bonded water molecules are observed within the grooves of all quadruplex structures

Role of molecular reorientation in the vitrification of molecular liquids: A Monte Carlo study of the liquid and glassy states of 2,2-dimethylbutane

The liquid and the glassy phases of 2,2-dimethylbutane have been investigated by isothermal isobaric ensemble Monte Carlo simulation. Thermodynamic Properties and radial distribution functions for both the liquid and the glass have been obtained. The radial distribution functions have been classified into three types based on the accessibility of the group. It has been shown that the structure of the Iiquid and the glass can be understood in terms of the above classification of the radial distribution functions. Molecular reorientation plays an important role in the structural rearrangement accompanying glass formation. As much as 35% of the contribution to the increase in the intermolecular interaction energy on vitrification is due to the reorientation of the neighbouring pairs of molecules. The observed changes in the dimerisation energy and the bonding energy distribution function are consistent with the observed structural changes.

Miscibility of Poly(aromatic acrylate)s and Methacrylates with Styrene/Acrylonitrile copolymers

Miscibilities of some poly[aromatic (meth)acrylatels namely, poly(pheny1 acrylate) (PPA), poly(pheny1 methacrylate) (PPMA), poly(benzy1 acrylate) (PBA), and poly(benzy1 methacrylate) (PBMAY polystyrene blends, have been studied through the so-called copolymer effect by incorporating acrylonitrile units in PS chains. In these systems, miscibility occurs on account of the strong repulsion between the acrylonitrile and styrene units in the copolymer. PBA and PBMA were blended with different styreneacrylonitrile (SAN) copolymers. A miscibility window has been identified for the latter system, and from these limits, the binary interaction energy density parameters (B,j.’sw) ere calculated. Using these values, the miscibilities in other homopolymer-copolymer and copolymer-copolymer systems containing benzyl methacrylate, acrylonitrile, and styrene monomer units have been predicted and subsequently verified experimentally. The miscibility window limits in poly[aromatic (meth)acrylate]s/SAN copolymer blends have been compared. PBA does not exhibit a miscibility window with SAN copolymers, which has been explained by the weak intramolecular hydrogen bonding in PBA. The miscibility window in the PBW SAN copolymer system, as observed by DSC, shows a considerable narrowing in nonradiative energy transfer (NRET) measurements, as this technique is more sensitive.