Molecular dynamics simulations of proteins: can the explicit water model be varied?

In molecular mechanics simulations of biological systems, the solvation water is typically represented by a default water model which is an integral part of the force field. Indeed, protein nonbonding parameters are chosen in order to obtain a balance between water-water and protein-water interactions and hence a reliable description of protein solvation. However, less attention has been paid to the question of whether the water model provides a reliable description of the water properties under the chosen simulation conditions, for which more accurate water models often exist. Here we consider the case of the CHARMM protein force field, which was parametrized for use with a modified TIP3P model. Using quantum mechanical and molecular mechanical calculations, we investigate whether the CHARMM force field can be used with other water models: TIP4P and TIP5P. Solvation properties of N-methylacetamide (NMA), other small solute molecules, and a small protein are examined. The results indicate differences in binding energies and minimum energy geometries, especially for TIP5P, but the overall description of solvation is found to be similar for all models tested. The results provide an indication that molecular mechanics simulations with the CHARMM force field can be performed with water models other than TIP3P, thus enabling an improved description of the solvent water properties.

Metal complexes of a tetraazacyclophane: solution and molecular modelling studies

An alternative synthetic approach to yield the compound 2,3,5,6,8,9,11,14-octahydrobenzo[1][ 1,4,7,10]tetraazacyclotetradecine (bz[14]N-4) is presented. The protonation constants of bz[14]N-4 and the stability constants of its complexes with Ni2+, Cu2+, Zn2+, Cd2+, Pb2+ were determined in H2O at 25degreesC with ionic strength 0.10 mol dm(-3) in KNO3 and they were compared with structurally related macrocycles cyclam (1,4,8,11-tetraazacyclotetradecane) and cyclen (1,4,7,10-tetraazacyclododecane). These studies indicate that only 1 : 1 ( M : L) species are formed in solution, and the ligand exhibits a high affinity for larger ions such as Cd2+ and Pb2+. The X-ray study of [bz[14]N4H3](3+) shows that an inclusion compound with a chloride counter-anion is formed through NH...Cl hydrogen bonds. Spectroscopic data in solution ( electronic and NMR spectra) showed that the macrocycle adopts a planar arrangement upon metal complexation. Molecular mechanics calculations reveal that in spite of the presence of the benzene ring in the macrocyclic framework this ligand can encapsulate metal ions with different stereo-electronic sizes in square planar arrangements. Our results indicate that the presence of the benzene ring in the backbone of the bz[14]N-4 confers a coordination behaviour intermediate between that of cyclam and cyclen.

Crystal and molecular structure of the thioether-analogue of PEEK from X-ray powder data and diffraction-modelling

Analysis of X-ray powder data for the melt-crystallisable aromatic poly(thioether thioether ketone) [-S-Ar-S-Ar-CO-Ar](n), ('PTTK', Ar= 1,4-phenylene), reveals that it adopts a crystal structure very different from that established for its ether-analogue PEEK. Molecular modelling and diffraction-simulation studies of PTTK show that the structure of this polymer is analogous to that of melt-crystallised poly(thioetherketone) [-SAr-CO-Ar](n) in which the carbonyl linkages in symmetry-related chains are aligned anti-parallel to one another. and that these bridging units are crystallographically interchangeable. The final model for the crystal structure of PTTK is thus disordered, in the monoclinic space group 121a (two chains per unit cell), with cell dimensions a = 7.83, b = 6.06, c = 10.35 angstrom, beta = 93.47 degrees. (c) 2005 Elsevier Ltd. All rights reserved.

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.

Applicability of molecular simulations for modelling the adsorption of the greenhouse gas CF4 on carbons

Tetrafluoromethane, CF4, is powerful greenhouse gas, and the possibility of storing it in microporous carbon has been widely studied. In this paper we show, for the first time, that the results of molecular simulations can be very helpful in the study of CF4 adsorption. Moreover, experimental data fit to the results collected from simulations. We explain the meaning of the empirical parameters of the supercritical Dubinin–Astakhov model proposed by Ozawa and finally the meaning of the parameter k of the empirical relation proposed by Amankwah and Schwarz.

Modelling habitat selection and distribution of the critically endangered Jerdon's courser Rhinoptilus bitorquatus in scrub jungle: an application of a new tracking method

