Investigating the molecular structures of solid dispersions by the simulated annealing method

Knowledge of the molecular structures of solid dispersions is vital, yet, despite thousands of reports in this area, it remains unclear. The aim of this research is to investigate the molecular structure of solid dispersions with hot melt preparation method by the simulated annealing method. Simulation results showed linear polymer chains form the random coils under heat and the drug molecules stick on the surface of polymer coils, while drug molecules are dispersed molecularly but irregularly within the amorphous low molecular weight carriers. This research presents more reasonable molecular images of solid dispersions than the existed theory.

Complex temporal patterns in molecular dynamics:a direct measure of the phase-space exploration by the trajectory at macroscopic time scales

Computer simulated trajectories of bulk water molecules form complex spatiotemporal structures at the picosecond time scale. This intrinsic complexity, which underlies the formation of molecular structures at longer time scales, has been quantified using a measure of statistical complexity. The method estimates the information contained in the molecular trajectory by detecting and quantifying temporal patterns present in the simulated data (velocity time series). Two types of temporal patterns are found. The first, defined by the short-time correlations corresponding to the velocity autocorrelation decay times (â‰0.1â€ps), remains asymptotically stable for time intervals longer than several tens of nanoseconds. The second is caused by previously unknown longer-time correlations (found at longer than the nanoseconds time scales) leading to a value of statistical complexity that slowly increases with time. A direct measure based on the notion of statistical complexity that describes how the trajectory explores the phase space and independent from the particular molecular signal used as the observed time series is introduced. © 2008 The American Physical Society.

Molecular phase space transport in water:non-stationary random walk model

Molecular transport in phase space is crucial for chemical reactions because it defines how pre-reactive molecular configurations are found during the time evolution of the system. Using Molecular Dynamics (MD) simulated atomistic trajectories we test the assumption of the normal diffusion in the phase space for bulk water at ambient conditions by checking the equivalence of the transport to the random walk model. Contrary to common expectations we have found that some statistical features of the transport in the phase space differ from those of the normal diffusion models. This implies a non-random character of the path search process by the reacting complexes in water solutions. Our further numerical experiments show that a significant long period of non-stationarity in the transition probabilities of the segments of molecular trajectories can account for the observed non-uniform filling of the phase space. Surprisingly, the characteristic periods in the model non-stationarity constitute hundreds of nanoseconds, that is much longer time scales compared to typical lifetime of known liquid water molecular structures (several picoseconds).

Study of neutral Fe(III) complexes of pyridoxal-N-substituted thiosemicarbazone with desolvation induced spin-state transformation above room temperature

The preparation and characterization of two new neutral ferric complexes with desolvation-induced discontinuous spin-state transformation above room temperature are reported. The compounds, Fe(Hthpy)(thpy).CH3OH.3H2O (1) and Fe(Hmthpy)(mthpy).2H2O (2), are low-spin (LS) at room temperature and below, whereas their nonsolvated forms are high-spin (HS), exhibiting zero-field splitting. In these complexes, Hthpy, Hmthpy, and thpy, mthpy are the deprotonated forms of pyridoxal thiosemicarbazone and pyridoxal methylthiosemicarbazone, respectively; each is an O,N,S-tridentate ligand. The molecular structures have been determined at 100(1) K using single-crystal X-ray diffraction techniques and resulted in a triclinic system (space group P1) and monoclinic unit cell (space group P21/c) for 1 and 2, respectively. Structures were refined to the final error indices, where RF = 0.0560 for 1 and RF = 0.0522 for 2. The chemical inequivalence of the ligands was clearly established, for the "extra" hydrogen atom on the monodeprotonated ligands (Hthpy, Hmthpy) was found to be bound to the nitrogen of the pyridine ring. The ligands are all of the thiol form; the doubly deprotonated chelates (thpy, mthpy) have C-S bond lengths slightly longer than those of the singly deprotonated forms. There is a three-dimensional network of hydrogen bonds in both compounds. The discontinuous spin-state transformation is accompanied with liberation of solvate molecules. This is evidenced also from DSC analysis. Heat capacity data for the LS and HS phases are tabulated at selected temperatures, the values of the enthalpy and entropy changes connected with the change of spin state were reckoned at DeltaH = 12.5 0.3 kJ mol-1 and DeltaS = 33.3 0.8 J mol-1 K-1, respectively, for 1 and DeltaH = 6.5 0.3 kJ mol-1 and DeltaS = 17.6 0.8 J mol-1 K-1, respectively, for 2

