Cluster glasses of ultrasoft particles

We present molecular dynamics (MD) simulations results for dense fluids of ultrasoft, fully penetrable particles. These are a binary mixture and a polydisperse system of particles interacting via the generalized exponential model, which is known to yield cluster crystal phases for the corresponding monodisperse systems. Because of the dispersity in the particle size, the systems investigated in this work do not crystallize and form disordered cluster phases. The clusteringtransition appears as a smooth crossover to a regime in which particles are mostly located in clusters, isolated particles being infrequent. The analysis of the internal cluster structure reveals microsegregation of the big and small particles, with a strong homo-coordination in the binary mixture. Upon further lowering the temperature below the clusteringtransition, the motion of the clusters" centers-of-mass slows down dramatically, giving way to a cluster glass transition. In the cluster glass, the diffusivities remain finite and display an activated temperature dependence, indicating that relaxation in the cluster glass occurs via particle hopping in a nearly arrested matrix of clusters. Finally we discuss the influence of the microscopic dynamics on the transport properties by comparing the MD results with Monte Carlo simulations.

Asymmetric stochastic switching driven by intrinsic molecular noise

Low-copy-number molecules are involved in many functions in cells. The intrinsic fluctuations of these numbers can enable stochastic switching between multiple steady states, inducing phenotypic variability. Herein we present a theoretical and computational study based on Master Equations and Fokker-Planck and Langevin descriptions of stochastic switching for a genetic circuit of autoactivation. We show that in this circuit the intrinsic fluctuations arising from low-copy numbers, which are inherently state-dependent, drive asymmetric switching. These theoretical results are consistent with experimental data that have been reported for the bistable system of the gallactose signaling network in yeast. Our study unravels that intrinsic fluctuations, while not required to describe bistability, are fundamental to understand stochastic switching and the dynamical relative stability of multiple states.

Cluster glasses of ultrasoft particles

We present molecular dynamics (MD) simulations results for dense fluids of ultrasoft, fully penetrable particles. These are a binary mixture and a polydisperse system of particles interacting via the generalized exponential model, which is known to yield cluster crystal phases for the corresponding monodisperse systems. Because of the dispersity in the particle size, the systems investigated in this work do not crystallize and form disordered cluster phases. The clusteringtransition appears as a smooth crossover to a regime in which particles are mostly located in clusters, isolated particles being infrequent. The analysis of the internal cluster structure reveals microsegregation of the big and small particles, with a strong homo-coordination in the binary mixture. Upon further lowering the temperature below the clusteringtransition, the motion of the clusters" centers-of-mass slows down dramatically, giving way to a cluster glass transition. In the cluster glass, the diffusivities remain finite and display an activated temperature dependence, indicating that relaxation in the cluster glass occurs via particle hopping in a nearly arrested matrix of clusters. Finally we discuss the influence of the microscopic dynamics on the transport properties by comparing the MD results with Monte Carlo simulations.

Asymmetric stochastic switching driven by intrinsic molecular noise

Low-copy-number molecules are involved in many functions in cells. The intrinsic fluctuations of these numbers can enable stochastic switching between multiple steady states, inducing phenotypic variability. Herein we present a theoretical and computational study based on Master Equations and Fokker-Planck and Langevin descriptions of stochastic switching for a genetic circuit of autoactivation. We show that in this circuit the intrinsic fluctuations arising from low-copy numbers, which are inherently state-dependent, drive asymmetric switching. These theoretical results are consistent with experimental data that have been reported for the bistable system of the gallactose signaling network in yeast. Our study unravels that intrinsic fluctuations, while not required to describe bistability, are fundamental to understand stochastic switching and the dynamical relative stability of multiple states.

Allosteric conversation in the androgen receptor ligand-binding domain surfaces

Androgen receptor (AR) is a major therapeutic target that plays pivotal roles in prostate cancer (PCa) and androgen insensitivity syndromes. We previously proposed that compounds recruited to ligand-binding domain (LBD) surfaces could regulate AR activity in hormone-refractory PCa and discovered several surface modulators of AR function. Surprisingly, the most effective compounds bound preferentially to a surface of unknown function [binding function 3 (BF-3)] instead of the coactivator-binding site [activation function 2 (AF-2)]. Different BF-3 mutations have been identified in PCa or androgen insensitivity syndrome patients, and they can strongly affect AR activity. Further, comparison of AR x-ray structures with and without bound ligands at BF-3 and AF-2 showed structural coupling between both pockets. Here, we combine experimental evidence and molecular dynamic simulations to investigate whether BF-3 mutations affect AR LBD function and dynamics possibly via allosteric conversation between surface sites. Our data indicate that AF-2 conformation is indeed closely coupled to BF-3 and provide mechanistic proof of their structural interconnection. BF-3 mutations may function as allosteric elicitors, probably shifting the AR LBD conformational ensemble toward conformations that alter AF-2 propensity to reorganize into subpockets that accommodate N-terminal domain and coactivator peptides. The induced conformation may result in either increased or decreased AR activity. Activating BF-3 mutations also favor the formation of another pocket (BF-4) in the vicinity of AF-2 and BF-3, which we also previously identified as a hot spot for a small compound. We discuss the possibility that BF-3 may be a protein-docking site that binds to the N-terminal domain and corepressors. AR surface sites are attractive pharmacological targets to develop allosteric modulators that might be alternative lead compounds for drug design.

