Dynamics of Water near a Protein Surface

Osmolyte induced protein stabilization: modulation of associated water dynamics might be a key factor

The mechanism of protein stabilization by osmolytes remains one of the most important and long-standing puzzles. The traditional explanation of osmolyte-induced stability through the preferential exclusion of osmolytes from the protein surface has been seriously challenged by the observations like the concentration-dependent reversal of osmolyte-induced stabilization/destabilization. The more modern explanation of protein stabilization/destabilization by osmolytes considers an indirect effect due to osmolyte-induced distortion of the water structure. It provides a general mechanism, but there are numerous examples of protein-specific effects, i.e., a particular osmolyte might stabilize one protein, but destabilize the other, that could not be rationalized through such an explanation. Herein, we hypothesized that osmolyte-induced modulation of associated water might be a critical factor in controlling protein stability in such a medium. Taking different osmolytes and papain as a protein, we proved that our proposal could explain protein stability in osmolyte media. Stabilizing osmolytes rigidify associated water structures around the protein, whereas destabilizing osmolytes make them flexible. The strong correlation between the stability and the associated water dynamics, and the fact that such dynamics are very much protein specific, established the importance of considering the modulation of associated water structures in explaining the osmolyte-induced stabilization/destabilization of proteins. More interestingly, we took another protein, bromelain, for which a traditionally stabilizing osmolyte, sucrose, acts as a stabilizer at higher concentrations but as a destabilizer at lower concentrations. Our proposal successfully explains such observations, which is probably impossible by any known mechanisms. We believe this report will trigger much research in this area.

Elucidating the Effect of Amine Charge State on Poly(β-amino ester) Degradation Using Permanently Charged Analogs

With synthetic ease and tunable degradation lifetimes, poly(β-amino ester)s (PBAEs) have found use in increasingly diverse applications, from gene therapy to thermosets. Protonatable amines in each repeating unit impart pH-dependent solution behavior and lifetimes, with acidic conditions favoring solubility, yet slowing hydrolysis. Due in part to these interconnected phenomena governing pH-dependent PBAE degradation, predictive degradation models, which would enable user-defined lifetimes, remain elusive. To separate the effects of charge state and solution pH on PBAE degradation, we synthesized poly(β-quaternary ammonium ester)s (PBQAEs), which differ from their parent PBAEs only by an additional methyl group, generating polymers with pH-independent cationic charge. Like PBAEs, PBQAE hydrolysis accelerates with increasing pH, although at a given pH, PBAE degradation outpaces PBQAE degradation. This difference is more pronounced in basic solutions, suggesting that deprotonated PBAE amines accelerate hydrolysis, providing an additional tuning parameter to PBAE lifetime and informing the degradation of PBAEs and other pH-responsive polymers.

Site-Specific Water Dynamics in the First Hydration Layer of an Anti-Freeze Glyco-Protein: A Simulation Study

Antifreeze glycoproteins (AFGPs) inhibit ice recrystallization by a mechanism remaining largely elusive. Dynamics of AFGPs’ hydration water and its involvement in the antifreeze activity, for instance, have not been identified...

Solvent-Dependent Excited-State Evolution of Prodan Dyes

Excited-state character and dynamics of two 6-(dimethylamino)-2-acylnaphthalene dyes (Prodan and Badan-SCH₂CH₂OH) were studied by picosecond time-resolved IR spectroscopy (TRIR) in solvents of different polarity and relaxation times: hexane, CD₃OD, and glycerol-d₈. In all these solvents, near-UV excitation initially produced the same S₁(ππ*) excited state characterized by a broad TRIR signal. A very fast decay (3, ∼100 ps) followed in hexane, whereas conversion to a distinct IR spectrum with a ν(C═O) band downshifted by 76 cm^-1 occurred in polar/H-bonding solvents, slowing down on going from CD₃OD (1, 23 ps) to glycerol-d₈ (5.5, 51, 330 ps). The final relaxed excited state was assigned as planar Me₂N → C═O intramolecular charge transfer S₁(ICT) by comparing experimental and TDDFT-calculated spectra. TRIR conversion kinetics are comparable to those of early stages of multiexponential fluorescence decay and dynamic fluorescence red-shift. This work presents a strong evidence that Prodan-type dyes undergo solvation-driven charge separation in their S₁ state, which is responsible for the dynamic fluorescence Stokes shift observed in polar/H-bonding solvents. The time evolution of the optically prepared S₁(ππ*) state to the S₁(ICT) final state reflects environment relaxation and solvation dynamics. This finding rationalizes the widespread use of Prodan-type dyes as probes of environment dynamics and polarity.

