A general method for chemogenetic control of peptide function.

Natural or engineered peptides serve important biological functions. A general approach to achieve chemical-dependent activation of short peptides will be valuable for spatial and temporal control of cellular processes. Here we present a pair of chemically activated protein domains (CAPs) for controlling the accessibility of both the N- and C-terminal portion of a peptide. CAPs were developed through directed evolution of an FK506-binding protein. By fusing a peptide to one or both CAPs, the function of the peptide is blocked until a small molecule displaces them from the FK506-binding protein ligand-binding site. We demonstrate that CAPs are generally applicable to a range of short peptides, including a protease cleavage site, a dimerization-inducing heptapeptide, a nuclear localization signal peptide, and an opioid peptide, with a chemical dependence up to 156-fold. We show that the CAPs system can be utilized in cell cultures and multiple organs in living animals.

Solvation free energy arithmetic for small organic molecules.

Solvent-mediated interactions contribute to ligand binding affinities in computational drug design and provide a challenge for theoretical predictions. In this study, we analyze the solvation free energy of benzene derivatives in water to guide the development of predictive models for solvation free energies and solvent-mediated interactions. We use a spatially resolved analysis of local solvation free energy contributions and define solvation free energy arithmetic, which enable us to construct additive models to describe the solvation of complex compounds. The substituents analyzed in this study are carboxyl and nitro-groups due to their similar sterical requirements but distinct interactions with water. We find that nonadditive solvation free energy contributions are primarily attributed to electrostatics, which are qualitatively reproduced with computationally efficient continuum models. This suggests a promising route for the development of efficient and accurate models for the solvation of complex molecules with varying substitution patterns using solvation arithmetic.

Electric-field induced entropic effects in liquid water.

Externally applied electric fields in liquid water can induce a plethora of effects with wide implications in electrochemistry and hydrogen-based technologies. Although some effort has been made to elucidate the thermodynamics associated with the application of electric fields in aqueous systems, to the best of our knowledge, field-induced effects on the total and local entropy of bulk water have never been presented so far. Here, we report on classical TIP4P/2005 and ab initio molecular dynamics simulations measuring entropic contributions carried by diverse field intensities in liquid water at room temperature. We find that strong fields are capable of aligning large fractions of molecular dipoles. Nevertheless, the order-maker action of the field leads to quite modest entropy reductions in classical simulations. Albeit more significant variations are recorded during first-principles simulations, the associated entropy modifications are small compared to the entropy change involved in the freezing phenomenon, even at intense fields slightly beneath the molecular dissociation threshold. This finding further corroborates the idea that electrofreezing (i.e., the electric-field-induced crystallization) cannot take place in bulk water at room temperature. In addition, here, we propose a molecular-dynamics-based analysis (3D-2PT) that spatially resolves the local entropy and the number density of bulk water under an electric field, which enables us to map their field-induced changes in the environment of reference H2O molecules. By returning detailed spatial maps of the local order, the proposed approach is capable of establishing a link between entropic and structural modifications with atomistic resolution.

Solvent-mediated forces in protein dielectrophoresis.

DEP is an established method to manipulate micrometer-sized particles, but standard continuum theories predict only negligible effects for nanometer-sized proteins despite contrary experimental evidence. A theoretical description of protein DEP needs to consider details on the molecular scale. Previous work toward this goal addressed the role of orientational polarization of static protein dipole moments for dielectrophoretic effects, which successfully predicts the general magnitude of dielectrophoretic forces on proteins but does not readily explain negative DEP forces observed for proteins in some experiments. However, contributions to the protein chemical potential due to protein-water interactions have not yet been considered in this context. Here, we utilize atomistic molecular dynamics simulations to evaluate polarization-induced changes in the protein solvation free energy, which result in a solvent-mediated contribution to dielectrophoretic forces. We quantify solvent-mediated dielectrophoretic forces for two proteins and a small peptide in water, which follow expectations for protein-water dipole-dipole interactions. The magnitude of solvent-mediated dielectrophoretic forces exceeds predictions of nonmolecular continuum theories, but plays a minor role for the total dielectrophoretic force for the simulated proteins due to dominant contributions from the orientational polarization of their static protein dipoles. However, we extrapolate that solvent-mediated contributions to negative protein DEP forces will become increasingly relevant for multidomain proteins, complexes and aggregates with large protein-water interfaces, as well as for high electric field frequencies, which provides a potential mechanism for corresponding experimental observations of negative protein DEP.

