A molecular dynamics model of the atomic structure of dysprosium alumino-phosphate glass

Molecular dynamics (MD) has been used to identify the relative distribution of dysprosium in the phosphate glass DyAl0.30P3.05O9.62. The MD model has been compared directly with experimental data obtained from neutron diffraction to enable a detailed comparison beyond the total structure factor level. The MD simulation gives Dy ... Dy correlations at 3.80(5) and 6.40(5) angstrom with relative coordination numbers of 0.8(1) and 7.3(5), thus providing evidence of minority rare-earth clustering within these glasses. The nearest neighbour Dy-O peak occurs at 2.30 angstrom with each Dy atom having on average 5.8 nearest neighbour oxygen atoms. The MD simulation is consistent with the phosphate network model based on interlinked PO4 tetrahedra where the addition of network modifiers Dy3+ depolymerizes the phosphate network through the breakage of P-(O)-P bonds whilst leaving the tetrahedral units intact. The role of aluminium within the network has been taken into explicit account, and A1 is found to be predominantly (78 tetrahedrally coordinated. In fact all four A1 bonds are found to be to P (via an oxygen atom) with negligible amounts of Al-O-Dy bonds present. This provides an important insight into the role of Al additives in improving the mechanical properties of these glasses.

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.

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).

Approximation of photonic crystal fibres with large air holes by the step index fibre model

An equivalent step index fibre with a silica core and air cladding is used to model photonic crystal fibres with large air holes. We model this fibre for linear polarisation (we focus on the lowest few transverse modes of the electromagnetic field). The equivalent step index radius is obtained by equating the lowest two eigenvalues of the model to those calculated numerically for the photonic crystal fibres. The step index parameters thus obtained can then be used to calculate nonlinear parameters like the nonlinear effective area of a photonic crystal fibre or to model nonlinear few-mode interactions using an existing model.

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.

Evidence-based nanoscopic and molecular framework for excipient functionality in compressed orally disintegrating tablets

The work investigates the adhesive/cohesive molecular and physical interactions together with nanoscopic features of commonly used orally disintegrating tablet (ODT) excipients microcrystalline cellulose (MCC) and D-mannitol. This helps to elucidate the underlying physico-chemical and mechanical mechanisms responsible for powder densification and optimum product functionality. Atomic force microscopy (AFM) contact mode analysis was performed to measure nano-adhesion forces and surface energies between excipient-drug particles (6-10 different particles per each pair). Moreover, surface topography images (100 nm2-10 μm2) and roughness data were acquired from AFM tapping mode. AFM data were related to ODT macro/microscopic properties obtained from SEM, FTIR, XRD, thermal analysis using DSC and TGA, disintegration testing, Heckel and tabletability profiles. The study results showed a good association between the adhesive molecular and physical forces of paired particles and the resultant densification mechanisms responsible for mechanical strength of tablets. MCC micro roughness was 3 times that of D-mannitol which explains the high hardness of MCC ODTs due to mechanical interlocking. Hydrogen bonding between MCC particles could not be established from both AFM and FTIR solid state investigation. On the contrary, D-mannitol produced fragile ODTs due to fragmentation of surface crystallites during compression attained from its weak crystal structure. Furthermore, AFM analysis has shown the presence of extensive micro fibril structures inhabiting nano pores which further supports the use of MCC as a disintegrant. Overall, excipients (and model drugs) showed mechanistic behaviour on the nano/micro scale that could be related to the functionality of materials on the macro scale. © 2014 Al-khattawi et al.

Molecular Dynamics implementation of BN2D or 'Mercedes Benz' water model

Two-dimensional 'Mercedes Benz' (MB) or BN2D water model (Naim, 1971) is implemented in Molecular Dynamics. It is known that the MB model can capture abnormal properties of real water (high heat capacity, minima of pressure and isothermal compressibility, negative thermal expansion coefficient) (Silverstein et al., 1998). In this work formulas for calculating the thermodynamic, structural and dynamic properties in microcanonical (NVE) and isothermal-isobaric (NPT) ensembles for the model from Molecular Dynamics simulation are derived and verified against known Monte Carlo results. The convergence of the thermodynamic properties and the system's numerical stability are investigated. The results qualitatively reproduce the peculiarities of real water making the model a visually convenient tool that also requires less computational resources, thus allowing simulations of large (hydrodynamic scale) molecular systems. We provide the open source code written in C/C++ for the BN2D water model implementation using Molecular Dynamics.

Multiscale molecular dynamics/hydrodynamics implementation of two dimensional “Mercedes Benz” water model

A multiscale Molecular Dynamics/Hydrodynamics implementation of the 2D Mercedes Benz (MB or BN2D) [1] water model is developed and investigated. The concept and the governing equations of multiscale coupling together with the results of the two-way coupling implementation are reported. The sensitivity of the multiscale model for obtaining macroscopic and microscopic parameters of the system, such as macroscopic density and velocity fluctuations, radial distribution and velocity autocorrelation functions of MB particles, is evaluated. Critical issues for extending the current model to large systems are discussed.