Water-protein interactions from high-resolution protein crystallography

The Effect of Lipopolysaccharide Core Oligosaccharide Size on the Electrostatic Binding of Antimicrobial Proteins to Models of the Gram Negative Bacterial Outer Membrane

Understanding the electrostatic interactions between bacterial membranes and exogenous proteins is crucial to designing effective antimicrobial agents against Gram-negative bacteria. Here we study, using neutron reflecometry under multiple isotopic contrast conditions, the role of the uncharged sugar groups in the outer core region of lipopolysaccharide (LPS) in protecting the phosphate-rich inner core region from electrostatic interactions with antimicrobial proteins. Models of the asymmetric Gram negative outer membrane on silicon were prepared with phopshatidylcholine (PC) in the inner leaflet (closest to the silicon), whereas rough LPS was used to form the outer leaflet (facing the bulk solution). We show how salt concentration can be used to reversibly alter the binding affinity of a protein antibiotic colicin N (ColN) to the anionic LPS confirming that the interaction is electrostatic in nature. By examining the interaction of ColN with two rough LPS types with different-sized core oligosaccharide regions we demonstrate the role of uncharged sugars in blocking short-range electrostatic interactions between the cationic antibiotics and the vulnerable anionic phosphate groups.

Lipid Bilayer Membrane Perturbation by Embedded Nanopores: A Simulation Study

A macromolecular nanopore inserted into a membrane may perturb the dynamic organization of the surrounding lipid bilayer. To better understand the nature of such perturbations, we have undertaken a systematic molecular dynamics simulation study of lipid bilayer structure and dynamics around three different classes of nanopore: a carbon nanotube, three related cyclic peptide nanotubes differing in the nature of their external surfaces, and a model of a β-barrel nanopore protein. Periodic spatial distributions of several lipid properties as a function of distance from the nanopore were observed. This was especially clear for the carbon nanotube system, for which the density of lipids, the bilayer thickness, the projection of lipid head-to-tail vectors onto the membrane plane, and lipid lateral diffusion coefficients exhibited undulatory behavior as a function of the distance from the surface of the channel. Overall, the differences in lipid behavior as a function of the nanopore structure reveal local adaptation of the bilayer structure and dynamics to different embedded nanopore structures. Both the local structure and dynamic behavior of lipids around membrane-embedded nanopores are sensitive to the geometry and nature of the outer surface of the macromolecule/molecular assembly forming the pore.

Interaction between bound water molecules and local protein structures: A statistical analysis of the hydrogen bond structures around bound water molecules

Water molecules play an important role in protein folding and protein interactions through their structural association with proteins. Examples of such structural association can be found in protein crystal structures, and can often explain protein functionality in the context of structure. We herein report the systematic analysis of the local structures of proteins interacting with water molecules, and the characterization of their geometric features. We first examined the interaction of water molecules with a large local interaction environment by comparing the preference of water molecules in three regions, namely, the protein-protein interaction (PPI) interfaces, the crystal contact (CC) interfaces, and the non-interfacial regions. High preference of water molecules to the PPI and CC interfaces was found. In addition, the bound water on the PPI interface was more favorably associated with the complex interaction structure, implying that such water-mediated structures may participate in the shaping of the PPI interface. The pairwise water-mediated interaction was then investigated, and the water-mediated residue-residue interaction potential was derived. Subsequently, the types of polar atoms surrounding the water molecules were analyzed, and the preference of the hydrogen bond acceptor was observed. Furthermore, the geometries of the structures interacting with water were analyzed, and it was found that the major structure on the protein surface exhibited planar geometry rather than tetrahedral geometry. Several previously undiscovered characteristics of water-protein interactions were unfolded in this study, and are expected to lead to a better understanding of protein structure and function.

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.

Structure of the extracellular domain of matrix protein 2 of influenza A virus in complex with a protective monoclonal antibody

The extracellular domain of influenza A virus matrix protein 2 (M2e) is conserved and is being evaluated as a quasiuniversal influenza A vaccine candidate. We describe the crystal structure at 1.6 Å resolution of M2e in complex with the Fab fragment of an M2e-specific monoclonal antibody that protects against influenza A virus challenge. This antibody binds M2 expressed on the surfaces of cells infected with influenza A virus. Five out of six complementary determining regions interact with M2e, and three highly conserved M2e residues are critical for this interaction. In this complex, M2e adopts a compact U-shaped conformation stabilized in the center by the highly conserved tryptophan residue in M2e. This is the first description of the three-dimensional structure of M2e.

IMPORTANCE

