Hydration structure of the α-chymotrypsin substrate binding pocket: the impact of constrained geometry

On the Nature of Guest Complexation in Water: Triggered Wetting-Water-Mediated Binding

The complexity of macromolecular surfaces means that there are still many open questions regarding how specific areas are solvated and how this might affect the complexation of guests. Contributing to the identification and classification of the different possible mechanisms of complexation events in aqueous solution, and as part of the recent SAMPL8 exercise, we report here on the synthesis and conformational properties of TEEtOA 2, a cavitand with conformationally flexible ethyl groups at its portal. Using a combination of ITC and NMR spectroscopy, we report the binding affinities of a series of carboxylates to 2 and compare it to a related cavitand TEMOA 1. Additionally, we report MD simulations revealing how the wetting of the pocket of 2 is controlled by the conformation of its rim ethyl groups and, correspondingly, a novel triggered wetting, guest complexation mechanism, whereby the approaching guest opens up the pocket of the host, inducing its wetting and ultimately allows the formation of a hydrated host-guest complex (H·G·H2O). A general classification of complexation mechanisms is also suggested.

Evidence for Entropically Controlled Interfacial Hydration in Mesoporous Organosilicas

At aqueous interfaces, the distribution and dynamics of adsorbates are modulated by the behavior of interfacial water. Hydration of a hydrophobic surface can store entropy via the ordering of interfacial water, which contributes to the Gibbs energy of solute binding. However, there is little experimental evidence for the existence of such entropic reservoirs, and virtually no precedent for their rational design in systems involving extended interfaces. In this study, two series of mesoporous silicas were modified in distinct ways: (1) progressively deeper thermal dehydroxylation, via condensation of surface silanols, and (2) increasing incorporation of nonpolar organic linkers into the silica framework. Both approaches result in decreasing average surface polarity, manifested in a blue-shift in the fluorescence of an adsorbed dye. For the inorganic silicas, hydrogen-bonding of water becomes less extensive as the number of surface silanols decreases. Overhauser dynamic nuclear polarization (ODNP) relaxometry indicates enhanced surface water diffusivity, reflecting a loss of enthalpic hydration. In contrast, organosilicas show a monotonic decrease in surface water diffusivity with decreasing polarity, reflecting enhanced hydrophobic hydration. Molecular dynamics simulations predict increased tetrahedrality of interfacial water for the organosilicas, implying increased ordering near the nm-size organic domains (relative to inorganic silicas, which necessarily lack such domains). These findings validate the prediction that hydrophobic hydration at interfaces is controlled by the microscopic length scale of the hydrophobic regions. They further suggest that the hydration thermodynamics of structurally heterogeneous silica surfaces can be tuned to promote adsorption, which in turn tunes the selectivity in catalytic reactions.

From Specific γ-CD/[Nb6 Cl12 (H2 O)6 ]2+ Recognition to Biological Activity Tuning

Specific molecular recognition of γ-cyclodextrin (γ-CD) by the cationic hexanuclear niobium [Nb₆ Cl₁₂ (H₂ O)₆ ]²⁺ cluster complex in aqueous solutions results in a 1:1 supramolecular assembly {[Nb₆ Cl₁₂ (H₂ O)₆ ]@γ-CD}²⁺ . NMR spectroscopy, isothermal titration calorimetry (ITC), and ESI-MS were used to study the interaction between the inorganic cluster and the organic macrocycle. Such molecular association affects the biological activity of [Nb₆ Cl₁₂ (H₂ O)₆ ]²⁺ , decreasing its cytotoxicity despite enhanced cellular uptake. The 1:1 stoichiometry is maintained in solution over a large window of the reagents' ratio, but crystallization by slow evaporation produces a 1:2 host-guest complex [Nb₆ Cl₁₂ (H₂ O)₆ @(γ-CD)₂ ]Cl₂ ⋅20 H₂ O featuring the cluster encapsulated between two molecules of γ-CD. The 1:2 complex was characterized by XRD, elemental analysis, IR spectroscopy, and thermogravimetric analysis (TGA). Quantum chemical calculations were performed to model host-guest interaction.

