Flexibility of the Bacterial Chaperone Trigger Factor in Microsecond-Timescale Molecular Dynamics Simulations

The bacterial chaperone trigger factor (TF) is the first chaperone to be encountered by a nascent protein chain as it emerges from the ribosome exit tunnel. Experimental results suggest that TF possesses considerable conformational flexibility, and in an attempt to provide an atomic-level view of this flexibility, we have performed independent 1.5-μs molecular dynamics simulations of TF in explicit solvent using two different simulation force fields (OPLS-AA/L and AMBER ff99SB-ILDN). Both simulations indicate that TF possesses tremendous flexibility, with huge excursions from the crystallographic conformation caused by reorientations of the protein’s constituent domains; both simulations also predict the formation of extensive contacts between TF’s PPIase domain and the Arm 1 domain that is involved in nascent-chain binding. In the OPLS simulation, however, TF rapidly settles into a very compact conformation that persists for at least 1 μs, whereas in the AMBER simulation, it remains highly dynamic; additional simulations in which the two force fields were swapped suggest that these differences are at least partly attributable to sampling issues. The simulation results provide potential rationalizations of a number of experimental observations regarding TF’s conformational behavior and have implications for using simulations to model TF’s function on translating ribosomes.     

Conformational plasticity and dynamics in the generic protein folding catalyst SlyD unraveled by single-molecule FRET

The relation between conformational dynamics and chemistry in enzyme catalysis recently has received increasing attention. While, in the past, the mechanochemical coupling was mainly attributed to molecular motors, nowadays, it seems that this linkage is far more general. Single-molecule fluorescence methods are perfectly suited to directly evidence conformational flexibility and dynamics. By labeling the enzyme SlyD, a member of peptidyl-prolyl cis-trans isomerases of the FK506 binding protein type with an inserted chaperone domain, with donor and acceptor fluorophores for single-molecule fluorescence resonance energy transfer, we directly monitor conformational flexibility and conformational dynamics between the chaperone domain and the FK506 binding protein domain. We find a broad distribution of distances between the labels with two main maxima, which we attribute to an open conformation and to a closed conformation of the enzyme. Correlation analysis demonstrates that the conformations exchange on a rate in the 100 Hz range. With the aid from Monte Carlo simulations, we show that there must be conformational flexibility beyond the two main conformational states. Interestingly, neither the conformational distribution nor the dynamics is significantly altered upon binding of substrates or other known binding partners. Based on these experimental findings, we propose a model where the conformational dynamics is used to search the conformation enabling the chemical step, which also explains the remarkable substrate promiscuity connected with a high efficiency of this class of peptidyl-prolyl cis-trans isomerases.     

Effects of force fields on the conformational and dynamic properties of amyloid β(1‐40) dimer explored by replica exchange molecular dynamics simulations

The conformational space and structural ensembles of amyloid beta (Aβ) peptides and their oligomers in solution are inherently disordered and proven to be challenging to study. Optimum force field selection for molecular dynamics (MD) simulations and the biophysical relevance of results are still unknown. We compared the conformational space of the Aβ(1‐40) dimers by 300 ns replica exchange MD simulations at physiological temperature (310 K) using: the AMBER‐ff99sb‐ILDN, AMBER‐ff99sb*‐ILDN, AMBER‐ff99sb‐NMR, and CHARMM22* force fields. Statistical comparisons of simulation results to experimental data and previously published simulations utilizing the CHARMM22* and CHARMM36 force fields were performed. All force fields yield sampled ensembles of conformations with collision cross sectional areas for the dimer that are statistically significantly larger than experimental results. All force fields, with the exception of AMBER‐ff99sb‐ILDN (8.8 ± 6.4%) and CHARMM36 (2.7 ± 4.2%), tend to overestimate the α‐helical content compared to experimental CD (5.3 ± 5.2%). Using the AMBER‐ff99sb‐NMR force field resulted in the greatest degree of variance (41.3 ± 12.9%). Except for the AMBER‐ff99sb‐NMR force field, the others tended to under estimate the expected amount of β‐sheet and over estimate the amount of turn/bend/random coil conformations. All force fields, with the exception AMBER‐ff99sb‐NMR, reproduce a theoretically expected β‐sheet‐turn‐β‐sheet conformational motif, however, only the CHARMM22* and CHARMM36 force fields yield results compatible with collapse of the central and C‐terminal hydrophobic cores from residues 17‐21 and 30‐36. Although analyses of essential subspace sampling showed only minor variations between force fields, secondary structures of lowest energy conformers are different.     