M2e of influenza A is under investigation as a universal influenza A vaccine, but its three-dimensional structure is unknown. We describe the structure of M2e stabilized with an M2e-specific monoclonal antibody that recognizes natural M2. We found that the conserved tryptophan is positioned in the center of the U-shaped structure of M2e and stabilizes its conformation. The structure also explains why previously reported in vivo escape viruses, selected with a similar monoclonal antibody, carried proline residue substitutions at position 10 in M2.

Hydration of protein-RNA recognition sites

We investigate the role of water molecules in 89 protein–RNA complexes taken from the Protein Data Bank. Those with tRNA and single-stranded RNA are less hydrated than with duplex or ribosomal proteins. Protein–RNA interfaces are hydrated less than protein–DNA interfaces, but more than protein–protein interfaces. Majority of the waters at protein–RNA interfaces makes multiple H-bonds; however, a fraction do not make any. Those making H-bonds have preferences for the polar groups of RNA than its partner protein. The spatial distribution of waters makes interfaces with ribosomal proteins and single-stranded RNA relatively ‘dry’ than interfaces with tRNA and duplex RNA. In contrast to protein–DNA interfaces, mainly due to the presence of the 2′OH, the ribose in protein–RNA interfaces is hydrated more than the phosphate or the bases. The minor groove in protein–RNA interfaces is hydrated more than the major groove, while in protein–DNA interfaces it is reverse. The strands make the highest number of water-mediated H-bonds per unit interface area followed by the helices and the non-regular structures. The preserved waters at protein–RNA interfaces make higher number of H-bonds than the other waters. Preserved waters contribute toward the affinity in protein–RNA recognition and should be carefully treated while engineering protein–RNA interfaces.

Role of water in protein folding, oligomerization, amyloidosis and miniprotein

The essential involvement of water in most fundamental extra-cellular and intracellular processes of proteins is critically reviewed and evaluated in this article. The role of water in protein behavior displays structural ambivalence; it can protect the disordered peptide-chain by hydration or helps the globular chain-folding, but promotes also the protein aggregation, as well (see: diseases). A variety of amyloid diseases begins as benign protein monomers but develops then into toxic amyloid aggregates of fibrils. Our incomplete knowledge of this process emphasizes the essential need to reveal the principles of governing this oligomerization. To understand the biophysical basis of the simpler in vitro amyloid formation may help to decipher also the in vivo way. Nevertheless, to ignore the central role of the water's effect among these events means to receive an uncompleted picture of the true phenomenon. Therefore this review represents a stopgap role, because the most published studies--with a few exceptions--have been neglected the crucial importance of water in the protein research. The following questions are discussed from the water's viewpoint: (i) interactions between water and proteins, (ii) protein hydration/dehydration, (iii) folding of proteins and miniproteins, (iv) peptide/protein oligomerization, and (v) amyloidosis.

Blind prediction of interfacial water positions in CAPRI

We report the first assessment of blind predictions of water positions at protein-protein interfaces, performed as part of the critical assessment of predicted interactions (CAPRI) community-wide experiment. Groups submitting docking predictions for the complex of the DNase domain of colicin E2 and Im2 immunity protein (CAPRI Target 47), were invited to predict the positions of interfacial water molecules using the method of their choice. The predictions-20 groups submitted a total of 195 models-were assessed by measuring the recall fraction of water-mediated protein contacts. Of the 176 high- or medium-quality docking models-a very good docking performance per se-only 44% had a recall fraction above 0.3, and a mere 6% above 0.5. The actual water positions were in general predicted to an accuracy level no better than 1.5 Å, and even in good models about half of the contacts represented false positives. This notwithstanding, three hotspot interface water positions were quite well predicted, and so was one of the water positions that is believed to stabilize the loop that confers specificity in these complexes. Overall the best interface water predictions was achieved by groups that also produced high-quality docking models, indicating that accurate modelling of the protein portion is a determinant factor. The use of established molecular mechanics force fields, coupled to sampling and optimization procedures also seemed to confer an advantage. Insights gained from this analysis should help improve the prediction of protein-water interactions and their role in stabilizing protein complexes.

An experimental point of view on hydration/solvation in halophilic proteins

Protein-solvent interactions govern the behaviors of proteins isolated from extreme halophiles. In this work, we compared the solvent envelopes of two orthologous tetrameric malate dehydrogenases (MalDHs) from halophilic and non-halophilic bacteria. The crystal structure of the MalDH from the non-halophilic bacterium Chloroflexus aurantiacus (Ca MalDH) solved, de novo, at 1.7 Å resolution exhibits numerous water molecules in its solvation shell. We observed that a large number of these water molecules are arranged in pentagonal polygons in the first hydration shell of Ca MalDH. Some of them are clustered in large networks, which cover non-polar amino acid surface. The crystal structure of MalDH from the extreme halophilic bacterium Salinibacter ruber (Sr) solved at 1.55 Å resolution shows that its surface is strongly enriched in acidic amino acids. The structural comparison of these two models is the first direct observation of the relative impact of acidic surface enrichment on the water structure organization between a halophilic protein and its non-adapted counterpart. The data show that surface acidic amino acids disrupt pentagonal water networks in the hydration shell. These crystallographic observations are discussed with respect to halophilic protein behaviors in solution

