A molecular dynamics study of solvent behavior around a protein

The solvent structure and behavior around a protein were examined by analyzing a trajectory of molecular dynamics simulation of thetrp-holorepressor in a periodic box of water. The calculated selfdiffusion coefficient indicated that the solvent within 10 Å of the protein had lower mobility. Examination of the solvent diffusion around different atoms of different kinds of residues showed no general tendency. Thisfact suggested that the solvent mobility is not influenced significantly bythe kind of the atom or residue they solvated. Distribution analysis aroundthe protein revealed two peaks of water oxygen: a sharp one at 2.8 Å around polar and charged atoms and a broad one at ～3.4 Å aroundapolar atoms. The former was stabilized by water–protein hydrogen bonds, and the latter was stabilized by water-lwater hydrogen bonds, suggesting the existence of a hydrophobic shell. An analysis of protein atom–water radial distribution functions confirmed these shell structures around polar or charged atoms and apolar ones. © 1993 Wiley-Liss, Inc.     

Water structure around dipeptides in aqueous solutions

The bulk water structure around small peptide fragments-glycyl-L: -alanine, glycyl-L: -proline and L: -alanyl-L: -proline-has been determined by a combination of neutron diffraction with isotopic substitution and empirical potential structural refinement techniques. The addition of each of the dipeptides to water gives rise to decreased water-water coordination in the surrounding water solvent. Additionally, both the O(w)-O(w) radial distribution functions and the water-water spatial density functions in all of the solutions indicate an electrostrictive effect in the second water coordination shell of the bulk water network. This effect is not observed in similar experiments on the amino acid L: -proline alone in solution, which is one component of two of the peptides measured here.     

Polyproline II helix is the preferred conformation for unfolded polyalanine in water

Does aqueous solvent discriminate among peptide conformers? To address this question, we computed the solvation free energy of a blocked, 12‐residue polyalanyl‐peptide in explicit water and analyzed its solvent structure. The peptide was modeled in each of 4 conformers: α‐helix, antiparallel β‐strand, parallel β‐strand, and polyproline II helix (P_II). Monte Carlo simulations in the canonical ensemble were performed at 300 K using the CHARMM 22 forcefield with TIP3P water. The simulations indicate that the solvation free energy of P_II is favored over that of other conformers for reasons that defy conventional explanation. Specifically, in these 4 conformers, an almost perfect correlation is found between a residue's solvent‐accessible surface area and the volume of its first solvent shell, but neither quantity is correlated with the observed differences in solvation free energy. Instead, solvation free energy tracks with the interaction energy between the peptide and its first‐shell water. An additional, previously unrecognized contribution involves the conformation‐dependent perturbation of first‐shell solvent organization. Unlike P_II, β‐strands induce formation of entropically disfavored peptide:water bridges that order vicinal water in a manner reminiscent of the hydrophobic effect. The use of explicit water allows us to capture and characterize these dynamic water bridges that form and dissolve during our simulations. Proteins 2004. © 2004 Wiley‐Liss, Inc.     

Myosin-induced volume increase of the hyper-mobile water surrounding actin filaments

Microwave dielectric spectroscopy can measure the rotational mobility of water molecules that hydrate proteins and the hydration-shell volume. Using this technique, we have recently shown that apart from typical hydrating water molecules with lowered mobility there are other water molecules around the actin filaments (F-actin) which have a much higher mobility than that of bulk water [Biophys. J. 85 (2003) 3154]. We report here that the volume of this water component (hyper-mobile water) markedly increases without significant change of the volume of the ordinary hydration shell when the myosin motor-domain (S1, myosin subfragment-1) binds to F-actin. No hyper-mobile component was found in the hydration shell of S1 itself. The present results strongly suggest that the solvent space around S1 bound to F-actin is diffusionally asymmetric, which supports our model of force generation by actomyosin proposed previously [op. cit.].     

Surface-directed structure formation of β-lactoglobulin inside droplets