1. Jerdon's courser Rhinoptilus bitorquatus is a nocturnally active cursorial bird that is only known to occur in a small area of scrub jungle in Andhra Pradesh, India, and is listed as critically endangered by the IUCN. Information on its habitat requirements is needed urgently to underpin conservation measures. We quantified the habitat features that correlated with the use of different areas of scrub jungle by Jerdon's coursers, and developed a model to map potentially suitable habitat over large areas from satellite imagery and facilitate the design of surveys of Jerdon's courser distribution. 2. We used 11 arrays of 5-m long tracking strips consisting of smoothed fine soil to detect the footprints of Jerdon's coursers, and measured tracking rates (tracking events per strip night). We counted the number of bushes and trees, and described other attributes of vegetation and substrate in a 10-m square plot centred on each strip. We obtained reflectance data from Landsat 7 satellite imagery for the pixel within which each strip lay. 3. We used logistic regression models to describe the relationship between tracking rate by Jerdon's coursers and characteristics of the habitat around the strips, using ground-based survey data and satellite imagery. 4. Jerdon's coursers were most likely to occur where the density of large (>2 m tall) bushes was in the range 300-700 ha(-1) and where the density of smaller bushes was less than 1000 ha(-1). This habitat was detectable using satellite imagery. 5. Synthesis and applications. The occurrence of Jerdon's courser is strongly correlated with the density of bushes and trees, and is in turn affected by grazing with domestic livestock, woodcutting and mechanical clearance of bushes to create pasture, orchards and farmland. It is likely that there is an optimal level of grazing and woodcutting that would maintain or create suitable conditions for the species. Knowledge of the species' distribution is incomplete and there is considerable pressure from human use of apparently suitable habitats. Hence, distribution mapping is a high conservation priority. A two-step procedure is proposed, involving the use of ground surveys of bush density to calibrate satellite image-based mapping of potential habitat. These maps could then be used to select priority areas for Jerdon's courser surveys. The use of tracking strips to study habitat selection and distribution has potential in studies of other scarce and secretive species.

Spatiotemporal modelling of CheY complexes in Escherichia coli chemotaxis

The chemotaxis pathway of Escherichia coli is one of the best studied and modelled biological signalling pathways. Here we extend existing modelling approaches by explicitly including a description of the formation and subcellular localization of intermediary complexes in the phosphotransfer pathway. The inclusion of these complexes shows that only about 60% of the total output response regulator (CheY) is uncomplexed at any moment and hence free to interact with its target, the flagellar motor. A clear strength of this model is its ability to predict the experimentally observable subcellular localization of CheY throughout a chemotactic response. We have found good agreement between the model output and experimentally determined CheY localization patterns. (C) 2009 Elsevier Ltd. All rights reserved.

Modelling heterotachy in phylogenetic inference by reversible-jump Markov chain Monte Carlo

The rate at which a given site in a gene sequence alignment evolves over time may vary. This phenomenon-known as heterotachy-can bias or distort phylogenetic trees inferred from models of sequence evolution that assume rates of evolution are constant. Here, we describe a phylogenetic mixture model designed to accommodate heterotachy. The method sums the likelihood of the data at each site over more than one set of branch lengths on the same tree topology. A branch-length set that is best for one site may differ from the branch-length set that is best for some other site, thereby allowing different sites to have different rates of change throughout the tree. Because rate variation may not be present in all branches, we use a reversible-jump Markov chain Monte Carlo algorithm to identify those branches in which reliable amounts of heterotachy occur. We implement the method in combination with our 'pattern-heterogeneity' mixture model, applying it to simulated data and five published datasets. We find that complex evolutionary signals of heterotachy are routinely present over and above variation in the rate or pattern of evolution across sites, that the reversible-jump method requires far fewer parameters than conventional mixture models to describe it, and serves to identify the regions of the tree in which heterotachy is most pronounced. The reversible-jump procedure also removes the need for a posteriori tests of 'significance' such as the Akaike or Bayesian information criterion tests, or Bayes factors. Heterotachy has important consequences for the correct reconstruction of phylogenies as well as for tests of hypotheses that rely on accurate branch-length information. These include molecular clocks, analyses of tempo and mode of evolution, comparative studies and ancestral state reconstruction. The model is available from the authors' website, and can be used for the analysis of both nucleotide and morphological data.

Modelling habitat selection and distribution of the critically endangered Jerdon's courser Rhinoptilus bitorquatus in scrub jungle: an application of a new tracking method