Real-time observation of photoionization-induced water migration dynamics in 4-methylformanilide-water by picosecond time-resolved infrared spectroscopy and ab initio molecular dynamics simulations

A novel time-resolved pump-probe spectroscopic approach that enables to keep high resolution in both the time and energy domain, nanosecond excitation-picosecond ionization-picosecond infrared probe (ns-ps-ps TRIR) spectroscopy, has been applied to the trans-4-methylformanilide-water (4MetFA-W) cluster. Water migration dynamics from the CO to the NH binding site in a peptide linkage triggered by photoionization of 4MetFA-W is directly monitored by the ps time evolution of IR spectra, and the presence of an intermediate state is revealed. The time evolution is analyzed by rate equations based on a four-state model of the migration dynamics. Time constants for the initial to the intermediate and hot product and to the final product are obtained. The acceleration of the dynamics by methyl substitution and the strong contribution of intracluster vibrational energy redistribution in the termination of the solvation dynamics is suggested. This picture is well confirmed by the ab initio on-the-fly molecular dynamics simulations. Vibrational assignments of 4MetFA and 4MetFA-W in the neutral (S₀ and S₁) and ionic (D₀) electronic states measured by ns IR dip and electron-impact IR photodissociation spectroscopy are also discussed prior to the results of time-resolved spectroscopy.

Structural and dynamical heterogeneity of water trapped inside Na+-pumping KR2 rhodopsin in the dark state

Photoisomerization in the retinal leads to a channel opening in rhodopsins that triggers translocation or pumping of ions/protons. Crystal structures of rhodopsins contain several structurally conserved water molecules. It has been suggested that water plays an active role in facilitating the ion pumping/translocation process by acting as a lubricant in these systems. In this paper, we systematically investigate the localization, structure, dynamics, and energetics of the water molecules along the channel for the resting/dark state of KR2 rhodopsin. By employing several microseconds long atomistic molecular dynamics simulation of this trans-membrane protein system, we demonstrate the presence of five distinct water containing pockets/cavities separated by gateways controlled by protein side-chains. There exists a strong hydrogen bonded network involving these buried water molecules and functionally important key residues. We present evidence of significant structural and dynamical heterogeneity in the water molecules present in these cavities, with very rare exchange between them. The exchange time scale of such buried water with the bulk has an extremely wide range, from tens of nanoseconds to >1.5 µs. The translational and rotational dynamics of buried water are found to be strongly dependent on the protein cavity size and local interactions with a classic signature of trapped diffusion and rotational anisotropy.

Ultrafast Fluorescence Spectroscopy via Upconversion and Its Applications in Biophysics

In this review, the experimental set-up and functional characteristics of single-wavelength and broad-band femtosecond upconversion spectrophotofluorometers developed in our laboratory are described. We discuss applications of this technique to biophysical problems, such as ultrafast fluorescence quenching and solvation dynamics of tryptophan, peptides, proteins, reduced nicotinamide adenine dinucleotide (NADH), and nucleic acids. In the tryptophan dynamics field, especially for proteins, two types of solvation dynamics on different time scales have been well explored: ~1 ps for bulk water, and tens of picoseconds for “biological water”, a term that combines effects of water and macromolecule dynamics. In addition, some proteins also show quasi-static self-quenching (QSSQ) phenomena. Interestingly, in our more recent work, we also find that similar mixtures of quenching and solvation dynamics occur for the metabolic cofactor NADH. In this review, we add a brief overview of the emerging development of fluorescent RNA aptamers and their potential application to live cell imaging, while noting how ultrafast measurement may speed their optimization.

Hydrogen-Bond Dynamics and Energetics of Biological Water

