Temperature and length scale dependence of hydrophobic effects and their possible implications for protein folding

Single polymer studies of hydrophobic hydration

Hydrophobic interactions guide protein folding, multidomain protein assembly, receptor-ligand binding, membrane formation, and cellular transportation. On the macroscale, hydrophobic interactions consist of the aggregation of "oil-like" objects in water by minimizing the interfacial energy. However, studies of the hydration behavior of small hydrophobic molecules have shown that the microscopic (~1 nm) hydration mechanism differs fundamentally from its macroscopic counterpart. Theoretical studies over the last two decades have pointed to an intricate dependence of molecular hydration mechanisms on the length scale. The microscopic-to-macroscopic crossover length scale is critically important to hydrophobic interactions in polymers, proteins, and other macromolecules. Accurate experimental determination of hydration mechanisms and interaction strengths directly influence our understanding of protein folding. In this Account, we discuss our recent measurements of the hydration energies of single hydrophobic homopolymers as they unfold. We describe in detail our single molecule force spectroscopy technique, the interpretation of the single polymer force curve, and how it relates to the hydration free energy of a hydrophobic polymer. Specifically, we show how temperature, side-chain sizes and solvent conditions, affect the driving force of hydrophobic collapse. The experiments reveal that the size of the nonpolar polymer side-chains changes the thermal signatures of hydration. The sizes of the polymer side-chains bridge the length scale where theories had predicted a transition between entropically driven microscopic hydration and enthalpically driven macroscopic hydrophobic hydration. Our experimental results revealed a crossover length scale of approximately 1 nm, similar to the results from recent theoretical studies. Experiments that probe solvent dependency show that the microscopic polymer hydration is correlated with macroscopic interfacial tension. Consistent with theoretical predictions, the solvent conditions affect the microscopic and macroscopic hydrophobic strengths in similar ways. Although the extended polymers and proteins span hundreds of nanometers, the experiments show that their hydration behavior is determined by the size of a single hydrophobic monomer. As the hydrophobic particle size decreases from the macroscopic to the microscopic regime, the scaling relationship changes from a dependence on interfacial area to a dependence on volume. Therefore, under these conditions, the driving force for the aggregation of hydrophobic molecules is reduced, which has significant implications for the strength of hydrophobic interactions in molecular systems, particularly in protein folding.

Physical origin underlying the entropy loss upon hydrophobic hydration

The hydrophobic effect (HE) is commonly associated with the demixing of oil and water at ambient conditions and plays the leading role in determining the structure and stability of biomolecular assembly in aqueous solutions. On the molecular scale HE has an entropic origin. It is believed that hydrophobic particles induce order in the surrounding water by reducing the volume of configuration space available for hydrogen bonding. Here we show with computer simulation results that this traditional picture, based on average structural features of hydration water, configurational properties of single water molecules, and up to pairwise correlations, is not correct. Analyzing collective fluctuations in water clusters we are able to provide a fundamentally new picture of HE based on pronounced many-body correlations affecting the switching of hydrogen bonds (HBs) between molecules. These correlations emerge as a nonlocal compensation of reduced fluctuations of local electrostatic fields in the presence of an apolar solute. We propose an alternative view which may also be formulated as a maximization principle: The electrostatic noise acting on water molecules is maximized under the constraint that each water molecule on average maintains as many HBs as possible. In the presence of the solute the maximized electrostatic noise is a result of nonlocal fluctuations in the labile HB network giving rise to strong correlations among at least up to four water molecules.

Evidence for water-tuned structural differences in proteins: an approach emphasizing variations in local hydrophilicity

We present experimental evidence for the significant effect that water can have on the functional structure of proteins in solution. Human (HSA) and Bovine Serum Albumin (BSA) have an amino acid sequence identity of 75.52% and are chosen as model proteins. We employ EPR-based nanoscale distance measurements using double electron-electron resonance (DEER) spectroscopy and both albumins loaded with long chain fatty acids (FAs) in solution to globally (yet indirectly) characterize the tertiary protein structures from the bound ligands’ points of view. The complete primary structures and crystal structures of HSA and as of recently also BSA are available. We complement the picture as we have recently determined the DEER-derived solution structure of HSA and here present the corresponding BSA solution structure. The characteristic asymmetric FA distribution in the crystal structure of HSA can surprisingly be observed by DEER in BSA in solution. This indicates that the BSA conformational ensemble in solution seems to be narrow and close to the crystal structure of HSA. In contrast, for HSA in solution a much more symmetric FA distribution was found. Thus, conformational adaptability and flexibility dominate in the HSA solution structure while BSA seems to lack these properties. We further show that differences in amino acid hydropathies of specific structural regions in both proteins can be used to correlate the observed difference in the global (tertiary) solution structures with the differences on the primary structure level.

