Chemical unfolding of chicken villin headpiece in aqueous dimethyl sulfoxide solution: cosolvent concentration dependence, pathway, and microscopic mechanism

Unfolding of a protein often proceeds through partial unfolded intermediate states (PUIS). PUIS have been detected in several experimental and simulation studies. However, complete analyses of transitions between different PUIS and the unfolding trajectory are sparse. To understand such dynamical processes, we study chemical unfolding of a small protein, chicken villin head piece (HP-36), in aqueous dimethyl sulfoxide (DMSO) solution. We carry out molecular dynamics simulations at various solution compositions under ambient conditions. In each concentration, the initial step of unfolding involves separation of two adjacent native contacts, between phenyl alanine residues (11-18 and 7-18). This first step induces, under appropriate conditions, subsequent separation among other hydrophobic contacts, signifying a high degree of cooperativity in the unfolding process. The observed sequence of structural changes in HP-36 on increasing DMSO concentration and the observed sequence of PUIS, are in approximate agreement with earlier simulation results (in pure water) and experimental observations on unfolding of HP-36. Peculiar to water-DMSO mixture, an intervening structural transformation (around 15% of DMSO) in the binary mixture solvent retards the progression of unfolding as composition is increased. This is reflected in a remarkable nonmonotonic composition dependence of RMSD, radius of gyration and the fraction of native contacts. At 30% mole fraction of DMSO, we find the extended randomly coiled structure of the unfolded protein. The molecular mechanism of DMSO induced unfolding process is attributed to the initial preferential solvation of the hydrophobic side chain atoms through the methyl groups of DMSO, followed by the hydrogen bonding of the oxygen atom of DMSO to the exposed backbone NH groups of HP-36.

Molecular Origin of the Intrinsic Bending Force for Helical Morphology Observed in Chiral Amphiphilic Assemblies: Concentration and Size Dependence

The formation of the helical morphology in monolayers and bilayers of chiral amphiphilic assemblies is believed to be driven at least partly by the interactions at the chiral centers of the amphiphiles. However, a detailed microscopic understanding of these interactions and their relation with the helix formation is still not clear. In this article a study of the molecular origin of the chirality-driven helix formation is presented by calculating, for the first time, the effective pair potential between a pair of chiral molecules. This effective potential depends on the relative sizes of the groups attached to the two chiral centers, on the orientation of the amphiphile molecules, and also on the distance between them. We find that for the mirror-image isomers (in the racemic modification) the minimum energy conformation is a nearly parallel alignment of the molecules. On the other hand, the same for a pair of molecules of one kind of enantiomer favors a tilt angle between them, thus leading to the formation of a helical morphology of the aggregate. The tilt angle is determined by the size of the groups attached to the chiral centers of the pair of molecules considered and in many cases predicted it to be close to 45 degrees. The present study, therefore, provides a molecular origin of the intrinsic bending force, suggested by Helfrich (J. Chem. Phys. 1986, 85, 1085-1087), to be responsible for the formation of helical structure. This effective potential may explain many of the existing experimental results, such as the size and the concentration dependence of the formation of helical morphology. It is further found that the elastic forces can significantly modify the pitch predicted by the chiral interactions alone and that the modified real pitch is close to the experimentally observed value. The present study is expected to provide a starting point for future microscopic studies.

Free Energy Gap Dependence of the Electron-Transfer Rate from the Inverted to the Normal Region

numerical study of the free energy gap (FEG) dependence of the electron-transfer rate in polar solvents is presented. This study is based on the generalized multidimensional hybrid model, which not only includes the solvent polarization and the molecular vibration modes, but also the biphasic polar response of the solvent. The free energy gap dependence is found to be sensitive to several factors, including the solvent relaxation rate, the electronic coupling between the surfaces, the frequency of the high-frequency quantum vibrational mode, and the magnitude of the solvent reorganization energy. It is shown that in some cases solvent relaxation can play an important role even in the Marcus normal regime. The minimal hybrid model involves a large number of parameters, giving rise to a diverse non-Marcus FEG behavior which is often determined collectively by these parameters. The model gives the linear free energy gap dependence of the logarithmic rate over a substantial range of FEG, spanning from the normal to the inverted regime. However, even for favorable values of the relevant parameters, a linear free energy gap dependence of the rate could be obtained only over a range of 5000-6000 cm(-1) (compared to the experimentally observed range of 10000 cm(-1) reported by Benniston et al.). The present work suggests several extensions/generalizations of the hybrid model which might be necessary to fully understand the observed free energy gap dependence.

Long range resonance energy transfer from a dye molecule to graphene has (distance)-4 dependence

In our previous report on resonance energy transfer from a dye molecule to graphene [J. Chem. Phys.129, 054703 (2008)], we had derived an expression for the rate of energy transfer from a dye to graphene. An integral in the expression for the rate was evaluated approximately. We found a Yuwaka-type dependence of the rate on the distance. We now present an exact evaluation of the integral involved, leading to very interesting results. For short distances (z < 20 A), the present rate and the previous rate are in good agreement. For larger distances, the rate is found to have a z(-4) dependence on the distance, exactly. Thus we predict that for the case of pyrene on graphene, it is possible to observe fluorescence quenching up to a distance of 300 A. This is in sharp contrast to the traditional fluorescence resonance energy transfer where the quenching is observable only up to 100 A.