How proteins modify water dynamics

Much of biology happens at the protein-water interface, so all dynamical processes in this region are of fundamental importance. Local structural fluctuations in the hydration layer can be probed by ¹⁷O magnetic relaxation dispersion (MRD), which, at high frequencies, measures the integral of a biaxial rotational time correlation function (TCF)-the integral rotational correlation time. Numerous ¹⁷O MRD studies have demonstrated that this correlation time, when averaged over the first hydration shell, is longer than in bulk water by a factor 3-5. This rotational perturbation factor (RPF) has been corroborated by molecular dynamics simulations, which can also reveal the underlying molecular mechanisms. Here, we address several outstanding problems in this area by analyzing an extensive set of molecular dynamics data, including four globular proteins and three water models. The vexed issue of polarity versus topography as the primary determinant of hydration water dynamics is resolved by establishing a protein-invariant exponential dependence of the RPF on a simple confinement index. We conclude that the previously observed correlation of the RPF with surface polarity is a secondary effect of the correlation between polarity and confinement. Water rotation interpolates between a perturbed but bulk-like collective mechanism at low confinement and an exchange-mediated orientational randomization (EMOR) mechanism at high confinement. The EMOR process, which accounts for about half of the RPF, was not recognized in previous simulation studies, where only the early part of the TCF was examined. Based on the analysis of the experimentally relevant TCF over its full time course, we compare simulated and measured RPFs, finding a 30% discrepancy attributable to force field imperfections. We also compute the full ¹⁷O MRD profile, including the low-frequency dispersion produced by buried water molecules. Computing a local RPF for each hydration shell, we find that the perturbation decays exponentially with a decay "length" of 0.3 shells and that the second and higher shells account for a mere 3% of the total perturbation measured by ¹⁷O MRD. The only long-range effect is a weak water alignment in the electric field produced by an electroneutral protein (not screened by counterions), but this effect is negligibly small for ¹⁷O MRD. By contrast, we find that the ¹⁷O TCF is significantly more sensitive to the important short-range perturbations than the other two TCFs examined here.

Natural Products and Their Mimics as Targets of Opportunity for Discovery

Diverse structural types of natural products and their mimics have served as targets of opportunity in our laboratory to inspire the discovery and development of new methods and strategies to assemble polyfunctional and polycyclic molecular architectures. Furthermore, our efforts toward identifying novel compounds having useful biological properties led to the creation of new targets, many of which posed synthetic challenges that required the invention of new methodology. In this Perspective, selected examples of how we have exploited a diverse range of natural products and their mimics to create, explore, and solve a variety of problems in chemistry and biology will be discussed. The journey was not without its twists and turns, but the unexpected often led to new revelations and insights. Indeed, in our recent excursion into applications of synthetic organic chemistry to neuroscience, avoiding the more-traveled paths was richly rewarding.

Recent developments in the theoretical, simulational, and experimental studies of the role of water hydrogen bonding in hydrophobic phenomena

