No Arabic abstract
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 simulations. 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).
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.
Dry lakes covered with a salt crust organised into beautifully patterned networks of narrow ridges are common in arid regions. Here, we consider the initial instability and the ultimate fate of buoyancy-driven convection that could lead to such patterns. Specifically, we look at convection in a deep porous medium with a constant through-flow boundary condition on a horizontal surface, which resembles the situation found below an evaporating salt lake. The system is scaled to have only one free parameter, the Rayleigh number, which characterises the relative driving force for convection. We then solve the resulting linear stability problem for the onset of convection. Further exploring the non-linear regime of this model with pseudo-spectral numerical methods, we demonstrate how the growth of small downwelling plumes is itself unstable to coarsening, as the system develops into a dynamic steady state. In this mature state we show how the typical speeds and length-scales of the convective plumes scale with forcing conditions, and the Rayleigh number. Interestingly, a robust length-scale emerges for the pattern wavelength, which is largely independent of the driving parameters. Finally, we introduce a spatially inhomogeneous boundary condition -- a modulated evaporation rate -- to mimic any feedback between a growing salt crust and the evaporation over the dry salt lake. We show how this boundary condition can introduce phase-locking of the downwelling plumes below sites of low evaporation, such as at the ridges of salt polygons.
Secondary structure formation of nucleic acids strongly depends on salt concentration and temperature. We develop a theory for RNA folding that correctly accounts for sequence effects, the entropic contributions associated with loop formation, and salt effects. Using an iterative expression for the partition function that neglects pseudoknots, we calculate folding free energies and minimum free energy configurations based on the experimentally derived base pairing free energies. The configurational entropy of loop formation is modeled by the asymptotic expression -c ln m, where m is the length of the loop and c the loop exponent, which is an adjustable constant. Salt effects enter in two ways: first, we derive salt induced modifications of the free energy parameters for describing base pairing and, second, we include the electrostatic free energy for loop formation. Both effects are modeled on the Debye-Hueckel level including counterion condensation. We validate our theory for two different RNA sequences: For tRNA-phe, the resultant heat capacity curves for thermal denaturation at various salt concentrations accurately reproduce experimental results. For the P5ab RNA hairpin, we derive the global phase diagram in the three-dimensional space spanned by temperature, stretching force, and salt concentration and obtain good agreement with the experimentally determined critical unfolding force. We show that for a proper description of RNA melting and stretching, both salt and loop entropy effects are needed.
Solid-state nanopores are promising tools for single molecule detection of both DNA and proteins. In this study, we investigate the patterns of ionic current blockades as DNA translocates into or out of the geometric confinement of such conically shaped pores. We studied how the geometry of a nanopore affects the detected ionic current signal of a translocating DNA molecule over a wide range of salt concentration. The blockade level in the ionic current depends on the translocation direction at a high salt concentration, and at lower salt concentrations we find a non-intuitive ionic current decrease and increase within each single event for the DNA translocations exiting from confinement. We use recently developed DNA rulers with markers and finite element calculations to explain our observations. Our calculations explain the shapes of the signals observed at all salt concentrations and show that the unexpected current decrease and increase are due to the competing effects of ion concentration polarization and geometric exclusion of ions. Our analysis shows that over a wide range of geometry, voltage and salt concentrations we are able to understand the ionic current signals of DNA in asymmetric nanopores enabling signal optimization in molecular sensing applications.
We review the background, theory and general equations for the analysis of equilibrium protein unfolding experiments, focusing on denaturant and heat-induced unfolding. The primary focus is on the thermodynamics of reversible folding/unfolding transitions and the experimental methods that are available for extracting thermodynamic parameters. We highlight the importance of modelling both how the folding equilibrium depends on a perturbing variable such as temperature or denaturant concentration, and the importance of modelling the baselines in the experimental observables.