Thermodynamic and structural studies of cavity formation in proteins suggest that loss of packing interactions rather than the hydrophobic effect dominates the observed energetics

The hydrophobic effect is widely believed to be an important determinant of protein stability. However, it is difficult to obtain unambiguous experimental estimates of the contribution of the hydrophobic driving force to the overall free energy of folding. Thermodynamic and structural studies of large to small substitutions in proteins are the most direct method of measuring this contribution. We have substituted the buried residue Phe8 in RNase S with alanine, methionine, and norleucine, Binding thermodynamics and structures were characterized by titration calorimetry and crystallography, respectively. The crystal structures of the RNase S F8A, F8M, and F8Nle mutants indicate that the protein tolerates the changes without any main chain adjustments, The correlation of structural and thermodynamic parameters associated with large to small substitutions was analyzed for nine mutants of RNase S as well as 32 additional cavity-containing mutants of T4 lysozyme, human lysozyme, and barnase. Such substitutions were typically found to result in negligible changes in Delta C-p and positive values of both Delta Delta H degrees and aas of folding. Enthalpic effects were dominant, and the sign of Delta Delta S is the opposite of that expected from the hydrophobic effect. Values of Delta Delta G degrees and Delta Delta H degrees correlated better with changes in packing parameters such as residue depth or occluded surface than with the change in accessible surface area upon folding. These results suggest that the loss of packing interactions rather than the hydrophobic effect is a dominant contributor to the observed energetics for large to small substitutions. Hence, estimates of the magnitude of the hydrophobic driving force derived from earlier mutational studies are likely to be significantly in excess of the actual value.

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.

Electronic structure of square planar CuO46− clusters

The results of spin-polarized MSXagr calculations show that the ground state of the CuO 4 6– cluster is essentially non-magnetic in spite of odd number of electrons in the system for short Cu–O distances (1.90 Å) as found in the highT c superconductors. This arises due to the fact that the unpaired electron resides in a molecular orbital with primarily oxygen 3s character. The stability of this molecular orbital is found to be sensitive to the cluster geometry and thus, increase in Cu–O distance (as well as other changes affecting oxygen-oxygen distance) tend to favour a magnetic state. From these calculations we have also estimated the Coulomb correlation strength within the Cu 3d to be about 5.3 eV.