Dependence of metal Auger intensity ratios on the effective charge on metal atoms

Metal Auger line intensity ratios were shown by Rao and others to be directly related to the occupancy of valence states. It is now shown that these intensity ratios are more generally related to the effective charge on the metal atom. The Auger intensity ratios are also directly proportional to valence band intensities of metals. Correlations of the intensity ratios with Auger parameters have also been examined.

Kinetics of Chemical Vapor Deposition

A generalized mass transport model is developed for predicting the rate ofdeposition in chemical vapor deposition (CVD) systems. This combines thegeneralized method of obtaining equilibrium composition, with elemental fluxbalance expressions. This procedure avoids the usual problems encountered incalculating the rates in multicomponent systems, like writing overall reactionschemes. The dependence of multicomponent diffusivities on the fluxes is accountedin this model using an iterative procedure. The model developed isapplied to the deposition of titanium carbide on cemented carbide tool bitsfrom a gas mixture of titanium tetrachloride, toluene, and hydrogen. Experimentaldeposition rates were obtained using a thermogravimetric assembly.Mass transport controlled rates give an order of magnitude estimates of theobserved rates.

Size dependence and dispersion behavior of the first hyperpolarizability of copper nanoparticles

We have probed the size dependency of the first hyperpolarizability (b) of copper nanoparticles by hyper-Rayleigh scattering (HRS). Our results indicate that second harmonic generation (SHG) originates predominantly at the surface of the nanoparticles as long as the size (d) remains small compared to the wavelength (k). However, volume contribution to the SH response due to the retardation effect becomes important when particle size grows beyond the `small particle limit'. There is a significant dispersion in the b values of copper nanoparticles owing tothe presence of the strong surface plasmon resonance (SPR) band.

Fractional power dependence of rate on activation energy for reactions with very low internal barriers

The dynamics of reactions with low internal barriers are studied both analytically and numerically for two different models. Exact expressions for the average rate,kI, are obtained by solving the associated first passage time problems. Both the average rate constant, kI, and the numerically calculated long-time rate constant, kL, show a fractional power law dependence on the barrier height for very low barriers. The crossover of the reaction dynamics from low to high barrier is investigated.

A chitotetrose specific lectin from Luffa acutangula :Physico-chemical properties and the assignment of orientation of sugars in the lectin binding site

A chitooligosaccharide specific lectin (Luffa acutangula agglutinin) has been purified from the exudate of ridge gourd fruits by affinity chromatography on soybean agglutininglycopeptides coupled to Sepharose-6B. The affinity purified lectin was found homogeneous by polyacrylamide gel electrophoresis, in sodium dodecyl sulphate-polyacrylamide gels, by gel filtration on Sephadex G-100 and by sedimentation velocity experiments. The relative molecular weight of this lectin is determined to be 48,000 ± 1,000 by gel chromatography and sedimentation equilibrium experiments. The sedimentation coefficient (S20, w) was obtained to be 4·06 S. The Stokes’ radius of the protein was found to be 2·9 nm by gel filtration. In sodium dodecyl sulphate-polyacrylamide gel electrophoresis the lectin gave a molecular weight of 24,000 in the presence as well as absence of 2-mercaptoethanol. The subunits in this dimeric lectin are therefore held by non-covalent interactions alone. The lectin is not a glycoprotein and circular dichroism spectral studies indicate that this lectin has 31% α-helix and no ß-sheet. The lectin is found to bind specifically to chitooligosaccharides and the affinity of the lectin increases with increasing oligosaccharide chain length as monitored by near ultra-violetcircular dichroism and intrinsic fluorescence titration. The values of ΔG, ΔΗ and ΔS for the binding process showed a pronounced dependence on the size of the oligosaccharide. The values for both ΔΗ and ΔS show a significant increase with increase in the oligosaccharide chain length showing that the binding of higher oligomers is progressively more favoured thermodynamically than chitobiose itself. The thermodynamic data is consistent with an extended binding site in the lectin which accommodates a tetrasaccharide. Based on the thermodynamic data, blue shifts and fluorescence enhancement, spatial orientation of chitooligosaccharides in the combining site of the lectin is assigned.

Nature of dependence of spin-orbit splittings on atomic number

It has been shown that Dirac equation employing a constant value of the screening constant Z0 does not explain the variation of spin-orbit splittings of 2p and 3p levels with atomic number Z. A model which takes into account the variation of Z0 withZ is shown to satisfactorily predict the dependence of spinorbit splittings onZ.

Distance dependence of fluorescence resonance energy transfer

Deviations from the usual R (-6) dependence of the rate of fluorescence resonance energy transfer (FRET) on the distance between the donor and the acceptor have been a common scenario in the recent times. In this paper, we present a critical analysis of the distance dependence of FRET, and try to illustrate the non R (-6) type behaviour of the rate for the case of transfer from a localized electronic excitation on the donor, a dye molecule to three different energy acceptors with delocalized electronic excitations namely, graphene,two-dimensional semiconducting sheet and the case of such a semiconducting sheet rolled to obtain a nanotube. We use simple analytic models to understand the distance dependence in each case.

Aggregation of apolar peptides in organic solvents. Concentration dependence of 1H-nmr parameters for peptide NH groups in 310 helical decapeptide fragment of suzukacillin

Peptide NH chemical shifts and their temperature dependences have been monitored as a function of concentration for the decapeptide, Boc-Aib-Pro-Val-Aib-Val-Ala-Aib-Ala-Aib-Aib-OMe in CDCl3 (0.001-0.06M) and (CD3)2SO (0.001-0.03M). The chemical shifts and temperature coefficients for all nine NH groups show no significant concentration dependence in (CD3)2SO. Seven NH groups yield low values of temperature coefficients over the entire range, while one yields an intermediate value. In CDCl3, the Aib(1) NH group shows a large concentration dependence of both chemical shift and temperature coefficient, in contrast to the other eight NH groups. The data suggest that in (CD3)2SO, the peptide adopts a 310 helical conformation and is monomeric over the entire concentration range. In CDCl3, the 310 helical peptide associates at a concentration of 0.01M, with the Aib(1) NH involved in an intermolecular hydrogen bond. Association does not disrupt the intramolecular hydrogen-bonding pattern in the decapeptide.

Diffusion of flexible, charged, nanoscopic molecules in solution: Size and pH dependence for PAMAM dendrimer

In order to understand self-diffusion (D) of a charged, flexible, and porous nanoscopic molecule in water, we carry out very long, fully atomistic molecular dynamics simulation of PAMAM dendrimer up to eight generations in explicit salt water under varying pH. We find that while the radius of gyration (R-g) varies as N-1/3, the self-diffusion constant (D) scales, surprisingly, as N-alpha, with alpha=0.39 at high pH and 0.5 at neutral pH, indicating a dramatic breakdown of Stokes-Einstein relation for diffusion of charged nanoscopic molecules. The variation in D as a function of radius of gyration demonstrates the importance of treating water and ions explicitly in the diffusion process of a flexible nanoscopic molecule. In agreement with recent experiments, the self-diffusion constant increases with pH, revealing the importance of dielectric friction in the diffusion process. The shape of a dendrimer is found to fluctuate on a nanosecond time scale. We argue that this flexibility (and also the porosity) of the dendrimer may play an important role in determining the mean square displacement of the dendrimer and the breakdown of the Stokes-Einstein relation between diffusion constant and the radius.

Dynamics of Barrierless and Activated Chemical Reactions in a Dispersive Medium within the Fractional Diffusion Equation Approach

Barrierless chemical reactions have often been modeled as a Brownian motion on a one-dimensional harmonic potential energy surface with a position-dependent reaction sink or window located near the minimum of the surface. This simple (but highly successful) description leads to a nonexponential survival probability only at small to intermediate times but exponential decay in the long-time limit. However, in several reactive events involving proteins and glasses, the reactions are found to exhibit a strongly nonexponential (power law) decay kinetics even in the long time. In order to address such reactions, here, we introduce a model of barrierless chemical reaction where the motion along the reaction coordinate sustains dispersive diffusion. A complete analytical solution of the model can be obtained only in the frequency domain, but an asymptotic solution is obtained in the limit of long time. In this case, the asymptotic long-time decay of the survival probability is a power law of the Mittag−Leffler functional form. When the barrier height is increased, the decay of the survival probability still remains nonexponential, in contrast to the ordinary Brownian motion case where the rate is given by the Smoluchowski limit of the well-known Kramers' expression. Interestingly, the reaction under dispersive diffusion is shown to exhibit strong dependence on the initial state of the system, thus predicting a strong dependence on the excitation wavelength for photoisomerization reactions in a dispersive medium. The theory also predicts a fractional viscosity dependence of the rate, which is often observed in the reactions occurring in complex environments.

Subquadratic wavenumber dependence of the structural relaxation of supercooled liquid in the crossover regime

As a liquid is progressively supercooled toward its glass transition temperature, an intriguing weakening of the wavenumber (q) dependence of the structural relaxation time tau(q) in the intermediate-to-large q limit is observed both in experiments and simulation studies. Neither continuous Brownian diffusive dynamics nor discontinuous activated events can alone explain the anomalous wavenumber dependence. Here we use our recently developed theory that unifies the mode coupling theory for continuous dynamics, with the random first order transition theory treatment of activated discontinuous motion as a nucleationlike instanton process to understand the wavenumber dependence of density relaxation. The predicted smooth change in mechanism of relaxation from diffusive to activated, in the crossover regime, is wavevector dependent and appears to be responsible for the observed subquadratic,nalmost linear, q dependence of the relaxation time.