Examination of the quality of various force fields and solvation models for the equilibrium simulations of GA88 and GB88

Elucidating the relationship between sequence and conformation is essential for the understanding of functions of proteins. While sharing 88 % sequence identity and differing by only seven residues, GA88 and GB88 have completely different structures and serve as ideal systems for investigating the relationship between sequence and function. Benefiting from the continuous advancement of the computational ability of modern computers, molecular dynamics (MD) simulation is now playing an increasingly important role in the study of proteins. However, the reliability of MD simulations is limited by the accuracy of the force fields and solvent model approximations. In this work, several AMBER force fields (AMBER03, AMBER99SB, AMBER12SB, AMBER14SB, AMBER96) and solvent models (TIP3P, IGB5, IGB7, IGB8) have been employed in the simulations of GA88 and GB88. The statistical results from 19 simulations show that GA88 and GB88 both adopt more compact structures than the native structures. GB88 is more stable than GA88 regardless of the force fields and solvent models utilized. Most of the simulations overestimated the salt bridge interaction. The combination of AMBER14SB force field and IGB8 solvent model shows the best overall performance in the simulations of both GA88 and GB88. AMBER03 and AMBER12SB also yield reasonable results but only in the TIP3P explicit solvent model.     

Effects of different force fields on the structural character of α synuclein β-hairpin peptide (35–56) in aqueous environment

The hallmark of Parkinson’s disease (PD) is the intracellular protein aggregation forming Lewy Bodies (LB) and Lewy neuritis which comprise mostly of a protein, alpha synuclein (α-syn). Molecular dynamics (MD) simulation methods can augment experimental techniques to understand misfolding and aggregation pathways with atomistic resolution. The quality of MD simulations for proteins and peptides depends greatly on the accuracy of empirical force fields. The aim of this work is to investigate the effects of different force fields on the structural character of β hairpin fragment of α-syn (residues 35–56) peptide in aqueous solution. Six independent MD simulations are done in explicit solvent using, AMBER03, AMBER99SB, GROMOS96 43A1, GROMOS96 53A6, OPLS-AA, and CHARMM27 force fields with CMAP corrections. The performance of each force field is assessed from several structural parameters such as root mean square deviation (RMSD), root mean square fluctuation (RMSF), radius of gyration (Rg), solvent accessible surface area (SASA), formation of β-turn, the stability of folded β-hairpin structure, and the favourable conformations obtained for different force fields. In this study, CMAP correction of CHARMM27 force field is found to overestimate the helical conformation, while GROMOS96 53A6 is found to most successfully capture the conformational dynamics of α-syn β-hairpin fragment as elicited from NMR.     

MMC and LD simulations of α-D-Manp-(1→2)-β-D-Glcp-OMe: comparison to long-range heteronuclear NMR coupling constants and to the crystal structure

The conformational flexibility and the dynamics of -D-Manp-(12)--D-Glcp-OMe have been investigated by Metropolis Monte Carlo (MMC) and Langevin dynamics (LD) simulations. The two simulation techniques employ different force fields, namely the HSEA force field and a CHARMm-based force field. The former shows less conformational flexibility than the latter, in which a multiple energy minima conformational space is sampled. Long-range heteronuclear nuclear magnetic resonance (NMR) coupling constants have been measured by selective excitations of the carbons at the glycosidic linkage. Calculated 3JC, H values from MMC and LD simulations show excellent agreement to those from NMR experiments. The X-ray crystal structure has a conformation within a region of the conformational space populated in both force fields.     