Hydrophobicity of proteins and interfaces: insights from density fluctuations

Macroscopic characterizations of hydrophobicity (e.g., contact angle measurements) do not extend to the surfaces of proteins and nanoparticles. Molecular measures of hydrophobicity of such surfaces need to account for the behavior of hydration water. Theory and state-of-the-art simulations suggest that water density fluctuations provide such a measure; fluctuations are enhanced near hydrophobic surfaces and quenched with increasing surface hydrophilicity. Fluctuations affect conformational equilibria and dynamics of molecules at interfaces. Enhanced fluctuations are reflected in enhanced cavity formation, more favorable binding of hydrophobic solutes, increased compressibility of hydration water, and enhanced water-water correlations at hydrophobic surfaces. These density fluctuation-based measures can be used to develop practical methods to map the hydrophobicity/philicity of heterogeneous surfaces including those of proteins. They highlight that the hydrophobicity of a group is context dependent and is significantly affected by its environment (e.g., chemistry and topography) and especially by confinement. The ability to include information about hydration water in mapping hydrophobicity is expected to significantly impact our understanding of protein-protein interactions as well as improve drug design and discovery methods and bioseparation processes.

Impact of chemical heterogeneity on protein self-assembly in water

Hydrophobicity is thought to underlie self-assembly in biological systems. However, the protein surface comprises hydrophobic and hydrophilic patches, and understanding the impact of such a chemical heterogeneity on protein self-assembly in water is of fundamental interest. Here, we report structural and thermodynamic investigations on the dimer formation of full-length amyloid-β proteins in water associated with Alzheimer's disease. Spontaneous dimerization process--from the individual diffusive regime at large separations, through the approach stage in which two proteins come close to each other, to the structural adjustment stage toward compact dimer formation--was captured in full atomic detail via unguided, explicit-water molecular dynamics simulations. The integral-equation theory of liquids was then applied to simulated protein structures to analyze hydration thermodynamic properties and the water-mediated interaction between proteins. We demonstrate that hydrophilic residues play a key role in initiating the dimerization process. A long-range hydration force of enthalpic origin acting on the hydrophilic residues provides the major thermodynamic force that drives two proteins to approach from a large separation to a contact distance. After two proteins make atomic contacts, the nature of the water-mediated interaction switches from a long-range enthalpic attraction to a short-range entropic one. The latter acts both on the hydrophobic and hydrophilic residues. Along with the direct protein-protein interactions that lead to the formation of intermonomer hydrogen bonds and van der Waals contacts, the water-mediated attraction of entropic origin brings about structural adjustment of constituent monomer proteins toward the formation of a compact dimer structure.

Temperature dependence of the structure of protein hydration water and the liquid-liquid transition

We study the temperature dependence of the structure and orientation of the first hydration layers of the protein lysozyme and compare it with the situation for a model homogeneous hydrophobic surface, a graphene sheet. We show that in both cases these layers are significantly better structured than bulk water. The geometrical constraint of the interface makes the water molecules adjacent to the surface lose one water-water hydrogen bond and expel the fourth neighbors away from the surface, lowering local density. We show that a decrease in temperature improves the ordering of the hydration water molecules, preserving such a geometrical effect. For the case of graphene, this favors an ice Ih-like local structuring, similar to the water-air interface but in the opposite way along the c axis of the basal plane (while the vicinal water molecules of the air interface orient a hydrogen atom toward the surface, the oxygens of the water molecules close to the graphene plane orient a lone pair in such a direction). In turn, the case of the first hydration layers of the lysozyme molecule is shown to be more complicated, but still displaying signs of both kinds of behavior, together with a tendency of the proximal water molecules to hydrogen bond to the protein both as donors and as acceptors. Additionally, we make evident the existence of signatures of a liquid-liquid transition (Widom line crossing) in different structural parameters at the temperature corresponding to the dynamic transition incorrectly referred to as "the protein glass transition."