The morphology of β-lactoglobulin structures inside droplets was studied during aggregation and gelation using confocal laser scanning microscopy (CLSM) equipped with a temperature stage and transmission electron microscopy (TEM). The results showed that there is a strong driving force for the protein to move to the interface between oil and water in the droplet, and the β-lactoglobulin formed a dense shell around the droplet built up from the inside of the droplets. Less protein was found inside the droplets. The longer the β-lactoglobulin was allowed to aggregate prior to gel formation, the larger the part of the protein went to the interface, resulting in a thicker shell and very little material being left inside the droplets. The droplets were easily deformed because no network stabilizes them. When 0.5% emulsifier, polyglycerol polyresinoleat (PGPR), was added to the oil phase, the β-lactoglobulin was situated both inside the droplets and at the interface between the droplets and the oil phase; when 2% PGPR was added, the β-lactoglobulin structure was concentrated to the inside of the droplets. The possibility to use the different morphological structures of β-lactoglobulin in droplets to control the diffusion rate through a β-lactoglobulin network was evaluated by fluorescence recovery after photobleaching (FRAP). The results show differences in the diffusion rate due to heterogeneities in the structure: the diffusion of a large water-soluble molecule, FITC-dextran, in a dense particulate gel was 1/4 of the diffusion rate in a more open particulate β-lactoglobulin gel in which the diffusion rate was similar to that in pure water.     

Analysis of peptide rotational diffusion by homonuclear NMR

The analysis of the rotational diffusion of a molecule using homonuclear NMR is investigated. The homonuclear longitudinal and transverse cross-relaxation rates, which can be quantitatively measured using off-Resonance Rotating frame nuclear Overhauser Effect Spectroscopy (ROESY), are used to build a distribution, which exhibits a solid-state-like pattern characteristic of the diffusion tensor. The distributions of the antimicrobial peptide ranalexin in water and in 30% of trifluoracetic acid (TFE) are compared, and the peptide rotational diffusion is shown to be more isotropic in water than in 30% TFE. This difference is further supported by the analysis of NMR ranalexin conformers in 30% TFE, and by the analysis of a molecular dynamics simulation of peptide in water.     

The molecular basis of the temperature- and pH-induced conformational transitions in elastin-based peptides

Elastin undergoes an inverse temperature transition and collapses at high temperatures in both simulation and experiment. We investigated a pH-dependent modification of this transition by simulating a glutamic acid (Glu)-substituted elastin at varying pHs and temperatures. The Glu-substituted peptide collapsed at higher temperature than the unsubstituted elastin when Glu was charged. The charge effects could be reversed by neutralization of the Glu carboxyl groups at low pH, and in that case the peptide collapsed at a lower temperature. The collapse was accompanied by the formation of beta-turns and short distorted beta-sheets. Formation of contacts between hydrophobic side chains drives the collapse at high temperature, but interactions between water and polar groups (Glu and main chain) can attenuate this effect at high pH. The overall competition and balance of the polar and nonpolar groups determined the conformational states of the peptide. Water hydration contributed to the conformational transition, and the peptide and its hydration shell must be considered. Structurally, waters near polar residues mainly formed hydrogen bonds with the protein atoms, while waters around the hydrophobic side chains tended to be parallel to the peptide groups to maximize water-water interactions.     

Preparation of tri(ethylene glycol) grafted core‐shell type polymer support for solid‐phase peptide synthesis

A core‐shell type polymer support for solid‐phase peptide synthesis has been developed for high coupling efficiency of peptides and versatile applications such as on‐bead bioassays. Although various kinds of polymer supports have been developed, they have their own drawbacks including poor accessibility of reagents and incompatibility in aqueous solution. In this paper, we prepared hydrophilic tri(ethylene glycol) (TEG) grafted core‐shell type polymer supports (TEG SURE) for efficient solid‐phase peptide synthesis and on‐bead bioassays. TEG SURE was prepared by grafting TEG derivative on the surface of AM PS resin via biphasic diffusion control method and subsequent acetylation of amine groups which are located at the core region of AM PS resin. The performance of TEG SURE was evaluated by synthesizing several peptides. Three points can be highlighted: (1) easy control of loading level of TEG, (2) improved efficiency of peptide synthesis compared with the conventional resins, and (3) applicability of on‐bead bioassays.     

