Browse
Submit Data
Databases
API
Help

Dataset Information

16 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Length Dependent Folding Kinetics of Alanine-Based Helical Peptides from Optimal Dimensionality Reduction.

ABSTRACT: We present a computer simulation study of helix folding in alanine homopeptides (ALA)n of length n = 5, 8, 15, and 21 residues. Based on multi-microsecond molecular dynamics simulations at room temperature, we found helix populations and relaxation times increasing from about 6% and ~2 ns for ALA5 to about 60% and ~500 ns for ALA21, and folding free energies decreasing linearly with the increasing number of residues. The helix folding was analyzed with the Optimal Dimensionality Reduction method, yielding coarse-grained kinetic models that provided a detailed representation of the folding process. The shorter peptides, ALA5 and ALA8, tended to convert directly from coil to helix, while ALA15 and ALA21 traveled through several intermediates. Coarse-grained aggregate states representing the helix, coil, and intermediates were heterogeneous, encompassing multiple peptide conformations. The folding involved multiple pathways and interesting intermediate states were present on the folding paths, with partially formed helices, turns, and compact coils. Statistically, helix initiation was favored at both termini, and the helix was most stable in the central region. Importantly, we found the presence of underlying universal local dynamics in helical peptides with correlated transitions for neighboring hydrogen bonds. Overall, the structural and dynamical parameters extracted from the trajectories are in good agreement with experimental observables, providing microscopic insights into the complex helix folding kinetics.

SUBMITTER: Kuczera K

PROVIDER: S-EPMC8170890 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Similar Datasets

Evaluation of Dimensionality-reduction Methods from Peptide Folding-unfolding Simulations.

Project description:Dimensionality reduction methods have been widely used to study the free energy landscapes and low-free energy pathways of molecular systems. It was shown that the non-linear dimensionality-reduction methods gave better embedding results than the linear methods, such as principal component analysis, in some simple systems. In this study, we have evaluated several non linear methods, locally linear embedding, Isomap, and diffusion maps, as well as principal component analysis from the equilibrium folding/unfolding trajectory of the second ?-hairpin of the B1 domain of streptococcal protein G. The CHARMM parm19 polar hydrogen potential function was used. A series of criteria which reflects different aspects of the embedding qualities were employed in the evaluation. Our results show that principal component analysis is not worse than the non-linear ones on this complex system. There is no clear winner in all aspects of the evaluation. Each dimensionality-reduction method has its limitations in a certain aspect. We emphasize that a fair, informative assessment of an embedding result requires a combination of multiple evaluation criteria rather than any single one. Caution should be used when dimensionality-reduction methods are employed, especially when only a few of top embedding dimensions are used to describe the free energy landscape.

| S-EPMC3678838 | biostudies-literature

Optimal Dimensionality Reduction of Multistate Kinetic and Markov-State Models.

Project description:We develop a systematic procedure for obtaining rate and transition matrices that optimally describe the dynamics of aggregated superstates formed by combining (clustering or lumping) microstates. These reduced dynamical models are constructed by matching the time-dependent occupancy-number correlation functions of the superstates in the full and aggregated systems. Identical results are obtained by using a projection operator formalism. The reduced dynamic models are exact for all times in their full non-Markovian formulation. In the approximate Markovian limit, we derive simple analytic expressions for the reduced rate or Markov transition matrices that lead to exact auto- and cross-relaxation times. These reduced Markovian models strike an optimal balance between matching the dynamics at short and long times. We also discuss how this approach can be used in a hierarchical procedure of constructing optimal superstates through aggregation of microstates. The results of the general reduced-matrix theory are illustrated with applications to simple model systems and a more complex master-equation model of peptide folding derived previously from atomistic molecular dynamics simulations. We find that the reduced models faithfully capture the dynamics of the full systems, producing substantial improvements over the common local-equilibrium approximation.

| S-EPMC4516310 | biostudies-other

Low-dimensional, free-energy landscapes of protein-folding reactions by nonlinear dimensionality reduction.

Project description:The definition of reaction coordinates for the characterization of a protein-folding reaction has long been a controversial issue, even for the "simple" case in which one single free-energy barrier separates the folded and unfolded ensemble. We propose a general approach to this problem to obtain a few collective coordinates by using nonlinear dimensionality reduction. We validate the usefulness of this method by characterizing the folding landscape associated with a coarse-grained protein model of src homology 3 as sampled by molecular dynamics simulations. The folding free-energy landscape projected on the few relevant coordinates emerging from the dimensionality reduction can correctly identify the transition-state ensemble of the reaction. The first embedding dimension efficiently captures the evolution of the folding process along the main folding route. These results clearly show that the proposed method can efficiently find a low-dimensional representation of a complex process such as protein folding.

