Towards an exact factorization of the molecular wave function

An exact single-product factorisation of the molecular wave function for the timedependent Schrodinger equation is investigated by using an ansatz involving a phasefactor. By using the Frenkel variational method, we obtain the Schrodinger equations for the electronic and nuclear wave functions. The concept of a potential energy surface (PES) is retained by introducing a modified Hamiltonian as suggested earlier by Cederbaum. The parameter in the phase factor is chosen such that the equations of motion retain the physically appealing Born- Oppenheimer-like form, and is therefore unique.

Crystal chemistry and photomechanical behavior of 3,4-dimethoxycinnamic acid: correlation between maximum yield in the solid-state topochemical reaction and cooperative molecular motion

A new monoclinic polymorph, form II (P2(1)/c, Z = 4), has been isolated for 3,4-dimethoxycinnamic acid (DMCA). Its solid-state 2 + 2 photoreaction to the corresponding alpha-truxillic acid is different from that of the first polymorph, the triclinic form I (P (1) over bar, Z = 4) that was reported in 1984. The crystal structures of the two forms are rather different. The two polymorphs also exhibit different photomechanical properties. Form I exhibits photosalient behavior but this effect is absent in form II. These properties can be explained on the basis of the crystal packing in the two forms. The nanoindentation technique is used to shed further insights into these structure-property relationships. A faster photoreaction in form I and a higher yield in form II are rationalized on the basis of the mechanical properties of the individual crystal forms. It is suggested that both Schmidt-type and Kaupp-type topochemistry are applicable for the solid-state trans-cinnamic acid photodimerization reaction. Form I of DMCA is more plastic and seems to react under Kaupp-type conditions with maximum molecular movements. Form II is more brittle, and its interlocked structure seems to favor Schmidt-type topochemistry with minimum molecular movement.

Rational molecular dynamics scheme for predicting optimum concentration loading of nano-additive in phase change materials

The present study deals with the diffusion and phase transition behaviour of paraffin reinforced with carbon nano-additives namely graphene oxide (GO) and surface functionalized single walled carbon nanotubes (SWCNT). Bulk disordered systems of paraffin hydrocarbons impregnated with carbon nano-additives have been generated in realistic equilibrium conformations for potential application as latent heat storage systems. Ab initio molecular dynamics(MD) in conjugation with COMPASS forcefield has been implemented using periodic boundary conditions. The proposed scheme allows determination of optimum nano-additive loading for improving thermo-physical properties through analysis of mass, thermal and transport properties; and assists in determination of composite behaviour and related performance from microscopic point of view. It was observed that nanocomposites containing 7.8% surface functionalised SWCNT and 55% GO loading corresponds to best latent heat storage system. The propounded methodology could serve as a by-pass route for economically taxing and iterative experimental procedures required to attain the optimum composition for best performance. The results also hint at the large unexplored potential of ab-initio classical MD techniques for predicting performance of new nanocomposites for potential phase change material applications. (C) 2015 Author(s).

Estimation of activation energy for electroporation and pore growth rate in liquid crystalline and gel phases of lipid bilayers using molecular dynamics simulations

Molecular dynamics simulations of electroporation in POPC and DPPC lipid bilayers have been carried out at different temperatures ranging from 230 K to 350 K for varying electric fields. The dynamics of pore formation, including threshold field, pore initiation time, pore growth rate, and pore closure rate after the field is switched off, was studied in both the gel and liquid crystalline (L-alpha) phases of the bilayers. Using an Arrhenius model of pore initiation kinetics, the activation energy for pore opening was estimated to be 25.6 kJ mol(-1) and 32.6 kJ mol(-1) in the L-alpha phase of POPC and DPPC lipids respectively at a field strength of 0.32 V nm(-1). The activation energy decreases to 24.2 kJ mol(-1) and 23.7 kJ mol(-1) respectively at a higher field strength of 1.1 V nm(-1). At temperatures below the melting point, the activation energy in the gel phase of POPC and DPPC increases to 28.8 kJ mol(-1) and 34.4 kJ mol(-1) respectively at the same field of 1.1 V nm(-1). The pore closing time was found to be higher in the gel than in the L-alpha phase. The pore growth rate increases linearly with temperature and quadratically with field, consistent with viscosity limited growth models.

H2O-CH4 and H2S-CH4 complexes: a direct comparison through molecular beam experiments and ab initio calculations

