Solvation dynamics in a Brownian dipolar lattice. Comparison between computer simulation and various molecular theories of solvation dynamics

Several recent theoretical and computer simulation studies have considered solvation dynamics in a Brownian dipolar lattice which provides a simple model solvent for which detailed calculations can be carried out. In this article a fully microscopic calculation of the solvation dynamics of an ion in a Brownian dipolar lattice is presented. The calculation is based on the non‐Markovian molecular hydrodynamic theory developed recently. The main assumption of the present calculation is that the two‐particle orientational correlation functions of the solid can be replaced by those of the liquid state. It is shown that such a calculation provides an excellent agreement with the computer simulation results. More importantly, the present calculations clearly demonstrate that the frequency‐dependent dielectric friction plays an important role in the long time decay of the solvation time correlation function. We also find that the present calculation provides somewhat better agreement than either the dynamic mean spherical approximation (DMSA) or the Fried–Mukamel theory which use the simulated frequency‐dependent dielectric function. It is found that the dissipative kernels used in the molecular hydrodynamic approach and in the Fried–Mukamel theory are vastly different, especially at short times. However, in spite of this disagreement, the two theories still lead to comparable results in good agreement with computer simulation, which suggests that even a semiquantitatively accurate dissipative kernel may be sufficient to obtain a reliable solvation time correlation function. A new wave vector and frequency‐dependent dissipative kernel (or memory function) is proposed which correctly goes over to the appropriate expressions in both the single particle and the collective limits. This form is expected to lead to better results than all the existing descriptions.

Dynamic coupling between the LID and NMP domain motions in the catalytic conversion of ATP and AMP to ADP by adenylate kinase

The catalytic conversion of adenosine triphosphate (ATP) and adenosine monophosphate (AMP) to adenosine diphosphate (ADP) by adenylate kinase (ADK) involves large amplitude, ligand induced domain motions, involving the opening and the closing of ATP binding domain (LID) and AMP binding domain (NMP) domains, during the repeated catalytic cycle. We discover and analyze an interesting dynamical coupling between the motion of the two domains during the opening, using large scale atomistic molecular dynamics trajectory analysis, covariance analysis, and multidimensional free energy calculations with explicit water. Initially, the LID domain must open by a certain amount before the NMP domain can begin to open. Dynamical correlation map shows interesting cross-peak between LID and NMP domain which suggests the presence of correlated motion between them. This is also reflected in our calculated two-dimensional free energy surface contour diagram which has an interesting elliptic shape, revealing a strong correlation between the opening of the LID domain and that of the NMP domain. Our free energy surface of the LID domain motion is rugged due to interaction with water and the signature of ruggedness is evident in the observed root mean square deviation variation and its fluctuation time correlation functions. We develop a correlated dynamical disorder-type theoretical model to explain the observed dynamic coupling between the motion of the two domains in ADK. Our model correctly reproduces several features of the cross-correlation observed in simulations. (C) 2011 American Institute of Physics. doi:10.1063/1.3516588]

Importance of orientational rearrangement during vitrification of hydrocarbons: dependence on molecular shape

On the basis of Monte Carlo calculations of 2,2-dimethylpropane (neopentane), n-pentane, and 2,2-dimethylbutane (neohexane) at several temperatures, thermodynamic properties and radial distribution functions as well as dimerization and bonding energy distribution functions are reported for both liquid and glassy states. Changes in the radial distribution functions on cooling depend on whether the groups are accessible (peripheral) or inaccessible. Peaks in the radial distribution functions corresponding to peripheral groups do not shift to lower distances on cooling and at times display a large increase in the intensity of the first peak. The peaks due to inaccessible groups, on the other hand, shift to lower distances on cooling. The magnitude of the reorientational contribution in determining the resulting structure of the glass is estimated for the different hydrocarbon molecules investigated. The reorientational contribution is highest for neopentane (26%) followed by isopentane (24%), neohexane (22%), and n-pentane (0%). It appears that molecular geometry has an important role in determining the magnitude of the reorientational contribution to the structure of the glass.

Using Correlated Parameters for Improved Ranking of Protein-Protein Docking Decoys

