On the relationship of thermodynamic parameters with the buried surface area in protein-ligand complex formation

Prediction of thermodynamic parameters of protein-protein and antigen-antibody complex formation from high resolution structural parameters has recently received much attention, since an understanding of the contributions of different fundamental processes like hydrophobic interactions, hydrogen bonding, salt bridge formation, solvent reorganization etc. to the overall thermodynamic parameters and their relations with the structural parameters would lead to rational drug design. Using the results of the dissolution of hydrocarbons and other model compounds the changes in heat capacity (DeltaCp), enthalpy (DeltaH) and entropy (DeltaS) have been empirically correlated with the polar and apolar surface areas buried during the process of protein folding/unfolding and protein-ligand complex formation. In this regard, the polar and apolar surfaces removed from the solvent in a protein-ligand complex have been calculated from the experimentally observed values of changes in heat capacity (DeltaCp) and enthalpy (DeltaH) for protein-ligand complexes for which accurate thermodynamic and high resolution structural data are available, and the results have been compared with the x-ray crystallographic observations. Analyses of the available results show poor correlation between the thermodynamic and structural parameters. Probable reasons for this discrepancy are mostly related with the reorganization of water accompanying the reaction which is indeed proven by the analyses of the energetics of the binding of the wheat germ agglutinin to oligosaccharides.

Theoretical Studies on Coupled Electron and Proton Transfer in Cytochrome c Oxidase

The complexity of life is based on an effective energy transduction machinery, which has evolved during the last 3.5 billion years. In aerobic life, the utilization of the high oxidizing potential of molecular oxygen powers this machinery. Oxygen is safely reduced by a membrane bound enzyme, cytochrome c oxidase (CcO), to produce an electrochemical proton gradient over the mitochondrial or bacterial membrane. This gradient is used for energy-requiring reactions such as synthesis of ATP by F0F1-ATPase and active transport. In this thesis, the molecular mechanism by which CcO couples the oxygen reduction chemistry to proton-pumping has been studied by theoretical computer simulations. By building both classical and quantum mechanical model systems based on the X-ray structure of CcO from Bos taurus, the dynamics and energetics of the system were studied in different intermediate states of the enzyme. As a result of this work, a mechanism was suggested by which CcO can prevent protons from leaking backwards in proton-pumping. The use and activation of two proton conducting channels were also enlightened together with a mechanism by which CcO sorts the chemical protons from pumped protons. The latter problem is referred to as the gating mechanism of CcO, and has remained a challenge in the bioenergetics field for more than three decades. Furthermore, a new method for deriving charge parameters for classical simulations of complex metalloenzymes was developed.

Description of accretion induced outflows from ultra-luminous sources to under-luminous AGNs

We study the energetics of the accretion-induced outflow and then plausible jet around black holes/compact objects using a newly developed disc-outflow coupled model. Inter-connecting dynamics of outflow and accretion essentially upholds the conservation laws. The energetics depend strongly on the viscosity parameter α and the cooling factor f which exhibit several interesting features. The bolometric luminosities of ultra-luminous X-ray binaries (e.g. SS433) and family of highly luminous AGNs and quasars can be reproduced by the model under the super-Eddington accretion flows. Under appropriate conditions, low-luminous AGNs (e.g. Sagittarius A*) also fit reasonably well with the luminosity corresponding to a sub-Eddington accretion flow with f→1.

Titration calorimetric studies to elucidate the specificity of the interactions of polymyxin B with lipopolysaccharides and lipid A

Lipopolysaccharide (LPS), the major cell wall constituent of Gram-negative bacteria, evokes a multitude of biological effects in mammals including pyrogenicity and toxic shock syndrome. Polymyxin B (PmB), a polycationic cyclic peptide, is known to neutralize most of its activities. The nature of the interaction of PmB with LPS and lipid A was investigated by isothermal titration calorimetry. PmB binds to LPS as well as lipid A stoichiometrically and non-co-operatively with micromolar affinity. These interactions are driven primarily by a favourable change in entropy (delta S) and are endothermic in nature. These positive changes in enthalpies decrease with increasing temperature, yielding a heat capacity change, delta Cp, of -2385 J.mol-1.degree-1 for PmB-LPS interactions while the binding of PmB to lipid A displays a delta Cp of -2259 J.mol-1.degree-1. The negative heat capacity changes provide strong evidence for the role of hydrophobic interactions as the driving force for the association of PmB with LPS and lipid A. A correlation of the energetics of these interactions with analyses of the molecular models of PmB suggests that a cluster of solvent-exposed non-polar amino acid side-chains that line one surface of the molecule, together with a ring of positively charged residues on its other surface, are responsible for its strong and stoichiometric binding to LPS.