Water molecules in the immediate vicinity of biomacromolecules and biomimetic organized assemblies often exhibit a markedly distinct behavior from that of their bulk counterparts. The overall sluggish behavior of biological water substantially affects the stability and integrity of biomolecules, as well as the successful execution of various crucial water-mediated biochemical phenomena. In this Minireview, insights are provided into the features of truncated hydrogen-bond networks that grant biological water its unique characteristics. In particular, experimental results and theoretical investigations, based on chemical kinetics, are presented that have shed light on the dynamics and energetics governing such characteristics. It is emphasized how such details help us to understand the energetics of biological water, an aspect relatively less explored than its dynamics. For instance, when biological water at hydrophilic or charged protein surfaces was explored, the free energy of H-bond breakage was found to be of the order of 0.4 kcal mol^-1 higher than that of bulk water.

Dynamic behaviors of interfacial water on the self-assembly monolayer (SAM) heterogeneous surface

Dynamic behaviors of water molecules near the surface with mixed hydrophobic and hydrophilic areas are studied by molecular dynamics simulation. More specifically, the diffusion coefficient and hydrogen bond lifetime of interfacial water on the self-assembly monolayer composed of hydrophobic and hydrophilic groups and their dependence on the mixing ratio are studied. The diffusion dramatically slows down, and the hydrogen bond lifetime considerably increases when a few hydrophilic groups are added to the hydrophobic surface. When the percentage of hydrophilic groups increases to 25%, the behavior of interfacial water is similar to the case of the pure hydrophilic surface. The sensitivity to the hydrophilic group can be attributed to the fact that the grafted hydrophilic groups can not only retard the directly bound water molecules but also affect indirectly bound water by stabilizing hydrogen bonds among interfacial water molecules.

Femtosecond Fluorescence Spectra of NADH in Solution: Ultrafast Solvation Dynamics

The ultrafast solvation dynamics of reduced nicotinamide adenine dinucleotide (NADH) free in solution has been investigated, using both a femtosecond upconversion spectrophotofluorometer and a picosecond time-correlated single-photon counting (TCSPC) apparatus. The familiar time constant of solvent relaxation originating in "bulk water" was found to be ∼1.4 ps, revealing ultrafast solvent reorientation upon excitation. We also found a slower spectral relaxation process with an apparent time of 27 ps, suggesting there could either be dissociable "biological water" hydration sites on the surface of NADH or internal dielectric rearrangements of the flexible solvated molecule on that timescale. In contrast, the femtosecond fluorescence anisotropy measurement revealed that rotational diffusion happened on two different timescales (3.6 ps (local) and 141 ps (tumbling)); thus, any dielectric rearrangement scenario for the 27 ps relaxation must occur without significant chromophore oscillator rotation. The coexistence of quasi-static self quenching (QSSQ) with the slower relaxation is also discussed.

Surface Topography Effects of Globular Biomolecules on Hydration Water

Hydration water serves as a microscopic manifestation of structural stability and functions of biomolecules. To develop bio-nanomaterials in applications, it is important to study how the surface topography and heterogeneity of biomolecules result in their diversity of the hydration dynamics and energetics. We here performed molecular dynamics simulations combined with the steered molecular dynamics and umbrella sampling to investigate the dynamics and escape process associated with the free energy change of water molecules close to a globular biomolecule, i.e., hemoglobin (Hb) and G-quadruplex DNA (GDNA). The residence time, power of long-time tail, and dipole relaxation time were found to display drastic changes within the averaged hydration shell of 3.0-5.0 Å. Compared with bulk water, in the inner hydration shell, the water dipole moment displays a slower relaxation process and is more oriented toward GDNA than toward Hb, forming a hedgehog-like structure when it surrounds GDNA. In particular, a spine water structure is observed in the GDNA narrow groove. The water isotope effect not only prolongs the dynamic time scales of libration motion in the inner hydration shell and the dipole relaxation processes in the bulk but also strengthens the DNA spine water structure. The potential of the mean force profile reflects the integrity of the hydration shell structure and enables us to obtain detailed insights into the structures formed by water, such as the caged H-bond network and the edge bridge structures; it also reveals that local hydration shell free energy (LHSFE) depends on H-bond rupture processes and ranges from 0.2 to 4.2 kcal/mol. Our results demonstrate that the surface topography of a biomolecule influences the integrity of the hydration shell structure and LHSFE. Our studies are able to identify various further applications in the areas of microfluid devices and nano-dewetting on bioinspired surfaces.