A successful protein-protein docking study culminates in identification of decoys at top ranks with near-native quaternary structures. However, this task remains enigmatic because no generalized scoring functions exist that effectively infer decoys according to the similarity to near-native quaternary structures. Difficulties arise because of the highly irregular nature of the protein surface and the significant variation of the nonbonding and solvation energies based on the chemical composition of the protein-protein interface. In this work, we describe a novel method combining an interface-size filter, a regression model for geometric compatibility (based on two correlated surface and packing parameters), and normalized interaction energy (calculated from correlated nonbonded and solvation energies), to effectively rank decoys from a set of 10,000 decoys. Tests on 30 unbound binary protein-protein complexes show that in 16 cases we can identify at least one decoy in top three ranks having <= 10 angstrom backbone root mean square deviation from true binding geometry. Comparisons with other state-of-art methods confirm the improved ranking power of our method without the use of any experiment-guided restraints, evolutionary information, statistical propensities, or modified interaction energy equations. Tests on 118 less-difficult bound binary protein-protein complexes with <= 35% sequence redundancy at the interface showed that in 77% cases, at least 1 in 10,000 decoys were identified with <= 5 angstrom backbone root mean square deviation from true geometry at first rank. The work will promote the use of new concepts where correlations among parameters provide more robust scoring models. It will facilitate studies involving molecular interactions, including modeling of large macromolecular assemblies and protein structure prediction. (C) 2010 Wiley Periodicals, Inc. J Comput Chem 32: 787-796, 2011.

Sorbate properties of xenon in cloverite: A molecular dynamics study

Molecular dynamics calculations on xenon adsorbed in the cubic cavity of cloverite, a gallophosphate, is presented. Guest-host radial distribution functions, guest-host energy distribution functions, power spectra, and diffusion coefficients for xenon are reported at 397, 494, and 716 K. Results suggest that xenon is highly mobile at 700 K. A shift in the peak position of the power spectra toward lower frequency is observed with increase in temperature. (C) 1994 Academic Press, Inc.

Correlated electronic states of C-60 fragments

We have studied electronic states of various fragments of C-60 within the Pariser-Parr-Pople (PPP) model and have obtained structural, magnetic and spectral properties of these molecules. The fragments studied include corannulene, fluoranthene and pyracylene. Pyracylene is studied using the exact valence bond (VB) approach while fluoranthene and corannulene are studied using a novel restricted CI technique which employs molecular orbitals for constructing the VB functions. The electronic excitations, bond order and ring currents are calculated for these systems. From these studies, the wide range of absorptions in C-60 can be viewed as those localized on pyracylene units or on the corannulene/fluoranthene units. The bond orders and ring currents show the hexagons to be similar to benzene rings. The pentagon-hexagon bonds are also found to be longer than the hexagon-hexagon bonds.

Electron donor-acceptor complex of ICl with diethyl ether - He I photoelectron spectroscopy and ab initio molecular orbital study

The He I photoelectron spectrum of the diethyl ether-ICl complex has been obtained. The oxygen orbitals are shifted to higher binding energies and that of ICl to lower binding energies owing to complex formation. Ab initio molecular orbital (MO) calculations of the complex molecule showed that the bonding is between the sigma-type lone pair of oxygen and the I atom and that the complex has C-2v symmetry. The binding energy of the complex is computed to be 8.06 kcal mol(-1) at the MP2/3-21G* level. The orbital energies obtained from the photoelectron spectra of the complex are compared and assigned with orbital energies obtained by MO calculations. Natural bond orbital analysis (NBO) shows that charge transfer is from the sigma-type oxygen lone pair to the iodine atom and the magnitude of charge transfer is 0.0744 e.

Solid-state photobehaviour and crystal packing of o-chlorobenzylidene-DL-piperitone: influence of molecular topology on photobehaviour

C17H19ClO, M(r) = 274.7, triclinic, P1BAR, a = 11.154 (3), b = 12.685 (2), c = 12.713 (2) angstrom, alpha = 100.68 (1), beta = 113.58 (1), gamma = 104.50 (2)-degrees, V = 1511.1 (6) angstrom3, Z = 4, D(m) = 1.22, D(x) = 1.215 Mg m-3, Cu K-alpha, lambda = 1.5418 angstrom, mu = 2.16 mm-1, F(000) = 584, T = 293 K, R = 0.057 for 3481 observed reflections. The title compound is photostable in the crystalline state and lattice-energy calculations have been employed to rationalize the photobehaviour. The well-known beta-steering ability of the chloro group is not operative in this system as there are no Cl...Cl interactions in the crystal lattice. All five benzylidene-DL-piperitone structures so far studied are alpha-packed and the molecular topology appears to be a deciding factor even in the presence of steering groups.

