We report on a combined atomistic molecular dynamics simulation and implicit solvent analysis of a generic hydrophobic pocket-ligand (host-guest) system. The approaching ligand induces complex wetting/dewetting transitions in the weakly solvated pock
et. The transitions lead to bimodal solvent fluctuations which govern magnitude and range of the pocket-ligand attraction. A recently developed implicit water model, based on the minimization of a geometric functional, captures the sensitive aqueous interface response to the concave-convex pocket-ligand configuration semi-quantitatively.
Approximately 1% of the known protein structures display knotted configurations in their native fold but their function is not understood. It has been speculated that the entanglement may inhibit mechanical protein unfolding or transport, e.g., as in
cellular threading or translocation processes through narrow biological pores. Here we investigate tigh peptide knot (TPK) characteristics in detail by pulling selected 3_1 and 4_1-knotted peptides using all-atom molecular dynamics computer simulations. We find that the 3_1 and 4_1-TPK lengths are typically Delta l~4.7 nm and 6.9 nm, respectively, for a wide range of tensions (F < 1.5 nN), pointing to a pore diameter of ~2 nm below which a translocated knotted protein might get stuck. The 4_1-knot length is in agreement with recent AFM pulling experiments. Detailed TPK characteristics however, may be sequence-specific: we find a different size and structural behavior in polyglycines, and, strikingly, a strong hydrogen bonding and water molecule trapping capability of hydrophobic TPKs due to side chain shielding of the polar TPK core. Water capturing and release is found to be controlla ble by the tightening force in a few cases. These mechanisms result into a sequence-specific locking and metastability of TPKs what might lead to a blocking of knotted peptide transport at designated sequence-positions. Intriguingly, macroscopic tight 4_1-knot structures are reproduced microscopically (figure-of-eight vs. the pretzel) and can be tuned by sequence in contrast to mathematical predictions. Our findings may explain a function of knots in native proteins, challenge previous studies on macromolecular knots, and may find use in bio- and nanotechnology.
A level-set method is developed for numerically capturing the equilibrium solute-solvent interface that is defined by the recently proposed variational implicit solvent model (Dzubiella, Swanson, and McCammon, Phys. Rev. Lett. {bf 104}, 527 (2006) an
d J. Chem.Phys. {bf 124}, 084905 (2006)). In the level-set method, a possible solute-solvent interface is represented by the zero level-set (i.e., the zero level surface) of a level-set function and is eventually evolved into the equilibrium solute-solvent interface. The evolution law is determined by minimization of a solvation free energy {it functional} that couples both the interfacial energy and the van der Waals type solute-solvent interaction energy. The surface evolution is thus an energy minimizing process, and the equilibrium solute-solvent interface is an output of this process. The method is implemented and applied to the solvation of nonpolar molecules such as two xenon atoms, two parallel paraffin plates, helical alkane chains, and a single fullerene $C_{60}$. The level-set solutions show good agreement for the solvation energies when compared to available molecular dynamics simulations. In particular, the method captures solvent dewetting (nanobubble formation) and quantitatively describes the interaction in the strongly hydrophobic plate system.
The structure of a single alanine-based Ace-AEAAAKEAAAKA-Nme peptide in explicit aqueous electrolyte solutions (NaCl, KCl, NaI, and KF) at large salt concentrations (3-4 M) is investigated using 1 microsecond molecular dynamics (MD) computer simulati
ons. The peptide displays 71 alpha-helical structure without salt and destabilizes with the addition of NaCl in agreement with experiments of a somewhat longer version. It is mainly stabilized by direct and indirect (i+4)EK salt bridges between the Lys and Glu side chains and a concomitant backbone shielding mechanism. NaI is found to be a stronger denaturant than NaCl, while the potassium salts hardly show influence. Investigation of the molecular structures reveals that consistent with recent experiments Na+ has a much stronger affinity to side chain carboxylates and backbone carbonyls than K+, thereby weakening salt bridges and secondary structure hydrogen bonds. At the same time the large I- has a considerable affinity to the nonpolar alanine in line with recent observations of a large propensity of I- to adsorb to simple hydrophobes, and thereby assists Na+ in its destabilizing action. In the denatured states of the peptide novel long-lived (10-20 ns) loop-configurations are observed in which single Na+ ions and water molecules are hydrogen-bonded to multiple backbone carbonyls. In an attempt to analyze the denaturation behavior within the preferential interaction formalism, we find indeed that for the strongest denaturant, NaI, the protein is least hydrated. Additionally, a possible indication for protein denaturation might be a preferential solvation of the first solvation shell of the peptide backbone by the destabilizing cosolute (sodium).