1. Jerdon's courser Rhinoptilus bitorquatus is a nocturnally active cursorial bird that is only known to occur in a small area of scrub jungle in Andhra Pradesh, India, and is listed as critically endangered by the IUCN. Information on its habitat requirements is needed urgently to underpin conservation measures. We quantified the habitat features that correlated with the use of different areas of scrub jungle by Jerdon's coursers, and developed a model to map potentially suitable habitat over large areas from satellite imagery and facilitate the design of surveys of Jerdon's courser distribution. 2. We used 11 arrays of 5-m long tracking strips consisting of smoothed fine soil to detect the footprints of Jerdon's coursers, and measured tracking rates (tracking events per strip night). We counted the number of bushes and trees, and described other attributes of vegetation and substrate in a 10-m square plot centred on each strip. We obtained reflectance data from Landsat 7 satellite imagery for the pixel within which each strip lay. 3. We used logistic regression models to describe the relationship between tracking rate by Jerdon's coursers and characteristics of the habitat around the strips, using ground-based survey data and satellite imagery. 4. Jerdon's coursers were most likely to occur where the density of large (>2 m tall) bushes was in the range 300-700 ha(-1) and where the density of smaller bushes was less than 1000 ha(-1). This habitat was detectable using satellite imagery. 5. Synthesis and applications. The occurrence of Jerdon's courser is strongly correlated with the density of bushes and trees, and is in turn affected by grazing with domestic livestock, woodcutting and mechanical clearance of bushes to create pasture, orchards and farmland. It is likely that there is an optimal level of grazing and woodcutting that would maintain or create suitable conditions for the species. Knowledge of the species' distribution is incomplete and there is considerable pressure from human use of apparently suitable habitats. Hence, distribution mapping is a high conservation priority. A two-step procedure is proposed, involving the use of ground surveys of bush density to calibrate satellite image-based mapping of potential habitat. These maps could then be used to select priority areas for Jerdon's courser surveys. The use of tracking strips to study habitat selection and distribution has potential in studies of other scarce and secretive species.

Use of essential and molecular dynamics to study gamma B-crystallin unfolding after non-enzymic post-translational modifications

Essential and Molecular Dynamics (ED/MD) have been used to model the conformational changes of a protein implicated in a conformational disease-cataract, the largest cause of blindness in the world-after non-enzymic post-translational modification. Cyanate modification did not significantly alter flexibility, while the Schiff's base adduct produced a more flexible N-terminal domain, and intra-secondary structure regions, than either the cyanate adduct or the native structure. Glycation also increased linker flexibility and disrupted the charge network. A number of post-translational adducts showed structural disruption around Cys15 and increased linker flexibility; this may be important in subsequent protein aggregation. Our modelling results are in accord with experimental evidence, and show that ED/MD is a useful tool in modelling conformational changes in proteins implicated in disease processes. (C) 2003 Published by Elsevier Ltd.

Metal complexes of a dipyridine octaazamacrocycle: stability constants, structural and modelling studies

Two 28-membered octaazamacrocycles, [28]py(2)N(6) and Me-2[28]py(2)N(6), have been synthesized. The protonation constants of the N-methyl. derivative and the stability constants of its complexes with Ni2+, Cu2+, Zn2+, Cd2+, and Pb2+ were determined at 25degreesC in 0.10 mol dm(-3) KNO3. The high overall basicity of Me-2[28]py(2)N(6) is ascribed to the weaker repulsion between protonated contiguous charged ammonium sites separated by propyl chains. These studies together with NMR, UV-vis and EPR spectroscopies indicated the presence of mono- and di-nuclear species, The single crystal structure of the complex [Ni-2([28]py(2)N(6))(H2O)(4)]Cl-4.3H(2)O was determined, and showed each nickel centre in a distorted octahedral co-ordination environment. The nickel centres are held within the macrocycle at a large distance of 6.991(g) Angstrom from each other. The formation of mononuclear complexes was evaluated theoretically via molecular mechanics (MM) and molecular dynamics (MD) calculations and showed that these large macrocycles have sufficient flexibility to encapsulate metal ions with different stereo-electronic sizes. Structures for small and large metal ions are proposed.

Selectivity of bis(calix[4]diquinone) ionophores towards metal ions in solvent dimethylsulfoxide: a molecular mechanics and molecular dynamics study