Protein flexibility reduces solvent-mediated friction barriers of ligand binding to a hydrophobic surface patch.

Solvent fluctuations have been explored in detail for idealized and rigid hydrophobic model systems, but so far it has remained unclear how internal protein motions and their coupling to the surrounding solvent affect the dynamics of ligand binding to biomolecular surfaces. Here, molecular dynamics simulations were used to elucidate the solvent-mediated binding of a model ligand to the hydrophobic surface patch of ubiquitin. The ligand's friction profiles reveal pronounced long-time correlations and enhanced friction in the vicinity of the protein, similar to idealized hydrophobic surfaces. Interestingly, these effects are shaped by internal protein motions. Protein flexibility modulates water density fluctuations near the hydrophobic surface patch and smooths out the friction profile of ligand binding.

Cooperativity and ion pairing in magnesium sulfate aqueous solutions from the dilute regime to the solubility limit.

We report a terahertz absorption spectroscopy study of MgSO4 aqueous solutions in the concentration range 0.1 mol dm-3 to 2.4 mol dm-3. Accompanying classical MD simulations were carried out that use a polarizable force field parameterized to reproduce the solution thermodynamics. Contrary to prior reports, we find no evidence of contact ion pairs, even close to the solubility limit. Only solvent separated and different types of solvent shared ion pairs are found, being abundant even at the lowest concentration investigated here. The structure of the solution is concentration-dependent: the number of both types of ion pairs grows with increasing salt concentration. The combined theoretical and experimental analysis of the spectra in the frequency region 50-640 cm-1 suggests that the dynamics of water directly between two ions in solvent shared configuration is very strongly perturbed, via a cooperative, supra-additive, effect arising from the two ions. At high concentrations, the results support a scenario, where the perturbations in the water dynamics extend up to the third hydration layer via a cooperative, but additive, effect involving multiple ions. The SO42- and its hydration shell are much more strongly perturbed by the presence of the counterions than the first hydration shell of Mg2+. It is further shown that our simulations and observations are in agreement with thermodynamic properties of aqueous MgSO4 solutions derived by other methods.

Pressure effects on collective density fluctuations in water and protein solutions.

Neutron Brillouin scattering and molecular dynamics simulations have been used to investigate protein hydration water density fluctuations as a function of pressure. Our results show significant differences between the pressure and density dependence of collective dynamics in bulk water and in concentrated protein solutions. Pressure-induced changes in the tetrahedral order of the water HB network have direct consequences for the high-frequency sound velocity and damping coefficients, which we find to be a sensitive probe for changes in the HB network structure as well as the wetting of biomolecular surfaces.

Stability Effect of Quinary Interactions Reversed by Single Point Mutations.

In cells, proteins are embedded in a crowded environment that controls their properties via manifold avenues including weak protein-macromolecule interactions. A molecular level understanding of these quinary interactions and their contribution to protein stability, function, and localization in the cell is central to modern structural biology. Using a mutational analysis to quantify the energetic contributions of single amino acids to the stability of the ALS related protein superoxide dismutase I (SOD1) in mammalian cells, we show that quinary interactions destabilize SOD1 by a similar energetic offset for most of the mutants, but there are notable exceptions: Mutants that alter its surface properties can even lead to a stabilization of the protein in the cell as compared to the test tube. In conclusion, quinary interactions can amplify and even reverse the mutational response of proteins, being a key aspect in pathogenic protein misfolding and aggregation.

Atomistic characterization of collective protein-water-membrane dynamics.

Correlated vibrational motion on the sub-picosecond timescale and associated collective dynamics in a protein-membrane environment are characterized using molecular dynamics simulations. We specifically analyze correlated motion of a membrane-associated protein and a lipid bilayer for distinct separation distances. Correlated vibrations persist up to distances of 25 Å between both biomolecular surfaces. These correlations are mediated by separating layers of water molecules, whose collective properties are altered by the simultaneous presence of protein and lipid bilayer interfaces.

