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.

Sorbate properties and cage-to-cage diffusion of argon in NaCaA: a molecular dynamics study

Detailed molecular dynamics simulations of argon in zeolite NaCaA are reported. Thermodynamic, structural, and dynamical properties of the sorbate as a function of temperature have been obtained. The properties calculated include various site-site radial distribution functions, different energy distribution functions, selfdiffusion coefficients, the power spectra, and properties relating to cage-to-cage diffusion. The results suggest that sorbate is delocalized above 300 K. Both modes of cage-to-cage diffusion-the surface-mediated and centralized diffusion-are associated with negative barrier heights. Surprisingly, rate of cage-to-cage diffusion is associated with negative and positive activation energies below and above 500 K. The observed differences in the behavior of the rate of cage-to-cage diffusion between Xe-NaY and Ar-NaCaA systems and the nature of the potential energy surface are discussed. Presence of sorbatezeolite interactions results in significant enhancement in the rate of cage-to-cage diffusion and rate of cage visits. It is shown that properties dependent on the long-time behavior such as the diffusion coefficient and the rate of cages visited exhibit the expected Arrhenius dependence on temperature.

Stereochemistry of 2'-5' linked nucleic acids: crystal and molecular structure of ammonium adenylyl-2',5'-adenosine tetrahydrate: a core fragment of 2'-5' oligo A's produced by interferon induced adenylate synthetase

The preponderance of 3'-5' phosphodiester links in nucleic acids is well known. Albeit less prevalent, the 2'-5' links are specifically utilised in the formation of 'lariat' in group II introns and in the msDNA-RNA junction in myxobacterium. As a sequel to our earlier study on cytidylyl-2',5'-adenosine we have now obtained the crystal structure of adenylyl-2',5'-adenosine (A2'p5'A) at atomic resolution. This dinucleoside monophosphate crystallizes in the orthorhombic space group P2(1)2(1)2(1) with a = 7.956(3) A, b = 12.212(3) A and c = 36.654(3) A. CuK alpha intensity data were collected on a diffractometer. The structure was sloved by direct methods and refined by full matrix least squares methods to R = 10.8%. The 2' terminal adenine is in the commonly observed anti (chi 2 = 161 degrees) conformation and the 5' terminal base has a syn (chi 1 = 55 degrees) conformation more often seen in purine nucleotides. A noteworthy feature of A2'p5'A is the intranucleotide hydrogen bond between N3 and O5' atoms of the 5' adenine base. The two furanose rings in A2'p5'A show different conformations - C2' endo, C3' endo puckering for the 5' and 2' ends respectively. In this structure too there is a stacking of the purine base on the ribose O4' just as in other 2'-5' dinucleoside structures, a feature characteristically seen in the left handed Z DNA. In having syn, anti conformation about the glycosyl bonds, C2' endo, C3' endo mixed sugar puckering and N3-O5' intramolecular hydrogen bond A2'p5'A resembles its 3'-5' analogue and several other 2'-5' dinucleoside monophosphate structures solved so far. Striking similarities between the 2'-5' dinucleoside monophosphate structures suggest that the conformation of the 5'-end nucleoside dictates the conformation of the 2' end nucleoside. Also, the 2'-5' dimers do not favour formation of miniature classical double helical structures like the 3'-5' dimers. It is conceivable, 2-5(A) could be using the stereochemical features of A2'p5'A which accounts for its higher activity.

Ab initio molecular orbital calculations on ion pair-water complexes of metal halides and oxides

Ab initio MO calculations are performed on a series of ion-molecular and ion pair-molecular complexes of H2O + MX (MX = LiF, LiCl, NaCl, BeO and MgO) systems. BSSE-corrected stabilization energies, optimized geometrical parameters, internal force constants and harmonic vibrational frequencies have been evaluated for all the structures of interest. The trends observed in the geometrical parameters and other properties calculated for the mono-hydrated contact ion pair complexes parallel those computed for the complexes of the individual ions. The bifurcated structures are found to be saddle points with an imaginary frequency corresponding to the rocking mode of water molecules. The solvent-shared ion pair complexes have high interaction energies. Trends in the internal force constant and harmonic frequency values are discussed in terms of ion-molecular and ion-pair molecular interactions.

A Molecular Dynamics Study of Atomic Correlations in Glassy BzS3t