Specific ions modulate diffusion dynamics of hydration water on lipid membrane surfaces

Effects of specific ions on the local translational diffusion of water near large hydrophilic lipid vesicle surfaces were measured by Overhauser dynamic nuclear polarization (ODNP). ODNP relies on an unpaired electron spin-containing probe located at molecular or surface sites to report on the dynamics of water protons within ∼10 Å from the spin probe, which give rise to spectral densities for electron–proton cross-relaxation processes in the 10 GHz regime. This pushes nuclear magnetic resonance relaxometry to more than an order of magnitude higher frequencies than conventionally feasible, permitting the measurement of water moving with picosecond to subnanosecond correlation times. Diffusion of water within ∼10 Å of, i.e., up to ∼3 water layers around the spin probes located on hydrophilic lipid vesicle surfaces is ∼5 times retarded compared to the bulk water translational diffusion. This directly reflects on the activation barrier for surface water diffusion, i.e., how tightly water is bound to the hydrophilic surface and surrounding waters. We find this value to be modulated by the presence of specific ions in solution, with its order following the known Hofmeister series. While a molecular description of how ions affect the hydration structure at the hydrophilic surface remains to be answered, the finding that Hofmeister ions directly modulate the surface water diffusivity implies that the strength of the hydrogen bond network of surface hydration water is directly modulated on hydrophilic surfaces.

KOTOBUKI-1 apparatus for cryogenic coherent X-ray diffraction imaging

We have developed an experimental apparatus named KOTOBUKI-1 for use in coherent X-ray diffraction imaging experiments of frozen-hydrated non-crystalline particles at cryogenic temperature. For cryogenic specimen stage with small positional fluctuation for a long exposure time of more than several minutes, we here use a cryogenic pot cooled by the evaporation cooling effect for liquid nitrogen. In addition, a loading device is developed to bring specimens stored in liquid nitrogen to the specimen stage in vacuum. The apparatus allows diffraction data collection for frozen-hydrated specimens at 66 K with a positional fluctuation of less than 0.4 μm and provides an experimental environment to easily exchange specimens from liquid nitrogen storage to the specimen stage. The apparatus was developed and utilized in diffraction data collection of non-crystalline particles with dimensions of μm from material and biological sciences, such as metal colloid particles and chloroplast, at BL29XU of SPring-8. Recently, it has been applied for single-shot diffraction data collection of non-crystalline particles with dimensions of sub-μm using X-ray free electron laser at BL3 of SACLA.

Hydration-state change of horse heart cytochrome c corresponding to trifluoroacetic-acid-induced unfolding

We investigate the hydration state of horse-heart cytochrome c (hh cyt c) in the unfolding process induced by trifluoroacetic acid (TFA). The conformation of hh cyt c changes from the native (N) state (2.9 < pH < 6.0) to the acid-unfolded (U(A)) state (1.7 < pH < 2.0) to the acid-induced molten globule (A) state (pH ∼1.2). Hydration properties of hh cyt c during this process are measured at 20°C by high-resolution dielectric relaxation (DR) spectroscopy, UV-vis absorbance, and circular dichroism spectroscopy. Constrained water of hh cyt c is observed at every pH as an ∼5-GHz Debye component (DC) (DR time, τ(D) ∼30 ps) and its DR amplitude (DRA) is increased by 77% upon N-to-U(A) transition, when pH changes from 6.0 to 2.0. Even in the N state, the DRA of the constrained-water component is found to be increased by 22% with decreasing pH from 6.0 to 2.9, suggesting an increase in the accessible surface area of native hh cyt c. Moreover, hypermobile water around native hh cyt c is detected at pH 6.0 as a 19-GHz DC (τ(D) ∼ 8.4 ps <τ(DW) = 9.4 ps), but is not found at other pH values. The DRA signal of constrained water is found to return to the pH 2.9 (N-state) level upon U(A)-to-A transition. Fast-response water (slightly slower than bulk) around A-state hh cyt c is detected at pH 1.2, and this suggests some accumulation of TFA(-) ions around the peptide chain. Thus, this high-resolution DR spectroscopy study reveals that hh cyt c exhibits significant hydration-state change in the TFA-unfolding process.