Fundamental limitations of the time-dependent Stokes shift for investigating protein hydration dynamics

In the TDSS measured in protein systems, large protein contributions fully obscure hydration dynamics.

Changes in protein hydration dynamics by encapsulation or crowding of ubiquitin: strong correlation between time-dependent Stokes shift and intermolecular nuclear Overhauser effect

The local changes in protein hydration dynamics upon encapsulation of the protein or macromolecular crowding are essential to understand protein function in cellular environments. We were able to obtain a spatially-resolved picture of the influence of confinement and crowding on the hydration dynamics of the protein ubiquitin by analyzing the time-dependent Stokes shift (TDSS), as well as the intermolecular Nuclear Overhauser Effect (NOE) at different sites of the protein by large-scale computer simulation of single and multiple proteins in water and confined in reverse micelles. Besides high advanced space resolved information on hydration dynamics we found a strong correlation of the change in NOE upon crowding or encapsulation and the change in the integral TDSS relaxation times in all investigated systems relative to the signals in a diluted protein solution.

Changes in local protein hydration dynamics caused by encapsulation or crowding are reflected in the TDSS and the intermolecular NOE alike.

Distinguishing dynamical features of water inside protein hydration layer: Distribution reveals what is hidden behind the average

Since the pioneering works of Pethig, Grant, and Wüthrich on a protein hydration layer, many studies have been devoted to find out if there are any "general and universal" characteristic features that can distinguish water molecules inside the protein hydration layer from bulk. Given that the surface itself varies from protein to protein, and that each surface facing the water is heterogeneous, search for universal features has been elusive. Here, we perform an atomistic molecular dynamics simulation in order to propose and demonstrate that such defining characteristics can emerge if we look not at average properties but the distribution of relaxation times. We present results of calculations of distributions of residence times and rotational relaxation times for four different protein-water systems and compare them with the same quantities in the bulk. The distributions in the hydration layer are unusually broad and log-normal in nature due to the simultaneous presence of peptide backbones that form weak hydrogen bonds, hydrophobic amino acid side chains that form no hydrogen bond, and charged polar groups that form a strong hydrogen bond with the surrounding water molecules. The broad distribution is responsible for the non-exponential dielectric response and also agrees with large specific heat of the hydration water. Our calculations reveal that while the average time constant is just about 2-3 times larger than that of bulk water, it provides a poor representation of the real behaviour. In particular, the average leads to the erroneous conclusion that water in the hydration layer is bulk-like. However, the observed and calculated lower value of static dielectric constant of hydration layer remained difficult to reconcile with the broad distribution observed in dynamical properties. We offer a plausible explanation of these unique properties.

Origin of diverse time scales in the protein hydration layer solvation dynamics: A simulation study

In order to inquire the microscopic origin of observed multiple time scales in solvation dynamics, we carry out several computer experiments. We perform atomistic molecular dynamics simulations on three protein-water systems, namely, lysozyme, myoglobin, and sweet protein monellin. In these experiments, we mutate the charges of the neighbouring amino acid side chains of certain natural probes (tryptophan) and also freeze the side chain motions. In order to distinguish between different contributions, we decompose the total solvation energy response in terms of various components present in the system. This allows us to capture the interplay among different self- and cross-energy correlation terms. Freezing the protein motions removes the slowest component that results from side chain fluctuations, but a part of slowness remains. This leads to the conclusion that the slow component approximately in the 20-80 ps range arises from slow water molecules present in the hydration layer. While the more than 100 ps component has multiple origins, namely, adjacent charges in amino acid side chains, hydrogen bonded water molecules and a dynamically coupled motion between side chain and water. In addition, the charges enforce a structural ordering of nearby water molecules and helps to form a local long-lived hydrogen bonded network. Further separation of the spatial and temporal responses in solvation dynamics reveals different roles of hydration and bulk water. We find that the hydration layer water molecules are largely responsible for the slow component, whereas the initial ultrafast decay arises predominantly (approximately 80%) due to the bulk. This agrees with earlier theoretical observations. We also attempt to rationalise our results with the help of a molecular hydrodynamic theory that was developed using classical time dependent density functional theory in a semi-quantitative manner.