Hydrophobic effects (hydrophobic hydration and hydrophobic interaction) constitute an important element of a wide variety of phenomena relevant to biological, physical, chemical, environmental, engineering, and pharmaceutical sciences, such as the immiscibility of oil and water, self-assembly of amphiphiles leading to micelle and membrane formation, folding and stability and unfolding of the native structure of a biologically active protein, gating of ion channels, wetting, froth floatation, and adhesion. On the other hand, the hydrogen bonding ability of water plays a major (if not crucial) role in hydrophobic phenomena. We present a review of most important and relatively recent experimental, simulational, and theoretical research on hydrophobic phenomena in various systems. With a particular interest we survey investigations clarifying the role of water hydrogen bonding therein, because it has been the main object of our own recent research. We have developed a probabilistic hydrogen bond (PHB) model that allows one to obtain an analytic expression for the number of bonds per water molecule as a function of its distance to a hydrophobe, hydrophobe radius, and temperature. Knowing that function, one can explicitly identify a water hydrogen bond contribution to the external potential whereto a water molecule is subjected near a hydrophobe. Combining the PHB model with the classical density functional theory (DFT), one can examine the contribution of water hydrogen bonding to the temperature and lengthscale effects on the hydration of particles and on their solvent-mediated interactions over the entire low-to-high temperature and small-to-large lengthscale ranges. We applied the combined DFT/PHB model to study a variety of hydrophobic phenomena such as (liquid) water in contact with a hydrophobic plate, solvation of spherical solutes of various radii in associated and non-associated liquids at various temperatures, the solvent-mediated interaction of spherical solutes and its temperature dependence, interaction of C60 fullerenes in water, temperature effect on the evaporation lengthscale of water confined between two hydrophobes, temperature dependence of the effective width of the solute-solvent transition layer and average density therein. These applications demonstrated that the DFT/PHB model can serve as a valuable tool in studying hydrophobic phenomena because it constitutes a balanced combination of simplicity, accuracy, and detail. The predictions of the combined DFT/PHB approach for the solvent density profiles and thermodynamic aspects of hydrophobic phenomena are generally in good agreement with experiments and simulations. For example, it predicts the small-to-large crossover lengthscale of its mechanism to be approximately in the range from 1nm to 4nm, and decreasing with increasing temperature. It also suggests that, in terms of the average fluid density in the solute-solvent transition layer, the transition layer for small hydrophobes (of radii ≲2 nm) becomes enriched with rather than depleted of fluid when both the solvent-solute affinity and hb-energy alteration ratio become large enough. The boundary values of these parameters, needed for the depletion-to-enrichment crossover, are predicted to decrease with increasing temperature.

Design of Catalytic Peptides and Proteins Through Rational and Combinatorial Approaches

This review focuses on recent progress in noncomputational methods to introduce catalytic function into proteins, peptides, and peptide assemblies. We discuss various approaches to creating catalytic activity and classification of noncomputational methods into rational and combinatorial classes. The section on rational design covers recent progress in the development of short peptides and oligomeric peptide assemblies for various natural and unnatural reactions. The section on combinatorial design describes recent advances in the discovery of catalytic peptides. We present the future prospects of these and other new approaches in a broader context, including implications for functional material design.

Molecular Self-Assembly Strategy for Generating Catalytic Hybrid Polypeptides

Recently, catalytic peptides were introduced that mimicked protease activities and showed promising selectivity of products even in organic solvents where protease cannot perform well. However, their catalytic efficiency was extremely low compared to natural enzyme counterparts presumably due to the lack of stable tertiary fold. We hypothesized that assembling these peptides along with simple hydrophobic pockets, mimicking enzyme active sites, could enhance the catalytic activity. Here we fused the sequence of catalytic peptide CP4, capable of protease and esterase-like activities, into a short amyloidogenic peptide fragment of Aβ. When the fused CP4-Aβ construct assembled into antiparallel β-sheets and amyloid fibrils, a 4.0-fold increase in the hydrolysis rate of p-nitrophenyl acetate (p-NPA) compared to neat CP4 peptide was observed. The enhanced catalytic activity of CP4-Aβ assembly could be explained both by pre-organization of a catalytically competent Ser-His-acid triad and hydrophobic stabilization of a bound substrate between the triad and p-NPA, indicating that a design strategy for self-assembled peptides is important to accomplish the desired functionality.

Molecular recognition in chemical and biological systems

Structure-based ligand design in medicinal chemistry and crop protection relies on the identification and quantification of weak noncovalent interactions and understanding the role of water. Small-molecule and protein structural database searches are important tools to retrieve existing knowledge. Thermodynamic profiling, combined with X-ray structural and computational studies, is the key to elucidate the energetics of the replacement of water by ligands. Biological receptor sites vary greatly in shape, conformational dynamics, and polarity, and require different ligand-design strategies, as shown for various case studies. Interactions between dipoles have become a central theme of molecular recognition. Orthogonal interactions, halogen bonding, and amide⋅⋅⋅π stacking provide new tools for innovative lead optimization. The combination of synthetic models and biological complexation studies is required to gather reliable information on weak noncovalent interactions and the role of water.