Water makes the difference: rearrangement of water solvation layer triggers non-additivity of functional group contributions in protein-ligand binding

The binding of four congeneric peptide-like thermolysin inhibitors has been studied by high-resolution crystal structure analysis and isothermal titration calorimetry. The ligands differ only by a terminal carboxylate and/or methyl group. A surprising non-additivity of functional group contributions for the carboxylate and/or methyl groups is detected. Adding the methyl first and then the carboxylate group results in a small Gibbs free energy increase and minor enthalpy/entropy partitioning for the first modification, whereas the second involves a strong affinity increase combined with large enthalpy/entropy changes. However, first adding the carboxylate and then the methyl group yields reverse effects: the acidic group attachment now causes minor effects, whereas the added methyl group provokes large changes. As all crystal structures show virtually identical binding modes, affinity changes are related to rearrangements of the first solvation layer next to the S(2)' pocket. About 20-25 water molecules are visible next to the studied complexes. The added COO(-) groups perturb the local water network in both carboxylated complexes, and the attached methyl groups provide favorable interaction sites for water molecules. Apart from one example, a contiguously connected water network between protein and ligand functional groups is observed in all complexes. In the complex with the carboxylated ligand, which still lacks the terminal methyl group, the water network is unfavorably ruptured. This results in a surprising thermodynamic signature showing only a minor affinity increase upon COO(-) group attachment. Because the further added methyl group provides a favorable interaction site for water, the network can be reestablished, and a strong affinity increase with a large enthalpy/entropy signature is then detected.

Humidity-controlled preparation of frozen-hydrated biological samples for cryogenic coherent x-ray diffraction microscopy

Coherent x-ray diffraction microscopy (CXDM) has the potential to visualize the structures of micro- to sub-micrometer-sized biological particles, such as cells and organelles, at high resolution. Toward advancing structural studies on the functional states of such particles, here, we developed a system for the preparation of frozen-hydrated biological samples for cryogenic CXDM experiments. The system, which comprised a moist air generator, microscope, micro-injector mounted on a micromanipulator, custom-made sample preparation chamber, and flash-cooling device, allowed for the manipulation of sample particles in the relative humidity range of 20%-94%rh at 293 K to maintain their hydrated and functional states. Here, we report the details of the system and the operation procedure, including its application to the preparation of a frozen-hydrated chloroplast sample. Sample quality was evaluated through a cryogenic CXDM experiment conducted at BL29XUL of SPring-8. Taking the performance of the system and the quality of the sample, the system was suitable to prepare frozen-hydrated biological samples for cryogenic CXDM experiments.

Progressive dry-core-wet-rim hydration trend in a nested-ring topology of protein binding interfaces

BACKGROUND

Water is an integral part of protein complexes. It shapes protein binding sites by filling cavities and it bridges local contacts by hydrogen bonds. However, water molecules are usually not included in protein interface models in the past, and few distribution profiles of water molecules in protein binding interfaces are known.

RESULTS

In this work, we use a tripartite protein-water-protein interface model and a nested-ring atom re-organization method to detect hydration trends and patterns from an interface data set which involves immobilized interfacial water molecules. This data set consists of 206 obligate interfaces, 160 non-obligate interfaces, and 522 crystal packing contacts. The two types of biological interfaces are found to be drier than the crystal packing interfaces in our data, agreeable to a hydration pattern reported earlier although the previous definition of immobilized water is pure distance-based. The biological interfaces in our data set are also found to be subject to stronger water exclusion in their formation. To study the overall hydration trend in protein binding interfaces, atoms at the same burial level in each tripartite protein-water-protein interface are organized into a ring. The rings of an interface are then ordered with the core atoms placed at the middle of the structure to form a nested-ring topology. We find that water molecules on the rings of an interface are generally configured in a dry-core-wet-rim pattern with a progressive level-wise solvation towards to the rim of the interface. This solvation trend becomes even sharper when counterexamples are separated.

CONCLUSIONS

Immobilized water molecules are regularly organized in protein binding interfaces and they should be carefully considered in the studies of protein hydration mechanisms.

Hydration shells of molecules in molecular association: A mechanism for biomolecular recognition

