Source author record

Nicholas M. Glykos

Nicholas M. Glykos appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Biomolecules Quantitative Methods

Catalog footprint

What is connected

6works

2topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2026arXiv

Quantifying the uncertainty of molecular dynamics simulations : Good-Turing statistics revisited

We have previously shown that Good-Turing statistics can be applied to molecular dynamics trajectories to estimate the probability of observing completely new (thus far unobserved) biomolecular structures, and showed that the method is stable, dependable and its predictions verifiable. The major problem with that initial algorithm was the requirement for calculating and storing in memory the two-dimensional RMSD matrix of the currently available trajectory. This requirement precluded the application of the method to very long simulations. Here we describe a new variant of the Good-Turing algorithm whose memory requirements scale linearly with the number of structures in the trajectory, making it suitable even for extremely long simulations. We show that the new method gives essentially identical results with the older implementation, and present results obtained from trajectories containing up to 22 million structures. A computer program implementing the new algorithm is available from standard repositories.

preprint2022arXiv

Folding Molecular Dynamics Simulations of the Transmembrane Peptides of Influenza A, B M2, and MERS-, SARS-CoV E Viral Proteins

Viroporins are small viral proteins that oligomerize in the membrane of host cells and induce the formation of hydrophilic pores in these membranes, thus altering the physiological properties of the host cells. Due to their significance for viral pathogenicity, they have become targets for pharmaceutical intervention, especially through compounds that block their pore-forming activity. Here we add to the growing literature concerning the structure and function of viroporins by studying and comparing -- through molecular dynamics simulations -- the folding of the transmembrane domain peptides of viroporins derived from four viruses : influenza A, influenza B, and the coronaviruses MERS-Cov-2 and SARS-CoV-2. Through a total of more than 50 μs of simulation time in explicit solvent (TFE) and with full electrostatics, we characterize the folding behavior, helical stability and helical propensity of these transmembrane peptides in their monomeric state and we identify common motifs that may reflect their quaternary organization and/or biological function. We show that the two influenza-derived peptides are significantly different in peptide sequence and secondary structure from the two coronavirus-derived peptides, and that they are organized in two structurally distinct parts : a significantly more stable N-terminal half, and a fast converting C-terminal half that continuously folds and unfolds between $α$-helical structures and non-canonical structures which are mostly turns. In contrast, the two coronavirus-derived transmembrane peptides are much more stable and fast helix formers when compared with the influenza ones. We discuss possible interpretations of these findings and their putative connection to the structural characteristics of the respective viroporins.

preprint2021arXiv

Folding molecular dynamics simulation of T-peptide, a HIV viral entry inhibitor : Structure, dynamics, and comparison with the experimental data

Peptide T is a synthetic octapeptide fragment, which corresponds to the region 185-192 of the gp120 HIV coat protein and functions as a viral entry inhibitor. In this work, a folding molecular dynamics simulation of peptide T in a membrane-mimicking (DMSO) solution was performed with the aim of characterizing the peptide's structural and dynamical properties. We show that peptide T is highly flexible and dynamic. The main structural characteristics observed were rapidly interconverting short helical stretches and turns, with a notable preference for the formation of $β$-turns. The simulation also indicated that the C-terminal part appears to be more stable than the rest of the peptide, with the most preferred conformation for residues 5-8 being a $β$-turn. In order to validate the accuracy of the simulations, we compared our results with the experimental NMR data obtained for the T-peptide in the same solvent. In agreement with the simulation, the NMR data indicated the presence of a preferred structure in solution that was consistent with a $β$-turn comprising the four C-terminal residues. An additional comparison between the experimental and simulation-derived chemical shifts also showed a reasonable agreement between experiment and simulation, further validating the simulation-derived structural characterization of the T-peptide.

preprint2020arXiv

On the propensity of Asn-Gly-containing heptapeptides to form $β$-turn structures : comparison between ab initio quantum mechanical calculations and Molecular Dynamics simulations

Both molecular mechanical and quantum mechanical calculations play an important role in describing the behavior and structure of molecules. In this work, we compare for the same peptide systems the results obtained from folding molecular dynamics simulations with previously reported results from quantum mechanical calculations. More specifically, three molecular dynamics simulations of 5 $μ$s each in explicit water solvent were carried out for three Asn-Gly-containing heptapeptides, in order to study their folding and dynamics. Previous data, based on quantum mechanical calculations and the DFT methods have shown that these peptides adopt $β$-turn structures in aqueous solution, with type I' $β$-turn being the most preferred motif. The results from our analyses indicate that for the given system the two methods diverge in their predictions. The possibility of a force field-dependent deficiency is examined as a possible source of the observed discrepancy.

preprint2015arXiv

Clustering of molecular dynamics trajectories via peak-picking in multidimensional PCA-derived distributions

We describe a robust, fast, and memory-efficient procedure that can cluster millions of structures derived from molecular dynamics simulations. The essence of the method is based on a peak-picking algorithm applied to three- and five-dimensional distributions of the principal components derived from the trajectory and automatically supports both Cartesian and dihedral PCA-based clustering. The density threshold required for identifying isolated peaks (which correspond to discrete clusters) is determined through the application of a variance-explained criterion which allows for a completely automated clustering procedure with no user intervention. In this communication we describe the algorithm and present some of the results obtained from the application of the method as implemented in the molecular dynamics analysis programs carma, grcarma. and cluster5D. We conclude with a discussion of the limitations and possible pitfalls of this method.

preprint2013arXiv

On the application of Good-Turing statistics to quantify convergence of biomolecular simulations

Quantifying convergence and sufficient sampling of macromolecular molecular dynamics simulations is more often than not a source of controversy (and of various ad hoc solutions) in the field. Clearly, the only reasonable, consistent and satisfying way to infer convergence (or otherwise) of a molecular dynamics trajectory must be based on probability theory. Ideally, the question we would wish to answer is the following : "What is the probability that a molecular configuration important for the analysis in hand has not yet been observed ?". Here we propose a method for answering a variant of this question by using the Good-Turing formalism for frequency estimation of unobserved species in a sample. Although several approaches may be followed in order to deal with the problem of discretizing the configurational space, for this work we use the classical RMSD matrix as a means to answering the following question: "What is the probability that a molecular configuration with an RMSD (from all other already observed configurations) higher than a given threshold has not actually been observed ?". We apply the proposed method to several different trajectories and show that the procedure appears to be both computationally stable and internally consistent. A free, open-source program implementing these ideas is immediately available for download via public repositories.