Flexibility of short-strand RNA in aqueous solution as revealed by molecular dynamics simulation:are A′-RNA and A-RNA distinct conformational structures?

We use molecular dynamics simulations to compare the conformational structure and dynamics of a 21-base pair RNA sequence initially constructed according to the canonical A-RNA and A'-RNA forms in the presence of counterions and explicit water. Our study aims to add a dynamical perspective to the solid-state structural information that has been derived from X-ray data for these two characteristic forms of RNA. Analysis of the three main structural descriptors commonly used to differentiate between the two forms of RNA namely major groove width, inclination and the number of base pairs in a helical twist over a 30 ns simulation period reveals a flexible structure in aqueous solution with fluctuations in the values of these structural parameters encompassing the range between the two crystal forms and more. This provides evidence to suggest that the identification of distinct A-RNA and A'-RNA structures, while relevant in the crystalline form, may not be generally relevant in the context of RNA in the aqueous phase. The apparent structural flexibility observed in our simulations is likely to bear ramifications for the interactions of RNA with biological molecules (e.g. proteins) and non-biological molecules (e.g. non-viral gene delivery vectors). © CSIRO 2009.

Nuclear magnetic resonance studies of molecular interactions and structures

This thesis is concerned with the investigation, by nuclear magnetic resonance spectroscopy, of the molecular interactions occurring in mixtures of benzene and cyclohexane to which either chloroform or deutero-chloroform has been added. The effect of the added polar molecule on the liquid structure has been studied using spin-lattice relaxation time, 1H chemical shift, and nuclear Overhauser effect measurements. The main purpose of the work has been to validate a model for molecular interaction involving local ordering of benzene around chloroform. A chemical method for removing dissolved oxygen from samples has been developed to encompass a number of types of sample, including quantitative mixtures, and its supremacy over conventional deoxygenation technique is shown. A set of spectrometer conditions, the use of which produces the minimal variation in peak height in the steady state, is presented. To separate the general diluting effects of deutero-chloroform from its effects due to the production of local order a series of mixtures involving carbon tetrachloride, instead of deutero-chloroform, have been used as non-interacting references. The effect of molecular interaction is shown to be explainable using a solvation model, whilst an approach involving 1:1 complex formation is shown not to account for the observations. It is calculated that each solvation shell, based on deutero-chloroform, contains about twelve molecules of benzene or cyclohexane. The equations produced to account for the T1 variations have been adapted to account for the 1H chemical shift variations in the same system. The shift measurements are shown to substantiate the solvent cage model with a cage capacity of twelve molecules around each chloroform molecule. Nuclear Overhauser effect data have been analysed quantitatively in a manner consistent with the solvation model. The results show that discrete shells only exist when the mole fraction of deutero-chloroform is below about 0.08.

Molecular dynamics simulations:bring biomolecular structures alive on a computer

The molecular dynamics (MD) simulations play a very important role in science today. They have been used successfully in binding free-energy calculations and rational design of drugs and vaccines. MD simulations can help visualize and understand structures and dynamics at an atomistic level when combined with molecular graphics programs. The molecular and atomistic properties can be displayed on a computer in a time-dependent way, which opens a road toward a better understanding of the relationship of structure, dynamics, and function. In this chapter, the basics of MD are explained, together with a step-by-step description of setup and running an MD simulation.

Characterization of the structure of RAMP1 by mutagenesis and molecular modeling

Receptor activity modifying proteins (RAMPs) are a family of single-pass transmembrane proteins that dimerize with G-protein-coupled receptors. They may alter the ligand recognition properties of the receptors (particularly for the calcitonin receptor-like receptor, CLR). Very little structural information is available about RAMPs. Here, an ab initio model has been generated for the extracellular domain of RAMP1. The disulfide bond arrangement (Cys 27-Cys82, Cys40-Cys72, and Cys 57-Cys104) was determined by site-directed mutagenesis. The secondary structure (a-helices from residues 29-51, 60-80, and 87-100) was established from a consensus of predictive routines. Using these constraints, an assemblage of 25,000 structures was constructed and these were ranked using an all-atom statistical potential. The best 1000 conformations were energy minimized. The lowest scoring model was refined by molecular dynamics simulation. To validate our strategy, the same methods were applied to three proteins of known structure; PDB:1HP8, PDB:1V54 chain H (residues 21-85), and PDB:1T0P. When compared to the crystal structures, the models had root mean-square deviations of 3.8 Å, 4.1 Å, and 4.0 Å, respectively. The model of RAMP1 suggested that Phe93, Tyr 100, and Phe101 form a binding interface for CLR, whereas Trp74 and Phe92 may interact with ligands that bind to the CLR/RAMP1 heterodimer. © 2006 by the Biophysical Society.

