A dielectric continuum molecular dynamics method

We introduce a novel method to simulate hydrated macromolecules with a dielectric continuum representation of the surrounding solvent. In our approach, the interaction between the solvent and the molecular degrees of freedom is described by means of a polarization density free energy functional which is minimum at electrostatic equilibrium. After a pseudospectral expansion of the polarization and a discretization of the functional, we construct the equations of motion for the system based on a Car-Parrinello technique. In the limit of the adiabatic evolution of the polarization field variables, our method provides the solution of the dielectric continuum problem "on the fly," while the molecular coordinates are propagated. In this first study, we show how our dielectric continuum molecular dynamics method can be successfully applied to hydrated biomolecules, with low cost compared to free energy simulations with explicit solvent. To our knowledge, this is the first time that stable and conservative molecular dynamic simulations of solutes can be performed for a dielectric continuum model of the solvent. (C) 2001 American Institute of Physics.

Thin Ta films: growth, stability, and diffusion studied by molecular-dynamics simulations

We report results of classical molecular-dynamics simulations of bcc and beta-Ta thin films. Thermal PVD film growth, surface roughness, argon ion bombardment, phase stability and transformation, vacancy and adatom diffusion, and thermal relaxation kinetics are discussed. Distinct differences between the two structures are observed, including a complex vacancy diffusion mechanism in beta-Ta. Embedded atom method potentials, which were fitted to bcc properties, have been used to model the Ta-Ta interactions. In order to verify the application of these potentials to the more complex beta-Ta structure, we have also performed density functional theory calculations. Results and implications of these calculations are discussed.

Liquid structure of the ionic liquid, 1-methyl-4-cyanopyridinium bis{(trifluoromethyl)sulfonyl}imide determined from neutron scattering and molecular dynamics simulations

The liquid structure of 1-methyl-4-cyanopyridinium bis {(trifluoromethyl)sulfonyl}imide, a prototypical ionic liquid containing an electron-withdrawing group on the cation, has been investigated at 368 K. Experimental neutron scattering combined with empirical potential structure refinement analysis of the data and classical molecular dynamics simulations have been used to probe the liquid structure in detail. Both techniques generated highly consistent results that provide valuable validation of the force fields and refinement approaches. A significant degree of apparent charge ordering is found in the liquid structure, although the nonspherical shape of the ions results in interpenetration of cations into the first shell of adjacent cations, with much shorter closest contact distances than the averaged center-of-mass cation-cation and cation-anion separations.

Using sketch-map coordinates to analyze and bias molecular dynamics simulations

When examining complex problems, such as the folding of proteins, coarse grained descriptions of the system drive our investigation and help us to rationalize the results. Oftentimes collective variables (CVs), derived through some chemical intuition about the process of interest, serve this purpose. Because finding these CVs is the most difficult part of any investigation, we recently developed a dimensionality reduction algorithm, sketch-map, that can be used to build a low-dimensional map of a phase space of high-dimensionality. In this paper we discuss how these machine-generated CVs can be used to accelerate the exploration of phase space and to reconstruct free-energy landscapes. To do so, we develop a formalism in which high-dimensional configurations are no longer represented by low-dimensional position vectors. Instead, for each configuration we calculate a probability distribution, which has a domain that encompasses the entirety of the low-dimensional space. To construct a biasing potential, we exploit an analogy with metadynamics and use the trajectory to adaptively construct a repulsive, history-dependent bias from the distributions that correspond to the previously visited configurations. This potential forces the system to explore more of phase space by making it desirable to adopt configurations whose distributions do not overlap with the bias. We apply this algorithm to a small model protein and succeed in reproducing the free-energy surface that we obtain from a parallel tempering calculation.

Interactions between pyrazole derived enantiomers and Chiralcel OJ: Prediction of enantiomer absolute configurations and elution order by molecular dynamics simulations

The separation of enantiomers and confirmation of their absolute configurations is significant in the development of chiral drugs. The interactions between the enantiomers of chiral pyrazole derivative and polysaccharide-based chiral stationary phase cellulose tris(4-methylbenzoate) (Chiralcel OJ) in seven solvents and under different temperature were studied using molecular dynamics simulations. The results show that solvent effect has remarkable influence on the interactions. Structure analysis discloses that the different interactions between two isomers and chiral stationary phase are dependent on the nature of solvents, which may invert the elution order. The computational method in the present study can be used to predict the elution order and the absolute configurations of enantiomers in HPLC separations and therefore would be valuable in development of chiral drugs.

Can the scaling behavior of electric conductivity be used to probe the self-organizational changes in solution with respect to the ionic liquid structure? The case of [C8MIM][NTf2]

Electrical conductivity of the supercooled ionic liquid [C8MIM][NTf2], determined as a function of temperature and pressure, highlights strong differences in its ionic transport behavior between low and high temperature regions. To date, the crossover effect which is very well known for low molecular van der Waals liquids has been rarely described for classical ionic liquids. This finding highlights that the thermal fluctuations could be dominant mechanisms driving the dramatic slowing down of ion motions near Tg. An alternative way to analyze separately low and high temperature dc-conductivity data using a density scaling approach was then proposed. Based on which a common value of the scaling exponent [gamma] = 2.4 was obtained, indicating that the applied density scaling is insensitive to the crossover effect. By comparing the scaling exponent [gamma] reported herein along with literature data for other ionic liquids, it appears that [gamma] decreases by increasing the alkyl chain length on the 1-alkyl-3-methylimidazolium-based ionic liquids. This observation may be related to changes in the interaction between ions in solution driven by an increase in the van der Waals type interaction by increasing the alkyl chain length on the cation. This effect may be related to changes in the ionic liquid nanostructural organization with the alkyl chain length on the cation as previously reported in the literature based on molecular dynamic simulations. In other words, the calculated scaling exponent [gamma] may be then used as a key parameter to probe the interaction and/or self-organizational changes in solution with respect to the ionic liquid structure.