Anomalous water dynamics at surfaces and interfaces: synergistic effects of confinement and surface interactions

In nature, water is often found in contact with surfaces that are extended on the scale of molecule size but small on a macroscopic scale. Examples include lipid bilayers and reverse micelles as well as biomolecules like proteins, DNA and zeolites, to name a few. While the presence of surfaces and interfaces interrupts the continuous hydrogen bond network of liquid water, confinement on a mesoscopic scale introduces new features. Even when extended on a molecular scale, natural and biological surfaces often have features (like charge, hydrophobicity) that vary on the scale of the molecular diameter of water. As a result, many new and exotic features, which are not seen in the bulk, appear in the dynamics of water close to the surface. These different behaviors bear the signature of both water-surface interactions and of confinement. In other words, the altered properties are the result of the synergistic effects of surface-water interactions and confinement. Ultrafast spectroscopy, theoretical modeling and computer simulations together form powerful synergistic approaches towards an understanding of the properties of confined water in such systems as nanocavities, reverse micelles (RMs), water inside and outside biomolecules like proteins and DNA, and also between two hydrophobic walls. We shall review the experimental results and place them in the context of theory and simulations. For water confined within RMs, we discuss the possible interference effects propagating from opposite surfaces. Similar interference is found to give rise to an effective attractive force between two hydrophobic surfaces immersed and kept fixed at a separation of d, with the force showing an exponential dependence on this distance. For protein and DNA hydration, we shall examine a multitude of timescales that arise from frustration effects due to the inherent heterogeneity of these surfaces. We pay particular attention to the role of orientational correlations and modification of the same due to interaction with the surfaces.

Interfacial water screens the protein-induced transmembrane voltage

Transmembrane proteins are crucial in cellular traffic, signal transduction, and energy storage in a form of transmembrane voltage. These proteins are stabilized by hydrophobic and hydrophilic interactions, in which cytoplasmic and exoplasmic water plays a special role. Water structural ordering generates the dipole potential that typically overcompensates for an intrinsic membrane-protein potential gradient, and thus it modifies and sustains an overall cellular electrostatics. Although the transmembrane voltage has been extensively studied, the dipole potential has attracted very little attention. Here, by using molecular dynamics, we examined water electrostatic response to the transmembrane charge, field, and potential asymmetry introduced by the presence of four integral membrane proteins: typical of inner (α-helix) and outer membrane (β-barrel). In all cases, the protein presence introduces electrostatic directionality in the transmembrane dipole field and voltage. In particular, water generates a deep potential sink if strongly polar residues are densely packed on one side of bilayer, as frequently occurs in a selectivity filter of the K(+) channel. We also found that protein secondary structure is less important than the polar residue distribution along the protein channel. Our findings are relevant for understanding the driving force behind biomembrane conductivity: the ability of biological water to electrostatically screen the transmembrane voltage.

Single water solvation dynamics in the 4-aminobenzonitrile–water cluster cation revealed by picosecond time-resolved infrared spectroscopy

The excess energy of photoionization can control the time scale of single water solvent orientation dynamics from picoseconds to infinitely long trapping in a local minimum.

Heterogeneous Hydration of p53/MDM2 Complex

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

Phase diagrams and kinetics of phase transitions in protein solutions

The phase behavior of proteins is of interest for fundamental and practical reasons. The nucleation of new phases is one of the last major unresolved problems of nature. The formation of protein condensed phases (crystals, polymers, and other solid aggregates, as well as dense liquids and gels) underlies pathological conditions, plays a crucial role in the biological function of the respective protein, or is an essential part of laboratory and industrial processes. In this review, we focus on phase transitions of proteins in their properly folded state. We first summarize the recently acquired understanding of physical processes underlying the phase diagrams of the protein solutions and the thermodynamics of protein phase transitions. Then we review recent findings on the kinetics of nucleation of dense liquid droplets and crystals. We explore the transition from nucleation to spinodal decomposition for liquid-liquid separation and introduce the new concept of solution-to-crystal spinodal. We review the two-step mechanism of protein crystal nucleation, in which mesoscopic metastable protein clusters serve as precursors to the ordered crystal nuclei. The concepts and mechanisms reviewed here provide powerful tools for control of the nucleation process by varying the solution thermodynamic parameters.

