Folding of a stable DNA motif involves a highly cooperative network of interactions

Hairpins are structural elements that play important roles in the folding and function of RNA and DNA. The extent of cooperativity in folding is an important aspect of the RNA folding problem. We reasoned that an investigation into the origin of cooperativity might be best carried out on a stable nucleic acid system with a limited number of interactions, such as a stable DNA hairpin loop. The stable d(cGNAg) hairpin loop motif (closing base pair in lower case; loop in upper case; N = A, C, G, or T) is stabilized through only three interactions: two loop-loop hydrogen bonds in a sheared GA base pair and a loop-closing base pair interaction. Herein, we investigate this network of interactions and test whether the loop-loop and loop-closing base pair interactions communicate. Thermodynamic measurements of nucleotide analogue substituted oligonucleotides were used to probe the additivity of the interactions. On the basis of double mutant cycles, all interactions were found to be nonadditive and interdependent, suggesting that loop-loop and loop-closing base pair interactions form in a highly cooperative manner. When double mutant cycles were repeated in the absence of the other interaction, nonadditivity was significantly reduced suggesting that coupling is indirect and requires all three interactions in order to be optimal. A cooperative network of interactions helps explain the structural and energetic bases of stability in certain DNA hairpins and paves the way for similar studies in more complex nucleic acid systems.     

Cation-pi interactions of a thiocarbonyl group and a carbonyl group with a pyridinium nucleus

Attractive interactions between a thiocarbonyl group and a pyridinium nucleus, and between a carbonyl group and a pyridinium nucleus have been proven by (1)H and (13)C NMR studies, UV-vis spectral analyses, and X-ray crystallographic analyses of nicotinic amides 1 and 3, and pyridinium salts 2 and 4. Comparison of the Deltadelta values, which are the differences in the chemical shifts with reference compounds 5 or 6, showed that the absolute Deltadelta values of 2 and 4 are much larger than those of 1 and 3. In the UV-vis spectra, the n-->pi absorption of the C=S group of 2a exhibited a significant blue shift in CHCl(3). X-ray crystallographic analysis of 1-4 clearly showed that the C=S group of 2a and the C=O group of 4 are very close to the pyridinium moiety compared to the case of 1 and 3. In addition, the X-ray crystal packing structure of 2a showed the C=S group is sandwiched between two pyridinium rings. These experimental results strongly suggested the existence of attractive (C=S)...Py(+) and (C=O)...Py(+) interactions in solution and in crystal. The optimized geometries of 1 and 2 calculated at the HF/6-311G level are in good agreement with their X-ray geometries. MP2/6-311G calculations for the model systems of pyridinium salts 2 and 4 predicted that the electrostatic and induction energies are the major source of the attractive interactions. Since the larger contribution of electrostatic and induction interactions are characteristic features of cation-pi interactions, the (C=S)...Py(+) and (C=O)...Py(+) interactions would be classified as a cation-pi interaction.     

Weak Intermolecular Interactions in a Series of Bioactive Oxazoles

The intermolecular interactions in a series of nine similar 4,5-phenyl-oxazoles were studied on the basis of crystal structures determined by X-ray diffraction. The crystal architectures were analyzed for the importance and hierarchies of different, weak intermolecular interactions using three approaches: the geometrical characteristics, topological analysis (for the model based on the transfer of multipolar parameters), and energetics of the molecule–molecule interactions. The geometries of the molecules were quite similar and close to the typical values. The results of the analysis of the interactions suggest that the number of nonspecific interactions is more important than the apparent strength of the specific interactions. The interactions involving covalently bound bromine and divalent sulfur atoms were classified as secondary, they certainly did not define the crystal packing, and they played a minor role in the overall crystal cohesion energies. Incidentally, another method for confirming the degree of isostructurality, according to the topologies of the interactions, is described.     

Interactions underpinning the plasticization of a polymer matrix: a dynamic and structural analysis of DMP-plasticized cellulose acetate

