Free energy calculations of the interactions of c-Jun-based synthetic peptides with the c-Fos protein

The c-Fos–c-Jun complex forms the activator protein 1 transcription factor, a therapeutic target in the treatment of cancer. Various synthetic peptides have been designed to try to selectively disrupt the interaction between c-Fos and c-Jun at its leucine zipper domain. To evaluate the binding affinity between these synthetic peptides and c-Fos, polarizable and nonpolarizable molecular dynamics (MD) simulations were conducted, and the resulting conformations were analyzed using the molecular mechanics generalized Born surface area (MM/GBSA) method to compute free energies of binding. In contrast to empirical and semiempirical approaches, the estimation of free energies of binding using a combination of MD simulations and the MM/GBSA approach takes into account dynamical properties such as conformational changes, as well as solvation effects and hydrophobic and hydrophilic interactions. The predicted binding affinities of the series of c-Jun-based peptides targeting the c-Fos peptide show good correlation with experimental melting temperatures. This provides the basis for the rational design of peptides based on internal, van der Waals, and electrostatic interactions.

Influence of Hydrophilic Surface Specificity on the Structural Properties of Confined Water

The influence of chemical specificity of hydrophilic surfaces on the structure of confined water in the subnanometer regime is investigated using grand canonical Monte Carlo Simulations. The structural variations for water confined between hydroxylated silica surfaces are contrasted with water confined between mica surfaces. Although both surfaces are hydrophilic, our Study shows that hydration of potassium ions on the mica surface has a strong influence on the water Structure and solvation force response of confined water. In contrast to the disrupted hydrogen bond network observed for water confined between Mica Surfaces, water between silica surfaces retains its hydrogen bond network displaying bulklike structural features down to surface separations as small as 0.45 nm. Hydrogen bonding of all invariant contact water layer with the surface silanol groups aids in maintaining a constant number of hydrogen bonds per water molecule for the silica surfaces. As a consequence water depletion and rearrangement upon decreasing confinement is a strong function of the hydrophilic surface specificity, particularly at smaller separations. An oscillatory solvation force response is only observed for water confined between Silica surfaces, and bulklike features are observed for both Surfaces above a surface separation of about 1.2 nm. We evaluate and contrast the water density, dipole moment distributions, pi pair correlation functions, and solvation forces as a function of the surface separation.

Understanding solvent effects on structure and reactivity of organic intermediates: a Raman study

The effect of solvent on chemical reactivity has generally been explained on the basis of the dielectric constant and viscosity. However a number of spectroscopic studies, including UV-VIS, IR and Raman, has led to numerous empirical parameters to define solvent effect based on either solvating ability or polarity scale. These parameters include solvent polarizability, dipolarity, Lewis acidity and Lewis basicity, E-T(30), pi*, alpha, beta etc. However, from a structural point of view, we can separate solvation as static and dynamic processes. The static solvation basically relates to stabilization of the molecular structure by the solvent to attain the equilibrium structure, both in the intermediate and ground state. Dynamic solvation relates to solvent reorganization-induced dynamics prior to the structural reorganization to reach the equilibrium state. In this paper, we present (a) structural distortions induced by the solvent due to preferential solvation of the triplet excited state, and (b) the importance of dynamic solvation induced by vibronic coupling (pseudo-Jahn-Teller coupling). The examples include the effect of solvent on structure and reactivity of excited states of 2,2,2-trifluoroacetophenone (TFA). Based on the comparison of time resolved resonance Raman (TR3) data of TFA and other substituted acetophenone systems, it was found that change in solvent polarity indeed results in electronic state switching and structural changes in the excited state, which explains the trend in reactivity. Further, a TR3 study of fluoranil (FA) in the triplet excited state in solvents of varying polarities indicates that the structure of FA in the triplet excited state is determined by vibronic coupling effects and thus distorted structure. These experimental results have been well supported by density functional theoretical computational studies.

Binding of nucleobases with single-walled carbon nanotubes: Theory and experiment

We report the binding energy of various nucleobases (guanine (G), adenine (A), thymine (T) and cytosine (C)) with (5,5) single-walled carbon nanotube (SWNT) calculated using first-principle Hartre–Fock method (HF) together with classical force field. The binding energy without including the solvation effects of water decreases in the order G>A>T>C. The inclusion of solvation energy changes the order of binding preference to be G>T>A>C. Using isothermal titration (micro) calorimetry experiments, we also show the relative binding affinity to be T>A>C, in agreement with our calculations.

Microcalorimetric indications for ligand binding as a function of the protein for galactoside-specific plant and avian lectins