Glassy B&, the parent compound of the superionic conductor LiI-Li&B& has been studied by the molecular dynamics technique using a new potential model. The results suggest that the glass is made up of local units of four-membered B2S2 rings bridged by sulfur atoms, leading to a chainlike structure. Various pair correlation functions have been analyzed, and the B2Sz rings have been found to be planar. The calculated neutron structure factor shows a peak at 1.4 A-' which has been attributed to B-B correlations at 5.6 A. The glass transition temperature of the simulated system has been calculated to be around 800 K.

Analysis of Hydrogen Bonding and Stability of Protein Secondary Structures in Molecular Dynamics Simulation

Molecular Dynamics (MD) simulations provide an atomic level account of the molecular motions and have proven to be immensely useful in the investigation of the dynamical structure of proteins. Once an MD trajectory is obtained, specific interactions at the molecular level can be directly studied by setting up appropriate combinations of distance and angle monitors. However, if a study of the dynamical behavior of secondary structures in proteins becomes important, this approach can become unwieldy. We present herein a method to study the dynamical stability of secondary structures in proteins, based on a relatively simple analysis of backbone hydrogen bonds. The method was developed for studying the thermal unfolding of beta-lactamases, but can be extended to other systems and adapted to study relevant properties.

1H NMR study of molecular dynamics and phase transition in (NH4)2ZnBr4

The proton second moment (M2) and spin-lattice relaxation time (T1) have been measured in (NH4)2ZnBr4 in the range 77-300 K. The room-temperature spectrum shows a structure which disappears around 243 K. The signal is strong and narrow even at 77 K. Proton T1 shows a maximum at 263 K, caused by spin rotation interaction and decreases with decreasing temperature till 235 K, where it shows a sudden increase. Below 235 K, again it decreases and shows a slope change around 216.5 K (reported Tc). From 216.5 K, T1 decreases continuously without exhibiting any minimum down to 77 K. The narrow line at 77 K, and absence of a T1 minimum down to 77 K indicate the possibility of quantum mechanical tunnelling in this system. Motional parameters such as activation energy and pre-exponential factor have been evaluated for the reorientational motion of the NH+4 ion.

Molecular systematics and conservation of the langurs and leaf monkeys of South Asia

Numerous morphology-based classification schemes have been proposed for langurs and leaf monkeys of South Asia but there is very little agreement between them. An incorrect classification scheme when used as a basis for biogeographic studies can support erroneous hypotheses. Further, lack of taxonomic resolution will also confound conservation efforts, given that conservation biologists use traditional morphology-based-classification schemes to prioritize species for conservation. Here, I have revisited recent molecular phylogenetic studies done on langurs and leaf monkeys of South Asia. Results from these studies are in turn used to derive a rational and scientific basis for prioritizing species for conservation. Molecular data support the classification of langurs of the Indian subcontinent-Hanuman, Nilgiri and purple-faced langurs-in the genus Semnopithecus, whereas Phayre's leaf monkey along with other Southeast Asian leaf monkeys form another distinct clade (Trachypithecus). The phylogenetic position of capped and golden langurs remains unresolved. Molecular data suggest that they are closely related to each other but this group might have evolved through past hybridization between Semnopithecus and Trachypithecus. Additionally, genetic data also support the splitting of the so-called Hanuman langurs into at least three species. The scores for taxonomic uniqueness of langurs and leaf monkeys of South Asia were revised using this molecular phylogeny-based classification. According to the revised scores, Phayres leaf monkey and golden langur are priority species for conservation followed by capped and Nilgiri langurs.

Self-assembly of neutral and cationic Pd-II organometallic molecular rectangles: synthesis, characterization and nitroaromatic sensing

Design and synthesis of three novel 2 + 2] self-assembled molecular rectangles 1-3 via coordination driven self-assembly of predesigned Pd(II) ligands is reported. 1,8-Diethynylanthracene was assembled with trans-Pd(PEt3)(2)Cl-2 in the presence of CuCl catalyst to yield a neutral rectangle 1 via Pd-C bond formation. Complex 1 represents the first example of a neutral molecular rectangle obtained via C-Pd coordination driven self-assembly. A new Pd-2(II) organometallic building block with 180 degrees bite-angle 1,4-bistrans-(ethynyl)Pd(PEt3)(2)(NO3)] benzene (M-2) containing ethynyl functionality was synthesized in reasonable yield by employing Sonagashira coupling reaction. Self-assembly of M-2 with two organic clip-type donors (L-2-L-3) afforded 2 + 2] self-assembled molecular rectangles 2 and 3, respectively L-2 = 1,8-bis(4-pyridylethynyl) anthracene; L-3 = 1,3-bis(3-pyridyl) isophthalamide]. The macrocycles 1-3 were fully characterized by multinuclear NMR and ESI-MS spectroscopic techniques, and in case of 1 the structure was unambiguously determined by single crystal X-ray diffraction analysis. Incorporation of Pd-ethynyl bonds helped to make the assemblies p-electron rich and fluorescent in nature. Complexes 1-2 showed quenching of fluorescence intensity in solution in presence of nitroaromatics, which are the chemical signatures of many commercially available explosives.