A smoothly decoupled particle interface: new methods for coupling explicit and implicit solvent

A common theme of studies using molecular simulation is a necessary compromise between computational efficiency and resolution of the forcefield that is used. Significant efforts have been directed at combining multiple levels of granularity within a single simulation in order to maintain the efficiency of coarse-grained models, while using finer resolution in regions where such details are expected to play an important role. A specific example of this paradigm is the development of hybrid solvent models, which explicitly sample the solvent degrees of freedom within a specified domain while utilizing a continuum description elsewhere. Unfortunately, these models are complicated by the presence of structural artifacts at or near the explicit/implicit boundary. The presence of these artifacts significantly complicates the use of such models, both undermining the accuracy obtained and necessitating the parameterization of effective potentials to counteract the artificial interactions. In this work, we introduce a novel hybrid solvent model that employs a smoothly decoupled particle interface (SDPI), a switching region that gradually transitions from fully interacting particles to a continuum solvent. The resulting SDPI model allows for the use of an implicit solvent model based on a simple theory that needs to only reproduce the behavior of bulk solvent rather than the more complex features of local interactions. In this study, the SDPI model is tested on spherical hybrid domains using a coarse-grained representation of water that includes only Lennard-Jones interactions. The results demonstrate that this model is capable of reproducing solvent configurations absent of boundary artifacts, as if they were taken from full explicit simulations.

Stability of Positively Charged Solutes in Water: A Transition from Hydrophobic to Hydrophilic

To improve the description of solvation thermodynamics of biomolecules, we report here the dependence of solvation on the curvature and surface charge of positively charged solutes in water based on extensive molecular dynamics simulations analyzed using the two-phase thermodynamic method. At a surface charge of +0.4e, the compensating forces of favorable electrostatic stabilization and entropic destabilization cancel almost exactly, representing a molecular crossover point from hydrophobic to hydrophilic behavior, independent of curvature. These results suggest that one should include charge-dependent entropic corrections to continuum models aimed at predicting the solvation free energies of large biomolecules.

The effect of temperature on single-polypeptide adsorption

The hydrophobic attraction (HA) is believed to be one of the main driving forces for protein folding. Understanding its temperature dependence promises a deeper understanding of protein folding. Herein, we present an approach to investigate the HA with a combined experimental and simulation approach, which is complementary to previous studies on the temperature dependence of the solvation of small hydrophobic spherical particles. We determine the temperature dependence of the free-energy change and detachment length upon desorption of single polypeptides from hydrophobic substrates in aqueous environment. Both the atomic force microscopy (AFM) based experiments and the molecular dynamics (MD) simulations show only a weak dependence of the free energy change on temperature. In fact, depending on the substrate, we find a maximum or a minimum in the temperature-dependent free energy change, meaning that the entropy increases or decreases with temperature for different substrates. These observations are in contrast to the solvation of small hydrophobic particles and can be rationalized by a compensation mechanism between the various contributions to the desorption force. On the one hand this is reminiscent of the protein folding process, where large entropic and enthalpic contributions compensate each other to result in a small free energy difference between the folded and unfolded state. On the other hand, the protein folding process shows much stronger temperature dependence, pointing to a fundamental difference between protein folding and adsorption. Nevertheless such temperature dependent single molecule desorption studies open large possibilities to study equilibrium and non-equilibrium processes dominated by the hydrophobic attraction.

Diversified nanoparticle assembly pathways: materials architecture control beyond the amphiphilicity paradigm

The functional versatility of a chemical system is ultimately dictated by the availability of distinctly accessible architectures. The generation of a diverse array of assembled constructs from a single type of nanoscale building block is a promising yet largely elusive goal. We report herein the utility of a monolayer-modified nanoparticle for the creation of a broad range of architectures. The versatile modes of assembly complement the conventionally used, amphiphilicity-driven strategy. We demonstrate that one can vary the nanoparticle assembly pathways within the confines of solvent media through the modulation of interactions and partitioning of nanoparticles. Merging of the molecular-scale design and higher-ordered arrangement enables diversified assembly through the manipulation of experimental parameters such as solvent, pH, affinity molecule, and temperature. Microfluidics provides an effective channel to control the monodispersity and size on all the architectures attainable in the bulk solution phase. These observations could be further explored for an understanding of diversified matter organization and order generation beyond the amphiphilicity paradigm.

