Hydration of proteins and nucleic acids: Advances in experiment and theory. A review

Carbonyl stretching band in amides as Lorentz oscillator. Insights into anharmonicity and local environment in the liquid phase from NIR and MIR spectra of N-methylformamide and di-N,N-methylformamide

We investigated the anharmonicity and intermolecular interactions of N-methylformamide (NMF) and di-N,N-methylformamide (DMF) in the neat liquid phase with particular interest in the amide bands. The vibrational spectra, complex refractive index, and complex electric permittivity were determined in in the mid- (MIR) and near-infrared (NIR) regions (11,500-560 cm-1; 870-17857 nm). Dispersion analysis was based on the Classical Damped Harmonic Oscillator (CDHO) and simultaneous modelling of the real and imaginary components of the spectra. This data delivered insights into the vibrational energy dissipation and self-association in liquid amides. Identification of the MIR and NIR bands was based on anharmonic GVPT2//B3LYP/6-311++G(d,p) calculations. DMF and NMF follow distinct self-association, evidenced in the MIR fingerprint by the two components of the νCO, the analog of the Amide I band. These conclusions are supported by the structural information derived from the NIR spectra. Furthermore, the contribution of overtones and combination bands in the MIR spectra of amides was examined. The conclusions on molecular interactions and structural dynamics of NMF and DMF contribute to a deeper understanding of the effects of changes in the local environment of the amide group.

Anomalously Large Heat Generation of Hydration Water under Microwave Irradiation

Much attention has been paid to the biological effects of microwave irradiation. The hydration water surrounding a biomolecule is crucial in its biological reactions and functions. Therefore, it is important to know the response of hydration water to microwaves to understand their biological effects; however, the scarcity of studies about it often leads to speculations and debates about that effect. In this study, we have made real-time temperature measurements of reverse micellar solutions with their water droplet size from ∼2.3 to ∼9.5 nm using a waveguide system combined with a microwave generator at 2.45 GHz. The heat generated by water in reverse micelles has been observed to depend on their size. It is about 10 times larger than that of liquid water at their small sizes (<∼3.5 nm) and diminishes with further enlarging the size, approaching the water's value at their large sizes (∼10 nm). These results indicate that the heat generation behavior has an interfacial effect; specifically, the hydration water on the surfactant layer produces heat 10 times larger than bulk water. Moreover, the hydration number per surfactant molecule decreases in a core-shell model with increasing the reverse micelle size. These features are also reflected in the heat generation rate. Our findings may offer a new and fundamental perspective for studies on the biological effects of microwave irradiation.

HydraProt: A New Deep Learning Tool for Fast and Accurate Prediction of Water Molecule Positions for Protein Structures

Water molecules are integral to the structural stability of proteins and vital for facilitating molecular interactions. However, accurately predicting their precise position around protein structures remains a significant challenge, making it a vibrant research area. In this paper, we introduce HydraProt (deep Hydration of Proteins), a novel methodology for predicting precise positions of water molecule oxygen atoms around protein structures, leveraging two interconnected deep learning architectures: a 3D U-net and a Multi-Layer Perceptron (MLP). Our approach starts by introducing a coarse voxel-based representation of the protein, which allows for rapid sampling of candidate water positions via the 3D U-net. These water positions are then assessed by embedding the water-protein relationship in the Euclidean space by means of an MLP. Finally, a postprocessing step is applied to further refine the MLP predictions. HydraProt surpasses existing state-of-the-art approaches in terms of precision and recall and has been validated on large data sets of protein structures. Notably, our method offers rapid inference runtime and should constitute the method of choice for protein structure studies and drug discovery applications. Our pretrained models, data, and the source code required to reproduce these results are accessible at https://github.com/azamanos/HydraProt.

On Water Arrangements in Right- and Left-Handed DNA Structures

DNA requires hydration to maintain its structural integrity. Crystallographic analyses have enabled patterns of water arrangements to be visualized. We survey these water motifs in this review, focusing on left- and right-handed duplex and quadruplex DNAs, together with the i-motif. Common patterns of linear spines of water organization in grooves have been identified and are widely prevalent in right-handed duplexes and quadruplexes. By contrast, a left-handed quadruplex has a distinctive wheel of hydration populating the almost completely circular single groove in this structure.

The Advances and Limitations of the Determination and Applications of Water Structure in Molecular Engineering