Atomistic Insights into Rhodopsin Activation from a Dynamic Model

Rhodopsin, the light sensitive receptor responsible for blue-green vision, serves as a prototypical G protein-coupled receptor (GPCR). Upon light absorption, it undergoes a series of conformational changes that lead to the active form, metarhodopsin II (META II), initiating a signaling cascade through binding to the G protein transducin (G(t)). Here, we first develop a structural model of META II by applying experimental distance restraints to the structure of lumi-rhodopsin (LUMI), an earlier intermediate. The restraints are imposed by using a combination of biased molecular dynamics simulations and perturbations to an elastic network model. We characterize the motions of the transmembrane helices in the LUMI-to-META II transition and the rearrangement of interhelical hydrogen bonds. We then simulate rhodopsin activation in a dynamic model to study the path leading from LUMI to our META II model for wild-type rhodopsin and a series of mutants. The simulations show a strong correlation between the transition dynamics and the pharmacological phenotypes of the mutants. These results help identify the molecular mechanisms of activation in both wild type and mutant rhodopsin. While static models can provide insights into the mechanisms of ligand recognition and predict ligand affinity, a dynamic model of activation could be applicable to study the pharmacology of other GPCRs and their ligands, offering a key to predictions of basal activity and ligand efficacy.

Addressing Selective Polypharmacology of Antipsychotic Drugs Targeting the Bioaminergic Receptors through the Receptor Dynamic Conformational Ensembles

Selective polypharmacology, where a drug acts on multiple rather than single molecular targets involved in a disease, emerges to develop a structure-based system biology approach to design drugs selectively targeting a disease-active protein network. We focus on the bioaminergic receptors that belong to the group of integral membrane signalling proteins coupled to the G protein and represent targets for therapeutic agents against schizophrenia and depression. Among them, it has been shown that the serotonin (5-HT2A and 5-HT6), dopamine (D2 and D3) receptors induce a cognition-enhancing effect (group 1), while the histamine (H1) and serotonin (5-HT2C) receptors lead to metabolic side effects and the 5-HT2B serotonin receptor causes pulmonary hypertension (group 2). Thus, the problem arises to develop an approach that allows identifying drugs targeting only the disease-active receptors, i.e. group 1. The recent release of several crystal structures of the bioaminergic receptors, involving the D3 and H1 receptors provides the possibility to model the structures of all receptors and initiate a study of the structural and dynamic context of selective polypharmacology. In this work, we use molecular dynamics simulations to generate a conformational space of the receptors and subsequently characterize its binding properties applying molecular probe mapping. All-against-all comparison of the generated probe maps of the selected diverse conformations of all receptors with the Tanimoto similarity coefficient (Tc) enable to separate the receptors of group 1 from group 2. The pharmacophore built based on the Tc-selected receptor conformations, using the multiple probe maps discovers structural features that can be used to design molecules selective towards the receptors of group 1. The importance of several predicted residues to ligand selectivity is supported by the available mutagenesis and ligand structure-activity relationships studies. In addition, the Tc-selected conformations of the receptors for group 1 show good performance in isolation of known ligands from a random decoy. Our computational structure-based protocol to tackle selective polypharmacology of antipsychotic drugs could be applied for other diseases involving multiple drug targets, such as oncologic and infectious disorders.

Ab initio molecular dynamics simulation of a room temperature ionic liquid

Ab initio molecular dynamics simulations have been performed for the first time on the room-temperature organic ionic liquid dimethyl imidazolium chloride [DMIM][Cl] using density functional theory. The aim is to compare the local liquid structure with both that obtained from two different classical force fields and from neutron scattering experiments. The local structure around the cation shows significant differences compared to both the classical calculations and the neutron results. In particular, and unlike in the gas-phase ion pair, chloride ions tend to be located near a ring C-H proton in a position suggesting hydrogen bonding. The results are used to suggest ways in which the classical potentials may be improved.

On the structure and dynamics of ionic liquids

The structure and dynamics of the ionic liquid 1-ethyl-3-methylimidazolium nitrate is studied by molecular dynamics simulations. We find long-range spatial correlations between the ions and a three-dimensional local structure that reflects the asymmetry of the cations. The main contribution to the configurational energy comes from the electrostatic interactions which leads to charge-ordering effects. Radial screening and threedimensional distribution of charge are also analyzed. The motion of a single ion is studied via velocity and reorientational correlation functions. It is found that ions "rattle" in a long-lived cage, while the orientational structure relaxes on a time scale longer than 200 ps. As in a supercooled liquid, the mean square displacements reveal a subdiffusive dynamics. In addition, the presence of dynamic heterogeneities can be detected by analyzing the non-Gaussian behavior of the van Hove correlation function and the spatial arrangement of the most mobile ions. The short-time collective dynamics is also studied through the electric current time correlation function.