Are Current Molecular Dynamics Force Fields too Helical?

Accurate force fields are essential for the success of molecular dynamics simulations. In apparent contrast to the conformational preferences of most force fields, recent NMR experiments suggest that short polyalanine peptides in water populate the polyproline II structure almost exclusively. To investigate this apparent contradiction, with its ramifications for the assessment of molecular force fields and the structure of unfolded proteins, we performed extensive simulations of Ala₅ in water (∼5 μs total time), using twelve different force fields and three different peptide terminal groups. Using either empirical or density-functional-based Karplus relations for the J-couplings, we find that most current force fields do overpopulate the α-region, with quantitative results depending on the choice of Karplus relation and on the peptide termini. Even after reweighting to match experiment, we find that Ala₅ retains significant α- and β-populations. In fact, several force fields match the experimental data well before reweighting and have a significant helical population. We conclude that radical changes to the best current force fields are not necessary, based on the NMR data. Nevertheless, experiments on short peptides open the way toward the systematic improvement of current simulation models.     

Folding simulations of three proteins having all α-helix,all β-strand and α/β-structures

We performed folding simulations of three proteins using four force fields, AMBER parm96, AMBER parm99, CHARMM 27 and OPLS-AA/L, in order to examine the features of these force fields. We studied three proteins, protein A (all α-helix), cold-shock protein (all β-strand) and protein G (α/β-structures), for the folding simulations. For the simulation, we used the simulated annealing molecular dynamics method, which was performed 50 times for each protein using the four force fields. The results showed that the secondary-structure-forming tendencies are largely different among the four force fields. AMBER parm96 favours β-bridge structures and extended β-strand structures, and AMBER parm99 favours α-helix structures and 3₁₀-helix structures. CHARMM 27 slightly favours α-helix structures, and there are also π-helix and β-bridge structures. OPLS-AA/L favours α-helix structures and 3₁₀-helix structures.     

Structural equilibrium of DNA represented with different force fields.