| S-EPMC1502548 | biostudies-literature

Incorporating Polarizability of Backbone Hydrogen Bonds Improved Folding of Short ?-Helical Peptides.

Project description:Reliability of force fields is an essential aspect of protein-folding simulation. In this work, we introduced a newly developed on-the-fly charge-updating scheme called the polarized structure-specific backbone charge (PSBC) model. The PSBC model was designed with the purpose of building the polarizability of backbone hydrogen bonds into the force field by updating the partial charges of backbone hydrogen-bond donor and acceptor atoms during folding simulation. This implementation was intended to mimic the heterogeneity of the protein surrounding during folding. Multiple single-trajectory molecular dynamics simulations were performed to fold a polyalanine peptide, namely ER (Ac-A(EAAAR)3A-NH2), using both polarizable (PSBC) and nonpolarizable (Amber03) force fields. Through the PSBC model, ER was folded into a helical peptide with helix content that agrees well with experiments. Comparison between simulations performed using the aforementioned force fields demonstrably showed the importance of electrostatic polarization effect in the folding of the short ?-helical peptide. The PSBC model was further validated by folding two other short peptides with different helicities.

| S-EPMC6895711 | biostudies-literature

Systematic Dimensionality Reduction for Quantum Walks: Optimal Spatial Search and Transport on Non-Regular Graphs.

Project description:Continuous time quantum walks provide an important framework for designing new algorithms and modelling quantum transport and state transfer problems. Often, the graph representing the structure of a problem contains certain symmetries that confine the dynamics to a smaller subspace of the full Hilbert space. In this work, we use invariant subspace methods, that can be computed systematically using the Lanczos algorithm, to obtain the reduced set of states that encompass the dynamics of the problem at hand without the specific knowledge of underlying symmetries. First, we apply this method to obtain new instances of graphs where the spatial quantum search algorithm is optimal: complete graphs with broken links and complete bipartite graphs, in particular, the star graph. These examples show that regularity and high-connectivity are not needed to achieve optimal spatial search. We also show that this method considerably simplifies the calculation of quantum transport efficiencies. Furthermore, we observe improved efficiencies by removing a few links from highly symmetric graphs. Finally, we show that this reduction method also allows us to obtain an upper bound for the fidelity of a single qubit transfer on an XY spin network.

| S-EPMC4557070 | biostudies-literature

Analysis of the kinetics of folding of proteins and peptides using circular dichroism.

Project description:Circular dichroism (CD) is a useful spectroscopic technique for studying the secondary structure, folding and binding properties of proteins. This protocol covers how to use the intrinsic circular dichroic properties of proteins to follow their folding and unfolding as a function of time. Included are methods of obtaining data and for analyzing the folding and unfolding data to determine the rate constants and the order of the folding and unfolding reactions. The protocol focuses on the use of CD to follow folding when it is relatively slow, on the order of minutes to days. The methods for analyzing the data, however, can also be applied to data collected with a CD machine equipped with stopped-flow accessories in the range of milliseconds to seconds and folding analyzed by other spectroscopic methods including changes in absorption or fluorescence spectra as a function of time.

| S-EPMC2752309 | biostudies-literature

Influence of Glu/Arg, Asp/Arg, and Glu/Lys Salt Bridges on ?-Helical Stability and Folding Kinetics.

Project description:Using a combination of ultraviolet circular dichroism, temperature-jump transient-infrared spectroscopy, and molecular dynamics simulations, we investigate the effect of salt bridges between different types of charged amino-acid residue pairs on ?-helix folding. We determine the stability and the folding and unfolding rates of 12 alanine-based ?-helical peptides, each of which has a nearly identical composition containing three pairs of positively and negatively charged residues (either Glu(-)/Arg(+), Asp(-)/Arg(+), or Glu(-)/Lys(+)). Within each set of peptides, the distance and order of the oppositely charged residues in the peptide sequence differ, such that they have different capabilities of forming salt bridges. Our results indicate that stabilizing salt bridges (in which the interacting residues are spaced and ordered such that they favor helix formation) speed up ?-helix formation by up to 50% and slow down the unfolding of the ?-helix, whereas salt bridges with an unfavorable geometry have the opposite effect. Comparing the peptides with different types of charge pairs, we observe that salt bridges between side chains of Glu(-) and Arg(+) are most favorable for the speed of folding, probably because of the larger conformational space of the salt-bridging Glu(-)/Arg(+) rotamer pairs compared to Asp(-)/Arg(+) and Glu(-)/Lys(+). We speculate that the observed impact of salt bridges on the folding kinetics might explain why some proteins contain salt bridges that do not stabilize the final, folded conformation.