Protein-ligand interactions: probing the energetics of a putative cation-π interaction

In order to probe the energetics associated with a putative cation-π interaction, thermodynamic parameters are determined for complex formation between the Grb2 SH2 domain and tripeptide derivatives of RCO-pTyr-Ac6c-Asn wherein the R group is varied to include different alkyl, cycloalkyl, and aryl groups. Although an indole ring is reputed to have the strongest interaction with a guanidinium ion, binding free energies, ΔG°, for derivatives of RCO-pTyr-Ac6c-Asn bearing cyclohexyl and phenyl groups were slightly more favorable than their indolyl analog. Crystallographic analysis of two complexes reveals that test ligands bind in similar poses with the notable exception of the relative orientation and proximity of the phenyl and indolyl rings relative to an arginine residue of the domain. These spatial orientations are consistent with those observed in other cation-π interactions, but there is no net energetic benefit to such an interaction in this biological system. Accordingly, although cation-π interactions are well documented as important noncovalent forces in molecular recognition, the energetics of such interactions may be mitigated by other nonbonded interactions and solvation effects in protein-ligand associations.

Heterogeneous Hydration of p53/MDM2 Complex

Water-mediated interactions play critical roles in biomolecular recognition processes. Explicit solvent molecular dynamics (MD) simulations and the variational implicit-solvent model (VISM) are used to study those hydration properties during binding for the biologically important p53/MDM2 complex. Unlike simple model solutes, in such a realistic and heterogeneous solute–solvent system with both geometrical and chemical complexity, the local water distribution sensitively depends on nearby amino acid properties and the geometric shape of the protein. We show that the VISM can accurately describe the locations of high and low density solvation shells identified by the MD simulations and can explain them by a local coupling balance of solvent–solute interaction potentials and curvature. In particular, capillary transitions between local dry and wet hydration states in the binding pocket are captured for interdomain distance between 4 to 6 Å, right at the onset of binding. The underlying physical connection between geometry and polarity is illustrated and quantified. Our study offers a microscopic and physical insight into the heterogeneous hydration behavior of the biologically highly relevant p53/MDM2 system and demonstrates the fundamental importance of hydrophobic effects for biological binding processes. We hope our study can help to establish new design rules for drugs and medical substances.

Quantitative study of state switching in proteins using a single probe combined with trilinear decomposition

Excitation–Emission Matrix fluorescence (EEM) signals from the state switching of α-chymotrypsin can be resolved and investigated quantitatively by trilinear decomposition.

Correlating structure and energetics in protein-ligand interactions: paradigms and paradoxes

Predicting protein-binding affinities of small molecules, even closely related ones, is a formidable challenge in biomolecular recognition and medicinal chemistry. A thermodynamic approach to optimizing affinity in protein-ligand interactions requires knowledge and understanding of how altering the structure of a small molecule will be manifested in protein-binding enthalpy and entropy changes; however, there is a relative paucity of such detailed information. In this review, we examine two strategies commonly used to increase ligand potency. The first of these involves introducing a cyclic constraint to preorganize a small molecule in its biologically active conformation, and the second entails adding nonpolar groups to a molecule to increase the amount of hydrophobic surface that is buried upon binding. Both of these approaches are motivated by paradigms suggesting that protein-binding entropy changes should become more favorable, but paradoxes can emerge that defy conventional wisdom.

A simple atomic-level hydrophobicity scale reveals protein interfacial structure

