Simultaneous computation of dynamical and equilibrium information using a weighted ensemble of trajectories

Equilibrium formally can be represented as an ensemble of uncoupled systems undergoing unbiased dynamics in which detailed balance is maintained. Many non-equilibrium processes can be described by suitable subsets of the equilibrium ensemble. Here, we employ the "weighted ensemble" (WE) simulation protocol [Huber and Kim, Biophys. J., 1996] to generate equilibrium trajectory ensembles and extract non-equilibrium subsets for computing kinetic quantities. States do not need to be chosen in advance. The procedure formally allows estimation of kinetic rates between arbitrary states chosen after the simulation, along with their equilibrium populations. We also describe a related history-dependent matrix procedure for estimating equilibrium and non-equilibrium observables when phase space has been divided into arbitrary non-Markovian regions, whether in WE or ordinary simulation. In this proof-of-principle study, these methods are successfully applied and validated on two molecular systems: explicitly solvated methane association and the implicitly solvated Ala4 peptide. We comment on challenges remaining in WE calculations

Large expert-curated database for benchmarking document similarity detection in biomedical literature search

Document recommendation systems for locating relevant literature have mostly relied on methods developed a decade ago. This is largely due to the lack of a large offline gold-standard benchmark of relevant documents that cover a variety of research fields such that newly developed literature search techniques can be compared, improved and translated into practice. To overcome this bottleneck, we have established the RElevant LIterature SearcH consortium consisting of more than 1500 scientists from 84 countries, who have collectively annotated the relevance of over 180 000 PubMed-listed articles with regard to their respective seed (input) article/s. The majority of annotations were contributed by highly experienced, original authors of the seed articles. The collected data cover 76% of all unique PubMed Medical Subject Headings descriptors. No systematic biases were observed across different experience levels, research fields or time spent on annotations. More importantly, annotations of the same document pairs contributed by different scientists were highly concordant. We further show that the three representative baseline methods used to generate recommended articles for evaluation (Okapi Best Matching 25, Term Frequency-Inverse Document Frequency and PubMed Related Articles) had similar overall performances. Additionally, we found that these methods each tend to produce distinct collections of recommended articles, suggesting that a hybrid method may be required to completely capture all relevant articles. The established database server located at https://relishdb.ict.griffith.edu.au is freely available for the downloading of annotation data and the blind testing of new methods. We expect that this benchmark will be useful for stimulating the development of new powerful techniques for title and title/abstract-based search engines for relevant articles in biomedical research.Peer reviewe

Role of length-dependent stability of collagen-like peptides

Understanding the structure, folding, and stability of collagen is complex because of its length and variations in the amino acid (AA) sequence composition. It is well known that the basic constituent of the collagen helix is the triplet repeating sequence of the form Gly-XAA-YAA. On the basis of previous models and with the frequency of occurrence of the triplets, the ((Gly-Pro-Hyp)n)3 (where n is the number of triplets) sequence replicate has been chosen as the model for the most stable form of the collagen-like sequence. With a view to understand the role of sequence length (or the number of triplets) on the stability of collagen, molecular dynamics simulations have been carried out by varying the number of triplet units on the model collagen-like peptides. The results reveal that five triplets are required to form the stable triple helix. Further analysis shows that the intermolecular structural rigidity of the imino acid residues, hydrogen bonding, and water structure around the three chains of the triple helix play the dominant roles on its structure, folding, and stabilization

Exploring the changes in the structure of α-helical peptides adsorbed onto a single walled carbon nanotube using classical molecular dynamics simulation

Classical molecular dynamics (MD) simulation has been carried out in an explicit solvent environment to understand the interaction between the single walled carbon nanotube (SWCNT) and α-helix. A polyalanine peptide consisting of 40 alanine residues has been chosen as the model for the α-helix (PA40). Results reveal that the SWCNT induces conformational changes in PA40. Furthermore, breakage of hydrogen bonds in the chosen model peptides has been observed, which leads to conformational transitions (α → turns) in different parts of the PA40. Owing to these transitions, regions of different structural and energetic stability are generated in PA40 which enable the PA40 to curl around the surface of the SWCNT. The overall observations obtained from the MD simulations are not significantly influenced by the starting geometry and the choice of the force field. Although the qualities of structural information obtained from the MD simulation using ff03 and OPLS are different, the overall observation derived from the ff03 is similar to that of OPLS. Results from the MD simulation on the interaction of the α-helical fragment of the SNARES protein with the SWCNT elicit that the amino acid composition influences the interaction pattern. The wrapping of the α-helical fragment of the SNARES onto the SWCNT is similar to that of PA40. Overall, there is a considerable decrease in the helical content of peptides upon interaction with SWCNTs, in agreement with the experimental findings

Role of aspartic acid in collagen structure and stability: a molecular dynamics investigation

A molecular dynamics (MD) simulation study has been carried out to understand the stability of the triple helical collagen models. The calculations show that the presence of the aspartic acid residue in different positions leads to the local variation in the structure. Analyses of root-mean-square deviation (RMSD), radial distribution function (RDF), puckering effect, dihedral angle variation, hydrogen bond (H-bond), and conformational changes during molecular dynamics simulation reveal that the local perturbation in the sequences, increase in chain flexibility due to removal of five membered rings in the collagen by aspartic acid, change of intermolecular H-bonding pattern, and differences in the association of water are mainly influencing the nature of stabilization of collagen by aspartic acid

Bader's electron density analysis of hydrogen bonding in secondary structural elements of protein

The hydrogen-bonding (H-bonding) interactions in α-helical and β-sheet model peptides have been studied by using the atoms-in-molecule (AIM) approach. The relative importance of NH···O and CH···O H-bonding interactions in the different secondary elements such as α-helix, parallel, and antiparallel β-sheets have been assessed. The electron density values at the NH···O bond are higher than those of the CH···O bonds in the α-helical conformation. The electron density values at the H-bonded critical points (HBCPs) corresponding to NH···O and CH···O interactions are nearly equal in the parallel β-sheet of the order of 10-3 au, whereas in the case of antiparallel β-sheets, ρ(rc) values for NH···O and CH···O interactions are of the order of 10-2 and 10-3 au, respectively. It is interesting to point out here that the weakening of NH···O interactions in the parallel β-sheet arrangement is evident from the AIM analysis. This is concomitant with the increase in the NH···O distance in the parallel β-sheet conformation. In addition to the clear description of H-bonding by electron density at the HBCP, possible good linear relationships between the electron density at ring critical points (RCP) and stabilization energy (SE) have been observed corresponding to the various β-sheet conformations