Extended surfaces modulate hydrophobic interactions of neighboring solutes

Interfaces are a most common motif in complex systems. To understand how the presence of interfaces affects hydrophobic phenomena, we use molecular simulations and theory to study hydration of solutes at interfaces. The solutes range in size from subnanometer to a few nanometers. The interfaces are self-assembled monolayers with a range of chemistries, from hydrophilic to hydrophobic. We show that the driving force for assembly in the vicinity of a hydrophobic surface is weaker than that in bulk water and decreases with increasing temperature, in contrast to that in the bulk. We explain these distinct features in terms of an interplay between interfacial fluctuations and excluded volume effects--the physics encoded in Lum-Chandler-Weeks theory [Lum K, Chandler D, Weeks JD (1999) J Phys Chem B 103:4570-4577]. Our results suggest a catalytic role for hydrophobic interfaces in the unfolding of proteins, for example, in the interior of chaperonins and in amyloid formation.

Behavior of water in contact with model hydrophobic cavities and tunnels and carbon nanotubes

By means of molecular dynamics simulations we analyze the behavior of water in contact with model hydrophobic cavities and tunnels. We study the hydration and filling propensity of cavities and tunnels carved in alkane monolayers and, for comparison, we also study single-walled carbon nanotubes of similar size. Our results will determine the dependence of the filling propensity as a function of cavity size while revealing the dynamical nature of the process with alternation of filled and dry states. Concerning the tunnels built across the monolayer, we shall show that the minimum diameter in order to get filled is about twice as large as that for the carbon nanotubes, thus evidencing a more hydrophobic behavior. The existence of water-water hydrogen bonds, a necessary condition for penetration, will also be made evident.

Signature of hydrophobic hydration in a single polymer

Hydrophobicity underpins self-assembly in many natural and synthetic molecular and nanoscale systems. A signature of hydrophobicity is its temperature dependence. The first experimental evaluation of the temperature and size dependence of hydration free energy in a single hydrophobic polymer is reported, which tests key assumptions in models of hydrophobic interactions in protein folding. Herein, the hydration free energy required to extend three hydrophobic polymers with differently sized aromatic side chains was directly measured by single molecule force spectroscopy. The results are threefold. First, the hydration free energy per monomer is found to be strongly dependent on temperature and does not follow interfacial thermodynamics. Second, the temperature dependence profiles are distinct among the three hydrophobic polymers as a result of a hydrophobic size effect at the subnanometer scale. Third, the hydration free energy of a monomer on a macromolecule is different from a free monomer; corrections for the reduced hydration free energy due to hydrophobic interaction from neighboring units are required.

Possible mechanism underlying high-pressure unfolding of proteins: formation of a short-period high-density hydration shell

Hydration effects on high-pressure unfolding of a hydrophobic polymer chain are investigated through a multiscale simulation based on density-functional theory. The results strongly suggest the following: a thermodynamic origin for high-pressure denaturation, i.e., the decrease in volume due to the unfolding can be explained by the formation of a short-period high-density hydration shell.

Quantitative investigation of the two-state picture for water in the normal liquid and the supercooled regime

Several evidences have helped to establish the two-state nature of liquid water. Thus, within the normal liquid and supercooled regimes water has been shown to consist of a mixture of well-structured, low-density molecules and unstructured, high-density ones. However, quantitative analyses have faced the burden of unambiguously determining both the presence and the fraction of each kind of water "species". A recent approach by combining a local structure index with potential-energy minimisations allows us to overcome this difficulty. Thus, in this work we extend such study and employ it to quantitatively determine the fraction of structured molecules as a function of temperature for different densities. This enables us to validate predictions of two-state models.

Temperature dependence of lysozyme hydration and the role of elastic energy