Replacement of water on a steel substrate by a surrounding oil reservoir

With an objective to replace a water droplet from a steel surface by oil we study here the impact of injecting a hydrophilic/lipophilic surfactant into the droplet or into the surrounding oil reservoir. Contact angle goniometery, Grazing angle FTIR spectroscopy and Atomic force microscopy are used to record the oil/water interfacial tension, surface energetics of the substrate under the oil and water phases as well as the corresponding physical states of the substrates. Such energetics reflect the rate at which the excess surfactant molecules accumulate at the water/oil interface and desorb into the phases. The molecules diffuse into the substrate from the phases and build up specific molecular configurations which, with the interfacial tension, control the non-equilibrium progress of and the equilibrium status of the contact line. The study shows that the most efficient replacement of water by the surrounding oil happens when a surfactant is sparingly soluble in the supplier oil phase and highly soluble in the recipient water phase.

Structure and Bonding in Cyclic Isomers of B2AlHnm (n = 3−6, m = −2 to +1): A Comparative Study with B3Hnm, BAl2Hnm and Al3Hnm†

The structure, bonding and energetics of B2AlHnm (n = 3−6, m = −2 to +1) are compared with corresponding homocyclic boron, aluminum analogues and BAl2Hnm using density functional theory (DFT). Divalent to hexacoordinated boron and aluminum atoms are found in these species. The geometrical and bonding pattern in B2AlH4− is similar to that for B2SiH4. Species with lone pairs on the divalent boron and aluminum atoms are found to be minima on the potential energy surface of B2AlH32−. A dramatic structural diversity is observed in going from B3Hnm to B2AlHnm, BAl2Hnm and Al3Hnm and this is attributable to the preference of lower coordination on aluminum, higher coordination on boron and the higher multicenter bonding capability of boron. The most stable structures of B3H6+, B2AlH5 and BAl2H4− and the trihydrogen bridged structure of Al3H32− show an isostructural relationship, indicating the isolobal analogy between trivalent boron and divalent aluminum anion.

Theory of the unusual doping and temperature dependence of photoemission spectra in manganites

A recent, major, puzzle in the core-level photoemission spectra of doped manganites is the observation of a 1–2 eV wide shoulder with intensity varying with temperature T as the square of the magnetization over a T scale of order 200 K, an order of magnitude less than electronic energies. This is addressed and resolved here, by extending a recently proposed two-fluid polaron–mobile electron model for these systems to include core-hole effects. The position of the shoulder is found to be determined by Coulomb and Jahn-Teller energies, while its spectral weight is determined by the mobile electron energetics which is strongly T and doping dependent, due to annealed disorder scattering from the polarons and the t2g core spins. Our theory accounts quantitatively for the observed T dependence of the difference spectra, and furthermore, explains the observed correspondence between spectral changes due to increasing doping and decreasing T.

Multiscale study on hydrogen mobility in metallic fusion divertor material

