ALBERT

Description: Coherence oscillations measured in two-dimensional (2D) electronic spectra of pigment-protein complexes may have electronic, vibrational, or mixed-character vibronic origins, which depend on the degree of electronic-vibrational mixing. Oscillations from intrapigment vibrations can obscure the inter-site coherence lifetime of interest in elucidating the mechanisms of energy transfer in photosynthetic light-harvesting. Huang-Rhys factors ( S ) for low-frequency vibrations in Chlorophyll and Bacteriochlorophyll are quite small ( S ≤ 0.05), so it is often assumed that these vibrations influence neither 2D spectra nor inter-site coherence dynamics. In this work, we explore the influence of S within this range on the oscillatory signatures in simulated 2D spectra of a pigment heterodimer. To visualize the inter-site coherence dynamics underlying the 2D spectra, we introduce a formalism which we call the “site-probe response.” By comparing the calculated 2D spectra with the site-probe response, we show that an on-resonance vibration with Huang-Rhys factor as small as S = 0.005 and the most strongly coupled off-resonance vibrations ( S = 0.05) give rise to long-lived, purely vibrational coherences at 77 K. We moreover calculate the correlation between optical pump interactions and subsequent entanglement between sites, as measured by the concurrence. At 77 K, greater long-lived inter-site coherence and entanglement appear with increasing S . This dependence all but vanishes at physiological temperature, as environmentally induced fluctuations destroy the vibronic mixing.

Description: We address the problem of simulating biochemical reaction networks with time-dependent rates and propose a new algorithm based on our rejection-based stochastic simulation algorithm (RSSA) [Thanh et al. , J. Chem. Phys. 141 (13), 134116 (2014)]. The computation for selecting next reaction firings by our time-dependent RSSA (tRSSA) is computationally efficient. Furthermore, the generated trajectory is exact by exploiting the rejection-based mechanism. We benchmark tRSSA on different biological systems with varying forms of reaction rates to demonstrate its applicability and efficiency. We reveal that for nontrivial cases, the selection of reaction firings in existing algorithms introduces approximations because the integration of reaction rates is very computationally demanding and simplifying assumptions are introduced. The selection of the next reaction firing by our approach is easier while preserving the exactness.

Description: Reference interaction site model (RISM) has recently become a popular approach in the study of thermodynamical and structural properties of the solvent around macromolecules. On the other hand, it was widely suggested that there exists water density depletion around large hydrophobic solutes (〉1 nm), and this may pose a great challenge to the RISM theory. In this paper, we develop a new analytical theory, the Reference Interaction Site Model with Hydrophobicity induced density Inhomogeneity (RISM-HI), to compute solvent radial distribution function (RDF) around large hydrophobic solute in water as well as its mixture with other polyatomic organic solvents. To achieve this, we have explicitly considered the density inhomogeneity at the solute-solvent interface using the framework of the Yvon-Born-Green hierarchy, and the RISM theory is used to obtain the solute-solvent pair correlation. In order to efficiently solve the relevant equations while maintaining reasonable accuracy, we have also developed a new closure called the D2 closure. With this new theory, the solvent RDFs around a large hydrophobic particle in water and different water-acetonitrile mixtures could be computed, which agree well with the results of the molecular dynamics simulations. Furthermore, we show that our RISM-HI theory can also efficiently compute the solvation free energy of solute with a wide range of hydrophobicity in various water-acetonitrile solvent mixtures with a reasonable accuracy. We anticipate that our theory could be widely applied to compute the thermodynamic and structural properties for the solvation of hydrophobic solute.

Description: Recently proposed methods for DNA sequencing involve the use of cleavage-based enzymes attached to the opening of a nanopore. The idea is that DNA interacting with either an exonuclease or polymerase protein will lead to a small molecule being cleaved near the mouth of the nanopore, and subsequent entry into the pore will yield information about the DNA sequence. The prospects for this approach seem promising, but it has been shown that diffusion related effects impose a limit on the capture probability of molecules by the pore, which limits the efficacy of the technique. Here, we revisit the problem with the goal of optimizing the capture probability via a step decrease in the nucleotide diffusion coefficient between the pore and bulk solutions. It is shown through random walk simulations and a simplified analytical model that decreasing the molecule’s diffusion coefficient in the bulk relative to its value in the pore increases the nucleotide capture probability. Specifically, we show that at sufficiently high applied transmembrane potentials (≥100 mV), increasing the potential by a factor f is equivalent to decreasing the diffusion coefficient ratio D bulk / D pore by the same factor f . This suggests a promising route toward implementation of cleavage-based sequencing protocols. We also discuss the feasibility of forming a step function in the diffusion coefficient across the pore-bulk interface.

