The infrared spectrum, equilibrium structure and harmonic and anharmonic force field of thioborine, HBS

Infrared spectra of the two stretching fundamentals of both HBS and DBS have been observed, using a continuous flow system through a multiple reflection long path cell at a pressure around 1 Torr and a Nicolet Fourier Transform spectrometer with a resolution of about 0•1 cm-1. The v3 BS stretching fundamental of DBS, near 1140 cm-1, is observed in strong Fermi resonance with the overtone of the bend 2v2. The bending fundamental v2 has not been observed and must be a very weak band. The analysis of the results in conjunction with earlier work gives the equilibrium structure (re(BH) = 1•1698(12) , re(BS) = 1•5978(3) ) and the harmonic and anharmonic force field.

Intensities in the v7, v8 and v9 bands of CH2NH, and the harmonic force field of methyleneimine

A high resolution Fourier transform infrared spectrum of methyleneimine, HN=CH2, has been obtained in the gas phase in the region 700 to 1300 cm−1. The rovibrational line intensities of the three lowest fundamentals ν7 (A′), ν8 (A″), and ν9 (A″) have been simulated including all Coriolis interactions between the three bands, and by fitting the observed spectrum the relative signs and magnitudes of the vibrational transition moments have been determined. All of the available spectroscopic data have been used to determine the harmonic force field of methyleneimine.

The general harmonic force field of methyl fluoride

The complete general harmonic force field of methyl flouride was recalculated using the most recent literature frequency, Coriolis ζ, and centrifugal distortion data for 12CH3F, 13CH3F, 12CD3F, 12CHD2F and 12CH2DF. The anharmonic corrections applied to the observed frequency data and the adopted molecular geometry are considered to be more realistic than those used hitherto. There is excellent overall agreement between the fitted force constants and the highest quality ab initio force field currently available.

Molecular modelling studies of binding of DACD derivatives into G-Quadruplex DNA: comparison of force field and quantum polarized ligand docking methods

The DNA G-qadruplexes are one of the targets being actively explored for anti-cancer therapy by inhibiting them through small molecules. This computational study was conducted to predict the binding strengths and orientations of a set of novel dimethyl-amino-ethyl-acridine (DACA) analogues that are designed and synthesized in our laboratory, but did not diffract in Synchrotron light.Thecrystal structure of DNA G-Quadruplex(TGGGGT)4(PDB: 1O0K) was used as target for their binding properties in our studies.We used both the force field (FF) and QM/MM derived atomic charge schemes simultaneously for comparing the predictions of drug binding modes and their energetics. This study evaluates the comparative performance of fixed point charge based Glide XP docking and the quantum polarized ligand docking schemes. These results will provide insights on the effects of including or ignoring the drug-receptor interfacial polarization events in molecular docking simulations, which in turn, will aid the rational selection of computational methods at different levels of theory in future drug design programs. Plenty of molecular modelling tools and methods currently exist for modelling drug-receptor or protein-protein, or DNA-protein interactionssat different levels of complexities.Yet, the capasity of such tools to describevarious physico-chemical propertiesmore accuratelyis the next step ahead in currentresearch.Especially, the usage of most accurate methods in quantum mechanics(QM) is severely restricted by theirtedious nature. Though the usage of massively parallel super computing environments resulted in a tremendous improvement in molecular mechanics (MM) calculations like molecular dynamics,they are still capable of dealing with only a couple of tens to hundreds of atoms for QM methods. One such efficient strategy that utilizes thepowers of both MM and QM are the QM/MM hybrid methods. Lately, attempts have been directed towards the goal of deploying several different QM methods for betterment of force field based simulations, but with practical restrictions in place. One of such methods utilizes the inclusion of charge polarization events at the drug-receptor interface, that is not explicitly present in the MM FF.

Optimal space manoeuvres in a non-Keplerian force field

In this work the problem of a spacecraft bi-impulsive transfer between two given non coplanar elliptical orbits, with minimum fuel consumption, is solved considering a non-Keplerian force field (the perturbing forces include Earth gravity harmonics and atmospheric drag). The problem is transformed in the Two Point Boundary Value Problem. It is developed and implemented a new algorithm, that uses the analytical expressions developed here. A dynamics that considered a Keplerian force field was used to produce an initial guess to solve the Two Point Boundary Value Problem. Several simulations were performed to observe the spacecraft orbital behaviour by different kind of perturbations and constraints, on a fuel consumption optimization point of view. (C) 2002 COSPAR. Published by Elsevier B.V. Ltd. All rights reserved.