Structure, dynamics, and energetics of siRNA-cationic vector complexation:a molecular dynamics study

The design and synthesis of safe and efficient nonviral vectors for gene delivery has attracted significant attention in recent years. Previous experiments have revealed that the charge density of a polycation (the carrier) plays a crucial role in complexation and the release of the gene from the complex in the cytosol. In this work, we adopt an atomistic molecular dynamics simulation approach to study the complexation of short strand duplex RNA with six cationic carrier systems of varying charge and surface topology. The simulations reveal detailed molecular-level pictures of the structures and dynamics of the RNA-polycation complexes. Estimates for the binding free energy indicate that electrostatic contributions are dominant followed by van der Waals interactions. The binding free energy between the 8(+)polymers and the RNA is found to be larger than that of the 4(+)polymers, in general agreement with previously published data. Because reliable binding free energies provide an effective index of the ability of the polycationic carrier to bind the nucleic acid and also carry implications for the process of gene release within the cytosol, these novel simulations have the potential to provide us with a much better understanding of key mechanistic aspects of gene-polycation complexation and thereby advance the rational design of nonviral gene delivery systems.

Structure and dynamics of multiple cationic vectors-siRNA complexation by all-atomic molecular dynamics simulations

Understanding the molecular mechanism of gene condensation is a key component to rationalizing gene delivery phenomena, including functional properties such as the stability of the gene-vector complex and the intracellular release of the gene. In this work, we adopt an atomistic molecular dynamics simulation approach to study the complexation of short strand duplex RNA with four cationic carrier systems of varying charge and surface topology at different charge ratios. At lower charge ratios, polymers bind quite effectively to siRNA, while at high charge ratios, the complexes are saturated and there are free polymers that are unable to associate with RNA. We also observed reduced fluctuations in RNA structures when complexed with multiple polymers in solution as compared to both free siRNA in water and the single polymer complexes. These novel simulations provide a much better understanding of key mechanistic aspects of gene-polycation complexation and thereby advance progress toward rational design of nonviral gene delivery systems.

Characterizing the hierarchical structures of bioactive sol-gel silicate glass and hybrid scaffolds for bone regeneration

Bone is the second most widely transplanted tissue after blood. Synthetic alternatives are needed that can reduce the need for transplants and regenerate bone by acting as active temporary templates for bone growth. Bioactive glasses are one of the most promising bone replacement/regeneration materials because they bond to existing bone, are degradable and stimulate new bone growth by the action of their dissolution products on cells. Sol-gel-derived bioactive glasses can be foamed to produce interconnected macropores suitable for tissue ingrowth, particularly cell migration and vascularization and cell penetration. The scaffolds fulfil many of the criteria of an ideal synthetic bone graft, but are not suitable for all bone defect sites because they are brittle. One strategy for improving toughness of the scaffolds without losing their other beneficial properties is to synthesize inorganic/organic hybrids. These hybrids have polymers introduced into the sol-gel process so that the organic and inorganic components interact at the molecular level, providing control over mechanical properties and degradation rates. However, a full understanding of how each feature or property of the glass and hybrid scaffolds affects cellular response is needed to optimize the materials and ensure long-term success and clinical products. This review focuses on the techniques that have been developed for characterizing the hierarchical structures of sol-gel glasses and hybrids, from atomicscale amorphous networks, through the covalent bonding between components in hybrids and nanoporosity, to quantifying open macroporous networks of the scaffolds. Methods for non-destructive in situ monitoring of degradation and bioactivity mechanisms of the materials are also included. © 2012 The Royal Society.

T-cell epitope prediction and immune complex simulation using molecular dynamics:state of the art and persisting challenges

Atomistic Molecular Dynamics provides powerful and flexible tools for the prediction and analysis of molecular and macromolecular systems. Specifically, it provides a means by which we can measure theoretically that which cannot be measured experimentally: the dynamic time-evolution of complex systems comprising atoms and molecules. It is particularly suitable for the simulation and analysis of the otherwise inaccessible details of MHC-peptide interaction and, on a larger scale, the simulation of the immune synapse. Progress has been relatively tentative yet the emergence of truly high-performance computing and the development of coarse-grained simulation now offers us the hope of accurately predicting thermodynamic parameters and of simulating not merely a handful of proteins but larger, longer simulations comprising thousands of protein molecules and the cellular scale structures they form. We exemplify this within the context of immunoinformatics.