Acid-induced movements in the glycoprotein shell of an alphavirus turn the spikes into membrane fusion mode

In the icosahedral (T = 4) Semliki Forest virus, the envelope protomers, i.e. E1-E2 heterodimers, make one-to-one interactions with capsid proteins below the viral lipid bilayer, transverse the membrane and form an external glycoprotein shell with projections. The shell is organized by protomer domains interacting as hexamers and pentamers around shell openings at icosahedral 2- and 5-fold axes, respectively, and the projections by other domains associating as trimers at 3- and quasi 3-fold axes. We show here, using cryo- electron microscopy, that low pH, as occurs in the endosomes during virus uptake, results in the relaxation of protomer interactions around the 2- and the 5-fold axes in the shell, and movement of protomers towards 3- and quasi 3-fold axes in a way that reciprocally relocates their putative E1 and E2 domains. This seemed to be facilitated by a trimerization of transmembrane segments at the same axes. The alterations observed help to explain several key features of the spike-mediated membrane fusion reaction, including shell dissolution, heterodimer dissociation, fusion peptide exposure and E1 homotrimerization.     

Observation of unfreezable water in aqueous solution of globule protein by microwave dielectric measurement

A freezing process analyzed by the dielectric method on aqueous solution of albumin has revealed water structure around protein molecule. A relaxation peak due to bound water attached on the protein surface around 100 MHz at room temperature was found. It could be seen commonly in globule proteins. Another peak due to a different kind of unfreezable water was found around 1 GHz at ?6°C. The amount of this water is estimated as 0.36 g water/g protein and in good agreement with that obtained by differential scanning calorimetry and nmr measurements. The water molecules form a shell layer around the protein molecule. © 1995 John Wiley & Sons, Inc.     

Calculations suggest a pathway for the transverse diffusion of a hydrophobic peptide across a lipid bilayer

Alamethicin is a hydrophobic antibiotic peptide 20 amino acids in length. It is predominantly helical and partitions into lipid bilayers mostly in transmembrane orientations. The rate of the peptide transverse diffusion (flip-flop) in palmitoyl-oleyl-phosphatidylcholine vesicles has been measured recently and the results suggest that it involves an energy barrier, presumably due to the free energy of transfer of the peptide termini across the bilayer. We used continuum-solvent model calculations, the known x-ray crystal structure of alamethicin and a simplified representation of the lipid bilayer as a slab of low dielectric constant to calculate the flip-flop rate. We assumed that the lipids adjust rapidly to each configuration of alamethicin in the bilayer because their motions are significantly faster than the average peptide flip-flop time. Thus, we considered the process as a sequence of discrete peptide-membrane configurations, representing critical steps in the diffusion, and estimated the transmembrane flip-flop rate from the calculated free energy of the system in each configuration. Our calculations indicate that the simplest possible pathway, i.e., the rotation of the helix around the bilayer midplane, involving the simultaneous burial of the two termini in the membrane, is energetically unfavorable. The most plausible alternative is a two-step process, comprised of a rotation of alamethicin around its C-terminus residue from the initial transmembrane orientation to a surface orientation, followed by a rotation around the N-terminus residue from the surface to the final reversed transmembrane orientation. This process involves the burial of one terminus at a time and is much more likely than the rotation of the helix around the bilayer midplane. Our calculations give flip-flop rates of approximately 10(-7)/s for this pathway, in accord with the measured value of 1.7 x 10(-6)/s.     