This work establishes that the plasticization effect of a classical petrochemical plasticizer, dimethyl phthalate (DMP), on a polymer matrix, cellulose acetate (CA), is due to the development of intermolecular interactions of dipolar type. Plasticized cellulose acetate films are studied with regard to the interactions between the polymer and plasticizer at the macroscopic scale by thermogravimetric analysis and differential scanning calorimetry. At the molecular level, Fourier transform infrared spectroscopy and dielectric relaxation spectroscopy are used to elucidate the nature of interactions that are responsible for the plasticizing effects. These static and dynamic complementary analyses evidenced that DMP does not establish H-bonding interactions with the polymer chains of cellulose acetate but rather weaker interactions of dipolar type. These dipole–dipole interactions that develop between acetyl side groups of CA and the ester phthalate moieties of DMP increase the overall mobility of CA chains and also locally influence the molecular mobility and the water uptake tendency.     

Charged particles at fluid interfaces as a probe into structural details of a double layer

Electrostatic interactions between charged, distant colloids in a bulk electrolyte solution do not depend on the inherent structure of ions and a solvent forming a double layer. For charged colloids trapped at an interface between an electrolyte and air this no longer holds; as the electrostatic interactions are mediated via air and the field lines determining the interactions originate at the charged surface, these details come into prominence. Using the Langevin-Poisson-Boltzmann equation we investigate how steric effects and the polarization saturation of a solvent effect the contact potential at the colloid surface and, in consequence, the long range interactions between colloids trapped at an interface. For a surface charge 0.4 C m(-2) the combination of these effects can increase the interactions by up to ～40 times when compared to Poisson-Boltzmann calculations. The validity of these enhancement mechanisms is supported by recent experimental data (K. Masschaele et al., Phys. Rev. Lett., 2010, 105, 048303).     

A 10,000-fold nuclear hyperpolarization of a membrane protein in the liquid phase via a solid-state mechanism

Several techniques rely on electron-nuclear interactions to boost the polarization of nuclear spins in the solid phase. Averaging out of anisotropic interactions as a result of molecular tumbling strongly reduces the applicability of such hyperpolarization approaches in liquids. Here we show for the first time that anisotropic electron-nuclear interactions in solution can survive sufficiently long to generate nuclear spin polarization by the solid-state photo-CIDNP mechanism. A 10,000-fold NMR signal increase in solution was observed for a giant biomolecular complex of a photosynthetic membrane protein with a tumbling correlation time in the submicrosecond regime, corresponding to a molecular weight close to 1 MDa.     

Interactions of polarizable media in water: a molecular dynamics approach

We investigate the interactions of polarizable solutes in water as a function of the solute permittivity. A generic and computationally efficient simulation methodology for the investigation of systems involving dielectric discontinuities is introduced. We report results for interactions between two polarizable cylindrical solutes of nanometer dimensions, which demonstrate that the interactions between the solutes strongly depend on the solute permittivity epsilon. For low permittivity, epsilon approximately 1-2, the interactions are dominated by surface tension forces whose origin lies in the formation of a vapor cavity between the two hydrophobic solutes. This effect leads to a drying transition, where the intersolute force changes discontinuously at a specific solute-solute separation. We find that a moderate permittivity, epsilon approximately 20, enhances the solvation of the polarizable objects inhibiting this drying transition. In the limit of moderately high permittivity, the interactions are dominated by solvation forces. These forces are much larger than those calculated using macroscopic models of dielectrics, which consider water as a continuum dielectric medium. Our results emphasize the importance of including the solvent explicitly to investigate dielectric discontinuities and interactions between polarizable media in water.     

Computational design of heterochiral peptides against a helical target

Polypeptides incorporating D-amino acids occasionally occur in nature and are an important class of pharmaceutical molecules. With the use of heterochiral Monte Carlo (HCMC), a method inspired by the de novo design of proteins, we develop peptide scaffolds for interacting with a molecular target, a left-handed alpha-helix. The HCMC approach concurrently seeks to optimize a peptide sequence, its internal conformation, and its docked conformation with a target surface. Several major classes of interactions are observed: (1) homochiral interactions between two alphaL helices, (2) heterochiral interactions between an alphaL and an alphaR helix, and (3) heterochiral interactions between the alphaL target and novel nonhelical structures. We explore the application of HCMC to simulating the preferential enantioselectivity of heterochiral complexes. Implications for biomimetic design in molecular recognition are discussed.     