Water plays a critical role in protein dynamics and functions. However, the most basic property of hydration--the water sorption isotherm--remains inadequately understood. Surface adsorption is the commonly adopted picture of hydration. Since it does not account for changes in the conformational entropy of proteins, it is difficult to explain why protein dynamics and activity change upon hydration. The solution picture of hydration provides an alternative approach to describe the thermodynamics of hydration. Here, the flexibility of proteins could influence the hydration level through the change of elastic energy upon hydration. Using nuclear magnetic resonance to measure the isotherms of lysozyme in situ between 18 and 2 °C, the present work provides evidence that the part of water uptake associated with the onset of protein function is significantly reduced below 8 °C. Quantitative analysis shows that such reduction is directly related to the reduction of protein flexibility and enhanced cost in elastic energy upon hydration at lower temperature. The elastic property derived from the water isotherm agrees with direct mechanical measurements, providing independent support for the solution model. This result also implies that water adsorption at charged and polar groups occurring at low vapor pressure, which is known for softening the protein, is crucial for the later stage of water uptake, leading to the activation of protein dynamics. The present work sheds light on the mutual influence of protein flexibility and hydration, providing the basis for understanding the role of hydration on protein dynamics.

An improved coarse-grained model of solvation and the hydrophobic effect

We present a coarse-grained lattice model of solvation thermodynamics and the hydrophobic effect that implements the ideas of Lum-Chandler-Weeks theory [J. Phys. Chem. B 134, 4570 (1999)] and improves upon previous lattice models based on it. Through comparison with molecular simulation, we show that our model captures the length-scale and curvature dependence of solvation free energies with near-quantitative accuracy and 2-3 orders of magnitude less computational effort, and further, correctly describes the large but rare solvent fluctuations that are involved in dewetting, vapor tube formation, and hydrophobic assembly. Our model is intermediate in detail and complexity between implicit-solvent models and explicit-water simulations.

Interplay between hydration water and headgroup dynamics in lipid bilayers

In this study, the interplay between water and lipid dynamics has been investigated by broadband dielectric spectroscopy and modulated differential scanning calorimetry (MDSC). The multilamellar lipid bilayer system 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) has been studied over a broad temperature range at three different water contents: about 3, 6, and 9 water molecules per lipid molecule. The results from the dielectric relaxation measurements show that at temperatures <250 K the lipid headgroup rotation is described by a super-Arrhenius temperature dependence at the lowest hydration level and by the Arrhenius law at the highest hydration level. This difference in the temperature dependence of the lipid headgroup rotation can be explained by the increasing interaction between the headgroups with decreasing water content, which causes their rotational motion to be more cooperative in character. The main water relaxation shows an anomalous dependence on the water content in the supercooled and glassy regime. In contrast to the general behavior of interfacial water, the water dynamics is fastest in the driest sample and its temperature dependence is best described by a super-Arrhenius temperature dependence. The best explanation for this anomalous behavior is that the water relaxation becomes more determined by fast local lipid motions than by the intrinsic water dynamics at low water contents. In support for this interpretation is the finding that the relaxation time of the main water process is faster than that in most other host systems at temperatures below 180 K. Thus, the dielectric relaxation data show clearly the strong interplay between water and lipid dynamics; the water influences the lipid dynamics and vice versa. In the MDSC data, we observe a weak enthalpy relaxation at 203 K for the driest sample and at 179 K for the most hydrated sample, attributed to the freezing-in of the lipid headgroup rotation observed in the dielectric data, since this motion reaches a time scale of about 100 s at about the same temperatures.

Role of electrostatics in modulating hydrophobic interactions and barriers to hydrophobic assembly

Hydrophobic effects continue to be an active area of research due to implications for a wide range of physicochemical phenomena. Molecular dynamics simulations have been used extensively in the study of such effects using various water potential models, with few studies addressing the differences between models. In particular, studies considering the explicit treatment of water polarizability are underrepresented in the literature. We present results from molecular dynamics simulations that systematically compare the dependence of large-scale hydrophobic effects on the water model. We consider three common nonpolarizable models (SPC/E, TIP3P, and TIP4P) and two common polarizable models (TIP4P-FQ and SWM4-NDP). Results highlight the similarities and differences of the different water models in the vicinity of two large hydrophobic plates. In particular, profiles of average density, density fluctuations, orientation, and hydrogen bonding show only minor differences among the water models studied. However, the potential of mean force for the hydrophobe dimerization is significantly reduced in the polarizable water systems. TIP4P-FQ shows the deepest minimum of approximately -54(+/-3) kcal/mol compared to -40(+/-3), -40(+/-2), -42(+/-3), and -45(+/-5) kcal/mol for TIP4P, TIP3P, SPC/E, and SWM4-NDP (all relative to the dissociated state). We discuss the relationship between hydrophobic association and the strength of water-water interactions in the liquid phase. Results suggest that models treating polarizability (both implicitly and explicitly) influence a stronger driving force toward hydrophobic assembly. Implications of these results, as well as prospectives on future work, are discussed.