Crystal structures of fluorinated aryl biscarbonates and a biscarbamate: a counterpoise between weak intermolecular interactions and molecular symmetry

Conformational features and supramolecular structural organization in three aryl biscarbonates and an aryl biscarbamate with rigid acetylenic unit providing variable spacer lengths have been probed to gain insights into the packing features associated with molecular symmetry and the intermolecular interactions involving `organic' fluorine. Four structures but-2-yne-1,4-diyl bis(2,3,4,5,6-pentafluorophenylcarbonate), 1; but-2-yne-1,4-diyl bis(4-fluorophenylcarbonate), 2; but-2-yne-1,4-diyl bis(2,3,4,5,6-pentafluorophenylcarbamate), 3 and hexa-2,4-diyne-1,6-diyl bis(2,3,4,5,6-pentafluorophenylcarbonate), 4 have been analyzed in this context. Compound 1 adopts a non-centrosymmetric ``twisted'' (syn) conformation, whereas 2, 3 and 4 acquire a centrosymmetric ``extended'' (anti) conformation. Weak intermolecular interactions and in particular those involving fluorine are found to dictate this conformational variation in the crystal structure of 1. A single-crystal neutron diffraction study at 90 K was performed on 1 to obtain further insights into these interactions involving `organic' fluorine.

Bile acid based semi-rigid molecular tweezers

The easily constructed bile acid-based semi-rigid molecular tweezer 2 binds guest 8 in chloroform with an association constant of 83 dm(3) mol(-1).

Binding of active site directed ligands to phospholipase A2: Implications on the molecular constraints and catalytic mechanism

Molecular constraints for the localization of active site directed ligands (competitive inhibitors and substrates) in the active site of phospholipase A2 (PLA2) are characterized. Structure activity relationships with known inhibitors suggest that the head : group interactions dominate the selectivity as well as a substantial part of the affinity. The ab initio fitting of the amide ligands in the active site was carried out to characterize the head group interactions. Based on a systematic coordinate space search, formamide is docked with known experimental constraints such as coordination of the carbonyl group to Ca2+ and hydrogen bond between amide nitrogen and ND1 of His48. An optimal position for a bound water molecule is identified and its significance for the catalytic mechanism is postulated. Unlike the traditional ''pseudo-triad'' mechanism, the ''Ca-coordinatedoxyanion'' mechanism proposed here invokes activation of the catalytic water to form the oxyanion in the coordination sphere of calcium. As it attacks the carbonyl carbon of the ester, a near-tetrahedral intermediate is formed. As the second proton of the catalytic water is abstracted by the ester oxygen, its reorientation and simultaneous cleavage form hydrogen bond with ND1 of His48. In this mechanism of esterolysis, a catalytic role for the water co-ordinated to Ca2+ is recognised.

Calcium ionophore, A23187 and its amino acid complexes: spectroscopic and molecular modeling studies.

The circular dichroism, fluorescence, Nuclear Magnetic Resonance and BLM conductance studies indicate that A23187 forms a stable complex with amino acids at low ionophore concentrations (<10(-4)M). However, A23187 prefers to be in a dimeric structure with no significant binding to amino acids, at concentrations higher than 10(-4)M. It was also observed that at lower concentrations, at which the amino acids bind to the ionophore, the affinity for calcium ions was several orders of magnitude lower than that at higher ionophore concentrations. We have also conducted molecular modeling studies to examine the structure of the A23187 dimer and its amino acid complexes. The results of these modeling studies strongly support our experimental results and validate the formation of a hydrogen bonded and energetically stable A23187 dimer and its amino acid complexes.

Control of lipid microstructures by molecular design and its implications

An overview of the current trends in the lipid design for specific applications has been presented. Lipids with different surface charge and hydrophobic backbone undergo aggregation to produce lamellae or bilayer and multilayer vesicles in aqueous media. Various aspects of present development of chiral superstructures and enzyme-mimics have been discussed. Utility of these molecules for potential applications in immunomodulation and sustained drug-delivery systems is also summarized.

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.