Synthesis and characterization of photo-crosslinkable main-chain liquid-crystalline polymers containing bis(benzylidene)cycloalkanone units

A series of new photo-crosslinkable main-chain liquid-crystalline polymers containing bis(benzylidene)cycloalkanone units have been studied. These units in the polymers function as mesogens as well as photoactive centres. Polyesters with three different bis(4-hydroxybenzylidene)cycloalkanones corresponding to three cycloalkanones, namely cyclopentanone, cyclohexanone and cycloheptanone, have been prepared. Three dicarboxylic acids with ether linkages, which were derived from oligoethylene oxides, namely triethylene glycol, tetraethylene glycol and pentaethylene glycol, have been used as spacers in these polymers. Polymerization was carried out by both solution and interfacial polycondensation; the latter method gave high-molecular-weight polymers. Structural characterizations were done by ultra-violet, infra-red and H-1 nuclear magnetic resonance spectroscopy. Liquid-crystalline properties were studied by differential scanning calorimetry and polarized-light optical microscopy. These polymers show a nematic mesophase. Liquid-crystalline transition temperatures were correlated with polymer structure. The decrease in transition temperature with increase in cycloalkanone ring size was explained in terms of the change in geometrical anisotropy of bis(benzylidene)cycloalkanone units. MNDO (modified neglect of differential overlap) calculations were performed on the model compounds, bis(4-acetyloxybenzylidene)cycloalkanone to elucidate the geometrical variation of the mesogenic units with cycloalkanone ring size. Studies of photolysis reveal the two kinds of photoreactions that proceed in these polymer systems, namely photoisomerization and photo-crosslinking. The former reaction disrupts the parallel stacking of the chromophores and is reflected as an increase in the ultra-violet spectral intensity. The favourability of these two reactions depends on the mobility of the polymer chains. When the photolysis was done below T-g, photo-crosslinking dominates over photoisomerization. Above T-g, photoisomerization is followed by photo-crosslinking. The photosensitivity of the polymers decreases with increase in size of the cycloalkanone ring.

UVPES and ab initio molecular orbital studies on the electron donor-acceptor complexes of bromine with methylamines

The He I photoelectron spectra of bromine, methylamine, and their complex have been obtained, and the spectra show that lone-pair orbital energy of nitrogen in methylamine is stabilized by 1.8 eV and the bromine orbital energies are destabilized by about 0.5 eV due to complexation. Ab initio calculations have been performed on the charge-transfer complexes of Br-2 with ammonia and methyl-, dimethyl-, and trimethylamines at the 3-21G*, 6-311G, and 6-311G* levels and also with effective core potentials. Calculations predict donor and acceptor orbital energy shifts upon complexation, and there is a reasonable agreement between the calculated and experimental results. Complexation energies have been corrected for BSSE. Frequency analysis has confirmed that ammonia and trimethylamine form complexes with C-3v symmetry and methylamine and dimethylamine with C-s symmetry. Calculations reveal that the lone-pair orbital of nitrogen in amine and the sigma* orbital of Br-2 are involved in the charge-transfer interaction. LANL1DZ basis seems to be consistent and give a reliable estimate of the complexation energy. The computed complexation energies, orbital energy shifts, and natural bond orbital analysis show that the strength of the complex gradually increases from ammonia to trimethylamine.

Molecular structure and reactions of 1,1 '-bi(acenaphthen-1-ylidene)-2,2 '-dione

The molecular structure of 1,1'-bi(acenaphthen-1-ylidene)-2,2'-dione 1, a potential building-block for the synthesis of fullerene fragments, has been investigated by X-ray crystallography and semi-empirical (AM1 and PM3) calculations. There is a good agreement between the calculated and crystal structure which is essentially planar and has E-configuration. In the solid state, molecules of 1 pack in an interesting manner as corrugated sheets sustained by a network of C-H ... O hydrogen bonds and resulting in the formation of tetrameric loops. While steric factors limit the reactivity of the carbonyl groups in 1, the ene double bond of the ene-dione moiety present in it exhibits propensity toward [4 + 2]-cycloadditions to furnish novel and highly compressed polycycles 8-10.

The molecular dynamics study of lithium ion conduction in phosphate glasses and the role of non-bridging oxygen