Structure and dynamics of the membrane attaching nitric oxide transporter nitrophorin 7

Nitrophorins represent a unique class of heme proteins that are able to perform the delicate transportation and release of the free-radical gaseous messenger nitric oxide (NO) in a pH-triggered manner. Besides its ability to bind to phospholipid membranes, the N-terminus contains an additional Leu-Pro-Gly stretch, which is a unique sequence trait, and the heme cavity is significantly altered with respect to other nitrophorins. These distinctive features encouraged us to solve the X-ray crystallographic structures of NP7 at low and high pH and bound with different heme ligands (nitric oxide, histamine, imidazole). The overall fold of the lipocalin motif is well preserved in the different X-ray structures and resembles the fold of other nitrophorins. However, a chain-like arrangement in the crystal lattice due to a number of head-to-tail electrostatic stabilizing interactions is found in NP7. Furthermore, the X-ray structures also reveal ligand-dependent changes in the orientation of the heme, as well as in specific interactions between the A-B and G-H loops, which are considered to be relevant for the biological function of nitrophorins. Fast and ultrafast laser triggered ligand rebinding experiments demonstrate the pH-dependent ligand migration within the cavities and the exit route. Finally, the topological distribution of pockets located around the heme as well as from inner cavities present at the rear of the protein provides a distinctive feature in NP7, so that while a loop gated exit mechanism to the solvent has been proposed for most nitrophorins, a more complex mechanism that involves several interconnected gas hosting cavities is proposed for NP7.

Speed of travelling fronts: Two-dimensional and ballistic dispersal probability distributions

The speed of traveling fronts for a two-dimensional model of a delayed reactiondispersal process is derived analytically and from simulations of molecular dynamics. We show that the one-dimensional (1D) and two-dimensional (2D) versions of a given kernel do not yield always the same speed. It is also shown that the speeds of time-delayed fronts may be higher than those predicted by the corresponding non-delayed models. This result is shown for systems with peaked dispersal kernels which lead to ballistic transport

A systematic and feasible method for computing nuclear contributions to electrical properties of polyatomic molecules

An analytic method to evaluate nuclear contributions to electrical properties of polyatomic molecules is presented. Such contributions control changes induced by an electric field on equilibrium geometry (nuclear relaxation contribution) and vibrational motion (vibrational contribution) of a molecular system. Expressions to compute the nuclear contributions have been derived from a power series expansion of the potential energy. These contributions to the electrical properties are given in terms of energy derivatives with respect to normal coordinates, electric field intensity or both. Only one calculation of such derivatives at the field-free equilibrium geometry is required. To show the useful efficiency of the analytical evaluation of electrical properties (the so-called AEEP method), results for calculations on water and pyridine at the SCF/TZ2P and the MP2/TZ2P levels of theory are reported. The results obtained are compared with previous theoretical calculations and with experimental values

Treatment of nonlinear optical properties due to large amplitude anharmonic vibrational motions: umbrella motion in NH3

A general reduced dimensionality finite field nuclear relaxation method for calculating vibrational nonlinear optical properties of molecules with large contributions due to anharmonic motions is introduced. In an initial application to the umbrella (inversion) motion of NH3 it is found that difficulties associated with a conventional single well treatment are overcome and that the particular definition of the inversion coordinate is not important. Future applications are described

Additional compact formulas for vibrational dynamic dipole polarizabilities and hyperpolarizabilities

Compact expressions, complete through second order in electrical and/or mechanical anharmonicity, are given for the dynamic dipole vibrational polarizability and dynamic first and second vibrational hyperpolarizabilities. Certain contributions not previously formulated are now included

Critical thoughts on computing atom condensed Fukui functions