For achieving efficient fusion energy production, the plasma-facing wall materials of the fusion reactor should ensure long time operation. In the next step fusion device, ITER, the first wall region facing the highest heat and particle load, i.e. the divertor area, will mainly consist of tiles based on tungsten. During the reactor operation, the tungsten material is slowly but inevitably saturated with tritium. Tritium is the relatively short-lived hydrogen isotope used in the fusion reaction. The amount of tritium retained in the wall materials should be minimized and its recycling back to the plasma must be unrestrained, otherwise it cannot be used for fueling the plasma. A very expensive and thus economically not viable solution is to replace the first walls quite often. A better solution is to heat the walls to temperatures where tritium is released. Unfortunately, the exact mechanisms of hydrogen release in tungsten are not known. In this thesis both experimental and computational methods have been used for studying the release and retention of hydrogen in tungsten. The experimental work consists of hydrogen implantations into pure polycrystalline tungsten, the determination of the hydrogen concentrations using ion beam analyses (IBA) and monitoring the out-diffused hydrogen gas with thermodesorption spectrometry (TDS) as the tungsten samples are heated at elevated temperatures. Combining IBA methods with TDS, the retained amount of hydrogen is obtained as well as the temperatures needed for the hydrogen release. With computational methods the hydrogen-defect interactions and implantation-induced irradiation damage can be examined at the atomic level. The method of multiscale modelling combines the results obtained from computational methodologies applicable at different length and time scales. Electron density functional theory calculations were used for determining the energetics of the elementary processes of hydrogen in tungsten, such as diffusivity and trapping to vacancies and surfaces. Results from the energetics of pure tungsten defects were used in the development of an classical bond-order potential for describing the tungsten defects to be used in molecular dynamics simulations. The developed potential was utilized in determination of the defect clustering and annihilation properties. These results were further employed in binary collision and rate theory calculations to determine the evolution of large defect clusters that trap hydrogen in the course of implantation. The computational results for the defect and trapped hydrogen concentrations were successfully compared with the experimental results. With the aforedescribed multiscale analysis the experimental results within this thesis and found in the literature were explained both quantitatively and qualitatively.

Preparation, characterization and thermal decomposition of praseodymium, terbium and neodymium carbonates

The preparation of three different types of carbonates of praseodymium, neodymium and terbium has been described. The carbonates have been characterized by potentiometry, chemical analysis, X-ray crystallography, infra-red spectroscopy and by their thermal behaviour. The thermal decomposition of several carbonates has been studied exhaustively under a variety of conditions and the stoicheiometry, thermodynamics and energetics of the reactions at various stages of decomposition have been examined. The stoicheiometry of the oxides obtained as final products of decomposition has been examined.

Studies on the brookite-rutile transformation

Brookite, the orthorhombic modification of titanium dioxide, transforms to the tetragonal modification, rutile, on heating. The kinetics and energetics of the transformation have been studied. Below 715±10°C, the rate of transformation is extremely slow. There appears to be little or no induction time. The kinetic data can be fitted reasonably well by the first-order equation. The energy of activation is about 60 kcal/mole and the frequency factor is of the order of 1013 h-1. The entropy of activation from Eyring's theory is about -18 cal/mole deg. at 800°C. The heat of this transformation is -100±75 cal/mole. The kinetic results may be explained qualitatively in terms of various analogies but more clearly by the application of the order-disorder theory to diffusionless transformation in solids. It has been shown that the ratio of propagation rate constant to the nucleation rate constant is small and that there is little or negligible phase aggregation.

InMnO3: A biferroic

InMnO3 which has a hexagonal structure similar to that of YMnO3 is found to show a canted antiferromagnetic behavior below 50 K (T-N) and a ferroelectric (FE) transition at 500 K accompanied by hystersis. We have determined the structure, polarization, and energetics of the FE and paraelectric (PE) phases of InMnO3 using first-principles density functional theory calculations based on pseudopotentials and a plane-wave basis, and find the polarization of the PE phase to be a half-integer quantum. The difference in polarization of the FE and PE phases calculated along a simple path is different from the absolute value of polarization of the FE phase. A weak piezoelectric response is exhibited by InMnO3 to uniaxial strain.

Elucidation of the mechanism of interaction of sheep spleen galectin-1 with splenocytes and its role in cell-matrix adhesion

The binding of a 14 kDa beta-galactoside animal lectin to splenocytes has been studied in detail. The binding data show that there are two classes of binding sites on the cells for the lectin: a high-affinity site with a K-a ranging from 1.1 x 10(6) to 5.1 x 10(5) M-1 and a low affinity binding site with a K-a ranging from 7.7 x 10(4) to 3.4 x 10(4) M-1 The number of receptors per cell for the high- and low-affinity sites is 9 +/- 3 x 10(6) and 2.5 +/- 0.5 x 10(6) respectively. The temperature dependence of the K value yielded the thermodynamic parameters. The energetics of this interaction shows that, although this interaction is essentially enthalpically driven (Delta H - 21 kJ lambda mol(-1)) for the high-affinity sites, there is a very favorable entropy contribution to the free energy of this interaction (-T Delta S - 17.5 Jmol(-1)), suggesting that hydrophobic interaction may also be playing a role in this interaction. Lactose brought about a 20% inhibition of this interaction, whereas the glycoprotein asialofetuin brought about a 75 % inhibition, suggesting that complex carbohydrate structures are involved in the binding of galectin-1 to splenocytes, Galectin-1 also mediated the binding and adhesion of splenocytes to the extracellular matrix glycoprotein laminin, suggesting a role for it in cell-matrix interactions. Copyright (C) 2000 John Wiley & Sons, Ltd.