The Permeation Mechanism of Cisplatin Through a Dioleoylphosphocholine Bilayer**

The investigation of the intermolecular interactions between platinum-based anticancer drugs and lipid bilayers is of special relevance to unveil the mechanisms involved in different steps of the anticancer mode of action of these drugs. We have simulated the permeation of cisplatin through a model membrane composed of 1,2-dioleoyl-sn-glycero-3-phosphocholine lipids by means of umbrella sampling classical molecular dynamics simulations. The initial physisorption of cisplatin into the polar region of the lipid membrane is controlled by long-range electrostatic interactions with the choline groups in a first step and, in a second step, by long-range electrostatic and hydrogen bond interactions with the phosphate groups. The second half of the permeation pathway, in which cisplatin diffuses through the nonpolar region of the bilayer, is characterized by the drop of the interactions with the polar heads and the rise of attractive interactions with the non-polar tails, which are dominated by van der Waals contributions. The permeation free-energy profile is explained by a complex balance between the drug/lipid interactions and the energy and entropy contributions associated with the dehydration of the drug along the permeation pathway and with the decrease and increase of the membrane ordering along the first and second half of the mechanism, respectively.     

Helical folding of alkanes in a self-assembled, cylindrical capsule

The reversible encapsulation of a series of normal alkane guests in a cylindrical host was studied by NMR methods. For small hydrocarbons such as n-pentane or n-hexane, two guests enter the host, and they move freely within. With n-heptane no encapsulation takes place. For longer alkanes such as n-decane, a single guest enters and the aromatic walls of the host are seen to twist to avoid empty spaces and increase favorable interactions with the hydrocarbon. The best guest (n-undecane) adopts a conformation with minimal gauche interactions. The longest alkane accommodated, n-tetradecane, adopts a helical conformation to fit in the cavity, a shape that maximizes CH/pi interactions with the aromatic walls of the receptor. These reciprocal conformational changes are discussed in terms of optimal host/guest interactions.     

Control of the stability of a protein-RNA complex by the position of fluorine in a base analogue

The effects of modifying the electronic characteristics of nonpolar base analogues substituted at positions involved in stacking interactions between SL2 RNA and the U1A protein are described. A surprisingly large difference in the stability between complexes formed with base analogues that differ only in the position of substitution of a single fluorine atom is observed. The results of high-level ab initio calculations of the interactions between the nonpolar base analogue and the amino acid side chain correlate with the experimentally observed trends in complex stability, which suggests that changes in stacking interactions that result from varying the position and degree of fluorine substitution contribute to the effects of fluorine substitution on the stability of the U1A-SL2 RNA complex.     

Efficient electron transfer in a protein network lacking specific interactions

In many biochemical processes, proteins need to bind partners amidst a sea of other molecules. Generally, partner selection is achieved by formation of a single-orientation complex with well-defined, short-range interactions. We describe a protein network that functions effectively in a metabolic electron transfer process but lacks such specific interactions. The soil bacterium Paracoccus denitrificans oxidizes a variety of compounds by channeling electrons into the main respiratory pathway. Upon conversion of methylamine by methylamine dehydrogenase, electrons are transported to the terminal oxidase to reduce molecular oxygen. Steady-state kinetic measurements and NMR experiments demonstrate a remarkable number of possibilities for the electron transfer, involving the cupredoxin amicyanin as well as four c-type cytochromes. The observed interactions appear to be governed exclusively by the electrostatic nature of each of the proteins. It is concluded that Paracoccus provides a pool of cytochromes for efficient electron transfer via weak, ill-defined interactions, in contrast with the view that functional biochemical interactions require well-defined molecular interactions. It is proposed that the lack of requirement for specificity in these interactions might facilitate the integration of new metabolic pathways.     

The dynamic behavior of a liquid ethanol-water mixture: a perspective from quantum chemical topology