| S-EPMC4906143 | biostudies-literature

Predicting unfolding thermodynamics and stable intermediates for alanine-rich helical peptides with the aid of coarse-grained molecular simulation.

Project description:This report focuses on the molecular-level processes and thermodynamics of unfolding of a series of helical peptides using a coarse-grained (CG) molecular model. The CG model was refined to capture thermodynamics and structural changes as a function of temperature for a set of published peptide sequences. Circular dichroism spectroscopy (CD) was used to experimentally monitor the temperature-dependent conformational changes and stability of published peptides and new sequences introduced here. The model predictions were quantitatively or semi-quantitatively accurate in all cases. The simulations and CD results showed that, as expected, in most cases the unfolding of helical peptides is well described by a simply 2-state model, and conformational stability increased with increased length of the helices. A notable exception in a 19-residue helix was when two Ala residues were each replaced with Phe. This stabilized a partly unfolded intermediate state via hydrophobic contacts, and also promoted aggregates at higher peptide concentrations.

| S-EPMC5011438 | biostudies-literature

Global tests of P-values for multifactor dimensionality reduction models in selection of optimal number of target genes.

Project description:BACKGROUND:Multifactor Dimensionality Reduction (MDR) is a popular and successful data mining method developed to characterize and detect nonlinear complex gene-gene interactions (epistasis) that are associated with disease susceptibility. Because MDR uses a combinatorial search strategy to detect interaction, several filtration techniques have been developed to remove genes (SNPs) that have no interactive effects prior to analysis. However, the cutoff values implemented for these filtration methods are arbitrary, therefore different choices of cutoff values will lead to different selections of genes (SNPs). METHODS:We suggest incorporating a global test of p-values to filtration procedures to identify the optimal number of genes/SNPs for further MDR analysis and demonstrate this approach using a ReliefF filter technique. We compare the performance of different global testing procedures in this context, including the Kolmogorov-Smirnov test, the inverse chi-square test, the inverse normal test, the logit test, the Wilcoxon test and Tippett's test. Additionally we demonstrate the approach on a real data application with a candidate gene study of drug response in Juvenile Idiopathic Arthritis. RESULTS:Extensive simulation of correlated p-values show that the inverse chi-square test is the most appropriate approach to be incorporated with the screening approach to determine the optimal number of SNPs for the final MDR analysis. The Kolmogorov-Smirnov test has high inflation of Type I errors when p-values are highly correlated or when p-values peak near the center of histogram. Tippett's test has very low power when the effect size of GxG interactions is small. CONCLUSIONS:The proposed global tests can serve as a screening approach prior to individual tests to prevent false discovery. Strong power in small sample sizes and well controlled Type I error in absence of GxG interactions make global tests highly recommended in epistasis studies.

| S-EPMC3508622 | biostudies-literature

Length-Dependent Formation of Transmembrane Pores by 310-Helical ?-Aminoisobutyric Acid Foldamers.

Project description:The synthetic biology toolbox lacks extendable and conformationally controllable yet easy-to-synthesize building blocks that are long enough to span membranes. To meet this need, an iterative synthesis of ?-aminoisobutyric acid (Aib) oligomers was used to create a library of homologous rigid-rod 310-helical foldamers, which have incrementally increasing lengths and functionalizable N- and C-termini. This library was used to probe the inter-relationship of foldamer length, self-association strength, and ionophoric ability, which is poorly understood. Although foldamer self-association in nonpolar chloroform increased with length, with a ? 14-fold increase in dimerization constant from Aib6 to Aib11, ionophoric activity in bilayers showed a stronger length dependence, with the observed rate constant for Aib11 ? 70-fold greater than that of Aib6. The strongest ionophoric activity was observed for foldamers with >10 Aib residues, which have end-to-end distances greater than the hydrophobic width of the bilayers used (? 2.8 nm); X-ray crystallography showed that Aib11 is 2.93 nm long. These studies suggest that being long enough to span the membrane is more important for good ionophoric activity than strong self-association in the bilayer. Planar bilayer conductance measurements showed that Aib11 and Aib13, but not Aib7, could form pores. This pore-forming behavior is strong evidence that Aibm (m ? 10) building blocks can span bilayers.

| S-EPMC4752191 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data