Is water structure around hydrophobic groups clathrate-like?

The term "clathrate structure" is quantified for solvation of nonpolar groups by enumerating hydrogen-bonded ring sizes both in the solvation shell and through the shell-bulk interface and comparing it to a bulk control using the ST4 water model. For clathrate-like structure to be evident, the distributions along the hydrophobic surface are expected to be dominated by pentagons, with significant depletion of hexagons and larger polygons. While the distribution in this region is indeed distinguished by a large number of pentagons, there are significant contributions from hexagons and larger rings as well. Calculated polygon distributions through the shell-bulk interface indicate that when water structure is highly cooperative along the hydrophobic surface, hydrogen-bonded pathways leading back into bulk are then reduced. These results are qualitatively consistent with the observation that hydrophobicity is proportional to the nonpolar solute surface area.     

A strategy for finding classes of minima on a hypersurface: implications for approaches to the protein folding problem.

Locating the native structure of a given protein is a task made difficult by the complexity of the potential energy hypersurface and by the huge number of local minima it contains. We have explored a strategy (the "antlion" method) for hypersurface modification that suppresses all minima but that of the native structure. Transferrable penalty functions with general applicability for modifying a hypersurface to retain the desired minimum are identified, and two blocked oligopeptides (alanine dipeptide and tetrapeptide) are used for specific numerical illustration of the dramatic simplification that ensues. In addition, an intermediary role for neural networks to manage some aspects of the antlion strategy applied to large polypeptides and proteins is introduced.     

Direct evidence for modified solvent structure within the hydration shell of a hydrophobic amino acid.

Neutron scattering experiments are used to determine scattering profiles for aqueous solutions of hydrophobic and hydrophilic amino acid analogs. Solutions of hydrophobic solutes show a shift in the main diffraction peak to smaller angle as compared with pure water, whereas solutions of hydrophilic solutes do not. The same difference for solutions of hydrophobic and hydrophilic side chains is also predicted by molecular dynamics simulations. The neutron scattering curves of aqueous solutions of hydrophobic amino acids at room temperature are qualitatively similar to differences between the liquid molecular structure functions measured for ambient and supercooled water. The nonpolar solute-induced expansion of water structure reported here is also complementary to recent neutron experiments where compression of aqueous solvent structure has been observed at high salt concentration.     

An isolated water droplet in the aqueous solution of a supramolecular tetrahedral cage

Water under nanoconfinement at ambient conditions has exhibited low-dimensional ice formation and liquid–solid phase transitions, but with structural and dynamical signatures that map onto known regions of water’s phase diagram. Using terahertz (THz) absorption spectroscopy and ab initio molecular dynamics, we have investigated the ambient water confined in a supramolecular tetrahedral assembly, and determined that a dynamically distinct network of 9 ± 1 water molecules is present within the nanocavity of the host. The low-frequency absorption spectrum and theoretical analysis of the water in the Ga₄L₆¹²⁻ host demonstrate that the structure and dynamics of the encapsulated droplet is distinct from any known phase of water. A further inference is that the release of the highly unusual encapsulated water droplet creates a strong thermodynamic driver for the high-affinity binding of guests in aqueous solution for the Ga₄L₆¹²⁻ supramolecular construct.

Supramolecular capsules create internal cavities that are thought to act like enzyme active sites (1). As aqueous enzymes provide inspiration for the design of supramolecular catalysts, one of the goals of supramolecular chemistry is the creation of synthetic “receptors” that have both a high affinity and a high selectivity for the binding of guests in water (2, 3). The Ga₄L₆¹²⁻ tetrahedral assembly formulated by Raymond and coworkers represents an excellent example of a water-soluble supramolecular cage that has provided host interactions that promotes guest encapsulation. Using steric interactions and electrostatic charge to chemically position the substrate while shielding the reaction from solvent, this host has been shown to provide enhanced reaction rates that approach the performance of natural biocatalysts (4–10). Moreover, aqueous solvation of the substrate, host, and encapsulated solvent also play an important role in the whole catalytic cycle. In particular, the driving forces that release water from the nanocage host to favor the direct binding with the substrate is thought to be a critical factor in successful catalysis, but is challenging to probe directly (7, 8, 11–14).In both natural and artificial nanometer-sized environments, confined water displays uniquely modified structure and dynamics with respect to the bulk liquid (15–18). Recently, these modified properties were also found to have significant implications for the mechanism and energetics of reactions taking place in confined water with respect to those observed in bulk aqueous solution (19–21). In a pioneering study on supramolecular assemblies, Cram and collaborators (22) concluded that the interior of those cages is a “new and unique phase of matter” for the incarcerated guests. In more recent studies, it was postulated that, similar to graphitic and zeolite nanopores (23, 24), confined water within supramolecular host cavities is organized in stable small clusters [(H₂O)_n, with n = 8 to 19] that are different from gas phase water clusters (25). In these studies, the hydrogen-bonded water clusters were reported to be mostly ice- or clathrate-like by X-ray and neutron diffraction in the solid state at both ambient and cryogenic temperatures (26–32). However, to the best of our knowledge, such investigations have not characterized the Ga₄L₆¹²⁻ supramolecular tetrahedral assembly in the liquid state near room temperature and pressure, where the [Ga₄L₆]¹²⁻ capsule can perform catalytic reactions (6, 8, 9).Here, we use terahertz (THz) absorption spectroscopy and ab initio molecular dynamics (AIMD) to characterize low-frequency vibrations and structural organization of water in the nanoconfined environment. THz is ideally suited to probe the intermolecular collective dynamics of the water hydrogen bond (HB) network with extremely high sensitivity, as illustrated for different phases of water (33–38), and for aqueous solutions of salts, osmolytes, alcohols, and amino acids (36, 39–42). The THz spectra of the water inside the nanocage has been quantitatively reproduced with AIMD, allowing us to confidently characterize the water network in the cage in order to provide a more complete dynamical, structural, and thermodynamic picture. We have determined that the spectroscopic signature of the confined water in the nanocage is a dynamically arrested state whose structure bears none of the features of water at any alternate thermodynamic state point such as pressurized liquid or ice. Our experimental and theoretical study provides insight into the role played by encapsulated water in supramolecular catalysis, creating a low entropy and low enthalpy water droplet readily displaced by a catalytic substrate.     