MC-PHS: A Monte Carlo Implementation of the Primary Hydration Shell for Protein Folding and Design

A primary hydration shell (PHS) approach is developed for Monte Carlo simulations of conformationally rich macromolecular systems in an environment that efficiently captures principal solvation effects. It has been previously demonstrated that molecular dynamics using PHS is an efficient method to study peptide structure and dynamics in aqueous solution. Here, we extend the PHS approach to Monte Carlo simulations, whereby a stable shell of water molecules is maintained with a flexible, nonspherical, half-harmonic potential, tuned to maintain a constant restraining energy, with the difference between the restraint and shell energies used to dynamically adjust the shell radius. Examination of the shell and system size dependence of the restraining potential reveals its robustness. Moreover, its suitability for biomolecular simulations is evaluated using small spheres of water, hydration properties of small biological molecules, and configurational sampling of β-hairpin pentapeptide YPGDV. This method, termed MC-PHS, appears to provide efficient representation of dominant solvation effects and should prove useful in the study of protein folding and design.     

Determination of the relative NH proton lifetimes of the peptide analogue viomycin in aqueous solution by NMR-based diffusion measurement

In aqueous solution, exchanging peptide NH protons experience two environments, that of the peptide itself with a relatively slow diffusion coefficient and that of the water solvent with a faster diffusion coefficient. Although in slow exchange on the NMR chemical shift timescale, the magnetic field gradient dependence of the NH peak intensities in an experiment used to measure diffusion coefficients reflects the relative time periods spent in the two environments and this allows the determination of the relative solvent accessibility of exchangeable protons in peptides or proteins. To test this approach, the magnetic field gradient dependent intensities of the chemically shifted amide and amine NH protons of the peptide antibiotic viomycin have been measured using the high resolution longitudinal-eddy-current-delay (LED) NMR method incorporating solvent water peak elimination by non-excitation. The NH resonances of viomycin have been assigned previously and their relative exchange rates determined. Here, the gradient dependence of each NH proton intensity is reported, and these, after a bi- exponential least squares fitting, yield the fractional lifetimes of the protons spent in the peptide and water environments during the diffusion period of the experiment.     

Quantitation of physical-chemical properties of the aqueous phase inside the phoE ionic channel.

The anion-specific channel of the phoE porine is a miniature body of water surrounded by peptide walls. The physical and chemical properties of the water in such a microscopic space were measured by monitoring the dynamics of a well-studied reaction--the protolytic dissociation of a strong acid. To attain this purpose, we allowed pyranine (8-hydroxypyrene-1,3,6-trisulfonate) to bind to the anion-specific channel. The dye is bound, with a 1:1 stoichiometry, with a delta G = -9.5 kcal/mol. Photoexcitation of the dye, to its first electronic singlet state (phi OH*), renders it very acidic and the hydroxyl proton dissociates to H+ and excited anion (phi O*-). We employed single photon-counting time-resolved fluorimetry, to monitor the reversible dissociation of pyranine as it proceeds within the channel and reconstructed the observed signal by a numerical integration of the differential diffusion equation pertinent for a proton within the channel. The most characteristic feature of the water-filled channel, is the intensified electrostatic interactions attained by the low dielectric constant of the diffusion space, epsilon eff = 24. For this reason, the electric field of a few positive charges is sufficient to ensure that an anion entering the channel will be effectively sucked in. The interaction of the water molecules with the peptide structure forming the channel affects the physical properties of the water. Their capacity to conduct proton, quantitated by the protons diffusion coefficient (4.5.10(-5) cm2/s), is reduced by 50% with respect to that of bulk water. The activity of the water in the channel is reduced to alpha H2O = 0.966. These observation are in accord with our previous studies of water in small defined cavities in proteins.     

Elucidating essential role of conserved carboxysomal protein CcmN reveals common feature of bacterial microcompartment assembly