Quantum Chemical Topology (QCT) is used to reveal the dynamics of atom-atom interactions in a liquid. A molecular dynamics simulation was carried out on an ethanol-water liquid mixture at its azeotropic concentration (X(ethanol)=0.899), using high-rank multipolar electrostatics. A thousand (ethanol)(9)-water heterodecamers, respecting the water-ethanol ratio of the azeotropic mixture, were extracted from the simulation. Ab initio electron densities were computed at the B3LYP/6-31+G(d) level for these molecular clusters. A video shows the dynamical behavior of a pattern of bond critical points and atomic interaction lines, fluctuating over 1 ns. A bond critical point distribution revealed the fluctuating behavior of water and ethanol molecules in terms of O-H···O, C-H···O and H···H interactions. Interestingly, the water molecule formed one to six C-H···O and one to four O-H···O interactions as a proton acceptor. We found that the more localized a dynamical bond critical point distribution, the higher the average electron density at its bond critical points. The formation of multiple C-H···O interactions affected the shape of the oxygen basin of the water molecule, which is shown in three dimensions. The hydrogen atoms of water strongly preferred to form H···H interactions with ethanol's alkyl hydrogen atoms over its hydroxyl hydrogen.     

Simple cation-pi interaction between a phenyl ring and a protonated amine stabilizes an alpha-helix in water

Cation-pi interactions have been proposed to be important contributors to protein structure and function. In particular, these interactions have been suggested to provide significant stability at the solvent-exposed surface of a protein. We have investigated the magnitude of cation-pi interactions between phenylalanine (Phe) and lysine (Lys), ornithine (Orn), and diaminobutanoic acid (Dab) in the context of an alpha-helix and have found that only the Phe...Orn interaction provides significant stability to the helix, stabilizing it by -0.4 kcal/mol. This interaction energy is in the same range as a salt bridge in an alpha-helix, and equivalent to the recently reported Trp...Arg interaction in an alpha-helix, despite the fact that Trp...guanidinium interactions have been proposed to be stronger than Phe...ammonium interactions. These results indicate that even the simplest cation-pi interaction can provide significant stability to protein structure and demonstrate the subtle factors that can influence the observed interaction energies in designed systems.     

Molecular mechanics and the deformation of macromolecules: The use of a very short cutoff combined with a quadratic approximation

An approximation to the molecular mechanical treatment of structural deformations of macromolecules is presented. The method is based on a partitioning of the conformational energy into three parts. The first part is covered by the condensed potential functions which absorb many short-range nonbonding interactions. The second part consists of a few nonbonded interactions below a very short cutoff radius of 4 Å. The third part, consisting of the vast majority of pairwise interactions, is approximated by a quadratic expression confined to a subspace of the conformational space. A detailed computational example on LH-RH, including an analysis of the errors resulting from other conventional approximation methods, is given. A comparison to the conventional cutoff approximation used in x-ray refinement delivers a speedup factor of at least two orders of magnitude.     

Density profiles and solvation forces for a Yukawa fluid in a slit pore

The effect of varying wall-particle and particle-particle interactions on the density profiles near a single wall and the solvation forces between two walls immersed in a fluid of particles is investigated by grand canonical Monte Carlo simulations. Attractive and repulsive particle-particle and particle-wall interactions are modeled by a versatile hard-core Yukawa form. These simulation results are compared to theoretical calculations using the hypernetted chain integral equation technique, as well as with fundamental measure density functional theory (DFT), where particle-particle interactions are either treated as a first order perturbation using the radial distribution function or else with a DFT based on the direct-correlation function. All three theoretical approaches reproduce the main trends fairly well, but exhibit inconsistent accuracy, particularly for attractive particle-particle interactions. We show that the wall-particle and particle-particle attractions can couple together to induce a nonlinear enhancement of the adsorption and a related "repulsion through attraction" effect for the effective wall-wall forces. We also investigate the phenomenon of bridging, where an attractive wall-particle interaction induces strongly attractive solvation forces.     

Model for reversible nanoparticle assembly in a polymer matrix