Variational calculation for the spin-(1/2 Heisenberg antiferromagnet on a square lattice

The ground-state properties of the spin-(1/2 Heisenberg antiferromagnet on a square lattice are studied by using a simple variational wave function that interpolates continuously between the Néel state and short-range resonating-valence-bond states. Exact calculations of the variational energy for small systems show that the state with the lowest energy has long-range antiferromagnetic order. The staggered magnetization in this state is approximately 70% of its maximum possible value. The variational estimate of the ground-state energy is substantially lower than the value obtained for the nearest-neighbor resonating-valence-bond wave function.

Mean square displacements of hydrogen interstitial and its neighbours in the intermetallic compound Fe0.5Ti0.5

A Green's function technique is used in the scattering matrix formalism to compute the mean square displacement of hydrogen and deuterium interstitials in the intermetallic compound Fe0.5Ti0.5 for low hydrogen/deuterium concentration. The mean square amplitudes of the metal atoms surrounding the interstitial are found to be smaller than those for the host crystal. This anomalous effect is due to the stiffening of the lattice by the dissolved hydrogen or deuterium at low concentration. This type of effect is experimentally observed in the case of NbHx at low hydrogen concentration.

Spiral states in the square-lattice Hubbard model

We present a variety of physical implications of a mean-field theory for spiral spin-density-wave states in the square-lattice Hubbard model for small deviations from half filling. The phase diagram with the paramagnetic metal, two spiral (semimetallic) states, and ferromagnet is calculated. The momentum distribution function and the (quasiparticle) density of states are discussed. There is a significant broadening of the quasiparticle bands when the antiferromagnetic insulator is doped. The evolution of the Fermi surface and the variation of the plasma frequency and a charge-stiffness constant with U/t and δ are calculated. The connection to results based on the Schwinger-boson-slave-fermion formalism is made.

A nanosecond molecular dynamics study of antiparallel d(G)(7) quadruplex structures: Effect of the coordinated cations

Nanosecond scale molecular dynamics simulations have been performed on antiparallel Greek key type d(G(7)) quadruplex structures with different coordinated ions, namely Na+ and K+ ion, water and Na+ counter ions, using the AMBER force field and Particle Mesh Ewald technique for electrostatic interactions. Antiparallel structures are stable during the simulation, with root mean square deviation values of similar to1.5 Angstrom from the initial structures. Hydrogen bonding patterns within the G-tetrads depend on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate different cations. However, alternating syn-anti arrangement of bases along a chain as well as in a quartet is maintained through out the MD simulation. Coordinated Na+ ions, within the quadruplex cavity are quite mobile within the central channel and can even enter or exit from the quadruplex core, whereas coordinated K+ ions are quite immobile. MD studies at 400 K indicate that K+ ion cannot come out from the quadruplex core without breaking the terminal G-tetrads. Smaller grooves in antiparallel structures are better binding sites for hydrated counter ions, while a string of hydrogen bonded water molecules are observed within both the small and large grooves. The hydration free energy for the K+ ion coordinated structure is more favourable than that for the Na+ ion coordinated antiparallel quadruplex structure.

Erosion and cavity characteristics in rotating components

An investigation of the initiation and growth of erosion and of the effect of velocity and pressure on erosion in a rotating disk is presented. Also, the role of an intervening noncavitating period on erosion is studied. The results indicate that at high intensities the peak rate of erosion decreases with increases in pressure. The erosion rate/time curves obtained for metallic materials are explained by the eroded particle distribution and the cavity size. The average size of the eroded particles decreased when pressure and tensile strength of the material were increased. The erosion rate peaked after an intervening noncavitating period. The use of the rate of erosion, defined as an average over the entire test duration, in the equation governing the theory of erosion resulted in reasonably good correlations. The correlations reveal that it is possible to predict the length, width, and area of a cavity when the cavitation parameter σ is known. The normalized width of a cavity may be estimated if its normalized length is known.

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

Mean-field results of the multiple-band extended Hubbard model for the square-planar CuO2 lattice

We obtain metal-insulator phase diagrams at half-filling for the five-band extended Hubbard model of the square-planar CuO2 lattice treated within a Hartree-Fock mean-field approximation, allowing for spiral spin-density waves. We indicate the existence of an insulating phase (covalent insulator) characterized by strong covalency effects, not identified in the earlier Zaanen-Sawatzky-Allen phase diagram. While the insulating phase is always antiferromagnetic, we also obtain an antiferromagnetic metallic phase for a certain range of interaction parameters. Performing a nonperturbative calculation of J(eff), the in-plane antiferromagnetic interaction is presented as a function of the parameters in the model. We also calculate the band gap and magnetic moments at various sites and discuss critically the contrasting interpretation of the electronic structure of high-T(c) materials arising from photoemission and neutron-scattering experiments.

Global solutions describing the collapse of a spherical or cylindrical cavity

The collapse of a spherical (cylindrical) cavity in air is studied analytically. The global solution for the entire domain between the sound front, separating the undisturbed and the disturbed gas, and the vacuum front is constructed in the form of infinite series in time with coefficients depending on an ldquoappropriaterdquo similarity variable. At timet=0+, the exact planar solution for a uniformly moving cavity is assumed to hold. The global analytic solution of this initial boundary value problem is found until the collapse time (=(gamma–1)/2) for gamma le 1+(2/(1+v)), wherev=1 for cylindrical geometry, andv=2 for spherical geometry. For higher values of gamma, the solution series diverge at timet — 2(beta–1)/ (v(1+beta)+(1–beta)2) where beta=2/(gamma–1). A close agreement is found in the prediction of qualitative features of analytic solution and numerical results of Thomaset al. [1].

Probing potential energy surfaces in confined systems: behavior of mean-square displacement in zeolites

Time evolution of mean-squared displacement based on molecular dynamics for a variety of adsorbate-zeolite systems is reported. Transition from ballistic to diffusive behavior is observed for all the systems. The transition times are found to be system dependent and show different types of dependence on temperature. Model calculations on a one-dimensional system are carried out which show that the characteristic length and transition times are dependent on the distance between the barriers, their heights, and temperature. In light of these findings, it is shown that it is possible to obtain valuable information about the average potential energy surface sampled under specific external conditions.

Solubilization of n-decane into n-decyl-polyoxyethylene surfactants: A cavity forming free energy based model

The present paper reports the results of a theoretical study of the forces and factors driving the solubilization of n-alkane solubilizates into the micellar core of some non-ionic surfactants, based on a micellar model which includes the cavity forming free energy as a component of micellization. The solubilizate is n-decane and the non-ionic surfactants considered are n-decyl-polyoxyethylene surfactants. The extent of solubilization, i.e. the mole fraction of the solubilizate within the core has been calculated. The results indicate that the incorporated solubilizate has more translational and rotational degrees of freedom as compared to those of the tail parts of the surfactants present in the core. This drives the total free energy of aggregation after solubilization into a more favourable direction. The results are in fair agreement with the experimental results.