Structural and dynamical aspects of water in contact with a hydrophobic surface

By means of molecular dynamics simulations we study the structure and dynamics of water molecules in contact with a model hydrophobic surface: a planar graphene-like layer. The analysis of the distributions of a local structural index indicates that the water molecules proximal to the graphene layer are considerably more structured than the rest and, thus, than the bulk. This structuring effect is lost in a few angstroms and is basically independent of temperature for a range studied comprising parts of both the normal liquid and supercooled states (240K to 320K). In turn, such structured water molecules present a dynamics that is slower than the bulk, as a consequence of their improved interactions with their first neighbors.

A probabilistic approach to the effect of water hydrogen bonds on the kinetics of protein folding and protein denaturation

Previously, we presented a review of our kinetic models for the nucleation mechanism of protein folding and for the protein thermal denaturation in a barrierless way. A protein was treated as a random heteropolymer consisting of hydrophobic, hydrophilic, and neutral beads. As a crucial idea of the model, an overall potential around the cluster of native residues wherein a residue performs a chaotic motion was considered as the combination of the average dihedral, effective pairwise, and confining potentials. The overall potential as a function of the distance from the cluster has a double well shape. This allowed one to develop kinetic models for the nucleation mechanism of protein folding (NMPF) and barrierless protein denaturation (BPD) by using the mean first passage time analysis. In the original models, however, hydrogen bonding effects were taken into account only indirectly which affected the accuracy of the models because hydrogen bonding does play a crucial role in the folding, stability, and denaturation of proteins. To improve the NMPF and BPD models and explicitly take into account the hydrogen bonding "water-water" and "water-protein residue", we have developed a probabilistic hydrogen bond (PHB) model for the effect of hydrogen bond networks of water molecules around two solute particles (immersed in water) on their interaction, and have then combined the PHB model with the NMPF and BPD models. In this paper, that can be regarded as sequel of our previous review, we analyze the modified NMPF and BPD models that explicitly take into account the effect of water-water hydrogen bonding on these processes. As expected, the application of the modified models to the folding/unfolding of two model proteins (one short, consisting of 124 residues and the other large, consisting of 2500 residues) demonstrate that the hydrogen bond networks play a very important role in the protein folding/unfolding phenomena.

Study of entropy-driven self-assembly of rigid macromolecules

A simple model composed of two rigid macromolecules (adsorbents) immersed in a large number of small molecules (adsorbates) is used to study entropy-driven association processes. The surfaces of the adsorbents are capable of adsorbing the smaller adsorbates. The partition function of the model is obtained analytically. The probability of dimerization and the number of adsorbed molecules are shown to depend on the enthalpy and the entropy differences between the assembled and the disassembled states. Under certain conditions, dimerization of the macromolecules occurs with increasing temperature. This entropy-driven self-assembly is originated from an overall entropy gain due to the release of the adsorbed small molecules, leading to a large peak in the heat capacity due to the dimer formation. The desorption of the adsorbates induces a sharp change in the first-order derivative of the free energy, resulting in another large peak in the heat capacity. A temperature-induced re-entrance into the dimer state is also contained in the model.

Force field bias in protein folding simulations

Long timescale (>1 micros) molecular dynamics simulations of protein folding offer a powerful tool for understanding the atomic-scale interactions that determine a protein's folding pathway and stabilize its native state. Unfortunately, when the simulated protein fails to fold, it is often unclear whether the failure is due to a deficiency in the underlying force fields or simply a lack of sufficient simulation time. We examine one such case, the human Pin1 WW domain, using the recently developed deactivated morphing method to calculate free energy differences between misfolded and folded states. We find that the force field we used favors the misfolded states, explaining the failure of the folding simulations. Possible further applications of deactivated morphing and implications for force field development are discussed.

Computational fragment-based binding site identification by ligand competitive saturation

