Volume-scaled common nearest neighbor clustering algorithm with free-energy hierarchy

CommonNNClustering─A Python Package for Generic Common-Nearest-Neighbor Clustering

Density-based clustering procedures are widely used in a variety of data science applications. Their advantage lies in the capability to find arbitrarily shaped and sized clusters and robustness against outliers. In particular, they proved effective in the analysis of molecular dynamics simulations, where they serve to identify relevant, low-energetic molecular conformations. As such, they can provide a convenient basis for the construction of kinetic (core-set) Markov-state models. Here we present the open-source Python project CommonNNClustering, which provides an easy-to-use and efficient reimplementation of the common-nearest-neighbor (CommonNN) method. The package provides functionalities for hierarchical clustering and an evaluation of the results. We put our emphasis on a generic API design to keep the implementation flexible and open for customization.

Effect of Flexibility, Lipophilicity, and the Location of Polar Residues on the Passive Membrane Permeability of a Series of Cyclic Decapeptides

Cyclic peptides have received increasing attention over the recent years as potential therapeutics for "undruggable" targets. One major obstacle is, however, their often relatively poor bioavailability. Here, we investigate the structure-permeability relationship of 24 cyclic decapeptides that share the same backbone N-methylation pattern but differ in their side chains. The peptides cover a large range of values for passive membrane permeability as well as lipophilicity and solubility. To rationalize the observed differences in permeability, we extracted for each peptide the population of the membrane-permeable conformation in water from extensive explicit-solvent molecular dynamics simulations and used this as a metric for conformational rigidity or "prefolding." The insights from the simulations together with lipophilicity measurements highlight the intricate interplay between polarity/lipophilicity and flexibility/rigidity and the possible compensating effects on permeability. The findings allow us to better understand the structure-permeability relationship of cyclic peptides and extract general guiding principles.

N-Glycosylation Enhances Conformational Flexibility of Protein Disulfide Isomerase Revealed by Microsecond Molecular Dynamics and Markov State Modeling

Secreted proteins of eukaryotes are decorated with branched carbohydrate oligomers called glycans. This fact is only starting to be considered for in silico investigations of protein dynamics. Using all-atom molecular dynamics (MD) simulations and Markov state modeling (MSM), we unveil the influence of glycans on the conformational flexibility of the multidomain protein disulfide isomerase (PDI), which is a ubiquitous chaperone in the endoplasmic reticulum (ER). Yeast PDI (yPDI) from Saccharomyces cerevisiae is glycosylated at asparagine side chains and the knowledge of its five modified sites enables a realistic computational modeling. We compare simulations of glycosylated and unglycosylated yPDI and find that the presence of glycan-glycan and glycan-protein interactions influences the flexibility of PDI in different ways. For example, glycosylation reduces interdomain interactions, shifting the conformational ensemble toward more open, extended structures. In addition, we compare our results on yPDI with structural information of homologous proteins such as human PDI (hPDI), which is natively unglycosylated. Interestingly, hPDI lacks a surface recess that is present in yPDI. We find that glycosylation of yPDI facilitates its catalytic site to reach close to this surface recess. Hence, this might point to a possible functional relevance of glycosylation in yeast to act on substrates, while glycosylation seems redundant for the human homologous protein. We conclude that glycosylation is fundamental for protein dynamics, making it a necessity for a truthful representation of the flexibility and function in in silico studies of glycoproteins.

Glycan-protein interactions determine kinetics of N-glycan remodeling

A hallmark of N-linked glycosylation in the secretory compartments of eukaryotic cells is the sequential remodeling of an initially uniform oligosaccharide to a site-specific, heterogeneous ensemble of glycostructures on mature proteins. To understand site-specific processing, we used protein disulfide isomerase (PDI), a model protein with five glycosylation sites, for molecular dynamics (MD) simulations and compared the result to a biochemical in vitro analysis with four different glycan processing enzymes. As predicted by an analysis of the accessibility of the N-glycans for their processing enzymes derived from the MD simulations, N-glycans at different glycosylation sites showed different kinetic properties for the processing enzymes. In addition, altering the tertiary structure of the glycoprotein PDI affected its N-glycan remodeling in a site-specific way. We propose that the observed differential N-glycan reactivities depend on the surrounding protein tertiary structure and lead to different glycan structures in the same protein through kinetically controlled processing pathways.

Atomistic glycoprotein simulations reveal a site-specific availability of glycan substrates in time-resolved mass spectrometry of maturating enzyme kinetics.

Modulation of the Passive Permeability of Semipeptidic Macrocycles: N- and C-Methylations Fine-Tune Conformation and Properties

Incorporating small modifications to peptidic macrocycles can have a major influence on their properties. For instance, N-methylation has been shown to impact permeability. A better understanding of the relationship between permeability and structure is of key importance as peptidic drugs are often associated with unfavorable pharmacokinetic profiles. Starting from a semipeptidic macrocycle backbone composed of a tripeptide tethered head-to-tail with an alkyl linker, we investigated two small changes: peptide-to-peptoid substitution and various methyl placements on the nonpeptidic linker. Implementing these changes in parallel, we created a collection of 36 compounds. Their permeability was then assessed in parallel artificial membrane permeability assay (PAMPA) and Caco-2 assays. Our results show a systematic improvement in permeability associated with one peptoid position in the cycle, while the influence of methyl substitution varies on a case-by-case basis. Using a combination of molecular dynamics simulations and NMR measurements, we offer hypotheses to explain such behavior.