Many amino acid residue hydrophobicity scales have been created in an effort to better understand and rapidly characterize water-protein interactions based only on protein structure and sequence. There is surprisingly low consistency in the ranking of residue hydrophobicity between scales, and their ability to provide insightful characterization varies substantially across subject proteins. All current scales characterize hydrophobicity based on entire amino acid residue units. We introduce a simple binary but atomic-level hydrophobicity scale that allows for the classification of polar and non-polar moieties within single residues, including backbone atoms. This simple scale is first shown to capture the anticipated hydrophobic character for those whole residues that align in classification among most scales. Examination of a set of protein binding interfaces establishes good agreement between residue-based and atomic-level descriptions of hydrophobicity for five residues, while the remaining residues produce discrepancies. We then show that the atomistic scale properly classifies the hydrophobicity of functionally important regions where residue-based scales fail. To illustrate the utility of the new approach, we show that the atomic-level scale rationalizes the hydration of two hydrophobic pockets and the presence of a void in a third pocket within a single protein and that it appropriately classifies all of the functionally important hydrophilic sites within two otherwise hydrophobic pores. We suggest that an atomic level of detail is, in general, necessary for the reliable depiction of hydrophobicity for all protein surfaces. The present formulation can be implemented simply in a manner no more complex than current residue-based approaches.

Picosecond-resolved solvent reorganization and energy transfer in biological and model cavities

Water molecules in hydrophobic biological cleft/cavities are of contemporary interest for the biomolecular structure and molecular recognition of hydrophobic ligands/drugs. Here, we have explored picosecond-resolved solvation dynamics of water molecules and associated polar amino acids in the hydrophobic cleft around Cys-34 position of Endogenous Serum Albumin (ESA). While site selective acrylodan labeling to Cys-34 allows us to probe solvation in the cleft, Förster resonance energy transfer (FRET) from intrinsic fluorescent amino acid Trp 214 to the extrinsic acrylodan probes structural integrity of the protein in our experimental condition. Temperature dependent solvation in the cleft clearly shows that the dynamics follows Arrhenius type behavior up to 60 °C, after which a major structural perturbation of the protein is evident. We have also monitored polarization gated dynamics of the acrylodan probe and FRET from Trp 214 to acrylodan at various temperatures. The dynamical behavior of the immediate environments around the probe acrylodan in the cleft has been compared with a model biomimetic cavity of a reverse micelle (w0 = 5). Using same fluorescent probe of acrylodan, we have checked the structural integrity of the model cavity at various temperatures using picosecond-resolved FRET from Trp to acrylodan in the cavity. We have also estimated possible distribution of donor-acceptor distances in the protein and reverse micelles. Our studies reveal that the energetics of the water molecules in the biological cleft is comparable to that in the model cavity indicating a transition from bound state to quasibound state, closely consistent with a recent MD simulation study.

Protein-ligand interactions: thermodynamic effects associated with increasing nonpolar surface area

Thermodynamic parameters were determined for complex formation between the Grb2 SH2 domain and Ac-pTyr-Xaa-Asn derived tripeptides in which the Xaa residue is an α,α-cycloaliphatic amino acid that varies in ring size from three- to seven-membered. Although the six- and seven-membered ring analogs are approximately equipotent, binding affinities of those having three- to six-membered rings increase incrementally with ring size because increasingly more favorable binding enthalpies dominate increasingly less favorable binding entropies, a finding consistent with an enthalpy-driven hydrophobic effect. Crystallographic analysis reveals that the only significant differences in structures of the complexes are in the number of van der Waals contacts between the domain and the methylene groups in the Xaa residues. There is a positive correlation between buried nonpolar surface area and binding free energy and enthalpy, but not with ΔC(p). Displacing a water molecule from a protein-ligand interface is not necessarily reflected in a favorable change in binding entropy. These findings highlight some of the fallibilities associated with commonly held views of relationships of structure and energetics in protein-ligand interactions and have significant implications for ligand design.

Structure of hydration water in proteins: a comparison of molecular dynamics simulations and database analysis

