Common folding processes of mini-proteins: Partial formations of secondary structures initiate the immediate protein folding

Protein Structure Validation Derives a Smart Conformational Search in a Physically Relevant Configurational Subspace

Since proteins perform biological functions through their dynamic properties, molecular dynamics (MD) simulation is a sophisticated strategy for investigating their functions. Analyses of trajectories provide statistical information about a specific protein as a free-energy landscape (FEL). However, the timescale of normal MD is shorter than that of biological functions, resulting in statistically insufficient conformational sampling, finally leading to unreliable FEL calculation. To search for a broad configurational subspace, an external bias is imposed on a target protein as biased sampling. However, its regulation is challenging because the optimal strength of the perturbation is unknown. Furthermore, a physically irrelevant configurational subspace was searched when imposing an inappropriate external bias. To address this issue, we newly proposed an external biased regulation scheme known as the G-factor external bias limiter (GERBIL). In GERBIL, protein configurations generated by external bias are structurally validated by an indicator (G-factor), enabling the search for a physically relevant subspace. In addition to biased sampling, nonbiased sampling might search for a physically irrelevant configurational subspace because repeating multiple MD simulations from several initial structures tends to search for an overly broad configurational subspace. For this issue, the structural qualities of configurations generated by nonbiased sampling have not been investigated. Therefore, we confirmed whether the G-factor screened the collapsed (low-quality) configurations generated by nonbiased sampling. To address this issue, the outlier flooding method (OFLOOD) was adopted in GERBIL as a nonbiased sampling method, which is referred to as OFLOOD-GERBIL. OFLOOD rapidly expands a configurational subspace by resampling the rarely occurring states of a given protein and tends to search an overly broad subspace. Thus, we considered that GERBIL might improve the excessive conformational search of OFLOOD for a physically irrelevant configurational subspace. As a demonstration, OFLOOD and OFLOOD-GERBIL were applied to a globular protein (T4 lysozyme) and their conformational search qualities were assessed. Based on our assessment, normal OFLOOD without the outlier validation frequently sampled low-quality configurations, whereas OFLOOD-GERBIL with the outlier validation intensively sampled high-quality configurations. In conclusion, OFLOOD-GERBIL derives a smart conformational search in a physically relevant configurational subspace, indicating that protein structure validation works in both nonbiased and biased sampling methods.

Residue Folding Degree-Relationship to Secondary Structure Categories and Use as Collective Variable

Recently, we have shown that the residue folding degree, a network-based measure of folded content in proteins, is able to capture backbone conformational transitions related to the formation of secondary structures in molecular dynamics (MD) simulations. In this work, we focus primarily on developing a collective variable (CV) for MD based on this residue-bound parameter to be able to trace the evolution of secondary structure in segments of the protein. We show that this CV can do just that and that the related energy profiles (potentials of mean force, PMF) and transition barriers are comparable to those found by others for particular events in the folding process of the model mini protein Trp-cage. Hence, we conclude that the relative segment folding degree (the newly proposed CV) is a computationally viable option to gain insight into the formation of secondary structures in protein dynamics. We also show that this CV can be directly used as a measure of the amount of α-helical content in a selected segment.

Selection Rules for Outliers in Outlier Flooding Method Regulate Its Conformational Sampling Efficiency

The outlier flooding method (OFLOOD) has been proposed as an enhanced conformational sampling method of proteins. In OFLOOD, rarely occurring states of proteins are detected as sparse conformational distributions (outliers) with a clustering algorithm. The detected outliers are intensively resampled with short-time molecular dynamics (MD) simulations. As a set of cycles, OFLOOD repeats selections of outliers and their conformational resampling. Herein, as an essential issue to be tackled to perform OFLOOD efficiently, a selection rule for outliers should be carefully specified. Generally, many outliers are detected from distributions on conformational subspaces with the clustering. Judging from its computational costs, it is unreasonable to select all the detected outliers upon the conformational resampling. Therefore, it is important to consider which outliers should be selected from the sparse distributions when restarting their short-time MD simulations with limited computational costs. In this sense, we investigated the conformational sampling efficiency of OFLOOD by changing the selection rules for outliers. To address the conformational sampling efficiency of OFLOOD depending on its selection rules, outliers to be resampled were selected by focusing their probability occurrences (populations of outliers). As a comparison, a random selection rule for outliers was also considered. Through the present assessment, the random selection of outliers showed the most efficient conformational sampling efficiency compared to the other OFLOOD trials using the biased selection rules, indicating that a variety of outliers should be selected and resampled during the OFLOOD cycles. In conclusion, the random outlier selection rule is the best strategy to perform OFLOOD efficiently.