An all-atom force field consistent with the general AMBER force field (GAFF) format for poly(ethylene glycol) dimethyl ether (diglyme or G2) was developed by fitting to experimental liquid densities and dielectric constants. Not surprisingly, the new force field gives excellent agreement with experimental liquid phase densities and dielectric constants over a wide temperature range. Other dynamic and thermodynamic properties of liquid G2 such as its self-diffusion coefficient, shear viscosity, and vaporization enthalpy were also calculated and compared to experimental data. For all of the properties studied, the performance of the proposed new force field is better than that of the standard GAFF force field. The force field parameters were transferred to model two other poly(ethylene glycol) ethers: monoglyme (G1) and tetraglyme (G4). The predictive ability of the modified force field for G1 and G4 was significantly better than that of the original GAFF force field. The proposed force field provides an alternative option for the simulation of mixtures containing glymes using GAFF-compatible force fields, particularly for electrochemical applications. The accuracy of a previously published force field based on the OPLS-AA format and the accuracies of two modified versions of that force field were also examined for G1, G2, and G4. It was found that the original OPLS-AA force field is superior to the modified versions of it, and that it has a similar accuracy to the proposed new GAFF-compatible force field.
A comprehensive molecular dynamics simulation study of n-alkanes using the optimized potential for liquid simulation with all-atoms (OPLS-AA) force field at ambient condition has been performed. Our results indicate that while simulations with the OPLS-AA force field accurately predict the liquid state mass density for n-alkanes with carbon number equal or less than 10, for n-alkanes with carbon number equal or exceeding 12, the OPLS-AA force field with the standard scaling factor for the 1–4 intramolecular Van der Waals and electrostatic interaction gives rise to a quasi-crystalline structure. We found that accurate predictions of the liquid state properties are obtained by successively reducing the aforementioned scaling factor for each increase of the carbon number beyond n-dodecane. To better understand the effects of reducing the scaling factor, its influence on the torsion potential profile, and the corresponding gauche-trans conformer distribution, heat of vaporization, melting point, and self-diffusion coefficient for n-dodecane were investigated. This relatively simple procedure enables more accurate predictions of the thermo-physical properties of longer n-alkanes.
In the absence of consolidated pipelines to archive biological data electronically, information dispersed in the literature must be captured by manual annotation. Unfortunately, manual annotation is time consuming and the coverage of published interaction data is therefore far from complete. The use of text-mining tools to identify relevant publications and to assist in the initial information extraction could help to improve the efficiency of the curation process and, as a consequence, the database coverage of data available in the literature. The 2006 BioCreative competition was aimed at evaluating text-mining procedures in comparison with manual annotation of protein-protein interactions.
Results
To aid the BioCreative protein-protein interaction task, IntAct and MINT (Molecular INTeraction) provided both the training and the test datasets. Data from both databases are comparable because they were curated according to the same standards. During the manual curation process, the major cause of data loss in mining the articles for information was ambiguity in the mapping of the gene names to stable UniProtKB database identifiers. It was also observed that most of the information about interactions was contained only within the full-text of the publication; hence, text mining of protein-protein interaction data will require the analysis of the full-text of the articles and cannot be restricted to the abstract.
Conclusion
The development of text-mining tools to extract protein-protein interaction information may increase the literature coverage achieved by manual curation. To support the text-mining community, databases will highlight those sentences within the articles that describe the interactions. These will supply data-miners with a high quality dataset for algorithm development. Furthermore, the dictionary of terms created by the BioCreative competitors could enrich the synonym list of the PSI-MI (Proteomics Standards Initiative-Molecular Interactions) controlled vocabulary, which is used by both databases to annotate their data content.
The structure of crystalline -cyclodextrin (-CD) hexahydrate, form I (C36H60O30·6H2O, space group P212121) is experimentally so well determined by X-ray and by neutron diffraction analyses that the positions of all the hydrogen atoms are available. This provides an opportunity for testing an empirical force field that is currently used in simulations of protein and nucleic acid structures by performing molecular dynamics studies employing the GROMOS program package on a system of 4 unit cells containing 16 -CD molecules and 96 water molecules.A detailed comparison of the simulated and experimentally determined crystal structures shows that the experimental positions of the -CD atoms are reproduced within 0.025 nm, well within the overall experimental accuracy of 0.036 nm; that the water molecules are on average within 0.072 nm from their experimental sites, with two thirds reproduced within experimental accuracy by the calculations; that high correlation is produced, between the occurrence of simulated and experimentally observed hydrogen bonds.The good agreement between simulated and experimental results suggests that the tested force field is reliable. 相似文献
Unknown force-field parameters for metal organic beryllium complexes used in emitting and electron transporting layers of OLED structures are determined. These parameters can be used for the predictive atomistic simulations of the structure and properties of amorphous organic layers containing beryllium complexes. The parameters are found for the AMBER force field using a relaxed scan procedure and quantum-mechanical DFT calculations of potential energy curves for specific internal (angular) coordinates in a series of three Be complexes (Bebq2; Be(4-mpp)2; Bepp2). The obtained parameters are verified in calculations of some molecular and crystal structures available from either quantum-mechanical DFT calculations or experimental data.
Understanding living systems is crucial for curing diseases. To achieve this task we have to understand biological networks based on protein-protein interactions. Bioinformatics has come up with a great amount of databases and tools that support analysts in exploring protein-protein interactions on an integrated level for knowledge discovery. They provide predictions and correlations, indicate possibilities for future experimental research and fill the gaps to complete the picture of biochemical processes. There are numerous and huge databases of protein-protein interactions used to gain insights into answering some of the many questions of systems biology. Many computational resources integrate interaction data with additional information on molecular background. However, the vast number of diverse Bioinformatics resources poses an obstacle to the goal of understanding. We present a survey of databases that enable the visual analysis of protein networks.
Results
We selected M =10 out of N =53 resources supporting visualization, and we tested against the following set of criteria: interoperability, data integration, quantity of possible interactions, data visualization quality and data coverage. The study reveals differences in usability, visualization features and quality as well as the quantity of interactions. StringDB is the recommended first choice. CPDB presents a comprehensive dataset and IntAct lets the user change the network layout. A comprehensive comparison table is available via web. The supplementary table can be accessed on http://tinyurl.com/PPI-DB-Comparison-2015.
Conclusions
Only some web resources featuring graph visualization can be successfully applied to interactive visual analysis of protein-protein interaction. Study results underline the necessity for further enhancements of visualization integration in biochemical analysis tools. Identified challenges are data comprehensiveness, confidence, interactive feature and visualization maturing. 相似文献
From the crystals of trans aquabis(N,N-dimethylglycinato-κN,κO)copper(II) dihydrate (compound 1, space group P212121) novel crystal structure of trans aquabis(N,N-dimethylglycinato-κN,κO)copper(II) (compound 2, space group Pbca) was obtained and analysed by X-ray diffraction. In the crystal structure 1, the O-H?O hydrogen bonds form three-dimensional network. In the crystal structure 2, two-dimensional layers stacking to each other are formed, with non-polar N,N-dimethyl groups placed on the opposite sides of the layers, and with the polar part in the middle forming CO?O-H and C-H?O hydrogen bonds. Different hydrogen bonding patterns in 1 and 2 do not pronouncedly affect molecular geometry of the title compound. Molecular mechanics force field suited for studying the properties of bis(amino acidato)copper(II) complexes in the solid state can follow the differences between the experimental molecular structures in the two diverse crystalline surroundings. To make possible direct comparison between crystal lattices, the force field was applied to predict unit cell packing of supposed anhydrous bis(N,N-dimethylglycinato)copper(II) in space group Pbca. Relative intermolecular energies of hypothetic anhydrous crystal and simulated 1 and 2 crystals are discussed. On the basis of experimental and theoretical results we conclude that the main effect of two water molecules of crystallisation in 1 is to stabilise the crystal packing via hydrogen bonding, whilst similar pyramidal copper(II) coordination geometry in 1 and 2 is due to axially coordinated water molecule and its intermolecular interactions. 相似文献
Prevention of disability (POD) is one of the key objectives of leprosy programmes. Recently, coverage and access have been identified as the priority issues in POD. Assessing the cost-effectiveness of POD interventions is highly relevant to understanding the barriers and opportunities to achieving universal coverage and access with limited resources. The purpose of this study was to systematically review the quality of existing cost-effectiveness evidence and discuss implications for future research and strategies to prevent disability in leprosy and other disabling conditions.
Methodology/Principal Findings
We searched electronic databases (NHS EED, MEDLINE, EMBASE, and LILACS) and databases of ongoing trials (www.controlled-trials.com/mrct/, www.who.int/trialsearch). We checked reference lists and contacted experts for further relevant studies. We included studies that reported both cost and effectiveness outcomes of two or more alternative interventions to prevent disability in leprosy. We assessed the quality of the identified studies using a standard checklist for critical appraisal of economic evaluations of health care programmes. We found 66 citations to potentially relevant studies and three met our criteria. Two were randomised controlled trials (footwear, management of neuritis) and one was a generic model-based study (cost per DALY). Generally, the studies were small in size, reported inadequately all relevant costs, uncertainties in estimates, and issues of concern and were based on limited data sources. No cost-effectiveness data on self-care, which is a key strategy in POD, was found.
Conclusion/Significance
Evidence for cost-effectiveness of POD interventions for leprosy is scarce. High quality research is needed to identify POD interventions that offer value for money where resources are very scarce, and to develop strategies aimed at available, affordable and sustainable quality POD services for leprosy. The findings are relevant for other chronically disabling conditions, such as lymphatic filariasis, Buruli ulcer and diabetes in developing countries. 相似文献
Diameter and wall thickness of self-assembled peptide nanotube of cyclo[(-d-Ala-l-Ala)4-] were characterised by molecular simulation. In order to verify the existence of peptide nanotube of cyclo[(-d-Ala-l-Ala)4-], cyclo[(-d-Ala-l-Ala)4-] was firstly synthesised through Fmoc solid-phase synthesis method and then self-assembled in trifluoroacetic acid. Based on the results of experiment, the single nanotube structure was further characterised by molecular dynamics (MD) employing the COMPASS force field. The results indicate that cyclo[(-d-Ala-l-Ala)4-] is self-assembled into nanotube bundles of about 0.5 μm in diameter and 10 μm in length; the inner and outer diameter of the single nanotube is 8.5 and 15.9 Å, respectively, and the nanotube wall thickness is 3.3 Å. 相似文献
The COMPASS tokamak (R = 0.56 m, a = 0.2 m, BT = 1.3 T, Ip ~ 300 kA, pulse duration 0.4 s) operates in ITER-like plasma shape in H-mode with Type-I ELMs. In 2019, we plan to install into the divertor a test target based on capillary porous system filled with liquid lithium/tin. This single target will be inclined toroidally in order to be exposed to ITER-relevant surface heat flux (20 MW/m2). Based on precisely measured actual heat fluxes, our simulations predict (for 45° inclination, without accounting for the lithium vapor shielding) the surface temperature rises up to 700°C within 120 ms of the standard ELMy H-mode heat flux with ELM filaments reaching hundreds MW/m2. Significant lithium vaporization is expected. The target surface will be observed by spectroscopy, fast visible and infrared cameras. The scientific program will be focused on operational issues (redeposition of the evaporated metal, ejection of droplets, if any) as well as on the effect on the plasma physics (improvement of plasma confinement, L–H power threshold, Zeff, etc.). After 2024, a closed liquid divertor may be installed into the planned COMPASS Upgrade tokamak (R = 0.84 m, a = 0.3 m, BT = 5 T, Ip = 2 MA, Pin = 8 MW, pulse duration ~2 s) with ITER-relevant heat fluxes loading the entire toroidal divertor. 相似文献
The conformational flexibility and the dynamics of -D-Manp-(12)--D-Glcp-OMe have been investigated by Metropolis Monte Carlo (MMC) and Langevin dynamics (LD) simulations. The two simulation techniques employ different force fields, namely the HSEA force field and a CHARMm-based force field. The former shows less conformational flexibility than the latter, in which a multiple energy minima conformational space is sampled. Long-range heteronuclear nuclear magnetic resonance (NMR) coupling constants have been measured by selective excitations of the carbons at the glycosidic linkage. Calculated 3JC, H values from MMC and LD simulations show excellent agreement to those from NMR experiments. The X-ray crystal structure has a conformation within a region of the conformational space populated in both force fields. 相似文献
Life cycle impact assessment (LCIA) results are used to assess potential environmental impacts of different products and services. As part of the UNEP-SETAC life cycle initiative flagship project that aims to harmonize indicators of potential environmental impacts, we provide a consensus viewpoint and recommendations for future developments in LCIA related to the ecosystem quality area of protection (AoP). Through our recommendations, we aim to encourage LCIA developments that improve the usefulness and global acceptability of LCIA results.
Methods
We analyze current ecosystem quality metrics and provide recommendations to the LCIA research community for achieving further developments towards comparable and more ecologically relevant metrics addressing ecosystem quality.
Results and discussion
We recommend that LCIA development for ecosystem quality should tend towards species-richness-related metrics, with efforts made towards improved inclusion of ecosystem complexity. Impact indicators—which result from a range of modeling approaches that differ, for example, according to spatial and temporal scale, taxonomic coverage, and whether the indicator produces a relative or absolute measure of loss—should be framed to facilitate their final expression in a single, aggregated metric. This would also improve comparability with other LCIA damage-level indicators. Furthermore, to allow for a broader inclusion of ecosystem quality perspectives, the development of an additional indicator related to ecosystem function is recommended. Having two complementary metrics would give a broader coverage of ecosystem attributes while remaining simple enough to enable an intuitive interpretation of the results.
Conclusions
We call for the LCIA research community to make progress towards enabling harmonization of damage-level indicators within the ecosystem quality AoP and, further, to improve the ecological relevance of impact indicators.
Mass spectrometric analysis of microbial metabolism provides a long list of possible compounds. Restricting the identification of the possible compounds to those produced by the specific organism would benefit the identification process. Currently, identification of mass spectrometry (MS) data is commonly done using empirically derived compound databases. Unfortunately, most databases contain relatively few compounds, leaving long lists of unidentified molecules. Incorporating genome-encoded metabolism enables MS output identification that may not be included in databases. Using an organism’s genome as a database restricts metabolite identification to only those compounds that the organism can produce.
Results
To address the challenge of metabolomic analysis from MS data, a web-based application to directly search genome-constructed metabolic databases was developed. The user query returns a genome-restricted list of possible compound identifications along with the putative metabolic pathways based on the name, formula, SMILES structure, and the compound mass as defined by the user. Multiple queries can be done simultaneously by submitting a text file created by the user or obtained from the MS analysis software. The user can also provide parameters specific to the experiment’s MS analysis conditions, such as mass deviation, adducts, and detection mode during the query so as to provide additional levels of evidence to produce the tentative identification. The query results are provided as an HTML page and downloadable text file of possible compounds that are restricted to a specific genome. Hyperlinks provided in the HTML file connect the user to the curated metabolic databases housed in ProCyc, a Pathway Tools platform, as well as the KEGG Pathway database for visualization and metabolic pathway analysis.
Conclusions
Metabolome Searcher, a web-based tool, facilitates putative compound identification of MS output based on genome-restricted metabolic capability. This enables researchers to rapidly extend the possible identifications of large data sets for metabolites that are not in compound databases. Putative compound names with their associated metabolic pathways from metabolomics data sets are returned to the user for additional biological interpretation and visualization. This novel approach enables compound identification by restricting the possible masses to those encoded in the genome. 相似文献
Next generation sequencing (NGS) platforms are currently being utilized for targeted sequencing of candidate genes or genomic intervals to perform sequence-based association studies. To evaluate these platforms for this application, we analyzed human sequence generated by the Roche 454, Illumina GA, and the ABI SOLiD technologies for the same 260 kb in four individuals.
Results
Local sequence characteristics contribute to systematic variability in sequence coverage (>100-fold difference in per-base coverage), resulting in patterns for each NGS technology that are highly correlated between samples. A comparison of the base calls to 88 kb of overlapping ABI 3730xL Sanger sequence generated for the same samples showed that the NGS platforms all have high sensitivity, identifying >95% of variant sites. At high coverage, depth base calling errors are systematic, resulting from local sequence contexts; as the coverage is lowered additional 'random sampling' errors in base calling occur.
Conclusions
Our study provides important insights into systematic biases and data variability that need to be considered when utilizing NGS platforms for population targeted sequencing studies. 相似文献
Male sexually selected vocalisations generally contain both individuality and quality cues that are crucial in intra- as well as inter-sexual communication. As individuality is a fixed feature whereas male phenotypic quality changes with age, individuality and quality cues may be subjected to different selection pressures over time. Individuality (for example, morphology of the vocal apparatus) and quality (for example, body size and dominance status) can both affect the vocal production mechanism, inducing the same components of vocalisations to convey both kinds of information. In this case, do quality-related changes to the acoustic structure of calls induce a modification of vocal cues to identity from year to year? We investigated this question in fallow deer (Dama dama), in which some acoustic parameters of vocalisations (groans) code for both individuality and quality.
Results
We carried out a longitudinal analysis of groan individuality, examining the effects of age and dominance rank on the acoustic structure of groans of the same males recorded during consecutive years. We found both age- and rank-related changes to groans; the minimum values of the highest formant frequencies and the fundamental frequency increased with the age of males and they decreased when males became more dominant. Both age- and rank-related acoustic parameters contributed to individuality. Male quality changed with age, inducing a change in quality-related parameters and thus, a modification of vocal cues to male individuality between years.
Conclusions
The encoding of individuality and quality information in the same components of vocalisations induces a tradeoff between these two kinds of signals over time. Fallow deer vocalisations are honest signals of quality that are not fixed over time but are modified dynamically according to male quality. As they are more reliable cues to quality than to individuality, they may not be used by conspecifics to recognize a given male from one year to another, but potentially used by both sexes to assess male quality during each breeding season.
Due to its protective properties of biological samples at low temperatures and under desiccation, dimethyl sulfoxide (DMSO) in aqueous solutions has been studied widely by many experimental approaches and molecular dynamics (MD) simulations. In the case of the latter, AMBER is among the most commonly used force fields for simulations of biomolecular systems; however, the parameters for DMSO published by Fox and Kollman in 1998 have only been tested for pure liquid DMSO. We have conducted an MD simulation study of DMSO in a water mixture and computed several structural and dynamical properties such as of the mean density, self-diffusion coefficient, hydrogen bonding and DMSO and water ordering. The AMBER force field of DMSO is seen to reproduce well most of the experimental properties of DMSO in water, with the mixture displaying strong and specific water ordering, as observed in experiments and multiple other MD simulations with other non-polarizable force fields.