It has become clear that water should not be treated as an inert environment, but rather as an integral and active component of molecules. Here, we consider molecules and their hydration shells together as single entities. We show that: (1) the rate of association of molecules should be determined by the energetic barriers arising from interactions between their hydration shells; (2) replacing non-polar atoms of molecular surfaces with polar atoms increases these barriers; (3) reduction of the hydration shells during molecular association is the driving force for association not only of non-polar, but of polar molecules as well; (4) in most cases the dehydration of polar atoms during molecular association thermodynamically counteracts association; (5) on balance the thermodynamic stability of associated complexes is basically determined by the action of these two opposing factors: reduction of the hydration shells and dehydration of polar atoms; (6) molecular crowding reduces the energetic barriers counteracting association and changes the thermodynamic stability of associated complexes. These results lead to a mechanism for biomolecular recognition in the context of which the formation of unique structures is provided by rapidly forming kinetic traps with a biologically necessary lifetime but with a marginal thermodynamic stability. The mechanism gives definitive answers to questions concerning the heart of specific interactions between biomolecules, their folding and intracellular organization. Predictions are given that can be subjected to direct experimental tests.

Immobilized Candida antarctica lipase B: Hydration, stripping off and application in ring opening polyester synthesis

This work reviews the stripping off, role of water molecules in activity, and flexibility of immobilized Candida antarctica lipase B (CALB). Employment of CALB in ring opening polyester synthesis emphasizing on a polylactide is discussed in detail. Execution of enzymes in place of inorganic catalysts is the most green alternative for sustainable and environment friendly synthesis of products on an industrial scale. Robust immobilization and consequently performance of enzyme is the essential objective of enzyme application in industry. Water bound to the surface of an enzyme (contact class of water molecules) is inevitable for enzyme performance; it controls enzyme dynamics via flexibility changes and has intensive influence on enzyme activity. The value of pH during immobilization of CALB plays a critical role in fixing the active conformation of an enzyme. Comprehensive selection of support and protocol can develop a robust immobilized enzyme thus enhancing its performance. Organic solvents with a log P value higher than four are more suitable for enzymatic catalysis as these solvents tend to strip away very little of the enzyme surface bound water molecules. Alternatively ionic liquid can work as a more promising reaction media. Covalent immobilization is an exclusively reliable technique to circumvent the leaching of enzymes and to enhance stability. Activated polystyrene nanoparticles can prove to be a practical and economical support for chemical immobilization of CALB. In order to reduce the E-factor for the synthesis of biodegradable polymers; enzymatic ring opening polyester synthesis (eROPS) of cyclic monomers is a more sensible route for polyester synthesis. Synergies obtained from ionic liquids and immobilized enzyme can be much effective eROPS.

Biological membranes: the importance of molecular detail

Are lipid interactions with membrane proteins best described in terms of the physical properties of the lipid bilayer or in terms of direct molecular interactions between particular lipid molecules and particular sites on a protein? A molecular interpretation is more challenging because it requires detailed knowledge of the 3D structure of a membrane protein, but recent studies have suggested that a molecular interpretation is necessary. Here, the idea is explored that lipid molecules modify the ways that transmembrane α-helices pack into bundles, by penetrating between the helices and by binding into clefts between the helices, and that these effects on helix packing will modulate the activity of a membrane protein.

Application of a real-space three-dimensional image reconstruction method in the structural analysis of noncrystalline biological macromolecules enveloped by water in coherent x-ray diffraction microscopy

Coherent x-ray diffraction microscopy is a novel technique in the structural analyses of particles that are difficult to crystallize, such as the biological particles composing living cells. As water is indispensable for maintaining particles in functional structures, sufficient hydration of targeted particles is required during sample preparation for diffraction microscopy experiments. However, the water enveloping particles also contributes significantly to the diffraction patterns and reduces the electron-density contrast of the sample particles. In this study, we propose a protocol for the structural analyses of particles in water by applying a three-dimensional reconstruction method in real space for the projection images phase-retrieved from diffraction patterns, together with a developed density modification technique. We examined the feasibility of the protocol through three simulations involving a protein molecule in a vacuum, and enveloped in either a droplet or a cube-shaped water. The simulations were carried out for the diffraction patterns in the reciprocal planes normal to the incident x-ray beam. This assumption and the simulation conditions corresponded to experiments using x-ray wavelengths of shorter than 0.03 Å. The analyses demonstrated that our protocol provided an interpretable electron-density map. Based on the results, we discuss the advantages and limitations of the proposed protocol and its practical application for experimental data. In particular, we examined the influence of Poisson noise in diffraction patterns on the reconstructed three-dimensional electron density in the proposed protocol.

A few low-frequency normal modes predominantly contribute to conformational responses of hen egg white lysozyme in the tetragonal crystal to variations of molecular packing controlled by environmental humidity