Molecular modelling studies have been carried out on two bis(calix[4]diqu(inone) ionophores, each created from two (calix[4]diquinone)arenes bridged at their bottom rims via alkyl chains (CH2)(n), 1: n = 3, 2; n = 4, in order to understand the reported selectivity of these ligands towards different sized metal ions such as Na+, K+, Rb+, and Cs+ in dmso solution. Conformational. analyses have been carried out which show that in the lowest energy conformations of the two macrocycles, the individual calix[4]diquinones exhibit a combination of partial cone, 1,3-alternate and cone conformations. The interactions of these alkali metals with the macrocycles have been studied in the gas phase and in a periodic box of solvent dmso by molecular mechanics and molecular dynamics calculations. Molecular mechanics calculations have been carried out on the mode of entry of the ions into the macrocycles and suggest that this is likely to occur from the side of the central cavity, rather than through the main axis of the calix[4]diquinones. There are energy barriers of ca. 19 kcal mol(-1) for this entry path in the gas phase, but in solution no energy barrier is found. Molecular dynamics simulations show that in both 1 and 2, though particularly in the latter macrocycle, one or two solvent molecules are bonded to the metal throughout the course of the simulation, often to the exclusion, of one or more of the ether oxygen atoms. By contrast the carbonyl oxygen atoms remain bonded to the metal atoms throughout with bond lengths that remain significantly less than those to the ether oxygen atoms. Free energy perturbation studies have been carried out in dmso and indicate that for 1, the selectivity follows the order Rb+ approximate to K+ > Cs+ >> Na+, which is partially in agreement with the experimental results. The energy differences are small and indeed the ratio between stability constants found for Cs+ and K+ complexes is only 0.60, showing that 1 has only a slight preference for K+. For the larger receptor 2, which is better suited to metal complexation, the binding affinity follows the pattern Cs+ >> Rb+ >> K+ >> Na+, with energy differences of 5.75, 2.61, 2.78 kcal mol(-1) which is perfectly consistent with experimental results.

Modelling for planning building design - a critique of available approaches

This paper reviews four approaches used to create rational tools to aid the planning and the management of the building design process and then proposes a fifth approach. The new approach that has been developed is based on the mechanical aspects of technology rather than subjective design issues. The knowledge base contains, for each construction technology, a generic model of the detailed design process. Each activity in the process is specified by its input and output information needs. By connecting the input demands of one technology with the output supply from another technology a map or network of design activity is formed. Thus, it is possible to structure a specific model from the generic knowledge base within a KBE system.

Modelling metal cutting using modern ductile fracture mechanics: Quantitative explanations for some longstanding problems

The assumption that negligible work is involved in the formation of new surfaces in the machining of ductile metals, is re-examined in the light of both current Finite Element Method (FEM) simulations of cutting and modern ductile fracture mechanics. The work associated with separation criteria in FEM models is shown to be in the kJ/m2 range rather than the few J/m2 of the surface energy (surface tension) employed by Shaw in his pioneering study of 1954 following which consideration of surface work has been omitted from analyses of metal cutting. The much greater values of surface specific work are not surprising in terms of ductile fracture mechanics where kJ/m2 values of fracture toughness are typical of the ductile metals involved in machining studies. This paper shows that when even the simple Ernst–Merchant analysis is generalised to include significant surface work, many of the experimental observations for which traditional ‘plasticity and friction only’ analyses seem to have no quantitative explanation, are now given meaning. In particular, the primary shear plane angle φ becomes material-dependent. The experimental increase of φ up to a saturated level, as the uncut chip thickness is increased, is predicted. The positive intercepts found in plots of cutting force vs. depth of cut, and in plots of force resolved along the primary shear plane vs. area of shear plane, are shown to be measures of the specific surface work. It is demonstrated that neglect of these intercepts in cutting analyses is the reason why anomalously high values of shear yield stress are derived at those very small uncut chip thicknesses at which the so-called size effect becomes evident. The material toughness/strength ratio, combined with the depth of cut to form a non-dimensional parameter, is shown to control ductile cutting mechanics. The toughness/strength ratio of a given material will change with rate, temperature, and thermomechanical treatment and the influence of such changes, together with changes in depth of cut, on the character of machining is discussed. Strength or hardness alone is insufficient to describe machining. The failure of the Ernst–Merchant theory seems less to do with problems of uniqueness and the validity of minimum work, and more to do with the problem not being properly posed. The new analysis compares favourably and consistently with the wide body of experimental results available in the literature. Why considerable progress in the understanding of metal cutting has been achieved without reference to significant surface work is also discussed.