Bacterial microcompartments are organelles composed of a protein shell that surrounds functionally related proteins. Bioinformatic analysis of sequenced genomes indicates that homologs to shell protein genes are widespread among bacteria and suggests that the shell proteins are capable of encapsulating diverse enzymes. The carboxysome is a bacterial microcompartment that enhances CO(2) fixation in cyanobacteria and some chemoautotrophs by sequestering ribulose-1,5-bisphosphate carboxylase/oxygenase and carbonic anhydrase in the microcompartment shell. Here, we report the in vitro and in vivo characterization of CcmN, a protein of previously unknown function that is absolutely conserved in β-carboxysomal gene clusters. We show that CcmN localizes to the carboxysome and is essential for carboxysome biogenesis. CcmN has two functionally distinct regions separated by a poorly conserved linker. The N-terminal portion of the protein is important for interaction with CcmM and, by extension, ribulose-1,5-bisphosphate carboxylase/oxygenase and the carbonic anhydrase CcaA, whereas the C-terminal peptide is essential for interaction with the carboxysome shell. Deletion of the peptide abolishes carboxysome formation, indicating that its interaction with the shell is an essential step in microcompartment formation. Peptides with similar length and sequence properties to those in CcmN can be bioinformatically detected in a large number of diverse proteins proposed to be encapsulated in functionally distinct microcompartments, suggesting that this peptide and its interaction with its cognate shell proteins are common features of microcompartment assembly.     

Decomposition of protein experimental compressibility into intrinsic and hydration shell contributions

The experimental determination of protein compressibility reflects both the protein intrinsic compressibility and the difference between the compressibility of water in the protein hydration shell and bulk water. We use molecular dynamics simulations to explore the dependence of the isothermal compressibility of the hydration shell surrounding globular proteins on differential contributions from charged, polar, and apolar protein-water interfaces. The compressibility of water in the protein hydration shell is accounted for by a linear combination of contributions from charged, polar, and apolar solvent-accessible surfaces. The results provide a formula for the deconvolution of experimental data into intrinsic and hydration contributions when a protein of known structure is investigated. The physical basis for the model is the variation in water density shown by the surface-specific radial distribution functions of water molecules around globular proteins. The compressibility of water hydrating charged atoms is lower than bulk water compressibility, the compressibility of water hydrating apolar atoms is somewhat larger than bulk water compressibility, and the compressibility of water around polar atoms is about the same as the compressibility of bulk water. We also assess whether hydration water compressibility determined from small compound data can be used to estimate the compressibility of hydration water surrounding proteins. The results, based on an analysis from four dipeptide solutions, indicate that small compound data cannot be used directly to estimate the compressibility of hydration water surrounding proteins.     

Influence of hydrophobic mismatch on structures and dynamics of gramicidin a and lipid bilayers

Gramicidin A (gA) is a 15-amino-acid antibiotic peptide with an alternating L-D sequence, which forms (dimeric) bilayer-spanning, monovalent cation channels in biological membranes and synthetic bilayers. We performed molecular dynamics simulations of gA dimers and monomers in all-atom, explicit dilauroylphosphatidylcholine (DLPC), dimyristoylphosphatidylcholine (DMPC), dioleoylphosphatidylcholine (DOPC), and 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) bilayers. The variation in acyl chain length among these different phospholipids provides a way to alter gA-bilayer interactions by varying the bilayer hydrophobic thickness, and to determine the influence of hydrophobic mismatch on the structure and dynamics of both gA channels (and monomeric subunits) and the host bilayers. The simulations show that the channel structure varied little with changes in hydrophobic mismatch, and that the lipid bilayer adapts to the bilayer-spanning channel to minimize the exposure of hydrophobic residues. The bilayer thickness, however, did not vary monotonically as a function of radial distance from the channel. In all simulations, there was an initial decrease in thickness within 4–5 Å from the channel, which was followed by an increase in DOPC and POPC or a further decrease in DLPC and DMPC bilayers. The bilayer thickness varied little in the monomer simulations—except one of three independent simulations for DMPC and all three DLPC simulations, where the bilayer thinned to allow a single subunit to form a bilayer-spanning water-permeable pore. The radial dependence of local lipid area and bilayer compressibility is also nonmonotonic in the first shell around gA dimers due to gA-phospholipid interactions and the hydrophobic mismatch. Order parameters, acyl chain dynamics, and diffusion constants also differ between the lipids in the first shell and the bulk. The lipid behaviors in the first shell around gA dimers are more complex than predicted from a simple mismatch model, which has implications for understanding the energetics of membrane protein-lipid interactions.     