Hydration-mediated stiffening of collective membrane dynamics by cholesterol.

The collective behaviour of individual lipid molecules determines the properties of phospholipid membranes. However, the collective molecular motions often remain challenging to characterise at the desired spatial and temporal resolution. Here we study collective vibrational motion on picosecond time scales in dioleoylphosphatidylcholine lipid bilayers with varying cholesterol content using all-atom molecular dynamics simulations. Cholesterol is found to not only laterally compact the lipid bilayer, but also to change the velocity of longitudinal density fluctuations propagating in the plane of the membrane. Cholesterol-induced reduction of the area per lipid alters the collective dynamics of the lipid headgroups, but not of the lipid tails. The introduction of cholesterol reduces the number of water molecules interacting with the lipid headgroups, leading to a decrease in the velocity of the laterally-propagating sound mode. Thus, the stiffening effect of cholesterol is found to be indirect: decreasing the area per lipid weakens the interactions between the lipid headgroups and water. The collective modes characterised in this work can enable the membrane to dissipate excess energy and thus maintain its structural integrity, e.g., under mechanical stress.

Heterogeneity of water structure and dynamics at the protein-water interface.

In this molecular dynamics simulation study, we analyze the local structural and dynamic properties of water hydrating the protein ubiquitin on a spatial grid with 1 Å resolution. This allows for insights into the spatial distribution of water number densities, molecular orientations, translations, and rotations as a function of distance from the protein surface. Water molecule orientations follow a heterogeneous distribution with preferred local orientations of water dipoles and O-H bond vectors up to 10-15 Å distances from the protein, while local variations of the water number density converge to homogeneous bulk-like values within less than 8 Å. Interestingly, we find that the long-ranged orientational structure of water does not impact either the translational or rotational dynamics of water. Instead, heterogeneous distributions of local dynamical parameters and averaged dynamical retardation factors are only found close to the protein surface and follow a distance dependence comparable to heterogeneities in the local water number density. This study shows that the formation of nanodomains of preferred water orientations far from the protein does not significantly impact dynamical processes probed as a non-local average in most experiments.

Anomalous behavior of water inside the SecY translocon.

The heterotrimeric SecY translocon complex is required for the cotranslational assembly of membrane proteins in bacteria and archaea. The insertion of transmembrane (TM) segments during nascent-chain passage through the translocon is generally viewed as a simple partitioning process between the water-filled translocon and membrane lipid bilayer, suggesting that partitioning is driven by the hydrophobic effect. Indeed, the apparent free energy of partitioning of unnatural aliphatic amino acids on TM segments is proportional to accessible surface area, which is a hallmark of the hydrophobic effect [Öjemalm K, et al. (2011) Proc Natl Acad Sci USA 108(31):E359-E364]. However, the apparent partitioning solvation parameter is less than one-half the value expected for simple bulk partitioning, suggesting that the water in the translocon departs from bulk behavior. To examine the state of water in a SecY translocon complex embedded in a lipid bilayer, we carried out all-atom molecular-dynamics simulations of the Pyrococcus furiosus SecYE, which was determined to be in a "primed" open state [Egea PF, Stroud RM (2010) Proc Natl Acad Sci USA 107(40):17182-17187]. Remarkably, SecYE remained in this state throughout our 450-ns simulation. Water molecules within SecY exhibited anomalous diffusion, had highly retarded rotational dynamics, and aligned their dipoles along the SecY transmembrane axis. The translocon is therefore not a simple water-filled pore, which raises the question of how anomalous water behavior affects the mechanism of translocon function and, more generally, the partitioning of hydrophobic molecules. Because large water-filled cavities are found in many membrane proteins, our findings may have broader implications.

Spatially Heterogeneous Surface Water Diffusivity around Structured Protein Surfaces at Equilibrium.