Asynchronous behaviour in some three-atom linear concerted radical-exchange reactions

A few simple three-atom thermoneutral radical exchange reactions (i.e. A + BC --> AB + C) are examined by ab initio SCF methods. Emphasis is laid on the detailed analysis of density matrices rather than on energetics. Results reveal that the sum of the bond orders of the breaking and forming bonds is not conserved to unity, due to development of free valence on the migrating atom 'B' in the transition state. Bond orders, free valence and spin densities on the atoms are calculated. The present analysis shows that the bond-cleavage process is always more advanced than the bond-formation process in the transition state. Further analysis shows a development of the negative spin density on the migrating atom 'B' in the transition state. The depletion of the alpha-spin density on the radical site "A" in the reactant during the reaction lags behind the growth of the alpha-spin density on the terminal atom "C" of the reactant bond, 'B-C' in the transition state. But all these processes are completed simultaneously at the end of the reaction. Hence, the reactions are asynchronous but kinetically concerted in most cases.

Reducibility of Ce1-xZrxO2: origin of enhanced oxygen storage capacity

We combine first-principles calculations with EXAFS studies to investigate the origin of high oxygen storage capacity in ceria-zirconia solid solution, prepared by solution combustion method. We find that nanocrystalline Ce0.5Zr0.5O2 can be reduced to Ce0.5Zr0.5O1.57 by H-2 upto 850 degrees C with an OSC of 65 cc/gm which is extremely high. Calculated local atomic-scale structure reveals the presence of long and short bonds resulting in four-fold coordination of the cations, confirmed by the EXAFS studies. Bond valence analysis of the microscopic structure and energetics is used to evaluate the strength of binding of different oxide ions and vacancies. We find the presence of strongly and weakly bound oxygens, of which the latter are responsible for the higher oxygen storage capacity in the mixed oxides than in the pure CeO2.

Identification of Critical Residues of the Mycobacterial Dephosphocoenzyme A Kinase by Site-Directed Mutagenesis

Dephosphocoenzyme A kinase performs the transfer of the c-phosphate of ATP to dephosphocoenzyme A, catalyzing the last step of coenzyme A biosynthesis. This enzyme belongs to the P-loop-containing NTP hydrolase superfamily, all members of which posses a three domain topology consisting of a CoA domain that binds the acceptor substrate, the nucleotide binding domain and the lid domain. Differences in the enzymatic organization and regulation between the human and mycobacterial counterparts, have pointed out the tubercular CoaE as a high confidence drug target (HAMAP database). Unfortunately the absence of a three-dimensional crystal structure of the enzyme, either alone or complexed with either of its substrates/regulators, leaves both the reaction mechanism unidentified and the chief players involved in substrate binding, stabilization and catalysis unknown. Based on homology modeling and sequence analysis, we chose residues in the three functional domains of the enzyme to assess their contributions to ligand binding and catalysis using site-directed mutagenesis. Systematically mutating the residues from the P-loop and the nucleotide-binding site identified Lys14 and Arg140 in ATP binding and the stabilization of the phosphoryl intermediate during the phosphotransfer reaction. Mutagenesis of Asp32 and Arg140 showed catalytic efficiencies less than 5-10% of the wild type, indicating the pivotal roles played by these residues in catalysis. Non-conservative substitution of the Leu114 residue identifies this leucine as the critical residue from the hydrophobic cleft involved in leading substrate, DCoA binding. We show that the mycobacterial enzyme requires the Mg2+ for its catalytic activity. The binding energetics of the interactions of the mutant enzymes with the substrates were characterized in terms of their enthalpic and entropic contributions by ITC, providing a complete picture of the effects of the mutations on activity. The properties of mutants defective in substrate recognition were consistent with the ordered sequential mechanism of substrate addition for CoaE.