The process cascade leading to the final accommodation of the carbohydrate ligand in the lectin’s binding site comprises enthalpic and entropic contributions of the binding partners and solvent molecules. With emphasis on lactose, N-acetyllactosamine, and thiodigalactoside as potent inhibitors of binding of galactoside-specific lectins, the question was addressed to what extent these parameters are affected as a function of the protein. The microcalorimetric study of carbohydrate association to the galectin from chicken liver (CG-16) and the agglutinin from Viscum album (VAA) revealed enthalpy–entropy compensation with evident protein type-dependent changes for N-acetyllactosamine. Reduction of the entropic penalty by differential flexibility of loops or side chains and/or solvation properties of the protein will have to be reckoned with to assign a molecular cause to protein type-dependent changes in thermodynamic parameters for lectins sharing the same monosaccharide specificity.

Interplay between multiple length and time scales in complex chemical systems

Processes in complex chemical systems, such as macromolecules, electrolytes, interfaces, micelles and enzymes, can span several orders of magnitude in length and time scales. The length and time scales of processes occurring over this broad time and space window are frequently coupled to give rise to the control necessary to ensure specificity and the uniqueness of the chemical phenomena. A combination of experimental, theoretical and computational techniques that can address a multiplicity of length and time scales is required in order to understand and predict structure and dynamics in such complex systems. This review highlights recent experimental developments that allow one to probe structure and dynamics at increasingly smaller length and time scales. The key theoretical approaches and computational strategies for integrating information across time-scales are discussed. The application of these ideas to understand phenomena in various areas, ranging from materials science to biology, is illustrated in the context of current developments in the areas of liquids and solvation, protein folding and aggregation and phase transitions, nucleation and self-assembly.

Influence of water and thermal history on ion transport in lithium salt-succinonitrile plastic crystalline electrolytes

Important issues of water and thermal history affecting ion transport in a representative plastic crystalline lithium salt electrolyte: succinonitrile (SN)-lithium perchlorate (LiClO4) are discussed here. Ionic conductivity of electrolytes with high lithium salt amounts (similar to 1 M) in SN at a particular temperature is known to be influenced both by the trans-gauche isomerism and ion association (solvation), the two most important intrinsic parameters of the plastic solvent. In the present study both water and thermal history influence SN and result in enhancement of ionic conductivity of 1 M LiClO4-SN electrolyte. Systematic observations reveal that the presence of water in varying amounts promote ion-pair dissociation in the electrolyte. While trace amounts (approximate to 1-15 ppm) do not affect the trans-gauche isomerism of SN, the presence of water in large amounts (approximate to 5500 ppm) submerges the plasticity of SN. Subjugating the electrolyte to different thermal protocol resulted in enhancement of trans concentration only. This is an interesting observation as it demonstrates a simple and effective procedure involving utilization of an optimized set of external parameters to decouple solvation from trans-gauche isomerism. Observations from the ionic conductivity of various samples were accounted by changes in signature isomer and ion-association bands in the mid-IR regime and also from plastic to normal crystal transition temperature peak obtained from thermal studies. (C) 2010 Elsevier B.V. All rights reserved.

Dielectric Relaxation Spectroscopy for Evaluation of the Influence of Solvent Dynamics on Ion Transport in Succinonitrile-Salt Plastic Crystalline Electrolytes

Influence of succinonitrile (SN) dynamics on ion transport in SN-lithium perchlorate (LiClO4) electrolytes is discussed here via dielectric relaxation spectroscopy. Dielectric relaxation spectroscopy (similar to 2 x 10(-3) Hz to 3 MHz) of SN and SN-LiClO4 was studied as a function of salt content (up to 7 mol % or 1 M) and temperature (-20 to +60 degrees C). Analyses of real and imaginary parts of permittivity convincingly reveal the influence Of trans gauche isomerism and solvent-salt association (solvation) effects on ion transport. The relaxation processes are highly dependent on the salt concentration and temperature. While pristine SN display only intrinsic dynamics (i.e., trans-gauche isomerism) which enhances with an increase in temperature, SN-LiClO4 electrolytes especially at high salt concentrations (similar to 0.04-1 M) show salt-induced relaxation processes. In the concentrated electrolytes, the intrinsic dynamics was observed to be a function of salt content, becoming faster with an increase in salt concentration. Deconvolution of the imaginary part of the permittivity spectra using Havriliak-Negami (HN) function show a relaxation process corresponding to the above phenomena. The permittivity data analyzed using HN and Kohlrausch-Williams-Watta (KWW) functions show non-Debye relaxation processes and enhancement in the trans phase (enhanced solvent dynamics) as a function of salt concentration and temperature.

Using Correlated Parameters for Improved Ranking of Protein-Protein Docking Decoys