Fragment-based drug discovery using NMR and x-ray crystallographic methods has proven utility but also non-trivial time, materials, and labor costs. Current computational fragment-based approaches circumvent these issues but suffer from limited representations of protein flexibility and solvation effects, leading to difficulties with rigorous ranking of fragment affinities. To overcome these limitations we describe an explicit solvent all-atom molecular dynamics methodology (SILCS: Site Identification by Ligand Competitive Saturation) that uses small aliphatic and aromatic molecules plus water molecules to map the affinity pattern of a protein for hydrophobic groups, aromatic groups, hydrogen bond donors, and hydrogen bond acceptors. By simultaneously incorporating ligands representative of all these functionalities, the method is an in silico free energy-based competition assay that generates three-dimensional probability maps of fragment binding (FragMaps) indicating favorable fragment∶protein interactions. Applied to the two-fold symmetric oncoprotein BCL-6, the SILCS method yields two-fold symmetric FragMaps that recapitulate the crystallographic binding modes of the SMRT and BCOR peptides. These FragMaps account both for important sequence and structure differences in the C-terminal halves of the two peptides and also the high mobility of the BCL-6 His116 sidechain in the peptide-binding groove. Such SILCS FragMaps can be used to qualitatively inform the design of small-molecule inhibitors or as scoring grids for high-throughput in silico docking that incorporate both an atomic-level description of solvation and protein flexibility.

Fragment-based drug discovery is based on a simple yet powerful principle: instead of trying to screen through the vast number of possible drug-like compounds during the drug discovery process, screen representative drug-like fragments, which are far fewer in number. Once a suitable fragment is discovered, it can then be built up or linked with other fragments to give a drug-like molecule. Because such fragments are small, even “good” fragments bind weakly to their targets, therefore requiring significant time, labor, and materials costs for experimental detection and characterization of binding. In the present work, we describe a computational approach to the problem of detecting and characterizing fragment binding. Importantly, the method provides atomic-resolution results and also explicitly takes into account the effect that molecular water has on binding and the inherent flexibility of protein targets. The methodology is demonstrated by application to the BCL-6 protein, which is implicated in a variety of cancers, is conceptually easy to understand, and can yield results in a matter of days using present-day commodity computers.

First passage time analysis of protein folding via nucleation and of barrierless protein denaturation

A review of the kinetic models, recently developed by the authors for the nucleation mechanism of protein folding and for the barrierless thermal denaturation, is presented. Both models are based on the mean first passage time analysis. A protein is treated as a random heteropolymer consisting of hydrophobic, hydrophilic, or neutral beads. As a crucial idea of the model, an overall potential around the cluster of native residues wherein a residue performs a chaotic motion is considered as the combination of the average dihedral, effective pairwise, and confining potentials. The overall potential as a function of the distance from the cluster center has a double well shape which allows one to determine its emission and absorption rates by the first passage time analysis. One can thus develop a theory for the nucleation mechanism of protein folding and calculate the temperature dependence of the folding time. A kinetic model for protein denaturation occurring in a barrierless way has been also developed by using the same approach. The numerical calculations for two model proteins (one consisting of 124 amino acids and the other of 2500 amino acids) demonstrate that the models can predict folding and unfolding times consistent with experimental data.

Novel behavior of heat of micellization of pluronics F68 and F88 in aqueous solutions

It is well understood that the heat of micellization for surfactants is monotonically decreased along with an increase in temperature. However, this behavior for polymeric surfactants has never been carefully examined. In this study, the heat of micellization of poly(ethylene oxide)-poly(propylene oxide)-poly(ethylene oxide) (PEO-PPO-PEO) triblock copolymers (Pluronics F68 and F88) in water as a function of temperature is carefully examined by using a high-sensitivity differential scanning calorimeter (HSDSC). The critical micelle temperature (CMT) decreases along with an increase in the concentration of Pluronic F68 (or F88). The heat of micellization decreases along with an increase in the temperature, as expected, when the CMT is higher than 55 and 42 degrees C for Pluronics F68 and F88, respectively. It is interesting to observe that the heat of micellization increases along with the temperature while the temperature is below 55 and 42 degrees C for Pluronics F68 and F88, respectively. The enthalpy-entropy compensation phenomenon for the micellization of Pluronics F68 and F88 in connection with the hydrophobicity is discussed.