Hydration water on the surface of a protein is thought to mediate the thermodynamics of protein-ligand interactions. For hydration water to play a role beyond modulating global protein solubility or stability, the thermodynamic properties of hydration water must reflect on the properties of the heterogeneous protein surface and thus spatially vary over the protein surface. A potent read-out of local variations in thermodynamic properties of hydration water is its equilibrium dynamics spanning picosecond to nanosecond time scales. In this study, we employ Overhauser dynamic nuclear polarization (ODNP) to probe the equilibrium hydration water dynamics at select sites on the surface of Chemotaxis Y (CheY) in dilute solution. ODNP reports on site-specific hydration water dynamics within 5-10 Å of a label tethered to the biomolecular surface on two separate time scales of motion, corresponding to diffusive water (DW) and protein-water coupled motions, referred to as bound water (BW). We find DW dynamics to be highly heterogeneous across the surface of CheY. We identify a significant correlation between DW dynamics and the local hydropathy of the CheY protein surface, empirically determined by molecular dynamics (MD) simulations, and find the more hydrophobic sites to be hydrated with slower diffusing water. Furthermore, we compare the hydration water dynamics on different polypeptides and liposome surfaces and find the DW dynamics on globular proteins to be significantly more heterogeneous than on intrinsically disordered proteins (IDPs), peptides, and liposomes. The heterogeneity in the hydration water dynamics suggests that structured proteins have the capacity to encode information into the surrounding hydration shell.

Hydration Dynamics of a Peripheral Membrane Protein.

Water dynamics in the hydration shell of the peripheral membrane protein annexin B12 were studied using MD simulations and Overhauser DNP-enhanced NMR. We show that retardation of water motions near phospholipid bilayers is extended by the presence of a membrane-bound protein, up to around 10 Å above that protein. Near the membrane surface, electrostatic interactions with the lipid head groups strongly slow down water dynamics, whereas protein-induced water retardation is weaker and dominates only at distances beyond 10 Å from the membrane surface. The results can be understood from a simple model based on additive contributions from the membrane and the protein to the activation free energy barriers of water diffusion next to the biomolecular surfaces. Furthermore, analysis of the intermolecular vibrations of the water network reveals that retarded water motions near the membrane shift the vibrational modes to higher frequencies, which we used to identify an entropy gradient from the membrane surface toward the bulk water. Our results have implications for processes that take place at lipid membrane surfaces, including molecular recognition, binding, and protein-protein interactions.

"Bind and Crawl" Association Mechanism of Leishmania major Peroxidase and Cytochrome c Revealed by Brownian and Molecular Dynamics Simulations.

Leishmania major, the parasitic causative agent of leishmaniasis, produces a heme peroxidase (LmP), which catalyzes the peroxidation of mitochondrial cytochrome c (LmCytc) for protection from reactive oxygen species produced by the host. The association of LmP and LmCytc, which is known from kinetics measurements to be very fast (∼10(8) M(-1) s(-1)), does not involve major conformational changes and has been suggested to be dominated by electrostatic interactions. We used Brownian dynamics simulations to investigate the mechanism of formation of the LmP-LmCytc complex. Our simulations confirm the importance of electrostatic interactions involving the negatively charged D211 residue at the LmP active site, and reveal a previously unrecognized role in complex formation for negatively charged residues in helix A of LmP. The crystal structure of the D211N mutant of LmP reported herein is essentially identical to that of wild-type LmP, reinforcing the notion that it is the loss of charge at the active site, and not a change in structure, that reduces the association rate of the D211N variant of LmP. The Brownian dynamics simulations further show that complex formation occurs via a "bind and crawl" mechanism, in which LmCytc first docks to a location on helix A that is far from the active site, forming an initial encounter complex, and then moves along helix A to the active site. An atomistic molecular dynamics simulation confirms the helix A binding site, and steady state activity assays and stopped-flow kinetics measurements confirm the role of helix A charges in the association mechanism.

Curvature dependence of hydrophobic hydration dynamics.

We investigate the solute curvature dependence of water dynamics in the vicinity of hydrophobic spherical solutes using molecular dynamics simulations. For both the lateral and perpendicular diffusivity, as well as for H-bond kinetics of water in the first hydration shell, we find a nonmonotonic solute-size dependence, exhibiting extrema close to the well-known structural crossover length scale for hydrophobic hydration. Additionally, we find an apparent anomalous diffusion for water moving parallel to the surface of small solutes, which, however, can be explained by topology effects. Our findings regarding the intimate connection between solute curvature and water dynamics has implications for our understanding of hydration dynamics at heterogeneous biomolecular surfaces.

Excluded-volume effects in living cells.