Vibrational spectra of proximal water in a thermo-sensitive polymer undergoing conformational transition across the lower critical solution temperature

The vibrational spectrum of water near a thermo-sensitive polymer poly(N-isopropylacrylamide) (PNIPAM) undergoing conformational transition through the lower critical solution temperature (LCST) is calculated using molecular dynamics simulations. The characteristic structural features observed at the atomic scale for these proximal water molecules in a solvated polymer chain while undergoing the conformational transition are strongly correlated to their vibrational densities of states. Comparison of the vibrational spectrum below LCST for the proximal water with the vibrational spectrum obtained for bulk water reveals a significant fraction of the hydrogen bonding between the proximal water molecules and the polymer side groups. Hydrogen-bonded bridges of water molecules are formed between two adjacent and alternate monomers. This network of hydrogen bonding results in formation of locally ordered water molecules at temperatures below the LCST. Analysis of the simulation trajectories confirms the presence of a quasi-stable solvation structure near the PNIPAM. The calculated vibrational spectra for proximal water above the LCST suggest significantly reduced hydrogen bonding with the polymer and indicate a reduction in the structural stability of proximal water around a collapsed polymer chain. Systematic trends in the observed peak intensities and frequency shifts at the low- and high-frequency ends of the spectrum can be correlated with the structural and dynamical changes of water molecules below and above the LCST transition, respectively, for various polymer chain lengths. The simulations reveal that, compared to bulk water, the libration bands are blue shifted and OH stretch bands red shifted for water in proximity to PNIPAM with 30 monomer units below the LCST. The simulations suggest that vibrational spectra can be used as a predictive tool for quantifying atomic-scale structural transitions in solvation of thermo-sensitive polymers such as PNIPAM.

Entropy-driven binding of opioid peptides induces a large domain motion in human dipeptidyl peptidase III

Opioid peptides are involved in various essential physiological processes, most notably nociception. Dipeptidyl peptidase III (DPP III) is one of the most important enkephalin-degrading enzymes associated with the mammalian pain modulatory system. Here we describe the X-ray structures of human DPP III and its complex with the opioid peptide tynorphin, which rationalize the enzyme's substrate specificity and reveal an exceptionally large domain motion upon ligand binding. Microcalorimetric analyses point at an entropy-dominated process, with the release of water molecules from the binding cleft ("entropy reservoir") as the major thermodynamic driving force. Our results provide the basis for the design of specific inhibitors that enable the elucidation of the exact role of DPP III and the exploration of its potential as a target of pain intervention strategies.

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.

Effect of the ordered interfacial water layer in protein complex formation: A nonlocal electrostatic approach

Using a nonlocal electrostatic approach that incorporates the short-range structure of the contacting media, we evaluated the electrostatic contribution to the energy of the complex formation of two model proteins. In this study, we have demonstrated that the existence of an ordered interfacial water layer at the protein-solvent interface reduces the charging energy of the proteins in the aqueous solvent, and consequently increases the electrostatic contribution to the protein binding (change in free energy upon the complex formation of two proteins). This is in contrast with the finding of the continuum electrostatic model, which suggests that electrostatic interactions are not strong enough to compensate for the unfavorable desolvation effects.

The implications of gene heterozygosity for protein folding and protein turnover

The offspring of closely related parents often suffer from inbreeding depression, sometimes resulting in a slower growth rate for inbred offspring relative to non-inbred offspring. Previous research has shown that some of the slower growth rate of inbred organisms can be attributed to the inbred organisms' increased levels of protein turnover. This paper attempts to show that the higher levels of protein turnover among inbred organisms can be attributed to accumulations of misfolded and aggregated proteins that require degradation by the inbred organisms' protein quality control systems. The accumulation of misfolded and aggregated proteins within inbred organisms are the result of more negative free energies of folding for proteins encoded at homozygous gene loci and higher concentrations of potentially aggregating non-native protein species within the cell. The theory presented here makes several quantitative predictions that suggest a connection between protein misfolding/aggregation and polyploidy that can be tested by future research.