Ramachandran-type Plots for Glycosidic Linkages: Examples from Molecular Dynamics Simulations using the Glycam06 Force Field

The goals of this article are to (1) provide further validation of the Glycam06 force field, specifically for its use in implicit solvent molecular dynamic (MD) simulations, and (2) to present the extension of G.N. Ramachandran's idea of plotting amino acid phi and psi angles to the glycosidic phi, psi, and omega angles formed between carbohydrates. As in traditional Ramachandran plots, these carbohydrate Ramachandran-type (carb-Rama) plots reveal the coupling between the glycosidic angles by displaying the allowed and disallowed conformational space. Considering two-bond glycosidic linkages, there are 18 possible conformational regions that can be defined by (α, ϕ, ψ) and (β, ϕ, ψ), whereas for three-bond linkages, there are 54 possible regions that can be defined by (α, ϕ, ψ, ω) and (β, ϕ, ψ, ω). Illustrating these ideas are molecular dynamic simulations on an implicitly hydrated oligosaccharide (700 ns) and its eight constituent disaccharides (50 ns/disaccharide). For each linkage, we compare and contrast the oligosaccharide and respective disaccharide carb-Rama plots, validate the simulations and the Glycam06 force field through comparison to experimental data, and discuss the general trends observed in the plots.

Energy-based de novo protein folding by conformational space annealing and an off-lattice united-residue force field: Application to the 10-55 fragment of staphylococcal protein A and to apo calbindin D9K

The conformational space annealing (CSA) method for global optimization has been applied to the 10-55 fragment of the B-domain of staphylococcal protein A (protein A) and to a 75-residue protein, apo calbindin D9K (PDB ID code 1CLB), by using the UNRES off-lattice united-residue force field. Although the potential was not calibrated with these two proteins, the native-like structures were found among the low-energy conformations, without the use of threading or secondary-structure predictions. This is because the CSA method can find many distinct families of low-energy conformations. Starting from random conformations, the CSA method found that there are two families of low-energy conformations for each of the two proteins, the native-like fold and its mirror image. The CSA method converged to the same low-energy folds in all cases studied, as opposed to other optimization methods. It appears that the CSA method with the UNRES force field, which is based on the thermodynamic hypothesis, can be used in prediction of protein structures in real time.

Cost-Effective Force Field Tailored for Solid-Phase Simulations of OLED Materials

A united atom force field is empirically derived by minimizing the difference between experimental and simulated crystal cells and melting temperatures for eight compounds representative of organic electronic materials used in OLEDs and other devices: biphenyl, carbazole, fluorene, 9,9′-(1,3-phenylene)bis(9H-carbazole)-1,3-bis(N-carbazolyl)benzene (mCP), 4,4′-bis(N-carbazolyl)-1,1′-biphenyl (pCBP), phenazine, phenylcarbazole, and triphenylamine. The force field is verified against dispersion-corrected DFT calculations and shown to also successfully reproduce the crystal structure for two larger compounds employed as hosts in phosphorescent and thermally activated delayed fluorescence OLEDs: N,N′-di(1-naphthyl)-N,N′-diphenyl-(1,1′-biphenyl)-4,4′-diamine (NPD), and 1,3,5-tri(1-phenyl-1H-benzo[d]imidazol-2-yl)phenyl (TPBI). The good performances of the force field coupled to the large computational savings granted by the united atom approximation make it an ideal choice for the simulation of the morphology of emissive layers for OLED materials in crystalline or glassy phases.

Optical force field mapping in microdevices

We present a method for characterizing microscopic optical force fields. Two dimensional vector force maps are generated by measuring the optical force applied to a probe particle for a grid of particle positions. The method is used to map Out the force field created by the beam from a lensed fiber inside a liquid filled microdevice. We find transverse gradient forces and axial scattering forces on the order of 2 pN per 10 mW laser power which are constant over a considerable axial range (> 35 mu m). These findings suggest Future useful applications of lensed fibers for particle guiding/sorting. The propulsion of a small particle at a constant velocity of 200 mu m s(-1) is shown.

Simulazioni atomiche di canali ionici del potassio: limiti del Force Field AMBER nelle stime di conduttanza