Water is a key actor of various processes of nature and, therefore, molecular engineering has to take the structural and energetic consequences of hydration into account. While the present review focuses on the target-ligand interactions in drug design, with a focus on biomolecules, these methods and applications can be easily adapted to other fields of the molecular engineering of molecular complexes, including solid hydrates. The review starts with the problems and solutions of the determination of water structures. The experimental approaches and theoretical calculations are summarized, including conceptual classifications. The implementations and applications of water models are featured for the calculation of the binding thermodynamics and computational ligand docking. It is concluded that theoretical approaches not only reproduce or complete experimental water structures, but also provide key information on the contribution of individual water molecules and are indispensable tools in molecular engineering.

Hydration and Structural Adaptations of the Human CYP1A1, CYP1A2, and CYP1B1 Active Sites by Molecular Dynamics Simulations

Cytochromes CYP1A1, CYP1A2, and CYP1B1, the members of the cytochrome P450 family 1, catalyze the metabolism of endogenous compounds, drugs, and non-drug xenobiotics which include substances involved in the process of carcinogenesis, cancer chemoprevention, and therapy. In the present study, the interactions of three selected polymethoxy-trans-stilbenes, analogs of a bioactive polyphenol trans-resveratrol (3,5,4'-trihydroxy-trans-stilbene) with the binding sites of CYP1 isozymes were investigated with molecular dynamics (MD) simulations. The most pronounced structural changes in the CYP1 binding sites were observed in two substrate recognition sites (SRS): SRS2 (helix F) and SRS3 (helix G). MD simulations show that the number and position of water molecules occurring in CYP1 APO and in the structures complexed with ligands are diverse. The presence of water in binding sites results in the formation of water-protein, water-ligand, and bridging ligand-water-protein hydrogen bonds. Analysis of the solvent and substrate channels opening during the MD simulation showed significant differences between cytochromes in relation to the solvent channel and the substrate channels 2c, 2ac, and 2f. The results of this investigation lead to a deeper understanding of the molecular processes that occur in the CYP1 binding sites and may be useful for further molecular studies of CYP1 functions.

Nonthermal acceleration of protein hydration by sub-terahertz irradiation

The collective intermolecular dynamics of protein and water molecules, which overlap in the sub-terahertz (THz) frequency region, are relevant for expressing protein functions but remain largely unknown. This study used dielectric relaxation (DR) measurements to investigate how externally applied sub-THz electromagnetic fields perturb the rapid collective dynamics and influence the considerably slower chemical processes in protein-water systems. We analyzed an aqueous lysozyme solution, whose hydration is not thermally equilibrated. By detecting time-lapse differences in microwave DR, we demonstrated that sub-THz irradiation gradually decreases the dielectric permittivity of the lysozyme solution by reducing the orientational polarization of water molecules. Comprehensive analysis combining THz and nuclear magnetic resonance spectroscopies suggested that the gradual decrease in the dielectric permittivity is not induced by heating but is due to a slow shift toward the hydrophobic hydration structure in lysozyme. Our findings can be used to investigate hydration-mediated protein functions based on sub-THz irradiation.

WaterKit: Thermodynamic Profiling of Protein Hydration Sites

Water desolvation is one of the key components of the free energy of binding of small molecules to their receptors. Thus, understanding the energetic balance of solvation and desolvation resulting from individual water molecules can be crucial when estimating ligand binding, especially when evaluating different molecules and poses as done in High-Throughput Virtual Screening (HTVS). Over the most recent decades, several methods were developed to tackle this problem, ranging from fast approximate methods (usually empirical functions using either discrete atom-atom pairwise interactions or continuum solvent models) to more computationally expensive and accurate ones, mostly based on Molecular Dynamics (MD) simulations, such as Grid Inhomogeneous Solvation Theory (GIST) or Double Decoupling. On one hand, MD-based methods are prohibitive to use in HTVS to estimate the role of waters on the fly for each ligand. On the other hand, fast and approximate methods show an unsatisfactory level of accuracy, with low agreement with results obtained with the more expensive methods. Here we introduce WaterKit, a new grid-based sampling method with explicit water molecules to calculate thermodynamic properties using the GIST method. Our results show that the discrete placement of water molecules is successful in reproducing the position of crystallographic waters with very high accuracy, as well as providing thermodynamic estimates with accuracy comparable to more expensive MD simulations. Unlike these methods, WaterKit can be used to analyze specific regions on the protein surface, (such as the binding site of a receptor), without having to hydrate and simulate the whole receptor structure. The results show the feasibility of a general and fast method to compute thermodynamic properties of water molecules, making it well-suited to be integrated in high-throughput pipelines such as molecular docking.

Knowledge-based prediction of DNA hydration using hydrated dinucleotides as building blocks

Water plays an important role in stabilizing the structure of DNA and mediating its interactions. Here, the hydration of DNA was analyzed in terms of dinucleotide fragments from an ensemble of 2727 nonredundant DNA chains containing 41 853 dinucleotides and 316 265 associated first-shell water molecules. The dinucleotides were classified into categories based on their 16 sequences and the previously determined structural classes known as nucleotide conformers (NtCs). The construction of hydrated dinucleotide building blocks allowed dinucleotide hydration to be calculated as the probability of water density distributions. Peaks in the water densities, known as hydration sites (HSs), uncovered the interplay between base and sugar-phosphate hydration in the context of sequence and structure. To demonstrate the predictive power of hydrated DNA building blocks, they were then used to predict hydration in an independent set of crystal and NMR structures. In ten tested crystal structures, the positions of predicted HSs and experimental waters were in good agreement (more than 40% were within 0.5 Å) and correctly reproduced the known features of DNA hydration, for example the `spine of hydration' in B-DNA. Therefore, it is proposed that hydrated building blocks can be used to predict DNA hydration in structures solved by NMR and cryo-EM, thus providing a guide to the interpretation of experimental data and computer models. The data for the hydrated building blocks and the predictions are available for browsing and visualization at the website https://watlas.datmos.org/watna/.

On the influence of water on THz vibrational spectral features of molecular crystals

The nanoscale structure of molecular assemblies plays a major role in many (μ)-biological mechanisms. Molecular crystals are one of the most simple of these assemblies and are widely used in a variety of applications from pharmaceuticals and agrochemicals, to nutraceuticals and cosmetics. The collective vibrations in such molecular crystals can be probed using terahertz spectroscopy, providing unique characteristic spectral fingerprints. However, the association of the spectral features to the crystal conformation, crystal phase and its environment is a difficult task. We present a combined computational-experimental study on the incorporation of water in lactose molecular crystals, and show how simulations can be used to associate spectral features in the THz region to crystal conformations and phases. Using periodic DFT simulations of lactose molecular crystals, the role of water in the observed lactose THz spectrum is clarified, presenting both direct and indirect contributions. A specific experimental setup is built to allow the controlled heating and corresponding dehydration of the sample, providing the monitoring of the crystal phase transformation dynamics. Besides the observation that lactose phases and phase transformation appear to be more complex than previously thought - including several crystal forms in a single phase and a non-negligible water content in the so-called anhydrous phase - we draw two main conclusions from this study. Firstly, THz modes are spread over more than one molecule and require periodic computation rather than a gas-phase one. Secondly, hydration water does not only play a perturbative role but also participates in the facilitation of the THz vibrations.

Observation of robust energy transfer in the photosynthetic protein allophycocyanin using single-molecule pump-probe spectroscopy

Photosynthetic organisms convert sunlight to electricity with near unity quantum efficiency. Absorbed photoenergy transfers through a network of chromophores positioned within protein scaffolds, which fluctuate due to thermal motion. The resultant variation in the individual energy transfer steps has not yet been measured, and so how the efficiency is robust to this variation has not been determined. Here, we describe single-molecule pump-probe spectroscopy with facile spectral tuning and its application to the ultrafast dynamics of single allophycocyanin, a light-harvesting protein from cyanobacteria. We disentangled the energy transfer and energetic relaxation from nuclear motion using the spectral dependence of the dynamics. We observed an asymmetric distribution of timescales for energy transfer and a slower and more heterogeneous distribution of timescales for energetic relaxation, which was due to the impact of the protein environment. Collectively, these results suggest that energy transfer is robust to protein fluctuations, a prerequisite for efficient light harvesting.

What Does Time-Dependent Fluorescence Shift (TDFS) in Biomembranes (and Proteins) Report on?

The organization of biomolecules and bioassemblies is highly governed by the nature and extent of their interactions with water. These interactions are of high intricacy and a broad range of methods based on various principles have been introduced to characterize them. As these methods view the hydration phenomena differently (e.g., in terms of time and length scales), a detailed insight in each particular technique is to promote the overall understanding of the stunning “hydration world.” In this prospective mini-review we therefore critically examine time-dependent fluorescence shift (TDFS)—an experimental method with a high potential for studying the hydration in the biological systems. We demonstrate that TDFS is very useful especially for phospholipid bilayers for mapping the interfacial region formed by the hydrated lipid headgroups. TDFS, when properly applied, reports on the degree of hydration and mobility of the hydrated phospholipid segments in the close vicinity of the fluorophore embedded in the bilayer. Here, the interpretation of the recorded TDFS parameters are thoroughly discussed, also in the context of the findings obtained by other experimental techniques addressing the hydration phenomena (e.g., molecular dynamics simulations, NMR spectroscopy, scattering techniques, etc.). The differences in the interpretations of TDFS outputs between phospholipid biomembranes and proteins are also addressed. Additionally, prerequisites for the successful TDFS application are presented (i.e., the proper choice of fluorescence dye for TDFS studies, and TDFS instrumentation). Finally, the effects of ions and oxidized phospholipids on the bilayer organization and headgroup packing viewed from TDFS perspective are presented as application examples.

Water molecules at protein-drug interfaces: computational prediction and analysis methods

The fundamental importance of water molecules at drug-protein interfaces is now widely recognised and a significant feature in structure-based drug design. Experimental methods for analysing the role of water in drug binding have many challenges, including the accurate location of bound water molecules in crystal structures, and problems in resolving specific water contributions to binding thermodynamics. Computational analyses of binding site water molecules provide an alternative, and in principle complete, structural and thermodynamic picture, and their use is now commonplace in the pharmaceutical industry. In this review, we describe the computational methodologies that are available and discuss their strengths and weaknesses. Additionally, we provide a critical analysis of the experimental data used to validate the methods, regarding the type and quality of experimental structural data. We also discuss some of the fundamental difficulties of each method and suggest directions for future study.

Nonthermal excitation effects mediated by sub-terahertz radiation on hydrogen exchange in ubiquitin

Water dynamics in the hydration layers of biomolecules play crucial roles in a wide range of biological functions. A hydrated protein contains multiple components of diffusional and vibrational dynamics of water and protein, which may be coupled at ∼0.1-THz frequency (10-ps timescale) at room temperature. However, the microscopic description of biomolecular functions based on various modes of protein-water-coupled motions remains elusive. A novel approach for perturbing the hydration dynamics in the subterahertz frequency range and probing them at the atomic level is therefore warranted. In this study, we investigated the effect of klystron-based, intense 0.1-THz excitation on the slow dynamics of ubiquitin using NMR-based measurements of hydrogen-deuterium exchange. We demonstrated that the subterahertz irradiation accelerated the hydrogen-deuterium exchange of the amides located in the interior of the protein and hydrophobic surfaces while decelerating this exchange in the amides located in the surface loop and short 3₁₀ helix regions. This subterahertz-radiation-induced effect was qualitatively contradictory to the increased-temperature-induced effect. Our results suggest that the heterogeneous water dynamics occurring at the protein-water interface include components that are nonthermally excited by the subterahertz radiation. Such subterahertz-excited components may be linked to the slow function-related dynamics of the protein.

Water spines and networks in G-quadruplex structures

Quadruplex DNAs can fold into a variety of distinct topologies, depending in part on loop types and orientations of individual strands, as shown by high-resolution crystal and NMR structures. Crystal structures also show associated water molecules. We report here on an analysis of the hydration arrangements around selected folded quadruplex DNAs, which has revealed several prominent features that re-occur in related structures. Many of the primary-sphere water molecules are found in the grooves and loop regions of these structures. At least one groove in anti-parallel and hybrid quadruplex structures is long and narrow and contains an extensive spine of linked primary-sphere water molecules. This spine is analogous to but fundamentally distinct from the well-characterized spine observed in the minor groove of A/T-rich duplex DNA, in that every water molecule in the continuous quadruplex spines makes a direct hydrogen bond contact with groove atoms, principally phosphate oxygen atoms lining groove walls and guanine base nitrogen atoms on the groove floor. By contrast, parallel quadruplexes do not have extended grooves, but primary-sphere water molecules still cluster in them and are especially associated with the loops, helping to stabilize loop conformations.

Motif orientation matters: Structural characterization of TEAD1 recognition of genomic DNA

TEAD transcription factors regulate gene expression through interactions with DNA and other proteins. They are crucial for the development of eukaryotic organisms and to control the expression of genes involved mostly in cell proliferation and differentiation; however, their deregulation can lead to tumorigenesis. To study the interactions of TEAD1 with M-CAT motifs and their inverted versions, the K_D of each complex was determined, and H/D exchange, quantitative chemical cross-linking, molecular docking, and smFRET were utilized for structural characterization. ChIP-qPCR was employed to correlate the results with a cell line model. The results obtained showed that although the inverted motif has 10× higher K_D, the same residues were affected by the presence of M-CAT in both orientations. Molecular docking and smFRET revealed that TEAD1 binds the inverted motif rotated 180°. In addition, the inverted motif was proven to be occupied by TEAD1 in Jurkat cells, suggesting that the low-affinity binding sites present in the human genome may possess biological relevance.

Structural alphabets for conformational analysis of nucleic acids available at dnatco.datmos.org

A detailed description of the dnatco.datmos.org web server implementing the universal structural alphabet of nucleic acids is presented. It is capable of processing any mmCIF- or PDB-formatted files containing DNA or RNA molecules; these can either be uploaded by the user or supplied as the wwPDB or PDB-REDO structural database access code. The web server performs an assignment of the nucleic acid conformations and presents the results for the intuitive annotation, validation, modeling and refinement of nucleic acids.

An assessment of water placement algorithms in quantum mechanics/molecular mechanics modeling: the case of rhodopsins' first spectral absorption band maxima

Quantum mechanics/molecular mechanics (QM/MM) models are a widely used tool to obtain detailed insight into the properties and functioning of proteins. The outcome of QM/MM studies heavily depends on the quality of the applied QM/MM model. Prediction and right placement of internal water molecules in protein cavities is one of the critical parts of any QM/MM model construction. Herein, we performed a systematic study of four protein hydration algorithms. We tested these algorithms for their ability to predict X-ray-resolved water molecules for a set of membrane photosensitive rhodopsin proteins, as well as the influence of the applied water placement algorithms on the QM/MM calculated absorption maxima (λ_max) of these proteins. We used 49 rhodopsins and their intermediates with available X-ray structures as the test set. We found that a proper choice of hydration algorithms and setups is needed to predict functionally important water molecules in the chromophore-binding cavity of rhodopsins, such as the water cluster in the N-H region of bacteriorhodopsin or two water molecules in the binding pocket of bovine visual rhodopsin. The QM/MM calculated λ_max of rhodopsins is also quite sensitive to the applied protein hydration protocols. The best methodology allows obtaining an 18.0 nm average value for the absolute deviation of the calculated λ_max from the experimental λ_max. Although the major effect of water molecules on λ_max originates from the water molecules located in the binding pocket, the water molecules outside the binding pocket also affect the calculated λ_max mainly by causing a reorganization of the protein structure. The results reported in this study can be used for the evaluation and further development of hydration methodologies, in general, and rhodopsin QM/MM models, in particular.

Applications of water molecules for analysis of macromolecule properties

Water molecules maintain proteins' structures, functions, stabilities and dynamics. They can occupy certain positions or pass quickly via a protein's interior. Regardless of their behaviour, water molecules can be used for the analysis of proteins' structural features and biochemical properties. Here, we present a list of several software programs that use the information provided by water molecules to: i) analyse protein structures and provide the optimal positions of water molecules for protein hydration, ii) identify high-occupancy water sites in order to analyse ligand binding modes, and iii) detect and describe tunnels and cavities. The analysis of water molecules' distribution and trajectories sheds a light on proteins' interactions with small molecules, on the dynamics of tunnels and cavities, on protein composition and also on the functionality, transportation network and location of functionally relevant residues. Finally, the correct placement of water molecules in protein crystal structures can significantly improve the reliability of molecular dynamics simulations.

Water in protein hydration and ligand recognition

This review describes selected basics of water in biomolecular recognition. We focus on a qualitative understanding of the most important physical aspects, how these change in magnitude between bulk water and protein environment, and how the roles that water plays for proteins arise from them. These roles include mechanical support, thermal coupling, dielectric screening, mass and charge transport, and the competition with a ligand for the occupation of a binding site. The presence or absence of water has ramifications that range from the thermodynamic binding signature of a single ligand up to cellular survival. The large inhomogeneity in water density, polarity and mobility around a solute is hard to assess in experiment. This is a source of many difficulties in the solvation of protein models and computational studies that attempt to elucidate or predict ligand recognition. The influence of water in a protein binding site on the experimental enthalpic and entropic signature of ligand binding is still a point of much debate. The strong water‐water interaction in enthalpic terms is counteracted by a water molecule's high mobility in entropic terms. The complete arrest of a water molecule's mobility sets a limit on the entropic contribution of a water displacement process, while the solvent environment sets limits on ligand reactivity.

Hydration Shell of Antifreeze Proteins: Unveiling the Role of Non-Ice-Binding Surfaces

Antifreeze proteins (AFPs) have the ability to inhibit ice growth by binding to ice nuclei. Their ice-binding mechanism is still unclear, yet the hydration layer is thought to play a fundamental role. Here, we use molecular dynamics simulations to characterize the hydration shell of two AFPs and two non-AFPs. The calculated shell thickness and density of the AFPs do not feature any relevant difference with respect to the non-AFPs. Moreover, the hydration shell density is always higher than the bulk density and, thus, no low-density, ice-like layer is detected at the ice-binding surface (IBS) of AFPs. Instead, we observe local water-density differences in AFPs between the IBS (lower density) and the non-IBS (higher density). The lower solvent density at the ice-binding site can pave the way to the protein binding to ice nuclei, while the higher solvent density at the non-ice-binding surfaces might provide protection against ice growth.

Aquaphotomics-From Innovative Knowledge to Integrative Platform in Science and Technology

Aquaphotomics is a young scientific discipline based on innovative knowledge of water molecular network, which as an intrinsic part of every aqueous system is being shaped by all of its components and the properties of the environment. With a high capacity for hydrogen bonding, water molecules are extremely sensitive to any changes the system undergoes. In highly aqueous systems-especially biological-water is the most abundant molecule. Minute changes in system elements or surroundings affect multitude of water molecules, causing rearrangements of water molecular network. Using light of various frequencies as a probe, the specifics of water structure can be extracted from the water spectrum, indirectly providing information about all the internal and external elements influencing the system. The water spectral pattern hence becomes an integrative descriptor of the system state. Aquaphotomics and the new knowledge of water originated from the field of near infrared spectroscopy. This technique resulted in significant findings about water structure-function relationships in various systems contributing to a better understanding of basic life phenomena. From this foundation, aquaphotomics started integration with other disciplines into systematized science from which a variety of applications ensued. This review will present the basics of this emerging science and its technological potential.

Structural insights into interactions between viral suppressor of RNA silencing protein p19 mutants and small RNAs

Viral suppressors of RNA silencing (VSRSs) are a diverse group of viral proteins that have evolved to disrupt eukaryotic RNA silencing pathways, thereby contributing to viral pathogenicity. The p19 protein is a VSRS that selectively binds to short interfering RNAs (siRNAs) over microRNAs (miRNAs). Mutational analysis has identified single amino acid substitutions that reverse this selectivity through new high-affinity interactions with human miR-122. Herein, we report crystal structures of complexed p19-T111S (2.6 Å), p19-T111H (2.3 Å) and wild-type p19 protein (2.2 Å) from the Carnation Italian ringspot virus with small interfering RNA (siRNA) ligands. Structural comparisons reveal that these mutations do not lead to major changes in p19 architecture, but instead promote subtle rearrangement of residues and solvent molecules along the p19 midline. These observations suggest p19 uses many small interactions to distinguish siRNAs from miRNAs and perturbing these interactions can create p19 variants with novel RNA-recognition properties. DATABASE: Model data are deposited in the PDB database under the accession numbers 6BJG, 6BJH and 6BJV.

The geometry of protein hydration

Based on molecular dynamics simulations of four globular proteins in dilute aqueous solution, with three different water models, we examine several, essentially geometrical, aspects of the protein-water interface that remain controversial or incompletely understood. First, we compare different hydration shell definitions, based on spatial or topological proximity criteria. We find that the best method for constructing monolayer shells with nearly complete coverage is to use a 5 Šwater-carbon cutoff and a 4 Šwater-water cutoff. Using this method, we determine a mean interfacial water area of 11.1 Ų which appears to be a universal property of the protein-water interface. We then analyze the local coordination and packing density of water molecules in the hydration shells and in subsets of the first shell. The mean polar water coordination number in the first shell remains within 1% of the bulk-water value, and it is 5% lower in the nonpolar part of the first shell. The local packing density is obtained from additively weighted Voronoi tessellation, arguably the most physically realistic method for allocating space between protein and water. We find that water in all parts of the first hydration shell, including the nonpolar part, is more densely packed than in the bulk, with a shell-averaged density excess of 6% for all four proteins. We suggest reasons why this value differs from previous experimental and computational results, emphasizing the importance of a realistic placement of the protein-water dividing surface and the distinction between spatial correlation and packing density. The protein-induced perturbation of water coordination and packing density is found to be short-ranged, with an exponential decay "length" of 0.6 shells. We also compute the protein partial volume, analyze its decomposition, and argue against the relevance of electrostriction.

Hydration Structure and Dynamics of Inhibitor-Bound HIV-1 Protease

Water is essential in many biological processes, and the hydration structure plays a critical role in facilitating protein folding, dynamics, and ligand binding. A variety of biophysical spectroscopic techniques have been used to probe the water solvating proteins, often complemented with molecular dynamics (MD) simulations to resolve the spatial and dynamic features of the hydration shell, but comparing relative water structure is challenging. In this study 1 μs MD simulations were performed to identify and characterize hydration sites around HIV-1 protease bound to an inhibitor, darunavir (DRV). The water density, hydration site occupancy, extent and anisotropy of fluctuations, coordinated water molecules, and hydrogen bonds were characterized and compared to the properties of bulk water. The water density of the principal hydration shell was found to be higher than bulk, dependent on the topology and physiochemical identity of the biomolecular surface. The dynamics of water molecules occupying principal hydration sites was highly dependent on the number of water-water interactions and inversely correlated with hydrogen bonds to the protein-inhibitor complex. While many waters were conserved following the symmetry of homodimeric HIV protease, the asymmetry induced by DRV resulted in asymmetric lower-occupancy hydration sites at the concave surface of the active site. Key interactions between water molecules and the protease, that stabilize the protein in the inhibited form, were altered in a drug resistant variant of the protease indicating that modulation of solvent-solute interactions might play a key role in conveying drug resistance. Our analysis provides insights into the interplay between an enzyme inhibitor complex and the hydration shell and has implications in elucidating water structure in a variety of biological processes and applications including ligand binding, inhibitor design, and resistance.

Small Angle Neutron Scattering Studies of R67 Dihydrofolate Reductase, a Tetrameric Protein with Intrinsically Disordered N-Termini

R67 dihydrofolate reductase (DHFR) is a homotetramer with a single active site pore and no sequence or structural homology with chromosomal DHFRs. The R67 enzyme provides resistance to trimethoprim, an active site-directed inhibitor of Escherichia coli DHFR. Sixteen to twenty N-terminal amino acids are intrinsically disordered in the R67 dimer crystal structure. Chymotrypsin cleavage of 16 N-terminal residues results in an active enzyme with a decreased stability. The space sampled by the disordered N-termini of R67 DHFR was investigated using small angle neutron scattering. From a combined analysis using molecular dynamics and the program SASSIE ( http://www.smallangles.net/sassie/SASSIE_HOME.html ), the apoenzyme displays a radius of gyration (Rg) of 21.46 ± 0.50 Å. Addition of glycine betaine, an osmolyte, does not result in folding of the termini as the Rg increases slightly to 22.78 ± 0.87 Å. SASSIE fits of the latter SANS data indicate that the disordered N-termini sample larger regions of space and remain disordered, suggesting they might function as entropic bristles. Pressure perturbation calorimetry also indicated that the volume of R67 DHFR increases upon addition of 10% betaine and decreased at 20% betaine because of the dehydration of the protein. Studies of the hydration of full-length R67 DHFR in the presence of the osmolytes betaine and dimethyl sulfoxide find around 1250 water molecules hydrating the protein. Similar studies with truncated R67 DHFR yield around 400 water molecules hydrating the protein in the presence of betaine. The difference of ∼900 waters indicates the N-termini are well-hydrated.

Structural study of the Fox-1 RRM protein hydration reveals a role for key water molecules in RRM-RNA recognition

The Fox-1 RNA recognition motif (RRM) domain is an important member of the RRM protein family. We report a 1.8 Å X-ray structure of the free Fox-1 containing six distinct monomers. We use this and the nuclear magnetic resonance (NMR) structure of the Fox-1 protein/RNA complex for molecular dynamics (MD) analyses of the structured hydration. The individual monomers of the X-ray structure show diverse hydration patterns, however, MD excellently reproduces the most occupied hydration sites. Simulations of the protein/RNA complex show hydration consistent with the isolated protein complemented by hydration sites specific to the protein/RNA interface. MD predicts intricate hydration sites with water-binding times extending up to hundreds of nanoseconds. We characterize two of them using NMR spectroscopy, RNA binding with switchSENSE and free-energy calculations of mutant proteins. Both hydration sites are experimentally confirmed and their abolishment reduces the binding free-energy. A quantitative agreement between theory and experiment is achieved for the S155A substitution but not for the S122A mutant. The S155 hydration site is evolutionarily conserved within the RRM domains. In conclusion, MD is an effective tool for predicting and interpreting the hydration patterns of protein/RNA complexes. Hydration is not easily detectable in NMR experiments but can affect stability of protein/RNA complexes.

WatAA: Atlas of Protein Hydration. Exploring synergies between data mining and ab initio calculations

Water molecules represent an integral part of proteins and a key determinant of protein structure, dynamics and function. WatAA is a newly developed, web-based atlas of amino-acid hydration in proteins. The atlas provides information about the ordered first hydration shell of the most populated amino-acid conformers in proteins. The data presented in the atlas are drawn from two sources: experimental data and ab initio quantum-mechanics calculations. The experimental part is based on a data-mining study of a large set of high-resolution protein crystal structures. The crystal-derived data include 3D maps of water distribution around amino-acids and probability of occurrence of each of the identified hydration sites. The quantum mechanics calculations validate and extend this primary description by optimizing the water position for each hydration site, by providing hydrogen atom positions and by quantifying the interaction energy that stabilizes the water molecule at the particular hydration site position. The calculations show that the majority of experimentally derived hydration sites are positioned near local energy minima for water, and the calculated interaction energies help to assess the preference of water for the individual hydration sites. We propose that the atlas can be used to validate water placement in electron density maps in crystallographic refinement, to locate water molecules mediating protein-ligand interactions in drug design, and to prepare and evaluate molecular dynamics simulations. WatAA: Atlas of Protein Hydration is freely available without login at .

Unifying the microscopic picture of His-containing turns: from gas phase model peptides to crystallized proteins

The presence in crystallized proteins of a local anchoring between the side chain of a His residue, located in the central position of a γ- or β-turn, and its local main chain environment, is assessed by the comparison of protein structures with relevant isolated model peptides.

Insulin Aggregation at a Dynamic Solid-Liquid-Air Triple Interface

Therapeutic proteins are privileged in drug development because of their exquisite specificity, which is due to their three-dimensional conformation in solution. During their manufacture, storage, and delivery, interactions with material surfaces and air interfaces are known to affect their stability. The growing use of automated devices for handling and injection of therapeutics increases their exposure to protocols involving intermittent wetting, during which the solid-liquid and liquid-air interfaces meet at a triple contact line, which is often dynamic. Using a microfluidic setup, we analyze the effect of a moving triple interface on insulin aggregation in real time over a hydrophobic surface. We combine thioflavin T fluorescence and reflection interference microscopy to concomitantly monitor insulin aggregation and the morphology of the liquid as it dewets the surface. We demonstrate that insulin aggregates in the region of a moving triple interface and not in regions submitted to hydrodynamic shear stress alone, induced by the moving liquid. During dewetting, liquid droplets form on the surface anchored by adsorbed proteins, and the accumulation of amyloid aggregates is observed exclusively as fluorescent rings growing eccentrically around these droplets. The fluorescent rings expand until the entire channel surface sweeped by the triple interface is covered by amyloid fibers. On the basis of our experimental results, we propose a model describing the growth mechanism of insulin amyloid fibers at a moving triple contact line, where proteins adsorbed at a hydrophobic surface are exposed to the liquid-air interface.