Relaxation dynamics of nucleosomal DNA

Recent experimental and theoretical evidence demonstrates that proteins and water in the hydration layer can follow complex stretched exponential or power law relaxation dynamics. Here, we report on a 50 ns all atom molecular dynamics (MD) simulation of the yeast nucleosome, where the interactions between DNA, histones, surrounding water and ions are explicitly included. DNA interacts with the histone core in 14 locations, approximately every 10.4 base pairs. We demonstrate that all sites of interaction exhibit anomalously slow power law relaxation, extending up to 10 ns, while fast exponential relaxation dynamics of hundreds of picoseconds applies to DNA regions outside these locations. The appearance of 1/f(alpha) noise or pink noise in DNA dynamics is ubiquitous. For histone-bound nucleotide dynamics alpha --> 1 and is a signature of complexity of the protein-DNA interactions. For control purposes two additional DNA simulations free of protein are conducted. Both utilize the same sequence of DNA, as found the in the nucleosome. In one simulation the initial conformation of the double helix is a straight B-form. In the other, the initial conformation is super helical. Neither of these simulations exhibits the variation of alpha as a function of position, the measure of power law for dynamical behavior, which we observe in the nucleosome simulation. The unique correspondence (high alpha to DNA-histone interaction sites, low alpha to free DNA sites), suggests that alpha may be an important and new quantification of protein-DNA interactions for future experiments.

Hydration profiles of amyloidogenic molecular structures

Hydration shells of normal proteins display regions of highly structured water as well as patches of less structured bulk-like water. Recent studies suggest that isomers with larger surface densities of patches of bulk-like water have an increased propensity to aggregate. These aggregates are toxic to the cellular environment. Hence, the early detection of these toxic deposits is of paramount medical importance. We show that various morphological states of association of such isomers can be differentiated from the normal protein background based on the characteristic partition between bulk, caged, and surface hydration water and the magnetic resonance (MR) signals of this water. We derive simple mathematical equations relating the compartmentalization of water to the local hydration fraction and the packing density of the newly formed molecular assemblies. Then, we employ these equations to predict the MR response of water constrained by protein aggregation. Our results indicate that single units and compact aggregates that contain no water between constituents induce a shift of the MR signal from normal protein background to values in the hyperintensity domain (bright spots), corresponding to bulk water. In contrast, large plaques that cage significant amounts of water between constituents are likely to generate MR responses in the hypointensity domain (dark spots), typical for strongly correlated water. The implication of these results is that amyloids can display both dark and bright spots when compared to the normal gray background tissue on MR images. In addition, our findings predict that the bright spots are more likely to correspond to amyloids in their early stage of development. The results help explain the MR contrast patterns of amyloids and suggest a new approach for identifying unusual protein aggregation related to disease.

Determination of the transition-state entropy for aggregation suggests how the growth of sickle cell hemoglobin polymers can be slowed

Sickle cell anemia is associated with the mutant hemoglobin HbS, which forms polymers in red blood cells of patients. The growth rate of the polymers is several micrometers per second, ensuring that a polymer fiber reaches the walls of an erythrocyte (which has a 7-microm diameter) within a few seconds after its nucleation. To understand the factors that determine this unusually fast rate, we analyze data on the growth rate of the polymer fibers. We show that the fiber growth follows a first-order Kramers-type kinetics model. The entropy of the transition state for incorporation into a fiber is 95 J mol(-1) K(-1), very close to the known entropy of polymerization. This agrees with a recent theoretical estimate for the hydrophobic interaction and suggests that the gain of entropy in the transition state is due to the release of the last layer of water molecules structured around contact sites on the surface of the HbS molecules. As a result of this entropy gain, the free-energy barrier for incorporation of HbS molecules into a fiber is negligible and fiber growth is unprecedentedly fast. This finding suggests that fiber growth can be slowed by components of the red cell cytosol, native or intentionally introduced, which restructure the hydration layer around the HbS molecules and thus lower the transition state entropy for incorporation of an incoming molecule into the growing fiber.