The structures of proteins in crystals are fixed by molecular interactions with neighboring molecules, except in non-contacting flexible regions. Thus, it is difficult to imagine what conformational changes occur in solution. However, if molecular interactions can be changed by manipulating molecular packing in crystals, it may be possible to visualize conformational responses of proteins at atomic resolution by diffraction experiments. For this purpose, it is suitable to control the molecular packing in protein crystals by changing the volume of solvent channels through variation of the environmental relative humidity. Here, we studied conformational responses of hen egg white lysozyme (HEWL) in the tetragonal crystal by X-ray diffraction experiments using a humidity-control apparatus, which provided air flow of 20-98%rh at 298 K. First, we monitored the lattice parameters and crystalline order during dehydration and rehydration of HEWL crystal between 61 and 94%rh at 300 K. Then two crystal structures at a resolution of 2.1 Å using diffraction data obtained at 84.2 and 71.9%rh were determined to discuss the conformational responses of HEWL against the external perturbation induced by changes in molecular packing. The structure at 71.9%rh displayed a closure movement that was likely induced by the molecular contacts formed during dehydration and could be approximated by ten low-frequency normal modes for the crystal structure obtained at 84.2%rh. In addition, we observed reorganization of hydration structures at the molecular interfaces between symmetry neighbors. These findings suggest that humidity-controlled X-ray crystallography is an effective tool to investigate the responses of inherent intramolecular motions of proteins to external perturbations.

Efficient synthesis of optically pure alcohols by asymmetric hydrogen-transfer biocatalysis: application of engineered enzymes in a 2-propanol-water medium

We describe an efficient method for producing both enantiomers of chiral alcohols by asymmetric hydrogen-transfer bioreduction of ketones in a 2-propanol (IPA)-water medium with E. coli biocatalysts expressing phenylacetaldehyde reductase (PAR: wild-type and mutant enzymes) from Rhodococcus sp. ST-10 and alcohol dehydrogenase from Leifsonia sp. S749 (LSADH). We also describe the detailed properties of mutant PARs, Sar268, and HAR1, which were engineered to have high activity and productivity in media composed of polar organic solvent and water, and the construction of three-dimensional structure of PAR by homology modeling. The K(m) and V(max) values for some substrates and the substrate specificity of mutant PARs were quite different from those of wild-type PAR. The results well explained the increased productivity of engineered PARs in IPA-water medium.

Site-resolved measurement of water-protein interactions by solution NMR

The interactions of biological macromolecules with water are fundamental to their structure, dynamics and function. Historically, characterization of the location and residence times of hydration waters of proteins in solution has been quite difficult. Confining proteins within the nanoscale interior of a reverse micelle slows water dynamics, allowing global protein-water interactions to be detected using nuclear magnetic resonance techniques. Complications that normally arise from hydrogen exchange and long-range dipolar coupling are overcome by the nature of the reverse micelle medium. Characterization of the hydration of ubiquitin demonstrates that encapsulation within a reverse micelle allows detection of dozens of hydration waters. Comparison of nuclear Overhauser effects obtained in the laboratory and rotating frames indicate a considerable range of hydration water dynamics is present on the protein surface. In addition, an unprecedented clustering of different hydration-dynamics classes of sites is evident.

The role of structured water in mediating general anesthetic action on alpha4beta2 nAChR

Water is an essential component for many biological processes. Pauling proposed that water might play a critical role in general anesthesia by forming water clathrates around anesthetic molecules. To examine potential involvement of water in general anesthesia, we analyzed water within alpha4beta2 nAChR, a putative protein target hypersensitive to volatile anesthetics. Experimental structure-derived closed- and open-channel nAChR systems in a fully hydrated lipid bilayer were examined using all-atom molecular dynamics simulations. At the majority of binding sites in alpha4beta2 nAChR, halothane replaced the slow-exchanging water molecules and caused a regional water population decrease. Only two binding sites had an increased quantity of water in the presence of halothane, where water arrangements resemble clathrate-like structures. The small number of such clathrate-like water clusters suggests that the formation of water clathrates is unlikely to be a primary cause for anesthesia. Despite the decrease in water population at most of the halothane binding sites, the number of sites that can be occupied transiently by water is actually increased in the presence of halothane. Many of these water sites were located between two subunits or in regions containing agonist binding sites or critical structural elements for transducing agonist binding to channel gating. Changes in water sites in the presence of halothane affected water-mediated protein-protein interactions and the protein dynamics, which can have direct impact on protein function. Collectively, water contributes to the action of anesthetics in proteins by mediating interactions between protein subunits and altering protein dynamics, instead of forming water clathrates around anesthetics.

Probability distributions of hydration water molecules around polar protein atoms obtained by a database analysis