We have recently indicated preliminary evidence of different equilibrium average structures with the CHARMM and AMBER force fields in explicit solvent molecular dynamics simulations on the DNA duplex d(C5T5) . d(A5G5) (Feig, M. and B.M. Pettitt, 1997, Experiment vs. Force Fields: DNA conformation from molecular dynamics simulations. J. Phys. Chem. B. (101:7361-7363). This paper presents a detailed comparison of DNA structure and dynamics for both force fields from extended simulation times of 10 ns each. Average structures display an A-DNA base geometry with the CHARMM force field and a base geometry that is intermediate between A- and B-DNA with the AMBER force field. The backbone assumes B form on both strands with the AMBER force field, while the CHARMM force field produces heterogeneous structures with the purine strand in A form and the pyrimidine strand in dynamical equilibrium between A and B conformations. The results compare well with experimental data for the cytosine/guanine part but fail to fully reproduce an overall B conformation in the thymine/adenine tract expected from crystallographic data, particularly with the CHARMM force field. Fluctuations between A and B conformations are observed on the nanosecond time scale in both simulations, particularly with the AMBER force field. Different dynamical behavior during the first 4 ns indicates that convergence times of several nanoseconds are necessary to fully establish a dynamical equilibrium in all structural quantities on the time scale of the simulations presented here.     

Comparison of force fields for Alzheimer's A : A case study for intrinsically disordered proteins

Intrinsically disordered proteins are essential for biological processes such as cell signalling, but are also associated to devastating diseases including Alzheimer's disease, Parkinson's disease or type II diabetes. Because of their lack of a stable three‐dimensional structure, molecular dynamics simulations are often used to obtain atomistic details that cannot be observed experimentally. The applicability of molecular dynamics simulations depends on the accuracy of the force field chosen to represent the underlying free energy surface of the system. Here, we use replica exchange molecular dynamics simulations to test five modern force fields, OPLS, AMBER99SB, AMBER99SB*ILDN, AMBER99SBILDN‐NMR and CHARMM22*, in their ability to model Aβ₄₂, an intrinsically disordered peptide associated with Alzheimer's disease, and compare our results to nuclear magnetic resonance (NMR) experimental data. We observe that all force fields except AMBER99SBILDN‐NMR successfully reproduce local NMR observables, with CHARMM22* being slightly better than the other force fields.     

Force field influences in beta-hairpin folding simulations

All-atom force fields are now routinely used for more detailed understanding of protein folding mechanisms. However, it has been pointed out that use of all-atom force fields does not guarantee more accurate representations of proteins; in fact, sometimes it even leads to biased structural distributions. Indeed, several issues remain to be solved in force field developments, such as accurate treatment of implicit solvation for efficient conformational sampling and proper treatment of backbone interactions for secondary structure propensities. In this study, we first investigate the quality of several recently improved backbone interaction schemes in AMBER for folding simulations of a beta-hairpin peptide, and further study their influences on the peptide's folding mechanism. Due to the significant number of simulations needed for a thorough analysis of tested force fields, the implicit Poisson-Boltzmann solvent was used in all simulations. The chosen implicit solvent was found to be reasonable for studies of secondary structures based on a set of simulations of both alpha-helical and beta-hairpin peptides with the TIP3P explicit solvent as benchmark. Replica exchange molecular dynamics was also utilized for further efficient conformational sampling. Among the tested AMBER force fields, ff03 and a revised ff99 force field were found to produce structural and thermodynamic data in comparably good agreement with the experiment. However, detailed folding pathways, such as the order of backbone hydrogen bond zipping and the existence of intermediate states, are different between the two force fields, leading to force field-dependent folding mechanisms.     

Bayesian Energy Landscape Tilting: Towards Concordant Models of Molecular Ensembles

Predicting biological structure has remained challenging for systems such as disordered proteins that take on myriad conformations. Hybrid simulation/experiment strategies have been undermined by difficulties in evaluating errors from computational model inaccuracies and data uncertainties. Building on recent proposals from maximum entropy theory and nonequilibrium thermodynamics, we address these issues through a Bayesian energy landscape tilting (BELT) scheme for computing Bayesian hyperensembles over conformational ensembles. BELT uses Markov chain Monte Carlo to directly sample maximum-entropy conformational ensembles consistent with a set of input experimental observables. To test this framework, we apply BELT to model trialanine, starting from disagreeing simulations with the force fields ff96, ff99, ff99sbnmr-ildn, CHARMM27, and OPLS-AA. BELT incorporation of limited chemical shift and ³J measurements gives convergent values of the peptide’s α, β, and PP_II conformational populations in all cases. As a test of predictive power, all five BELT hyperensembles recover set-aside measurements not used in the fitting and report accurate errors, even when starting from highly inaccurate simulations. BELT’s principled framework thus enables practical predictions for complex biomolecular systems from discordant simulations and sparse data.     

Conformational dynamics of Peb4 exhibit “mother’s arms” chain model: a molecular dynamics study

Peb4 from Campylobacter jejuni is an intertwined dimeric, periplasmic holdase, which also exhibits peptidyl prolyl cis/trans isomerase (PPIase) activity. Peb4 gene deletion alters the outer membrane protein profile and impairs cellular adhesion and biofilm formation for C. jejuni. Earlier crystallographic study has proposed that the PPIase domains are flexible and might form a cradle for holding the substrate and these aspects of Peb4 were explored using sub-microsecond molecular dynamics simulations in solution environment. Our simulations have revealed that PPIase domains are highly flexible and undergo a large structural change where they move apart from each other by 8 nm starting at .5 nm. Further, this large conformational change renders Peb4 as a compact protein with crossed-over conformation, forms a central cavity, which can “cradle” the target substrate. As reported for other chaperone proteins, flexibility of linker region connecting the chaperone and PPIase domains is key to forming the “crossed-over” conformation. The conformational transition of the Peb4 protein from the X-ray structure to the crossed-over conformation follows the “mother’s arms” chain model proposed for the FkpA chaperone protein. Our results offer insights into how Peb4 and similar chaperones can use the conformational heterogeneity at their disposal to perform its much-revered biological function.     

Conformational flexibility of the MHC class I alpha1-alpha2 domain in peptide bound and free states: a molecular dynamics simulation study

Major histocompatibility complex class I proteins play a key role in the recognition and presentation of peptide antigens to the host immune system. The structure of various major histocompatibility complex class I proteins has been determined experimentally in complex with several antigenic peptides. However, the structure in the unbound (empty) form is not known. To study the conformational dynamics of the empty major histocompatibility complex class I molecule comparative molecular dynamics simulations have been performed starting from the crystal structure of a peptide bound class I peptide-binding domain in the presence and absence of a peptide ligand. Simulations including the bound peptide stayed close to the experimental start structure at both simulation temperatures (300 and 355 K) during the entire simulation of 26 ns. Several independent simulations in the absence of peptide indicate that the empty domain may not adopt a single defined conformation but is conformationally significantly more heterogeneous in particular within the alpha-helices that flank the peptide binding cleft. The calculated conformational dynamics along the protein chain correlate well with available spectroscopic data and with the observed site-specific sensitivity of the empty class I protein to proteolytic digestion. During the simulations at 300 K the binding region for the peptide N-terminus stayed close to the conformation in the bound state, whereas the anchor region for the C-terminus showed significantly larger conformational fluctuations. This included a segment at the beginning of the second alpha-helix in the domain that is likely to be involved in the interaction with the chaperone protein tapasin during the peptide-loading process. The simulation studies further indicate that peptide binding at the C- and N-terminus may follow different mechanisms that involve different degrees of induced conformational changes in the peptide-binding domain. In particular binding of the peptide C-terminus may require conformational stabilization by chaperone proteins during peptide loading.     

Force Field Bias in Protein Folding Simulations

Long timescale (>1 μs) molecular dynamics simulations of protein folding offer a powerful tool for understanding the atomic-scale interactions that determine a protein's folding pathway and stabilize its native state. Unfortunately, when the simulated protein fails to fold, it is often unclear whether the failure is due to a deficiency in the underlying force fields or simply a lack of sufficient simulation time. We examine one such case, the human Pin1 WW domain, using the recently developed deactivated morphing method to calculate free energy differences between misfolded and folded states. We find that the force field we used favors the misfolded states, explaining the failure of the folding simulations. Possible further applications of deactivated morphing and implications for force field development are discussed.     

Glycan flexibility: insights into nanosecond dynamics from a microsecond molecular dynamics simulation explaining an unusual nuclear Overhauser effect

An atomistic all-atom molecular dynamics simulation of the trisaccharide β-d-ManpNAc-(1→4)[α-d-Glcp-(1→3)]-α-l-Rhap-OMe with explicit solvent molecules has been carried out. The trisaccharide represents a model for the branching region of the O-chain polysaccharide of a strain from Aeromonas salmonicida. The extensive MD simulations having a 1-μs duration revealed a conformational dynamics process on the nanosecond time scale, that is, a ‘time window’ not extensively investigated for carbohydrates to date. The results obtained from the MD simulation underscore the predictive power of molecular simulations in studies of biomolecular systems and also explain an unusual nuclear Overhauser effect originating from conformational exchange.     

CD4 binding partially locks the bridging sheet in gp120 but leaves the beta2/3 strands flexible

The structure of the free form HIV gp120, critical for therapeutic agent development, is unavailable due to its high flexibility. Previous thermodynamic data, structural analysis and simulation results have suggested a large conformational change in the core domain upon CD4 binding. The bridging sheet, which consists of four beta-strands with beta20/21 nestling against the inner/outer domains and beta2/3 facing outward, more exposed to the solvent, was proposed to be unfolded in the native state. In order to test this proposition and to characterize the native conformations, we performed potential mean force (PMF) molecular dynamics (MD) simulations on the CD4-bound crystal structure. We pushed the bridging sheet away from the inner and outer domain to explore the accessible conformational space for the bridging sheet. In addition, we performed conventional MD simulations on structures with the bridging sheet partially unfolded to investigate the stability of the association between the inner and outer domains. Based on the free energy profiles, we find that the whole bridging sheet is unlikely to unfold without other concurrent conformational changes. On the other hand, the partial bridging sheet, beta strands 2/3, can switch its conformation from the folded to the unfolded state. Furthermore, relaxation of conformation with partially unfolded bridging sheet through MD simulations leads to a conformation with beta strands 20/21 quickly re-anchoring against the inner and outer domains. Such a conformation, although lacking some of the hydrophobic interactions present in the CD4-bound structure, displayed high stability as further indicated by other restrained MD simulations. The relevance of this conformation to the free form structure and the pathway for conformational change from the free form to the CD4-bound structure is discussed in detail in light of the available unliganded SIV gp120 crystal structure.     

Derivation of original RESP atomic partial charges for MD simulations of the LDAO surfactant with AMBER: applications to a model of micelle and a fragment of the lipid kinase PI4KA

In this paper, we describe the derivation and the validation of original RESP atomic partial charges for the N, N-dimethyl-dodecylamine oxide (LDAO) surfactant. These charges, designed to be fully compatible with all the AMBER force fields, are at first tested against molecular dynamics simulations of pure LDAO micelles and with a fragment of the lipid kinase PIK4A (DI) modeled with the QUARK molecular modeling server. To model the micelle, we used two distinct AMBER force fields (i.e. Amber99SB and Lipid14) and a variety of starting conditions. We find that the micelle structural properties (such as the shape, size, the LDAO headgroup hydration, and alkyl chain conformation) slightly depend on the force field but not on the starting conditions and more importantly are in good agreement with experiments and previous simulations. We also show that the Lipid14 force field should be used instead of the Amber99SB one to better reproduce the C(sp3)C(sp3)C(sp3)C(sp3) conformation in the surfactant alkyl chain. Concerning the simulations with LDAO-DI protein, we carried out different runs at two NaCl concentrations (i.e. 0 and 300 mM) to mimic, in the latter case, the experimental conditions. We notice a small dependence of the simulation results with the LDAO parameters and the salt concentration. However, we find that in the simulations, three out of four tryptophans of the DI protein are not accessible to water in agreement with our fluorescence spectroscopy experiments reported in the paper.     

Molecular dynamics re-refinement of two different small RNA loop structures using the original NMR data suggest a common structure

Restrained molecular dynamics simulations are a robust, though perhaps underused, tool for the end-stage refinement of biomolecular structures. We demonstrate their utility-using modern simulation protocols, optimized force fields, and inclusion of explicit solvent and mobile counterions-by re-investigating the solution structures of two RNA hairpins that had previously been refined using conventional techniques. The structures, both domain 5 group II intron ribozymes from yeast ai5γ and Pylaiella littoralis, share a nearly identical primary sequence yet the published 3D structures appear quite different. Relatively long restrained MD simulations using the original NMR restraint data identified the presence of a small set of violated distance restraints in one structure and a possibly incorrect trapped bulge nucleotide conformation in the other structure. The removal of problematic distance restraints and the addition of a heating step yielded representative ensembles with very similar 3D structures and much lower pairwise RMSD values. Analysis of ion density during the restrained simulations helped to explain chemical shift perturbation data published previously. These results suggest that restrained MD simulations, with proper caution, can be used to "update" older structures or aid in the refinement of new structures that lack sufficient experimental data to produce a high quality result. Notable cautions include the need for sufficient sampling, awareness of potential force field bias (such as small angle deviations with the current AMBER force fields), and a proper balance between the various restraint weights.