Hydration layer water molecules play important structural and functional roles in proteins. Despite being a critical component in biomolecular systems, characterizing the properties of hydration water poses a challenge for both experiments and simulations. In this context we investigate the local structure of hydration water molecules as a function of the distance from the protein and water molecules respectively in 188 high resolution protein structures and compare it with those obtained from molecular dynamics simulations. Tetrahedral order parameter of water in proteins calculated from previous and present simulation studies show that the potential of bulk water overestimates the average tetrahedral order parameter compared to those calculated from crystal structures. Hydration waters are found to be more ordered at a distance between the first and second solvation shell from the protein surface. The values of the order parameter decrease sharply when the water molecules are located very near or far away from the protein surface. At small water-water distance, the values of order parameter of water are very low. The average order parameter records a maximum value at a distance equivalent to the first solvation layer with respect to the water-water radial distribution and asymptotically approaches a constant value at large distances. Results from present analysis will help to get a better insight into structure of hydration water around proteins. The analysis will also help to improve the accuracy of water models on the protein surface.

Interfaces and hydrophobic interactions in receptor-ligand systems: A level-set variational implicit solvent approach

A model nanometer-sized hydrophobic receptor-ligand system in aqueous solution is studied by the recently developed level-set variational implicit solvent model (VISM). This approach is compared to all-atom computer simulations. The simulations reveal complex hydration effects within the (concave) receptor pocket, sensitive to the distance of the (convex) approaching ligand. The ligand induces and controls an intermittent switching between dry and wet states of the hosting pocket, which determines the range and magnitude of the pocket-ligand attraction. In the level-set VISM, a geometric free-energy functional of all possible solute-solvent interfaces coupled to the local dispersion potential is minimized numerically. This approach captures the distinct metastable states that correspond to topologically different solute-solvent interfaces, and thereby reproduces the bimodal hydration behavior observed in the all-atom simulation. Geometrical singularities formed during the interface relaxation are found to contribute significantly to the energy barrier between different metastable states. While the hydration phenomena can thus be explained by capillary effects, the explicit inclusion of dispersion and curvature corrections seems to be essential for a quantitative description of hydrophobically confined systems on nanoscales. This study may shed more light onto the tight connection between geometric and energetic aspects of biomolecular hydration and may represent a valuable step toward the proper interpretation of experimental receptor-ligand binding rates.

Dewetting-controlled binding of ligands to hydrophobic pockets

We report on a combined atomistic molecular dynamics simulation and implicit solvent analysis of a generic hydrophobic pocket-ligand (host-guest) system. The approaching ligand induces complex wetting-dewetting transitions in the weakly solvated pocket. The transitions lead to bimodal solvent fluctuations which govern magnitude and range of the pocket-ligand attraction. A recently developed implicit water model, based on the minimization of a geometric functional, captures the sensitive aqueous interface response to the concave-convex pocket-ligand configuration semiquantitatively.

Enhanced surface hydrophobicity by coupling of surface polarity and topography

We use atomistic computer simulation to explore the relationship between mesoscopic (liquid drop contact angle) and microscopic (surface atomic polarity) characteristics for water in contact with a model solid surface based on the structure of silica. We vary both the magnitude and direction of the solid surface polarity at the atomic scale and characterize the response of an aqueous interface in terms of the solvent molecular organization and contact angle. We show that when the topography and polarity of the surface act in concert with the asymmetric charge distribution of water, the hydrophobicity varies substantially and, further, can be maximal for a surface with significant polarity. The results suggest that patterning of a surface on several length scales, from atomic to mum lengths, can make important independent contributions to macroscopic hydrophobicity.

The hydrophobic effect: molecular dynamics simulations of water confined between extended hydrophobic and hydrophilic surfaces