A successful protein-protein docking study culminates in identification of decoys at top ranks with near-native quaternary structures. However, this task remains enigmatic because no generalized scoring functions exist that effectively infer decoys according to the similarity to near-native quaternary structures. Difficulties arise because of the highly irregular nature of the protein surface and the significant variation of the nonbonding and solvation energies based on the chemical composition of the protein-protein interface. In this work, we describe a novel method combining an interface-size filter, a regression model for geometric compatibility (based on two correlated surface and packing parameters), and normalized interaction energy (calculated from correlated nonbonded and solvation energies), to effectively rank decoys from a set of 10,000 decoys. Tests on 30 unbound binary protein-protein complexes show that in 16 cases we can identify at least one decoy in top three ranks having <= 10 angstrom backbone root mean square deviation from true binding geometry. Comparisons with other state-of-art methods confirm the improved ranking power of our method without the use of any experiment-guided restraints, evolutionary information, statistical propensities, or modified interaction energy equations. Tests on 118 less-difficult bound binary protein-protein complexes with <= 35% sequence redundancy at the interface showed that in 77% cases, at least 1 in 10,000 decoys were identified with <= 5 angstrom backbone root mean square deviation from true geometry at first rank. The work will promote the use of new concepts where correlations among parameters provide more robust scoring models. It will facilitate studies involving molecular interactions, including modeling of large macromolecular assemblies and protein structure prediction. (C) 2010 Wiley Periodicals, Inc. J Comput Chem 32: 787-796, 2011.

Electron transfer reaction in the marcus inverted region: role of high frequency vibrational modes

A theoretical study of the dynamics of photo-electron transfer reactions in the Marcus inverted regime is presented. This study is motivated partly by the recent proposal of Barbara et al. (J. Phys. Chem. 96, 3728, 1991) that a minimal model of an electron transfer reaction should consist of a polar solvent mode (X), a low frequency vibrational mode (Q) and one high frequency mode (q). Interplay between these modes may be responsible for the crossover observed in the dynamics from a solvent controlled to a vibrational controlled electron transfer. The following results have been obtained. (i) In the case of slowly relaxing solvents, the proximity of the point of excitation to an effective sink on the excited surface is critical in determining the decay of the reactant population. This is because the Franck-Condon overlap between the reactant ground and the product excited states decreases rapidly with increase in the quantum number of the product vibrational state. (ii) Non-exponential solvation dynamics has an important effect in determining the rates of electron transfer. Especially, a biphasic solvation and a large coupling between the reactant and the product states both may be needed to explain the experimental results. ©1996 American Institute of Physics

Quantum and non-Markovian effects in the electron transfer reaction dynamics in the Marcus inverted region

Electron transfer reactions in large molecules may often be coupled to both the polar solvent modes and the intramolecular vibrational modes of the molecule. This can give rise to a complex dynamics which may in some systems, like betaine, be controlled more by vibrational rather than by solvent effects. Additionally, a significant contribution from an ultrafast relaxation component in the solvation dynamics may enhance the complexity. To explain the wide range of behavior that has been observed experimentally, Barbara et al. recently proposed that a model of an electron transfer reaction should minimally consist of a low-frequency classical solvent mode (X), a low-frequency vibrational mode (Q), and a high-frequency quantum mode (q) (J. Phys. Chem. 1991, 96, 3728). In the present work, a theoretical study of this model is described. This study generalizes earlier work by including the biphasic solvent response and the dynamics of the low-frequency vibrational mode in the presence of a delocalized, extended reaction zone. A novel Green's function technique has been developed which allowed us to study the non-Markovian dynamics on a multidimensional surface. The contributions from the high-frequency vibrational mode and the ultrafast component in the non-Markovian solvent dynamics are found to be primarily responsible for the dramatic increase in charge transfer rate over the prediction of the classical theories that neglect both these factors. These, along with a large coupling between the reactant and the product states, may combine to render the electron transfer rate both very large and constant over a wide range of solvent relaxation rates. A study on the free energy gap dependence of the electron transfer rate reveals that the rates are sensitive to changes in the quantum frequency particularly when the free energy gap is very large.

The intercalation reaction of pyridine with manganese thiophosphate, MnPS3

The intercalation of pyridine in the layered manganese thiophosphate, MnPS3, has been examined in detail by a variety of techniques. The reaction is interesting since none of the anticipated changes in optical and electrical properties associated with intercalation of electron donating molecules is observed. The only notable change in the properties of the host lattice is in the nature of the low-temperature magnetic ordering; while MnPS3 orders antiferromagnetically below 78 K, the intercalated compound shows weak ferromagnetism probably due to a canted spin structure. Vibrational spectra clearly show that the intercalated species are pyridinium ions solvated by neutral pyridine molecules. The corresponding reduced sites of the host lattice, however, were never observed. The molecules in the solvation shell are exchangeable. Although the reaction appears to be topotactic and reversible, from XRD, a more detailed analysis of the products of deintercalation reveal that it is not so. The intercalation proceeds by an ion exchange/intercalation mechanism wherein the intercalated species are pyridinium ions solvated by neutral molecules with charge neutrality being preserved not by electron transfer but by a loss of Mn2+ ions from the lattice. The experimental evidence leading to this conclusion is discussed and it is shown that this model can account satisfactorily for the observed changes (or lack of it) in the optical, electrical, vibrational, and magnetic properties.

