REMD and umbrella sampling simulations to probe the energy barrier of the folding pathways of engrailed homeodomain

Effects of number of parallel runs and frequency of bias-strength replacement in generalized ensemble molecular dynamics simulations

Background

The generalized ensemble approach with the molecular dynamics (MD) method has been widely utilized. This approach usually has two features. (i) A bias potential, whose strength is replaced during a simulation, is applied. (ii) Sampling can be performed by many parallel runs of simulations. Although the frequency of the bias-strength replacement and the number of parallel runs can be adjusted, the effects of these settings on the resultant ensemble remain unclear.

Method

In this study, we performed multicanonical MD simulations for a foldable mini-protein (Trp-cage) and two unstructured peptides (8- and 20-residue poly-glutamic acids) with various settings.

Results

As a result, running many short simulations yielded robust results for the Trp-cage model. Regarding the frequency of the bias-potential replacement, although using a high frequency enhanced the traversals in the potential energy space, it did not promote conformational changes in all the systems.

Influence of various parameters in the replica-exchange molecular dynamics method: Number of replicas, replica-exchange frequency, and thermostat coupling time constant

The replica-exchange molecular dynamics (REMD) method has been used for conformational sampling of various biomolecular systems. To maximize sampling efficiency, some adjustable parameters must be optimized. Although it is agreed that shorter intervals between the replica-exchange attempts enhance traversals in the temperature space, details regarding the artifacts caused by these short intervals are controversial. In this study, we revisit this problem by performing REMD simulations on an alanine octapeptide in an implicit solvent. Fifty different sets of conditions, which are a combination of five replica-exchange periods, five different numbers of replicas, and two thermostat coupling time constants, were investigated. As a result, although short replica-exchange intervals enhanced the traversals in the temperature space, they led to artifacts in the ensemble average of the temperature, potential energy, and helix content. With extremely short replica-exchange intervals, i.e., attempted at every time step, the ensemble average of the temperature deviated from the thermostat temperature by ca. 7 K. Differences in the ensembles were observed even for larger replica-exchange intervals (between 100 and 1,000 steps). In addition, the shorter thermostat coupling time constant reduced the artifacts found when short replica-exchange intervals were used, implying that these artifacts are caused by insufficient thermal relaxation between the replica-exchange events. Our results will be useful to reduce the artifacts found in REMD simulations by adjusting some key parameters.

Significance of NS5B Substitutions in Genotype 1b Hepatitis C Virus Evaluated by Bioinformatics Analysis

To evaluate the effects of HCV NS5B amino acid substitutions on treatment outcome in Ledipasvir (LDV)/Sofosbuvir (SOF) for Japanese patients with genotype 1b HCV infection, NS5B sequences were examined in i) seven patients experiencing virologic failure after LDV/SOF in real-world practice, ii) 109 SOF-naïve patients, iii) 165 patients enrolled in Phase-3 LDV/SOF trial. A218S and C316N were detected in all patients with viral relapse; the percentages of these substitutions in SOF-naïve patients were 64.2% and 55.0%, respectively. Genotype 1b HCV strains with NS5B-C316N mutation were located in the leaves different from those in which HCV strains without such substitutions were present on the phylogenetic tree. Structural modeling revealed that amino acid 218 was located on the surface of the NTP tunnel. Free energy analysis based on molecular dynamics simulations demonstrated that the free energy required to pass through the tunnel was larger for triphosphate SOF than for UTP in NS5B polymerase carrying A218S, but not in wild-type. However, no susceptibility change was observed for these substitutions to SOF in replicon assay. Furthermore, the SVR rate was 100% in patients enrolled the Phase-3 trial. In conclusion, NS5B A218S and C316N were detected in all patients who relapsed following LDV/SOF in real-world practice. These substitutions did not impact the overall SVR rate after LDV/SOF, however, further studies are needed to elucidate the impact of these substitutions.

Coupled Mutations-Enabled Glycerol Transportation in an Aquaporin Z Mutant

Aquaporins are transmembrane channel proteins with key function being transportation of water or other small substrates. Escherichia coli Aqp Z transports water molecules only, whereas Glp F is permeable to glycerol. It is intriguing to explore the possibility to induce glycerol permeability in Aqp Z by targeted mutations. The Aqp Z mutants with mutated selectivity filter (SF) residues exhibit poor permeability for both glycerol and water. For addressing the complexity of protein systems, pair correlation information in protein sequence analyses is instructive to identify residues that are coupled by coevolution and motion. In this study, we analyze the correlation between residues and unravel the clustering patterns of coupled residues, beyond SF residues, in aquaglyceroporins (AQGPs). The identified coupled motifs are proposed to be sequenced into aquaporin (Aqp Z) to introduce glycerol permeability. These residues are located in the vicinity of SF region, C-loop, and M6-M7 linkage domain. Significant enlargement of SF pore size of the proposed Aqp Z mutant is observed by an all-atom replica exchange molecular dynamics simulation, which is critical to facilitate considerable glycerol passage as characterized in calculated free-energy landscapes. Clearly, the hidden connections among residues play crucial roles in water/glycerol selectivity. In contrast, single-site mutation-based scheme may even lead to undesirable effects in AQGPs, such as the blocking of water transportation by aromatic π-stacked gate. As demonstrated in this work, the pair correlation analysis guided rational mutagenesis provides a feasible strategy to modulate proteins' functions.

Exploring the Origin of Differential Binding Affinities of Human Tubulin Isotypes αβII, αβIII and αβIV for DAMA-Colchicine Using Homology Modelling, Molecular Docking and Molecular Dynamics Simulations

Tubulin isotypes are found to play an important role in regulating microtubule dynamics. The isotype composition is also thought to contribute in the development of drug resistance as tubulin isotypes show differential binding affinities for various anti-cancer agents. Tubulin isotypes αβII, αβIII and αβIV show differential binding affinity for colchicine. However, the origin of differential binding affinity is not well understood at the molecular level. Here, we investigate the origin of differential binding affinity of a colchicine analogue N-deacetyl-N-(2-mercaptoacetyl)-colchicine (DAMA-colchicine) for human αβII, αβIII and αβIV isotypes, employing sequence analysis, homology modeling, molecular docking, molecular dynamics simulation and MM-GBSA binding free energy calculations. The sequence analysis study shows that the residue compositions are different in the colchicine binding pocket of αβII and αβIII, whereas no such difference is present in αβIV tubulin isotypes. Further, the molecular docking and molecular dynamics simulations results show that residue differences present at the colchicine binding pocket weaken the bonding interactions and the correct binding of DAMA-colchicine at the interface of αβII and αβIII tubulin isotypes. Post molecular dynamics simulation analysis suggests that these residue variations affect the structure and dynamics of αβII and αβIII tubulin isotypes, which in turn affect the binding of DAMA-colchicine. Further, the binding free-energy calculation shows that αβIV tubulin isotype has the highest binding free-energy and αβIII has the lowest binding free-energy for DAMA-colchicine. The order of binding free-energy for DAMA-colchicine is αβIV ≃ αβII >> αβIII. Thus, our computational approaches provide an insight into the effect of residue variations on differential binding of αβII, αβIII and αβIV tubulin isotypes with DAMA-colchicine and may help to design new analogues with higher binding affinities for tubulin isotypes.

Important roles of hydrophobic interactions in folding and charge interactions in misfolding of α-helix bundle protein

An enhanced-sampling molecular dynamics simulation is presented to quantitatively demonstrate the important roles of hydrophobic and charge interactions in the folding and misfolding of α-helix bundle protein, respectively.