Structural and dynamic properties of water confined between two parallel, extended, either hydrophobic or hydrophilic crystalline surfaces of n-alkane C(36)H(74) or n-alcohol C(35)H(71)OH, are studied by molecular dynamics simulations. Electron density profiles, directly compared with corresponding experimental data from x-ray reflectivity measurements, reveal a uniform weak de-wetting characteristic for the extended hydrophobic surface, while the hydrophilic surface is weakly wetted. These microscopic data are consistent with macroscopic contact angle measurements. Specific water orientation is present at both surfaces. The ordering is characteristically different between the surfaces and of longer range at the hydrophilic surface. Furthermore, the dynamic properties of water are different at the two surfaces and different from the bulk behavior. In particular, at the hydrophobic surface, time-correlation functions reveal that water molecules have characteristic diffusive behavior and orientational ordering due to the lack of hydrogen bonding interactions with the surface. These observations suggest that the altered dynamical properties of water in contact with extended hydrophobic surfaces together with a partial drying of the surfaces are more indicative of the hydrophobic effect than structural ordering, which we suggest to be independent of surface topology.

Simulation studies of the protein-water interface. I. Properties at the molecular resolution

We report molecular dynamics simulations of three globular proteins: ubiquitin, apo-calbindin D(9K), and the C-terminal SH2 domain of phospholipase C-gamma1 in explicit water. The proteins differ in their overall charge and fold type and were chosen to represent to some degree the structural variability found in medium-sized proteins. The length of each simulation was at least 15 ns, and larger than usual solvent boxes were used. We computed radial distribution functions, as well as orientational correlation functions about the surface residues. Two solvent shells could be clearly discerned about charged and polar amino acids. Near apolar amino acids the water density near such residues was almost devoid of structure. The mean residence time of water molecules was determined for water shells about the full protein, as well as for water layers about individual amino acids. In the dynamic properties, two solvent shells could be characterized as well. However, by comparison to simulations of pure water it could be shown that the influence of the protein reaches beyond 6 A, i.e., beyond the first two shells. In the first shell (r < or =3.5 A), the structural and dynamical properties of solvent waters varied considerably and depended primarily on the physicochemical properties of the closest amino acid side chain, with which the waters interact. By contrast, the solvent properties seem not to depend on the specifics of the protein studied (such as the net charge) or on the secondary structure element in which an amino acid is located. While differing considerably from the neat liquid, the properties of waters in the second solvation shell (3.5< r < or =6 A) are rather uniform; a direct influence from surface amino acids are already mostly shielded.

Structural effects of drying and rehydration for enzymes in soils: a kinetics-FTIR analysis of α-chymotrypsin adsorbed on montmorillonite

The effects of desiccation and rehydration cycles encountered by extracellular enzymes in soils are studied on -chymotrypsin adsorbed on montmorillonite. The controlled hygrometric FTIR cell used in this study enables to monitor drying and rehydration steps undergone by the -chymotrypsin-montmorillonite suspension or by the enzyme alone. Relative humidity (RH) determines the amount of deuterated water in the FTIR cell atmosphere. The molar water/protein ratio (W/P) as well as the conformational and solvation states of the enzyme have been determined using H/D exchange monitored by FTIR-transmission spectroscopy. When the W/P ratio decreases from 3500 to approximately 400, unfolding of beta-secondary structure in three different domains involves about 8% of the polypeptide backbone with respect to the most solvated states. Desiccation induces beta-unfolding, which opens channels allowing free vapor water molecules to diffuse into the enzyme at 15% RH. On drying to 0% RH, displacements of internal water (H2O) in the enzyme are demonstrated by reverse peptide isotopic exchanges (COND ==> CONH). Specific beta-structures, only formed in highly solvated states, sequester around 20 internal H2O molecules. Indeed, most of the unfolded secondary structures during the drying step are refolded at W/P approximately 1000 during rehydration. However, self-association hinders the recovery of the initial closed tertiary structure. The pD-dependent structural changes controlling inward and outward water diffusion are suppressed, whether the protein is initially in an adsorbed state or in solution. Changes in secondary structures encountered during desiccation/rehydration cycle are similar for the protein either free or in the adsorbed state. Thus domains that are unfolded by adsorption are not concerned by the desiccation/rehydration cycle.