Toward the atomistic simulation of T cell epitopes automated construction of MHC:peptide structures for free energy calculations

Epitopes mediated by T cells lie at the heart of the adaptive immune response and form the essential nucleus of anti-tumour peptide or epitope-based vaccines. Antigenic T cell epitopes are mediated by major histocompatibility complex (MHC) molecules, which present them to T cell receptors. Calculating the affinity between a given MHC molecule and an antigenic peptide using experimental approaches is both difficult and time consuming, thus various computational methods have been developed for this purpose. A server has been developed to allow a structural approach to the problem by generating specific MHC:peptide complex structures and providing configuration files to run molecular modelling simulations upon them. A system has been produced which allows the automated construction of MHC:peptide structure files and the corresponding configuration files required to execute a molecular dynamics simulation using NAMD. The system has been made available through a web-based front end and stand-alone scripts. Previous attempts at structural prediction of MHC:peptide affinity have been limited due to the paucity of structures and the computational expense in running large scale molecular dynamics simulations. The MHCsim server (http://igrid-ext.cryst.bbk.ac.uk/MHCsim) allows the user to rapidly generate any desired MHC:peptide complex and will facilitate molecular modelling simulation of MHC complexes on an unprecedented scale.

Innate recognition of apoptotic cells:novel apoptotic cell-associated molecular patterns revealed by crossreactivity of anti-LPS antibodies

Cells dying by apoptosis are normally cleared by phagocytes through mechanisms that can suppress inflammation and immunity. Molecules of the innate immune system, the pattern recognition receptors (PRRs), are able to interact not only with conserved structures on microbes (pathogen-associated molecular patterns, PAMPs) but also with ligands displayed by apoptotic cells. We reasoned that PRRs might therefore interact with structures on apoptotic cells-apoptotic cell-associated molecular patterns (ACAMPs)-that are analogous to PAMPs. Here we show that certain monoclonal antibodies raised against the prototypic PAMP, lipopolysaccharide (LPS), can crossreact with apoptotic cells. We demonstrate that one such antibody interacts with a constitutively expressed intracellular protein, laminin-binding protein, which translocates to the cell surface during apoptosis and can interact with cells expressing the prototypic PRR, mCD14 as well as with CD14-negative cells. Anti-LPS cross reactive epitopes on apoptotic cells colocalised with annexin V-and C1q-binding sites on vesicular regions of apoptotic cell surfaces and were released associated with apoptotic cell-derived microvesicles (MVs). These results confirm that apoptotic cells and microbes can interact with the immune system through common elements and suggest that anti-PAMP antibodies could be used strategically to characterise novel ACAMPs associated not only with apoptotic cells but also with derived MVs. © 2013 Macmillan Publishers Limited All rights reserved.

Copper(II) complexes with unsymmetrical pentadentate ed3a-type diamino-tricarboxylate ligands. Crystal structures, configurational analysis and DFT study of complexes

The O–O–N–N–O-type pentadentate ligands H3ed3a, H3pd3a and H3pd3p (H3ed3a stands ethylenediamine-N,N,N′-triacetic acid; H3pd3a stands 1,3-propanediamine-N,N,N′-triacetic acid and H3pd3p stands 1,3-propanediamine-N,N,N′-tri-3-propionic acid) and the corresponding novel octahedral or square-planar/trigonal-bipyramidal copper(II) complexes have been prepared and characterized. H3ed3a, H3pd3a and H3pd3p ligands coordinate to copper(II) ion via five donor atoms (three deprotonated carboxylate atoms and two amine nitrogens) affording octahedral in case of ed3a3− and intermediate square-pyramidal/trigonal-bipyramidal structure in case of pd3a3− and pd3p3−. A six coordinate, octahedral geometry has been established crystallographically for the [Mg(H2O)6][Cu(ed3a)(H2O)]2 · 2H2O complex and five coordinate square-pyramidal for the [Mg(H2O)5Cu(pd3a)][Cu(pd3a)] · 2H2O. Structural data correlating similar chelate Cu(II) complexes have been used for the better understanding the pathway: octahedral → square-pyramidal ↔ trigonal- bipyramid geometry. An extensive configuration analysis is discussed in relation to information obtained for similar complexes. The infra-red and electronic absorption spectra of the complexes are discussed in comparison with related complexes of known geometries. Molecular mechanics and density functional theory (DFT) programs have been used to model the most stable geometric isomer yielding, at the same time, significant structural data. The results from density functional studies have been compared with X-ray data.