Computational methods for biomolecular electrostatics

An understanding of intermolecular interactions is essential for insight into how cells develop, operate, communicate, and control their activities. Such interactions include several components: contributions from linear, angular, and torsional forces in covalent bonds, van der waals forces, as well as electrostatics. Among the various components of molecular interactions, electrostatics are of special importance because of their long range and their influence on polar or charged molecules, including water, aqueous ions, and amino or nucleic acids, which are some of the primary components of living systems. Electrostatics, therefore, play important roles in determining the structure, motion, and function of a wide range of biological molecules. This chapter presents a brief overview of electrostatic interactions in cellular systems, with a particular focus on how computational tools can be used to investigate these types of interactions.

Dynamic Stokes shift in green fluorescent protein variants

Solvent reorganization around the excited state of a chromophore leads to an emission shift to longer wavelengths during the excited-state lifetime. This solvation response is absent in wild-type green fluorescent protein, and this has been attributed to rigidity in the chromophore's environment necessary to exclude nonradiative transitions to the ground state. The fluorescent protein mPlum was developed via directed evolution by selection for red emission, and we use time-resolved fluorescence to study the dynamic Stokes shift through its evolutionary history. The far-red emission of mPlum is attributed to a picosecond solvation response that is observed at all temperatures above the glass transition. This time-dependent shift in emission is not observed in its evolutionary ancestors, suggesting that selective pressure has produced a chromophore environment that allows solvent reorganization. The evolutionary pathway and structures of related fluorescent proteins suggest the role of a single residue in close proximity to the chromophore as the primary cause of the solvation response.

Evaluation of the influence of the internal aqueous solvent structure on electrostatic interactions at the protein-solvent interface by nonlocal continuum electrostatic approach

The dielectric properties of the polar solvent on the protein-solvent interface at small intercharge distances are still poorly explored. To deconvolute this problem and to evaluate the pair-wise electrostatic interaction (PEI) energies of the point charges located at the protein-solvent interface we used a nonlocal (NL) electrostatic approach along with a static NL dielectric response function of water. The influence of the aqueous solvent microstructure (determined by a strong nonelectrostatic correlation effect between water dipoles within the orientational Debye polarization mode) on electrostatic interactions at the interface was studied in our work. It was shown that the PEI energies can be significantly higher than the energies evaluated by the classical (local) consideration, treating water molecules as belonging to the bulk solvent with a high dielectric constant. Our analysis points to the existence of a rather extended, effective low-dielectric interfacial water shell on the protein surface. The main dielectric properties of this shell (effective thickness together with distance- and orientation-dependent dielectric permittivity function) were evaluated. The dramatic role of this shell was demonstrated when estimating the protein association rate constants.

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

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

Hydration dynamics and time scales of coupled water-protein fluctuations

We report experimental and theoretical studies on water and protein dynamics following photoexcitation of apomyoglobin. Using site-directed mutation and with femtosecond resolution, we experimentally observed relaxation dynamics with a biphasic distribution of time scales, 5 and 87 ps, around the site Trp7. Theoretical studies using both linear response and direct nonequilibrium molecular dynamics (MD) calculations reproduced the biphasic behavior. Further constrained MD simulations with either frozen protein or frozen water revealed the molecular mechanism of slow hydration processes and elucidated the role of protein fluctuations. Observation of slow water dynamics in MD simulations requires protein flexibility, regardless of whether the slow Stokes shift component results from the water or protein contribution. The initial dynamics in a few picoseconds represents fast local motions such as reorientations and translations of hydrating water molecules, followed by slow relaxation involving strongly coupled water-protein motions. We observed a transition from one isomeric protein configuration to another after 10 ns during our 30 ns ground-state simulation. For one isomer, the surface hydration energy dominates the slow component of the total relaxation energy. For the other isomer, the slow component is dominated by protein interactions with the chromophore. In both cases, coupled water-protein motion is shown to be necessary for observation of the slow dynamics. Such biologically important water-protein motions occur on tens of picoseconds. One significant discrepancy exists between theory and experiment, the large inertial relaxation predicted by simulations but clearly absent in experiment. Further improvements required in the theoretical model are discussed.