New molecular beam scattering experiments have been performed to measure the total ( elastic plus inelastic) cross sections as a function of the velocity in collisions between water and hydrogen sulfide projectile molecules and the methane target. Measured data have been exploited to characterize the range and strength of the intermolecular interaction in such systems, which are of relevance as they drive the gas phase molecular dynamics and the clathrate formation. Complementary information has been obtained by rotational spectra, recorded for the hydrogen sulfide-methane complex, with a pulsed nozzle Fourier transform microwave spectrometer. Extensive ab initio calculations have been performed to rationalize all the experimental findings. The combination of experimental and theoretical information has established the ground for the understanding of the nature of the interaction and allows for its basic components to be modelled, including charge transfer, in these weakly bound systems. The intermolecular potential for H2S-CH4 is significantly less anisotropic than for H2O-CH4, although both of them have potential minima that can be characterized as `hydrogen bonded'.

Synthesis, ABTS-Radical Scavenging Activity, and Antiproliferative and Molecular Docking Studies of Novel Pyrrolo1,2-a]quinoline Derivatives

A new 2,2-azinobis-(3-ethylbenzothiazoline-6-sulfonic acid) (ABTS)-radical scavenging and antiproliferative agents of pyrrolo1,2-a]quinoline derivatives have been synthesized. An efficient method for the synthesis of 14 novel diversified pyrrolo1,2-a]quinoline derivatives has been described using 4-(1,3-dioxolan-2-yl)quinoline and different phenacyl bromides in acetone and followed by reacting with different acetylenes in dimethylformamide/K2CO3. The structure of the newly synthesized compounds was determined by infrared, H-1 NMR, C-13 NMR, mass spectrometry, and elemental analysis. The in vitro antioxidant activity revealed that among all the tested compounds 5n exhibited maximum scavenging activity with ABTS. Compound 5b has showed good antiproliferative activity as an inhibitor of epidermal growth factor receptor (EGFR) tyrosine kinase.

Structural studies on a non-toxic homologue of type II RIPs from bitter gourd: Molecular basis of non-toxicity, conformational selection and glycan structure

The structures of nine independent crystals of bitter gourd seed lectin (BGSL), a non-toxic homologue of type II RIPs, and its sugar complexes have been determined. The four-chain, two-fold symmetric, protein is made up of two identical two-chain modules, each consisting of a catalytic chain and a lectin chain, connected by a disulphide bridge. The lectin chain is made up of two domains. Each domain carries a carbohydrate binding site in type II RIPs of known structure. BGSL has a sugar binding site only on one domain, thus impairing its interaction at the cell surface. The adenine binding site in the catalytic chain is defective. Thus, defects in sugar binding as well as adenine binding appear to contribute to the non-toxicity of the lectin. The plasticity of the molecule is mainly caused by the presence of two possible well defined conformations of a surface loop in the lectin chain. One of them is chosen in the sugar complexes, in a case of conformational selection, as the chosen conformation facilitates an additional interaction with the sugar, involving an arginyl residue in the loop. The N-glycosylation of the lectin involves a plant-specific glycan while that in toxic type II RIPs of known structure involves a glycan which is animal as well as plant specific.

Molecular profiling of sepsis in mice using Fourier Transform Infrared Microspectroscopy

Sepsis is a life threatening condition resulting from a high burden of infection. It is a major health care problem and associated with inflammation, organ dysfunction and significant mortality. However, proper understanding and delineating the changes that occur during this complex condition remains a challenge. A comparative study involving intra-peritoneal injection of BALB/c mice with Salmonella Typhimurium (infection), lipopolysaccharide (endotoxic shock) or thioglycollate (sterile peritonitis) was performed. The changes in organs and sera were profiled using immunological assays and Fourier Transform Infrared (FTIR) micro-spectroscopy. There is a rapid rise in inflammatory cytokines accompanied with lowering of temperature, respiratory rate and glucose amounts in mice injected with S. Typhimurium or lipopolysaccharide. FTIR identifies distinct changes in liver and sera: decrease in glycogen and protein/lipid ratio and increase in DNA and cholesteryl esters. These changes were distinct from the pattern observed in mice treated with thioglycollate and the differences in the data obtained between the three models are discussed. The combination of FTIR spectroscopy and other biomarkers will be valuable in monitoring molecular changes during sepsis. GRAPHICS] Intra-peritoneal infection with high dose of Salmonella Typhimurium leads to rapid increase in inflammatory cytokines, e.g. Tnf alpha (A). FTIR analysis of liver (B) and sera (C) identifies several metabolic changes: glycogen, protein/lipid, cholesteryl esters and DNA.

MS3ALIGN: an efficient molecular surface aligner using the topology of surface curvature

Background: Aligning similar molecular structures is an important step in the process of bio-molecular structure and function analysis. Molecular surfaces are simple representations of molecular structure that are easily constructed from various forms of molecular data such as 3D atomic coordinates (PDB) and Electron Microscopy (EM) data. Methods: We present a Multi-Scale Morse-Smale Molecular-Surface Alignment tool, MS3ALIGN, which aligns molecular surfaces based on significant protrusions on the molecular surface. The input is a pair of molecular surfaces represented as triangle meshes. A key advantage of MS3ALIGN is computational efficiency that is achieved because it processes only a few carefully chosen protrusions on the molecular surface. Furthermore, the alignments are partial in nature and therefore allows for inexact surfaces to be aligned. Results: The method is evaluated in four settings. First, we establish performance using known alignments with varying overlap and noise values. Second, we compare the method with SurfComp, an existing surface alignment method. We show that we are able to determine alignments reported by SurfComp, as well as report relevant alignments not found by SurfComp. Third, we validate the ability of MS3ALIGN to determine alignments in the case of structurally dissimilar binding sites. Fourth, we demonstrate the ability of MS3ALIGN to align iso-surfaces derived from cryo-electron microscopy scans. Conclusions: We have presented an algorithm that aligns Molecular Surfaces based on the topology of surface curvature. Awebserver and standalone software implementation of the algorithm available at http://vgl.serc.iisc.ernet. in/ms3align.

On the critical thickness for non-localized to localized plastic flow transition in metallic glasses: A molecular dynamics study

Molecular dynamics simulations were employed to investigate the specimen thickness-dependent tensile behavior of a series of Cu(x)Z(100-x) (x = 20, 40, 50, 64 and 80 at%) metallic glass (MG) films, with a particular focus on the critical thickness, tc, below which non-localized plastic flow takes place. The simulation results reveal that while the transition occurs in all the alloys examined, t(c) is sensitive to the composition. We rationalize t(c) by postulating that the strain energy stored in the sample at the onset of plastic deformation has to be sufficient for the formation of shear bands. The composition-dependence of t(c) was found to correlate with the average activation energy of the atomic level plastic deformation events. (C) 2015 Elsevier Ltd. All rights reserved.

Docking, molecular dynamics and QM/MM studies to delineate the mode of binding of CucurbitacinE to F-actin

CucurbitacinE (CurE) has been known to bind covalently to F-actin and inhibit depolymerization. However, the mode of binding of CurE to F-actin and the consequent changes in the F-actin dynamics have not been studied. Through quantum mechanical/molecular mechanical (QM/MM) and density function theory (DFT) simulations after the molecular dynamics (MD) simulations of the docked complex of F-actin and CurE, a detailed transition state (TS) model for the Michael reaction is proposed. The TS model shows nucleophilic attack of the sulphur of Cys257 at the beta-carbon of Michael Acceptor of CurE producing an enol intermediate that forms a covalent bond with CurE. The MD results show a clear difference between the structure of the F-actin in free form and F-actin complexed with CurE. CurE affects the conformation of the nucleotide binding pocket increasing the binding affinity between F-actin and ADP, which in turn could affect the nucleotide exchange. CurE binding also limits the correlated displacement of the relatively flexible domain 1 of F-actin causing the protein to retain a flat structure and to transform into a stable ``tense'' state. This structural transition could inhibit depolymerization of F-actin. In conclusion, CurE allosterically modulates ADP and stabilizes F-actin structure, thereby affecting nucleotide exchange and depolymerization of F-actin. (C) 2015 Elsevier Inc. All rights reserved.

Differential dynamics of the serotonin(1A) receptor in membrane bilayers of varying cholesterol content revealed by all atom molecular dynamics simulation

The serotonin(1A) receptor belongs to the superfamily of G protein-coupled receptors (GPCRs) and is a potential drug target in neuropsychiatric disorders. The receptor has been shown to require membrane cholesterol for its organization, dynamics and function. Although recent work suggests a close interaction of cholesterol with the receptor, the structural integrity of the serotonin(1A) receptor in the presence of cholesterol has not been explored. In this work, we have carried out all atom molecular dynamics simulations, totaling to 3s, to analyze the effect of cholesterol on the structure and dynamics of the serotonin(1A) receptor. Our results show that the presence of physiologically relevant concentration of membrane cholesterol alters conformational dynamics of the serotonin(1A) receptor and, on an average lowers conformational fluctuations. Our results show that, in general, transmembrane helix VII is most affected by the absence of membrane cholesterol. These results are in overall agreement with experimental data showing enhancement of GPCR stability in the presence of membrane cholesterol. Our results constitute a molecular level understanding of GPCR-cholesterol interaction, and represent an important step in our overall understanding of GPCR function in health and disease.

Targeting human telomeric G-quadruplex DNA with curcumin and its synthesized analogues under molecular crowding conditions

The formation of telomeric G-quadruplexes has been shown to inhibit telomerase activity. Indeed, a number of small molecules capable of p-stacking with G-tetrads have shown the ability to inhibit telomerase activity through the stabilization of G-quadruplexes. Curcumin displays a wide spectrum of medicinal properties ranging from anti-bacterial, anti-viral, anti-protozoal, anti-fungal and anti-inflammatory to anti-cancer activity. We have investigated the interactions of curcumin and its structural analogues with the human telomeric sequence AG(3)(T(2)AG(3))(3) under molecular crowding conditions. Experimental studies indicated the existence of a AG(3)(T(2)AG(3))(3)/curcumin complex with binding affinity of 0.72 x 10(6) M-1 under molecular crowding conditions. The results from UV-visible absorption spectroscopy, a fluorescent TO displacement assay, circular dichroism and molecular docking studies, imply that curcumin and their analogues interact with G-quadruplex DNA via groove binding. While other analogs of curcumin studied here bind to G-quadruplexes in a qualitatively similar manner their affinities are relatively lower in comparison to curcumin. The Knoevenagel condensate, a methoxy-benzylidene derivative of curcumin, also exhibited significant binding to G-quadruplex DNA, although with two times decreased affinity. Our study establishes the potential of curcumin as a promising natural product for G-quadruplex specific ligands.