Biomolecules evolve and function in densely crowded and highly heterogeneous cellular environments. Such conditions are often mimicked in the test tube by the addition of artificial macromolecular crowding agents. Still, it is unclear if such cosolutes indeed reflect the physicochemical properties of the cellular environment as the in-cell crowding effect has not yet been quantified. We have developed a macromolecular crowding sensor based on a FRET-labeled polymer to probe the macromolecular crowding effect inside single living cells. Surprisingly, we find that excluded-volume effects, although observed in the presence of artificial crowding agents, do not lead to a compression of the sensor in the cell. The average conformation of the sensor is similar to that in aqueous buffer solution and cell lysate. However, the in-cell crowding effect is distributed heterogeneously and changes significantly upon cell stress. We present a tool to systematically study the in-cell crowding effect as a modulator of biomolecular reactions.

Resolving anisotropic distributions of correlated vibrational motion in protein hydration water.

In this study, we analyze correlations of vibrational motion on the surface of a small globular protein and in its hydration shell. In contrast to single particle hydration water dynamics, which are perturbed by interactions with the protein solute only in the first few hydration layers, we find that correlated, collective motions extend into the surrounding solvent on a 10 Å length scale, specifically at far-infrared frequencies below 100 cm(-1). As a function of frequency, we analyze the distribution of correlated longitudinal motions in the three-dimensional environment of the protein solute, as well as in the vicinity of different protein-water interfaces. An anisotropic distribution of these correlations is observed, which is related to specific protein-water vibrations and interactions at the interfaces, as well as flexibilities of solvent exposed sites. Our results show that coupling of protein and water dynamics leaves a three-dimensional imprint in the collective dynamics of its hydration shell, and we discuss potential implications for biomolecular function, e.g., molecular recognition and binding, and the dynamical coupling of proteins to their native solvation environment.

Exploring Conformational Landscapes Along Anharmonic Low-Frequency Vibrations.

We aim to automatize the identification of collective variables to simplify and speed up enhanced sampling simulations of conformational dynamics in biomolecules. We focus on anharmonic low-frequency vibrations that exhibit fluctuations on time scales faster than conformational transitions but describe a path of least resistance toward structural change. A key challenge is that harmonic approximations are ill-suited to characterize these vibrations, which are observed at far-infrared frequencies and are easily excited by thermal collisions at room temperature. Here, we approached this problem with a frequency-selective anharmonic (FRESEAN) mode analysis that does not rely on harmonic approximations and successfully isolates anharmonic low-frequency vibrations from short molecular dynamics simulation trajectories. We applied FRESEAN mode analysis to simulations of alanine dipeptide, a common test system for enhanced sampling simulation protocols, and compared the performance of isolated low-frequency vibrations to conventional user-defined collective variables (here backbone dihedral angles) in enhanced sampling simulations. The comparison shows that enhanced sampling along anharmonic low-frequency vibrations not only reproduces known conformational dynamics but can even further improve the sampling of slow transitions compared to user-defined collective variables. Notably, free energy surfaces spanned by low-frequency anharmonic vibrational modes exhibit lower barriers associated with conformational transitions relative to representations in backbone dihedral space. We thus conclude that anharmonic low-frequency vibrations provide a promising path for highly effective and fully automated enhanced sampling simulations of conformational dynamics in biomolecules.

Linear and Nonlinear Dielectric Response of Intrinsically Disordered Proteins.

Linear and nonlinear dielectric responses of solutions of intrinsically disordered proteins (IDPs) were analyzed by combining molecular dynamics simulations with formal theories. A large increment of the linear dielectric function over that of the solvent is found and related to large dipole moments of IDPs. The nonlinear dielectric effect (NDE) of the IDP far exceeds that of the bulk electrolyte, offering a route to interrogate protein conformational and rotational statistics and dynamics. Conformational flexibility of the IDP makes the dipole moment statistics consistent with the gamma/log-normal distributions and contributes to the NDE through the dipole moment's non-Gaussian parameter. The intrinsic non-Gaussian parameter of the dipole moment combines with the protein osmotic compressibility in the nonlinear dielectric susceptibility when dipolar correlations are screened by the electrolyte. The NDE is dominated by dipolar correlations when electrolyte screening is reduced.