Molecular dynamics (MD) simulation of lithium phosphate (Li2O-P2O5) glasses with varying Li2O content has been carried out. Two different P-O distances corresponding to phosphorus coordination with bridging oxygen (BO) and non-bridging oxygen (NBO) were identified in the simulated glasses. NBO-BO interconversion or bond switching was noted, which results in a dynamic equilibration of the tetrahedral phosphate units (P-n, n = 1,3 indicates the number of bridging oxygen atoms in the coordination of phosphorus). The NBO-BO bond switching is mildly activated with an effective activation barrier of 0.03-0.05 eV. Lithium ion jumps do not appear to be strongly coupled to bond switching. But the number of Li+ ions coordinated to an optimum number of NBOs and the number of Li+ ions jumping out of their sites appear to be correlated. Detailed analysis was made of the dynamics of P-n species and new insights have been obtained regarding ion migration in network-modified phosphate glasses.

Quaternary association in beta-prism I fold plant lectins: Insights from X-ray crystallography, modelling and molecular dynamics

Dimeric banana lectin and calsepa, tetrameric artocarpin and octameric heltuba are mannose-specific beta-prism I fold lectins of nearly the same tertiary structure. MD simulations on individual subunits and the oligomers provide insights into the changes in the structure brought about in the protomers on oligomerization, including swapping of the N-terminal stretch in one instance. The regions that undergo changes also tend to exhibit dynamic flexibility during MD simulations. The internal symmetries of individual oligomers are substantially retained during the calculations. Energy minimization and simulations were also carried out on models using all possible oligomers by employing the four different protomers. The unique dimerization pattern observed in calsepa could be traced to unique substitutions in a peptide stretch involved in dimerization. The impossibility of a specific mode of oligomerization involving a particular protomer is often expressed in terms of unacceptable steric contacts or dissociation of the oligomer during simulations. The calculations also led to a rationale for the observation of a heltuba tetramer in solution although the lectin exists as an octamer in the crystal, in addition to providing insights into relations among evolution, oligomerization and ligand binding.

Dependence of diffusivity on density and solute diameter in liquid phase: A molecular dynamics study of Lennard-Jones system

Investigations into the variation of self-diffusivity with solute radius, density, and degree of disorder of the host medium is explored. The system consists of a binary mixture of a relatively smaller sized solute, whose size is varied and a larger sized solvent interacting via Lennard-Jones potential. Calculations have been performed at three different reduced densities of 0.7, 0.8, and 0.933. These simulations show that diffusivity exhibits a maximum for some intermediate size of the solute when the solute diameter is varied. The maximum is found at the same size of the solute at all densities which is at variance with the prediction of the levitation effect. In order to understand this anomaly, additional simulations were carried out in which the degree of disorder has been varied while keeping the density constant. The results show that the diffusivity maximum gradually disappears with increase in disorder. Disorder has been characterized by means of the minimal spanning tree. Simulations have also been carried out in which the degree of disorder is constant and only the density is altered. The results from these simulations show that the maximum in diffusivity now shifts to larger distances with decrease in density. This is in agreement with the changes in void and neck distribution with density of the host medium. These results are in excellent agreement with the predictions of the levitation effect. They suggest that the effect of disorder is to shift the maximum in diffusivity towards smaller solute radius while that of the decrease in density is to shift it towards larger solute radius. Thus, in real systems where the degree of disorder is lower at higher density and vice versa, the effect due to density and disorder have opposing influences. These are confirmed by the changes seen in the velocity autocorrelation function, self part of the intermediate scattering function and activation energy. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.3701619]

Probing acid-amide intermolecular hydrogen bonding by NMR spectroscopy and DFT calculations

Benzene carboxylic acids and Benzamide act as their self-complement in molecular recognition to form inter-molecular hydrogen bonded dimers between amide and carboxylic acid groups, which have been investigated by H-1, C-13 and N-15 NMR spectroscopy. Extensive NMR studies using diffusion ordered spectroscopy (DOSY), variable temperature 1D, 2D NMR, established the formation of heterodimers of benzamide with benzoic acid, salicylic acid and phenyl acetic acid in deuterated chloroform solution. Association constants for the complex formation in the solution state have been determined. The results are ascertained by X-ray diffraction in the solid state. Intermolecular interactions in solution and in solid state were found to be similar. The structural parameters obtained by X-ray diffraction studies are compared with those obtained by DFT calculations. (C) 2012 Elsevier B.V. All rights reserved.