Hydration structures on protein surfaces are visualized by high-resolution cryogenic X-ray crystallography. We calculated the probability distributions of 4,831,570 hydration water molecules found around the 4,214,227 polar atoms in main chains and hydrophilic side chains from the 17,984 crystal structures in the Protein Data Bank. The structures are refined using the diffraction data collected below 150 K and at resolutions of better than 2.2 A. The calculated distributions were nonrandom but condensed into a few clusters. The clusters were decomposed into the distance and angular distributions by viewing from the polar coordinate system. The major peaks in the clusters were almost located along the directions of the N-H and O-H bonds or the lone pairs of oxygen atoms. The Gaussian fitting method was applied for the distribution profiles to evaluate quantitatively the peak positions and the widths. The parameters characterizing the distributions apparently depended on the hydrogen-bond partners of water molecules and on the modes whether the water molecules acted as donors or acceptors of protons. This led to propose the different roles of NH(n) (n = 1, 3), OH, and CO groups in protein hydration and possible in protein-ligand and protein-protein interaction: While C horizontal lineO groups appear to control the H-bond distances, NH(n) groups likely limit the angular range of H-bonds. The OH groups have both characteristics. In addition, it was also demonstrated that polar protein atoms were arranged to satisfy the tetrahedral hydrogen-bond geometry of water molecules, suggesting essential roles of water molecules in the folding process and in the stabilization of protein structures. These probability distributions are probably one of fundamental data to better understand the roles of hydration water molecules in the folding process and the stability of proteins in solution.

Charge density-dependent modifications of hydration shell waters by Hofmeister ions

Gadolinium (Gd(3+)) vibronic sideband luminescence spectroscopy (GVSBLS) is used to probe, as a function of added Hofmeister series salts, changes in the OH stretching frequency derived from first-shell waters of aqueous Gd(3+) and of Gd(3+) coordinated to three different types of molecules: (i) a chelate (EDTA), (ii) structured peptides (mSE3/SE2) of the lanthanide-binding tags (LBTs) family with a single high-affinity binding site, and (iii) a calcium-binding protein (calmodulin) with four binding sites. The vibronic sideband (VSB) corresponding to the OH stretching mode of waters coordinated to Gd(3+), whose frequency is inversely correlated with the strength of the hydrogen bonding to neighboring waters, exhibits an increase in frequency when Gd(3+) becomes coordinated to either EDTA, calmodulin, or mSE3 peptide. In all of these cases, the addition of cation chloride or acetate salts to the solution increases the frequency of the vibronic band originating from the OH stretching mode of the coordinated waters in a cation- and concentration-dependent fashion. The cation dependence of the frequency increase scales with charge density of the cations, giving rise to an ordering consistent with the Hofmeister ordering. On the other hand, water Raman spectroscopy shows no significant change upon addition of these salts. Additionally, it is shown that the cation effect is modulated by the specific anion used. The results indicate a mechanism of action for Hofmeister series ions in which hydrogen bonding among hydration shell waters is modulated by several factors. High charge density cations sequester waters in a configuration that precludes strong hydrogen bonding to neighboring waters. Under such conditions, anion effects emerge as anions compete for hydrogen-bonding sites with the remaining free waters on the surface of the hydration shell. The magnitude of the anion effect is both cation and Gd(3+)-binding site specific.

Water polygons in high-resolution protein crystal structures

We have analyzed the interstitial water (ISW) structures in 1500 protein crystal structures deposited in the Protein Data Bank that have greater than 1.5 A resolution with less than 90% sequence similarity with each other. We observed varieties of polygonal water structures composed of three to eight water molecules. These polygons may represent the time- and space-averaged structures of "stable" water oligomers present in liquid water, and their presence as well as relative population may be relevant in understanding physical properties of liquid water at a given temperature. On an average, 13% of ISWs are localized enough to be visible by X-ray diffraction. Of those, averages of 78% are water molecules in the first water layer on the protein surface. Of the localized ISWs beyond the first layer, almost half of them form water polygons such as trigons, tetragons, as well as expected pentagons, hexagons, higher polygons, partial dodecahedrons, and disordered networks. Most of the octagons and nanogons are formed by fusion of smaller polygons. The trigons are most commonly observed. We suggest that our observation provides an experimental basis for including these water polygon structures in correlating and predicting various water properties in liquid state.

Strategies to model the near-solute solvent molecular density/polarization