Orientational relaxation in dipolar systems: How much do we understand the role of correlations

Dipolar systems, both liquids and solids, constitute a class of naturally abundant systems that are important in all branches of natural science. The study of orientational relaxation provides a powerful method to understand the microscopic properties of these systems and, fortunately, there are many experimental tools to study orientational relaxation in the condensed phases. However, even after many years of intense research, our understanding of orientational relaxation in dipolar systems has remained largely imperfect. A major hurdle towards achieving a comprehensive understanding is the long range and complex nature of dipolar interactions which also made reliable theoretical study extremely difficult. These difficulties have led to the development of continuum model based theories, which although they provide simple, elegant expressions for quantities of interest, are mostly unsatisfactory as they totally neglect the molecularity of inter-molecular interactions. The situation has improved in recent years because of renewed studies, led by computer simulations. In this review, we shall address some of the recent advances, with emphasis on the work done in our laboratory at Bangalore. The reasons for the failure of the continuum model, as revealed by the recent Brownian dynamics simulations of the dipolar lattice, are discussed. The main reason is that the continuum model predicts too fast a decay of the torque-torque correlation function. On the other hand, a perturbative calculation, based on Zwanzig's projection operator technique, provides a fairly satisfactory description of the single particle orientational dynamics for not too strongly polar dipolar systems. A recently developed molecular hydrodynamic theory that properly includes the effects of intermolecular orientational pair correlations provides an even better description of the single-particle orientational dynamics. We also discuss the rank dependence of the dielectric friction. The other topics reviewed here includes dielectric relaxation and solvation dynamics, as they are intimately connected with orientational relaxation. Recent molecular dynamics simulations of the dipolar lattice are also discussed. The main theme of the present review is to understand the effects of intermolecular interactions on orientational relaxation. The presence of strong orientational pair correlation leads to a strong coupling between the single particle and the collective dynamics. This coupling can lead to rich dynamical properties, some of which are detailed here, while a major part remains yet unexplored.

Study of the reactions of pyridinium poly(hydrogen fluoride) with some ternary oxides

The products of the reaction of pyridinium poly(hydrogen fluoride), PPHF, with KIO3, Na2SnO3, NaBiO3, K2CrO4, Na2MoO4 and Na2WO4 were KIO2F2; Na2SnF6; NaHF2, BiF3; K3CrF6, KHF2, (PyH)(3)CrF6; NaHF2, (PyH)(2)MoO2F4.2NaHF(2); and (PyH)(2)WO2F4.2NaHF(2), respectively, while KClO3, KBrO3 and KlO(4) react with complete decomposition to form KHF2 as the fluorinated product. This differential reactivity and mode of reaction has been discussed in terms of the oxidation state of the central atom, the nature and strength of the bonds and the complex behaviour of the formed intermediate or fluorinated products that undergo complexation or solvation with pyridine and/or hydrogen fluoride.

Biological water: Femtosecond dynamics of macromolecular hydration

The unique features of a macromolecule and water as a solvent make the issue of solvation unconventional, with questions about the static versus dynamic nature of hydration and the, physics of orientational and translational diffusion at the boundary. For proteins, the hydration shell that covers the surface is critical to the stability of its structure and function. Dynamically speaking, the residence time of water at the surface is a signature of its mobility and binding. With femtosecond time resolution it is possible to unravel the shortest residence times which are key for the description of the hydration layer, static or dynamic. In this article we review these issues guided by experimental studies, from this laboratory, of polar hydration dynamics at the surfaces of two proteins (Subtilisin Carlsberg (SC) and Monellin). The natural probe tryptophan amino acid was used for the interrogation of the dynamics, and for direct comparison we also studied the behavior in bulk water - a complete hydration in 1 ps. We develop a theoretical description of solvation and relate the theory to the experimental observations. In this - theoretical approach, we consider the dynamical equilibrium in the hydration shell, defining the rate processes for breaking and making the transient hydrogen bonds, and the effective friction in the layer which is defined by the translational and orientational motions of water molecules. The relationship between the residence time of water molecules and the observed slow component in solvation dynamics is a direct one. For the two proteins studied, we observed a "bimodal decay" for the hydration correlation function, with two primary relaxation times: ultrafast, typically 1 ps or less, and longer, typically 15-40 ps, and both are related to the residence time at the protein surface, depending on the binding energies. We end by making extensions to studies of the denatured state of the protein, random coils, and the biomimetic micelles, and conclude with our thoughts on the relevance of the dynamics of native structures to their functions.