The clustering of nanoparticles (NPs) in solutions and polymer melts depends sensitively on the strength and directionality of the NP interactions involved, as well as the molecular geometry and interactions of the dispersing fluids. Since clustering can strongly influence the properties of polymer-NP materials, we aim to better elucidate the mechanism of reversible self-assembly of highly symmetric NPs into clusters under equilibrium conditions. Our results are based on molecular dynamics simulations of icosahedral NP with a long-ranged interaction intended to mimic the polymer-mediated interactions of a polymer-melt matrix. To distinguish effects of polymer-mediated interactions from bare NP interactions, we compare the NP assembly in our coarse-grained model to the case where the NP interactions are purely short ranged. For the "control" case of NPs with short-ranged interactions and no polymer matrix, we find that the particles exhibit ordinary phase separation. By incorporating physically plausible long-ranged interactions, we suppress phase separation and qualitatively reproduce the thermally reversible cluster formation found previously in computations for NPs with short-ranged interactions in an explicit polymer-melt matrix. We further characterize the assembly process by evaluating the cluster properties and the location of the self-assembly transition. Our findings are consistent with a theoretical model for equilibrium clustering when the particle association is subject to a constraint. In particular, the density dependence of the average cluster mass exhibits a linear concentration dependence, in contrast to the square root dependence found in freely associating systems. The coarse-grained model we use to simulate NP in a polymer matrix shares many features of potentials used to model colloidal systems. The model should be practically valuable for exploring factors that control the dispersion of NP in polymer matrices where explicit simulation of the polymer matrix is too time consuming.     

The adsorption of polar molecules on the surface of graphite modified with a monomolecular layer of a mesogen

An equation for the potential energy of the interaction of a quasi-solid polar adsorbate molecule with a semiinfinite graphite crystal containing a monomolecular layer of polar anisometric particles on its surface was obtained. Adsorbate-modifier electrostatic interactions were considered in the dipole-dipole approximation. Dispersion interactions were described using the Lennard-Jones potentials (6, 12). The suggestion was made that a modifier monolayer impermeable to adsorbate molecules could be modeled by force centers that formed a two-dimensional translationally periodic lattice and by force centers uniformly distributed in a plane parallel to the outer face of graphite crystals.     

Equilibrium sampling of self-associating polymer solutions: a parallel selective tempering approach

We present a novel simulation algorithm based on tempering a fraction of relaxation-limiting interactions to accelerate the process of obtaining uncorrelated equilibrium configurations of self-associating polymer solutions. This approach consists of tempering (turning off) the attractive interactions for a fraction of self-associating groups determined by a biasing field h. A number of independent configurations (replicas) with overlapping Hamiltonian distributions in the expanded (NVTh) ensemble with constant NVT but different biasing fields, forming a chain of Hamiltonians, were simulated in parallel with occasional attempts to exchange the replicas associated with adjacent fields. Each field had an associated distribution of tempered interactions, average fraction of tempered interactions, and structural decorrelation time. Tempering parameters (number of replicas, fields, and exchange frequencies) were chosen to obtain the highest efficiency in sampling equilibrium configurations of a self-association polymer solution based on short serial simulation runs and a statistical model. Depending on the strength of the relaxation-limiting interactions, system size, and thermodynamic conditions, the algorithm can be orders of magnitude more efficient than conventional canonical simulation and is superior to conventional temperature parallel tempering.     

Intramolecular rearrangement for regioselective complexation by intramolecular CH/pi interaction in a hydrophobic cavity of a ruthenium coordination sphere

A Ru(II) complex with a hydrophobic cavity formed from two 1-naphthoylamide groups was prepared. Its reactions with beta-diketones gave beta-diketonato complexes in which hydrophobic pi-pi or CH/pi interactions were confirmed by NMR spectroscopy and X-ray crystallography. In the case of the asymmetric beta-diketone benzoylacetone, an isomer with a CH/pi interaction was afforded as the sole product owing to thermodynamic control. The reaction was found to involve a novel intramolecular rearrangement from the phenyl-included isomer to the methyl-included one without rupture of Ru-beta-diketonato coordination bonds (activation energy 52 kJ mol(-1)). This indicates that CH/pi interactions can be more favored thermodynamically than pi-pi interactions in a suitable hydrophobic environment.