Bioluminescent signals spatially amplified by wavelength-specific diffusion through the shell of a marine snail

Some living organisms produce visible light (bioluminescence) for intra- or interspecific visual communication. Here, we describe a remarkable bioluminescent adaptation in the marine snail Hinea brasiliana. This species produces a luminous display in response to mechanical stimulation caused by encounters with other motile organisms. The light is produced from discrete areas on the snail's body beneath the snail's shell, and must thus overcome this structural barrier to be viewed by an external receiver. The diffusion and transmission efficiency of the shell is greater than a commercial diffuser reference material. Most strikingly, the shell, although opaque and pigmented, selectively diffuses the blue-green wavelength of the species bioluminescence. This diffusion generates a luminous display that is enlarged relative to the original light source. This unusual shell thus allows spatially amplified outward transmission of light communication signals from the snail, while allowing the animal to remain safely inside its hard protective shell.     

Effect of water and enzyme concentration on thermolysin-catalyzed solid-to-solid peptide synthesis

We have studied a thermolysin-catalyzed solid-to-solid dipeptide synthesis using equimolar amounts of Z-Gln-OH and H-Leu-NH2 as model substrates. The high substrate concentrations make this an effective alternative to enzymatic peptide synthesis in organic solvents. Water content was varied in the range of 0 to 600 mL water per mol substrate and enzyme concentration in the range of 0.5 to 10 g/mol of substrates. High yields around 80% conversion and initial rates from 5 to 20 mmol s-1 kg-1 were achieved. The initial rate increases 10-fold on reducing the water content, to reach a pronounced optimum at 40 mL water per mol substrate. Below this, the rate falls to much lower values in a system with no added water, and to zero in a rigorously dried system. This behavior is discussed in terms of two factors: At higher water contents the system is mass transfer limited (as shown by varying enzyme content), and the diffusion distances required vary. At low water levels, effects reflect the stimulation of the enzymatic activity by water.     

Deposition of CdS,CdS/ZnSe and CdS/ZnSe/ZnS shells around CdSeTe alloyed core quantum dots: effects on optical properties

In this work, we synthesized water‐soluble L ‐cysteine‐capped alloyed CdSeTe core quantum dots (QDs) and investigated the structural and optical properties of deposition of each of CdS, CdS/ZnSe and CdS/ZnSe/ZnS shell layers. Photophysical results showed that the overcoating of a CdS shell around the alloyed CdSeTe core [quantum yield (QY) = 8.4%] resulted in effective confinement of the radiative exciton with an improved QY value of 93.5%. Subsequent deposition of a ZnSe shell around the CdSeTe/CdS surface decreased the QY value to 24.7%, but an increase in the QY value of up to 49.5% was observed when a ZnS shell was overcoated around the CdSeTe/CdS/ZnSe surface. QDs with shell layers showed improved stability relative to the core. Data obtained from time‐resolved fluorescence measurements provided useful insight into variations in the photophysical properties of the QDs upon the formation of each shell layer. Our study suggests that the formation of CdSeTe/CdS core/shell QDs meets the requirements of quality QDs in terms of high photoluminescence QY and stability, hence further deposition of additional shells are not necessary in improving the optical properties of the core/shell QDs. Copyright © 2015 John Wiley & Sons, Ltd.