Le simulazioni atomiche dei canali ionici del K+ mettono in evidenza una significativa inaccuratezza: i valori di conduttanza del canale stimati sono molto distanti dai valori sperimentali. Per tentare di migliorare le simulazioni e i loro risultati, è necessario modificare il modello di Force Field impiegato. Si è scelto di adottare il modello di Force Field AMBER, uno dei più utilizzati per simulazioni di molecole biologiche in dinamica molecolare, e di agire sui valori dei parametri di Lennard-Jones cercando un set di valori che permetta di minimizzare la distanza tra la conduttanza stimata e la conduttanza sperimentale delle soluzioni saline. L’analisi ha dimostrato che non è possibile riprodurre in modo quantitativo i valori sperimentali di conduttanza di soluzioni saline di Na+ e K+ in simulazioni atomiche basate su un modello di Force Field additivo.

Molecular Dynamics Investigations of PRODAN in a DLPC Bilayer

Molecular dynamics computer simulations have been performed to identify preferred positions of the fluorescent probe PRODAN in a fully hydrated DLPC bilayer in the fluid phase. In addition to the intramolecular charge-transfer first vertical excited state, we considered different charge distributions for the electronic ground state of the PRODAN molecule by distinct atomic charge models corresponding to the probe molecule in vacuum as well as polarized in a weak and a strong dielectric solvent (cyclohexane and water). Independent on the charge distribution model of PRODAN, we observed a preferential orientation of this molecule in the bilayer with the dimethylamino group pointing toward the membrane's center and the carbonyl oxygen toward the membrane's interface. However, changing the charge distribution model of PRODAN, independent of its initial position in the equilibrated DLPC membrane, we observed different preferential positions. For the ground state representation without polarization and the in-cyclohexane polarization, the probe maintains its position close to the membrane's center. Considering the in-water polarization model, the probe approaches more of the polar headgroup region of the bilayer, with a strong structural correlation with the choline group, exposing its oxygen atom to water molecules. PRODAN's representation of the first vertical excited state with the in-water polarization also approaches the polar region of the membrane with the oxygen atom exposed to the bilayer's hydration shell. However, this model presents a stronger structural correlation with the phosphate groups than the ground state. Therefore, we conclude that the orientation of the PRODAN molecule inside the DLPC membrane is well-defined, but its position is very sensitive to the effect of the medium polarization included here by different models for the atomic charge distribution of the probe.

Modéliser la polarisation électronique par un continuum diélectrique intramoléculaire vers un champ de force polarisable pour la chimie bioorganique

Thèse numérisée par la Division de la gestion de documents et des archives de l'Université de Montréal

Studies on the dynamics and rheology of dilute suspensions of Brownian particles in simple shear flow under the action of an external force

We present a novel approach to computing the orientation moments and rheological properties of a dilute suspension of spheroids in a simple shear flow at arbitrary Peclct number based on a generalised Langevin equation method. This method differs from the diffusion equation method which is commonly used to model similar systems in that the actual equations of motion for the orientations of the individual particles are used in the computations, instead of a solution of the diffusion equation of the system. It also differs from the method of 'Brownian dynamics simulations' in that the equations used for the simulations are deterministic differential equations even in the presence of noise, and not stochastic differential equations as in Brownian dynamics simulations. One advantage of the present approach over the Fokker-Planck equation formalism is that it employs a common strategy that can be applied across a wide range of shear and diffusion parameters. Also, since deterministic differential equations are easier to simulate than stochastic differential equations, the Langevin equation method presented in this work is more efficient and less computationally intensive than Brownian dynamics simulations.We derive the Langevin equations governing the orientations of the particles in the suspension and evolve a procedure for obtaining the equation of motion for any orientation moment. A computational technique is described for simulating the orientation moments dynamically from a set of time-averaged Langevin equations, which can be used to obtain the moments when the governing equations are harder to solve analytically. The results obtained using this method are in good agreement with those available in the literature.The above computational method is also used to investigate the effect of rotational Brownian motion on the rheology of the suspension under the action of an external force field. The force field is assumed to be either constant or periodic. In the case of con- I stant external fields earlier results in the literature are reproduced, while for the case of periodic forcing certain parametric regimes corresponding to weak Brownian diffusion are identified where the rheological parameters evolve chaotically and settle onto a low dimensional attractor. The response of the system to variations in the magnitude and orientation of the force field and strength of diffusion is also analyzed through numerical experiments. It is also demonstrated that the aperiodic behaviour exhibited by the system could not have been picked up by the diffusion equation approach as presently used in the literature.The main contributions of this work include the preparation of the basic framework for applying the Langevin method to standard flow problems, quantification of rotary Brownian effects by using the new method, the paired-moment scheme for computing the moments and its use in solving an otherwise intractable problem especially in the limit of small Brownian motion where the problem becomes singular, and a demonstration of how systems governed by a Fokker-Planck equation can be explored for possible chaotic behaviour.