Description: Large-scale domain motions in alcohol dehydrogenase (ADH) have been observed previously by neutron spin-echo spectroscopy (NSE). We have extended the investigation on the dynamics of ADH in solution by using high-resolution neutron time-of-flight (TOF) and neutron backscattering (BS) spectroscopy in the incoherent scattering range. The observed hydrogen dynamics were interpreted in terms of three mobility classes, which allowed a simultaneous description of the measured TOF and BS spectra. In addition to the slow global protein diffusion and domain motions observed by NSE, a fast internal process could be identified. Around one third of the protons in ADH participate in the fast localized diffusive motion. The diffusion coefficient of the fast internal motions is around two third of the value of the surrounding D 2 O solvent. It is tempting to associate the fast internal process with solvent exposed amino acid residues with dangling side chains.

Description: We analyse the effect of intrinsic fluctuations on the properties of bistable stochastic systems with time scale separation operating under quasi-steady state conditions. We first formulate a stochastic generalisation of the quasi-steady state approximation based on the semi-classical approximation of the partial differential equation for the generating function associated with the chemical master equation. Such approximation proceeds by optimising an action functional whose associated set of Euler-Lagrange (Hamilton) equations provides the most likely fluctuation path. We show that, under appropriate conditions granting time scale separation, the Hamiltonian can be re-scaled so that the set of Hamilton equations splits up into slow and fast variables, whereby the quasi-steady state approximation can be applied. We analyse two particular examples of systems whose mean-field limit has been shown to exhibit bi-stability: an enzyme-catalysed system of two mutually inhibitory proteins and a gene regulatory circuit with self-activation. Our theory establishes that the number of molecules of the conserved species is order parameters whose variation regulates bistable behaviour in the associated systems beyond the predictions of the mean-field theory. This prediction is fully confirmed by direct numerical simulations using the stochastic simulation algorithm. This result allows us to propose strategies whereby, by varying the number of molecules of the three conserved chemical species, cell properties associated to bistable behaviour (phenotype, cell-cycle status, etc.) can be controlled.

Description: We report a facile method of doping cations from an electrospray ionization (ESI) source into superfluid helium droplets. By decelerating and stopping the ion pulse of reserpine and substance P from an ESI source in the path of the droplet beam, about 10 4 ion-doped droplets (one ion per droplet) can be recorded, corresponding to a pickup efficiency of nearly 1 out of 1000 ions. We attribute the success of this simple approach to the long residence time of the cations in the droplet beam. The resulting size of the doped droplets, on the order of 10 5 /droplet, is measured using deflection and retardation methods. Our method does not require an ion trap in the doping region, which significantly simplifies the experimental setup and procedure for future spectroscopic and diffraction studies.

Description: A number of important cellular functions rely on the formation of intracellular protein concentration gradients. Experimental studies discovered a number of mechanisms for the formation of such gradients. One of the mechanisms relies on the intracellular shuttling of a protein that interconverts between the two states with different diffusivities, under the action of two enzymes, one of which is localized to the plasma membrane, whereas the second is uniformly distributed in the cytoplasm. Recent work reported an analytical solution for the steady state gradient in this mechanism, obtained in the framework of a one-dimensional reaction-diffusion model. Here, we study the dynamics in this model and derive analytical expressions for the Laplace transforms of the time-dependent concentration profiles in terms of elementary transcendental functions. Inverting these transforms numerically, one can obtain time-dependent concentration profiles of the two forms of the protein.

Description: Using the probabilistic language of conditional expectations, we reformulate the force matching method for coarse-graining of molecular systems as a projection onto spaces of coarse observables. A practical outcome of this probabilistic description is the link of the force matching method with thermodynamic integration. This connection provides a way to systematically construct a local mean force and to optimally approximate the potential of mean force through force matching. We introduce a generalized force matching condition for the local mean force in the sense that allows the approximation of the potential of mean force under both linear and non-linear coarse graining mappings (e.g., reaction coordinates, end-to-end length of chains). Furthermore, we study the equivalence of force matching with relative entropy minimization which we derive for general non-linear coarse graining maps. We present in detail the generalized force matching condition through applications to specific examples in molecular systems.

Description: We introduce a set of charged patchy particle models (CPPMs) in order to systematically study the influence of electrostatic charge patchiness and multipolarity on macromolecular interactions by means of implicit-solvent, explicit-ion Langevin dynamics simulations employing the Gromacs software. We consider well-defined zero-, one-, and two-patched spherical globules each of the same net charge and (nanometer) size which are composed of discrete atoms. The studied mono- and multipole moments of the CPPMs are comparable to those of globular proteins with similar size. We first characterize ion distributions and electrostatic potentials around a single CPPM. Although angle-resolved radial distribution functions reveal the expected local accumulation and depletion of counter- and co-ions around the patches, respectively, the orientation-averaged electrostatic potential shows only a small variation among the various CPPMs due to space charge cancellations. Furthermore, we study the orientation-averaged potential of mean force (PMF), the number of accumulated ions on the patches, as well as the CPPM orientations along the center-to-center distance of a pair of CPPMs. We compare the PMFs to the classical Derjaguin-Verwey-Landau-Overbeek theory and previously introduced orientation-averaged Debye-Hückel pair potentials including dipolar interactions. Our simulations confirm the adequacy of the theories in their respective regimes of validity, while low salt concentrations and large multipolar interactions remain a challenge for tractable theoretical descriptions.