Ultrafast surface hydration dynamics and expression of protein functionality: alpha -Chymotrypsin

We report studies of hydration dynamics at the surface of the enzyme protein bovine pancreatic alpha-chymotrypsin. The probe is the well known 1-anilinonaphthalene-8-sulfonate, which binds selectively in the native state of the protein, not the molten globule, as shown by x-ray crystallography. With femtosecond time resolution, we examined the hydration dynamics at two pHs, when the protein is physiologically in the inactive state (pH 3.6) or the active state (pH 6.7); the global structure and the binding site remain the same. The hydration correlation function, C(t), whose decay is governed by the rotational and translational motions of water molecules at the site, shows the behavior observed in this laboratory for other proteins, Subtilisin Carlsberg and Monellin, using the intrinsic amino acid tryptophan as a probe for surface hydration. However, the time scales and amplitudes vary drastically at the two pHs. For the inactive protein state, C(t) decays with an ultrafast component, close to bulk-type behavior, but 50% of the C(t) decays at a much slower rate, tau = 43 ps. In contrast, for the active state, the ultrafast component becomes dominant (90%) and the slow component changes to a faster decay, tau = 28 ps. These results indicate that in the active state water molecules in the hydration layer around the site have a high degree of mobility, whereas in the inactive state the water is more rigidly structured. For the substrate-enzyme complex, the function and dynamics at the probe site are correlated, and the relevance to the enzymatic action is clear.

Connecting local structure to interface formation: a molecular scale van der Waals theory of nonuniform liquids

This article reviews a new and general theory of nonuniform fluids that naturally incorporates molecular scale information into the classical van der Waals theory of slowly varying interfaces. The method optimally combines two standard approximations, molecular (mean) field theory to describe interface formation and linear response (or Gaussian fluctuation) theory to describe local structure. Accurate results have been found in many different applications in nonuniform simple fluids and these ideas may have important implications for the theory of hydrophobic interactions in water.

Solvation and hydration of proteins and nucleic acids: a theoretical view of simulation and experiment

Many theoretical, computational, and experimental techniques recently have been successfully used for description of the solvent distribution around macromolecules. In this Account, we consider recent developments in the areas of protein and nucleic acid solvation and hydration as seen by experiment, theory, and simulations. We find that in most cases not only the general phenomena of solvation but even local hydration patterns are more accurately discussed in the context of water distributions rather than individual molecules of water. While a few localized or high-residency waters are often associated with macromolecules in solution (or crystals from aqueous liquors), these are readily and accurately included in this more general description. The goal of this Account is to review the theoretical models used for the description of the interfacial solvent structure on the border near DNA and protein molecules. In particular, we hope to highlight the progress in this field over the past five years with a focus on comparison of simulation and experimental results.

Probing the origins of increased activity of the E22Q "Dutch" mutant Alzheimer's beta-amyloid peptide

The amyloid peptide congener A beta(10--35)-NH(2) is simulated in an aqueous environment in both the wild type (WT) and E22Q "Dutch" mutant forms. The origin of the noted increase in deposition activity resulting from the Dutch mutation is investigated. Multiple nanosecond time scale molecular dynamics trajectories were performed and analyzed using a variety of measures of the peptide's average structure, hydration, conformational fluctuations, and dynamics. The results of the study support the conclusions that 1) the E22Q mutant and WT peptide are both stable in "collapsed coil" conformations consistent with the WT structure of, J. Struct. Biol. 130:130--141); 2) the E22Q peptide is more flexible in solution, supporting early claims that its equilibrium structural fluctuations are larger than those of the WT peptide; and 3) the local E22Q mutation leads to a change in the first solvation layer in the region of the peptide's "hydrophobic patch," resulting in a more hydrophobic solvation of the mutant peptide. The simulation results support the view that the noted increase in activity due to the Dutch mutation results from an enhancement of the desolvation process that is an essential step in the aggregation of the peptide.