Theoretical study of the rhenium-alkane interaction in transition metal-alkane sigma-complexes

Metal-alkane binding energies have been calculated for [CpRe(CO)2](alkane) and [(CO)2M(C5H4)C[triple bond]C(C5H4)M(CO)2](alkane), where M = Re or Mn. Calculated binding energies were found to increase with the number of metal-alkane interaction sites. In all cases examined, the manganese-alkane binding energies were predicted to be significantly lower than those for the analogous rhenium-alkane complexes. The metal (Mn or Re)-alkane interaction was predicted to be primarily one of charge transfer, both from the alkane to the metal complex (70-80% of total charge transfer) and from the metal complex to the alkane (20-30% of the total charge transfer).     

Solution X-ray scattering as a probe of hydration-dependent structuring of aqueous solutions

We report on new X-ray solution scattering experiments and molecular dynamics simulations conducted for increasing solute concentrations of N-acetyl-amino acid-amides and -methylamides in water, for the amino acids leucine, glutamine, and glycine. As the concentration increases, the main diffraction peak of pure water at Q = 2.0 Å^-1 shifts to smaller angle for the larger leucine and glutamine amino acids, and a new diffraction peak grows in at Q 0.8 Å^-1 for only the hydrophobic amino acid leucine. The unaltered value of the peak position at Q 0.8 Å^-1 over a large concentration range suggests that a stable and ordered leucine solute–solute distribution is sustained. Simulations of the distributions of leucines in water that reproduce the experimental observable show that mono-dispersed to small molecular aggregates of two to six hydrophobic amino acids are formed, as opposed to complete segregation of the hydrophobic solutes into one large cluster. The scattering results for the hydrophobic leucine amino acid are contrasted with experiments and simulations of the model hydrophilic side chain glutamine and the model backbone glycine. The self-assembly process of protein folding modeled with these experiments, in particular the condensation to a hydrophobic core, shares similar issues with the desolvation phenomena that are important in drug discovery.     

Tetrahedral structure or chains for liquid water

It has been suggested, based on x-ray absorption spectroscopy (XAS) experiments on liquid water [Wernet, Ph., et al. (2004) Science 304, 995-999], that a condensed-phase water molecule's asymmetric electron density results in only two hydrogen bonds per water molecule on average. The larger implication of the XAS interpretation is that the conventional view of liquid water being a tetrahedrally coordinated random network is now replaced by a structural organization that instead strongly favors hydrogen-bonded water chains or large rings embedded in a weakly hydrogen-bonded disordered network. This work reports that the asymmetry of the hydrogen density exhibited in the XAS experiments agrees with reported x-ray scattering structure factors and intensities for Q > 6.5 A(-1). However, the assumption that the asymmetry in the hydrogen electron density does not fluctuate and is persistent in all local molecular liquid water environments is inconsistent with longer-ranged tetrahedral network signatures present in experimental x-ray scattering intensity and structure factor data for Q < 6.5 A(-1).     

First-principles, quantum-mechanical simulations of electron solvation by a water cluster

Despite numerous experiments and static electronic structure calculations, the nature of hydrated-electron clusters, (H2O)(n)(-), remains poorly understood. Here, we introduce a hybrid ab initio molecular dynamics scheme, balancing accuracy against feasibility, to simulate vibrational and photoelectron spectra of (H(2)O)(n)(-), treating all electrons quantum-mechanically. This methodology provides a computational tool for understanding the spectra of weakly bound and supramolecular anions and for elucidating the fingerprint of dynamics in these spectra. Simulations of (H2O)(4)(-) provide quantitative agreement with experimental spectra and furnish direct evidence of the nonequilibrium nature of the cluster ensemble that is probed experimentally. The simulations also provide an estimate of the cluster temperature (T approximately 150-200 K) that is not available from experiment alone. The "double acceptor" electron-binding motif is found to be highly stable with respect to thermal fluctuations, even at T = 300 K, whereas the extra electron stabilizes what would otherwise be unfavorable water configurations.     

Poly(L-alanine) as a universal reference material for understanding protein energies and structures.

We present a proposition, the "poly(L-alanine) hypothesis," which asserts that the native backbone geometry for any polypeptide or protein of M residues has a closely mimicking, mechanically stable, image in poly(L-alanine) of the same number of residues. Using a molecular mechanics force field to represent the relevant potential energy hypersurfaces, we have carried out calculations over a wide range of M values to show that poly(L-alanine) possesses the structural versatility necessary to satisfy the proposition. These include poly(L-alanine) representatives of minima corresponding to secondary and supersecondary structures, as well as poly(L-alanine) images for tertiary structures of the naturally occurring proteins bovine pancreatic trypsin inhibitor, crambin, ribonuclease A, and superoxide dismutase. The successful validation of the hypothesis presented in this paper indicates that poly(L-alanine) will serve as a good reference material in thermodynamic perturbation theory and calculations aimed at evaluating relative free energies for competing candidate tertiary structures in real polypeptides and proteins.