Different procedures to obtain atom condensed Fukui functions are described. It is shown how the resulting values may differ depending on the exact approach to atom condensed Fukui functions. The condensed Fukui function can be computed using either the fragment of molecular response approach or the response of molecular fragment approach. The two approaches are nonequivalent; only the latter approach corresponds in general with a population difference expression. The Mulliken approach does not depend on the approach taken but has some computational drawbacks. The different resulting expressions are tested for a wide set of molecules. In practice one must make seemingly arbitrary choices about how to compute condensed Fukui functions, which suggests questioning the role of these indicators in conceptual density-functional theory

Critical analysis and extension of the Hirshfeld atoms in molecules

The computational approach to the Hirshfeld [Theor. Chim. Acta 44, 129 (1977)] atom in a molecule is critically investigated, and several difficulties are highlighted. It is shown that these difficulties are mitigated by an alternative, iterative version, of the Hirshfeld partitioning procedure. The iterative scheme ensures that the Hirshfeld definition represents a mathematically proper information entropy, allows the Hirshfeld approach to be used for charged molecules, eliminates arbitrariness in the choice of the promolecule, and increases the magnitudes of the charges. The resulting "Hirshfeld-I charges" correlate well with electrostatic potential derived atomic charges

Selected configuration interaction with truncation energy error and application to the Ne atom

Selected configuration interaction (SCI) for atomic and molecular electronic structure calculations is reformulated in a general framework encompassing all CI methods. The linked cluster expansion is used as an intermediate device to approximate CI coefficients BK of disconnected configurations (those that can be expressed as products of combinations of singly and doubly excited ones) in terms of CI coefficients of lower-excited configurations where each K is a linear combination of configuration-state-functions (CSFs) over all degenerate elements of K. Disconnected configurations up to sextuply excited ones are selected by Brown's energy formula, ΔEK=(E-HKK)BK2/(1-BK2), with BK determined from coefficients of singly and doubly excited configurations. The truncation energy error from disconnected configurations, Δdis, is approximated by the sum of ΔEKS of all discarded Ks. The remaining (connected) configurations are selected by thresholds based on natural orbital concepts. Given a model CI space M, a usual upper bound ES is computed by CI in a selected space S, and EM=E S+ΔEdis+δE, where δE is a residual error which can be calculated by well-defined sensitivity analyses. An SCI calculation on Ne ground state featuring 1077 orbitals is presented. Convergence to within near spectroscopic accuracy (0.5 cm-1) is achieved in a model space M of 1.4× 109 CSFs (1.1 × 1012 determinants) containing up to quadruply excited CSFs. Accurate energy contributions of quintuples and sextuples in a model space of 6.5 × 1012 CSFs are obtained. The impact of SCI on various orbital methods is discussed. Since ΔEdis can readily be calculated for very large basis sets without the need of a CI calculation, it can be used to estimate the orbital basis incompleteness error. A method for precise and efficient evaluation of ES is taken up in a companion paper

Anharmonicity contributions to the vibrational second hyperpolarizability of conjugated oligomers

Restricted Hartree-Fock 6-31G calculations of electrical and mechanical anharmonicity contributions to the longitudinal vibrational second hyperpolarizability have been carried out for eight homologous series of conjugated oligomers - polyacetylene, polyyne, polydiacetylene, polybutatriene, polycumulene, polysilane, polymethineimine, and polypyrrole. To draw conclusions about the limiting infinite polymer behavior, chains containing up to 12 heavy atoms along the conjugated backbone were considered. In general, the vibrational hyperpolarizabilities are substantial in comparison with their static electronic counterparts for the dc-Kerr and degenerate four-wave mixing processes (as well as for static fields) but not for electric field-induced second harmonic generation or third harmonic generation. Anharmonicity terms due to nuclear relaxation are important for the dc-Kerr effect (and for the static hyperpolarizability) in the σ-conjugated polymer, polysilane, as well as the nonplanar π systems polymethineimine and polypyrrole. Restricting polypyrrole to be planar, as it is in the crystal phase, causes these anharmonic terms to become negligible. When the same restriction is applied to polymethineimine the effect is reduced but remains quantitatively significant due to the first-order contribution. We conclude that anharmonicity associated with nuclear relaxation can be ignored, for semiquantitative purposes, in planar π-conjugated polymers. The role of zero-point vibrational averaging remains to be evaluated

Linear response functions for a vibrational configuration interaction state

Linear response functions are implemented for a vibrational configuration interaction state allowing accurate analytical calculations of pure vibrational contributions to dynamical polarizabilities. Sample calculations are presented for the pure vibrational contributions to the polarizabilities of water and formaldehyde. We discuss the convergence of the results with respect to various details of the vibrational wave function description as well as the potential and property surfaces. We also analyze the frequency dependence of the linear response function and the effect of accounting phenomenologically for the finite lifetime of the excited vibrational states. Finally, we compare the analytical response approach to a sum-over-states approach