The solvent molecular distribution significantly affects the behavior of the solute molecules and is thus important in studying many biological phenomena. It can be described by the solvent molecular density distribution, g, and the solvent electric dipole distribution, p. The g and p can be computed directly by counting the number of solvent molecules/dipoles in a microscopic volume centered at r during a simulation or indirectly from the mean force F and electrostatic field E acting on the solvent molecule at r, respectively. However, it is not clear how the g and p derived from simulations depend on the solvent molecular center or the solute charge and if the g(F) and p(E) computed from the mean force and electric field acting on the solvent molecule, respectively, could reproduce the corresponding g and p obtained by direct counting. Hence, we have computed g, p, g(F), and p(E) using different water centers from simulations of a solute atom of varying charge solvated in TIP3P water. The results show that g(F) and p(E) can reproduce the g and p obtained using a given count center. This implies that rather than solving the coordinates of each water molecule by MD simulations, the distribution of water molecules could be indirectly obtained from analytical formulas for the mean force F and electrostatic field E acting on the solvent molecule at r. Furthermore, the dependence of the g and p distributions on the solute charge revealed provides an estimate of the change in g and p surrounding a biomolecule upon a change in its conformation.

Molecular recognition in biomolecules studied by statistical-mechanical integral-equation theory of liquids

Recent progress in the theory of molecular recognition in biomolecules is reviewed, which has been made based on the statistical mechanics of liquids or the RISM/3D-RISM theory during the last five years in the authors' group. The method requires just the structure of protein and the potential energy parameters for the biomolecules and solutions as inputs. The calculation is carried out in two steps. The first step is to obtain the pair correlation functions for solutions consisting of water and ligands based on the RISM theory. Then, given the pair correlation functions prepared in the first step, we calculate the 3D-distribution functions of water and ligands around and inside protein based on the 3D-RISM theory. The molecular recognition of a ligand by the protein is realized by the 3D-distribution functions: if one finds some conspicuous peaks in the distribution of a ligand inside protein, then the ligand is regarded as "recognized" by the protein. Some biochemical processes are investigated, which are intimately related to the molecular recognition of small ligands including water, noble gases, and ions by a protein.

Incorporating dipolar solvents with variable density in Poisson-Boltzmann electrostatics

We describe a new way to calculate the electrostatic properties of macromolecules that goes beyond the classical Poisson-Boltzmann treatment with only a small extra CPU cost. The solvent region is no longer modeled as a homogeneous dielectric media but rather as an assembly of self-orienting interacting dipoles of variable density. The method effectively unifies both the Poisson-centric view and the Langevin Dipole model. The model results in a variable dielectric constant epsilon(r) in the solvent region and also in a variable solvent density rho(r) that depends on the nature of the closest exposed solute atoms. The model was calibrated using small molecules and ions solvation data with only two adjustable parameters, namely the size and dipolar moment of the solvent. Hydrophobicity scales derived from the solvent density profiles agree very well with independently derived hydrophobicity scales, both at the atomic or residue level. Dimerization interfaces in homodimeric proteins or lipid-binding regions in membrane proteins clearly appear as poorly solvated patches on the solute accessible surface. Comparison of the thermally averaged solvent density of this model with the one derived from molecular dynamics simulations shows qualitative agreement on a coarse-grained level. Because this calculation is much more rapid than that from molecular dynamics, applications of a density-profile-based solvation energy to the identification of the true structure among a set of decoys become computationally feasible. Various possible improvements of the model are discussed, as well as extensions of the formalism to treat mixtures of dipolar solvents of different sizes.

Hydration and conformational mechanics of single, end-tethered elastin-like polypeptides

We investigated the effect of temperature, ionic strength, solvent polarity, and type of guest residue on the force-extension behavior of single, end-tethered elastin-like polypeptides (ELPs), using single molecule force spectroscopy (SMFS). ELPs are stimulus-responsive polypeptides that contain repeats of the five amino acids Val-Pro-Gly-Xaa-Gly (VPGXG), where Xaa is a guest residue that can be any amino acid with the exception of proline. We fitted the force-extension data with a freely jointed chain (FJC) model which allowed us to resolve small differences in the effective Kuhn segment length distributions that largely arise from differences in the hydrophobic hydration behavior of ELP. Our results agree qualitatively with predictions from recent molecular dynamics simulations and demonstrate that hydrophobic hydration modulates the molecular elasticity for ELPs. Furthermore, our results show that SMFS, when combined with our approach for data analysis, can be used to study the subtleties of polypeptide-water interactions and thus provides a basis for the study of hydrophobic hydration in intrinsically unstructured biomacromolecules.

Structure-based analysis reveals hydration changes induced by arginine hydrochloride

Arginine hydrochloride has been used to suppress protein aggregation during refolding and in various other applications. We investigated the structure of hen egg-white lysozyme (HEL) and solvent molecules in arginine hydrochloride solution by X-ray crystallography. Neither the backbone nor side-chain structure of HEL was altered by the presence of arginine hydrochloride. In addition, no stably bound arginine molecules were observed. The number of hydration water molecules, however, changed with the arginine hydrochloride concentration. We suggest that arginine hydrochloride suppresses protein aggregation by altering the hydration structure and the transient binding of arginine molecules that could not be observed.