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.     

Force Distribution Reveals Signal Transduction in E.?coli Hsp90

Heat-shock protein 90 (Hsp90) is an ubiquitous chaperone that is essential for cell function in that it promotes client-protein folding and stabilization. Its function is tightly controlled by an ATP-dependent large conformational transition between the open and closed states of the Hsp90 dimer. The underlying allosteric pathway has remained largely unknown, but it is revealed here in atomistic detail for the Escherichia coli homolog HtpG. Using force-distribution analysis based on molecular-dynamics simulations (>1 μs in total), we identify an internal signaling pathway that spans from the nucleotide-binding site to an ∼2.3-nm-distant region in the HtpG middle domain, that serves as a dynamic hinge region, and to a putative client-protein-binding site in the middle domain. The force transmission is triggered by ATP capturing a magnesium ion and thereby rotating and bending a proximal long α-helix, which represents the major force channel into the middle domain. This allosteric mechanism is, with statistical significance, distinct from the dynamics in the ADP and apo states. Tracking the distribution of forces is likely to be a promising tool for understanding and guiding experiments of complex allosteric proteins in general.     

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.     

In Silico Adoption of an Orphan Nuclear Receptor NR4A1

A 4.1μs molecular dynamics simulation of the NR4A1 (hNur77) apo-protein has been undertaken and a previously undetected druggable pocket has become apparent that is located remotely from the ‘traditional’ nuclear receptor ligand-binding site. A NR4A1/bis-indole ligand complex at this novel site has been found to be stable over 1 μs of simulation and to result in an interesting conformational transmission to a remote loop that has the capacity to communicate with a NBRE within a RXR-α/NR4A1 heterodimer. Several features of the simulations undertaken indicate how NR4A1 can be affected by alternate-site modulators.     

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.     

Extended conformational states dominate the Hsp90 chaperone dynamics

The heat shock protein 90 (Hsp90) is a molecular chaperone central to client protein folding and maturation in eukaryotic cells. During its chaperone cycle, Hsp90 undergoes ATPase-coupled large-scale conformational changes between open and closed states, where the N-terminal and middle domains of the protein form a compact dimerized conformation. However, the molecular principles of the switching motion between the open and closed states remain poorly understood. Here we show by integrating atomistic and coarse-grained molecular simulations with small-angle X-ray scattering experiments and NMR spectroscopy data that Hsp90 exhibits rich conformational dynamics modulated by the charged linker, which connects the N-terminal with the middle domain of the protein. We show that the dissociation of these domains is crucial for the conformational flexibility of the open state, with the separation distance controlled by a β-sheet motif next to the linker region. Taken together, our results suggest that the conformational ensemble of Hsp90 comprises highly extended states, which could be functionally crucial for client processing.     

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.     

Structural Determinants of the Mechanical Stability of α-Catenin

α-Catenin plays a crucial role in cadherin-mediated adhesion by binding to β-catenin, F-actin, and vinculin, and its dysfunction is linked to a variety of cancers and developmental disorders. As a mechanotransducer in the cadherin complex at intercellular adhesions, mechanical and force-sensing properties of α-catenin are critical to its proper function. Biochemical data suggest that α-catenin adopts an autoinhibitory conformation, in the absence of junctional tension, and biophysical studies have shown that α-catenin is activated in a tension-dependent manner that in turn results in the recruitment of vinculin to strengthen the cadherin complex/F-actin linkage. However, the molecular switch mechanism from autoinhibited to the activated state remains unknown for α-catenin. Here, based on the results of an aggregate of 3 μs of molecular dynamics simulations, we have identified a dynamic salt-bridge network within the core M region of α-catenin that may be the structural determinant of the stability of the autoinhibitory conformation. According to our constant-force steered molecular dynamics simulations, the reorientation of the MII/MIII subdomains under force may constitute an initial step along the transition pathway. The simulations also suggest that the vinculin-binding domain (subdomain MI) is intrinsically much less stable than the other two subdomains in the M region (MII and MIII). Our findings reveal several key insights toward a complete understanding of the multistaged, force-induced conformational transition of α-catenin to the activated conformation.     

Molecular Dynamics of Class A β-lactamases—Effects of Substrate Binding

The effects of substrate binding on class A β-lactamase dynamics were studied using molecular dynamics simulations of two model enzymes; 40 100-ns trajectories of the free and substrate-bound forms of TEM-1 (with benzylpenicillin) and PSE-4 (with carbenicillin) were recorded (totaling 4.0 μs). Substrates were parameterized with the CHARMM General Force Field. In both enzymes, the Ω loop exhibits a marked flexibility increase upon substrate binding, supporting the hypothesis of substrate gating. However, specific interactions that are formed or broken in the Ω loop upon binding differ between the two enzymes: dynamics are conserved, but not specific interactions. Substrate binding also has a global structuring effect on TEM-1, but not on PSE-4. Changes in TEM-1’s normal modes show long-range effects of substrate binding on enzyme dynamics. Hydrogen bonds observed in the active site are mostly preserved upon substrate binding, and new, transient interactions are also formed. Agreement between NMR relaxation parameters and our theoretical results highlights the dynamic duality of class A β-lactamases: enzymes that are highly structured on the ps-ns timescale, with important flexibility on the μs-ms timescale in regions such as the Ω loop.     

Structural Model of the Proline-Rich Domain of Huntingtin Exon-1 Fibrils

Huntington’s disease is a heritable neurodegenerative disease that is caused by a CAG expansion in the first exon of the huntingtin gene. This expansion results in an elongated polyglutamine domain that increases the propensity of huntingtin exon-1 to form cross-β fibrils. Although the polyglutamine domain is important for fibril formation, the dynamic, C-terminal proline-rich domain (PRD) of huntingtin exon-1 makes up a large fraction of the fibril surface. Because potential fibril toxicity has to be mediated by interactions of the fibril surface with its cellular environment, we wanted to model the conformational space adopted by the PRD. We ran 800-ns long molecular dynamics simulations of the PRD using an explicit water model optimized for intrinsically disordered proteins. These simulations accurately predicted our previous solid-state NMR data and newly acquired electron paramagnetic resonance double electron-electron resonance distances, lending confidence in their accuracy. The simulations show that the PRD generally forms an imperfect polyproline (polyP) II helical conformation. The two polyP regions within the PRD stay in a polyP II helix for most of the simulation, whereas occasional kinks in the proline-rich linker region cause an overall bend in the PRD structure. The dihedral angles of the glycine at the end of the second polyP region are very variable, effectively decoupling the highly dynamic 12 C-terminal residues from the rest of the PRD.     

The Free Energy Profile of Tubulin Straight-Bent Conformational Changes,with Implications for Microtubule Assembly and Drug Discovery

αβ-tubulin dimers need to convert between a ‘bent’ conformation observed for free dimers in solution and a ‘straight’ conformation required for incorporation into the microtubule lattice. Here, we investigate the free energy landscape of αβ-tubulin using molecular dynamics simulations, emphasizing implications for models of assembly, and modulation of the conformational landscape by colchicine, a tubulin-binding drug that inhibits microtubule polymerization. Specifically, we performed molecular dynamics, potential-of-mean force simulations to obtain the free energy profile for unpolymerized GDP-bound tubulin as a function of the ∼12° intradimer rotation differentiating the straight and bent conformers. Our results predict that the unassembled GDP-tubulin heterodimer exists in a continuum of conformations ranging between straight and bent, but, in agreement with existing structural data, suggests that an intermediate bent state has a lower free energy (by ∼1 kcal/mol) and thus dominates in solution. In agreement with predictions of the lattice model of microtubule assembly, lateral binding of two αβ-tubulins strongly shifts the conformational equilibrium towards the straight state, which is then ∼1 kcal/mol lower in free energy than the bent state. Finally, calculations of colchicine binding to a single αβ-tubulin dimer strongly shifts the equilibrium toward the bent states, and disfavors the straight state to the extent that it is no longer thermodynamically populated.     

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.     

Dynamic Regulation of α-Actinin’s Calponin Homology Domains on F-Actin

α-Actinin is an essential actin cross-linker involved in cytoskeletal organization and dynamics. The molecular conformation of α-actinin’s actin-binding domain (ABD) regulates its association with actin and thus mutations in this domain can lead to severe pathogenic conditions. A point mutation at lysine 255 in human α-actinin-4 to glutamate increases the binding affinity resulting in stiffer cytoskeletal structures. The role of different ABD conformations and the effect of K255E mutation on ABD conformations remain elusive. To evaluate the impact of K255E mutation on ABD binding to actin we use all-atom molecular dynamics and free energy calculation methods and study the molecular mechanism of actin association in both wild-type α-actinin and in the K225E mutant. Our models illustrate that the strength of actin association is indeed sensitive to the ABD conformation, predict the effect of K255E mutation—based on simulations with the K237E mutant chicken α-actinin—and evaluate the mechanism of α-actinin binding to actin. Furthermore, our simulations showed that the calmodulin domain binding to the linker region was important for regulating the distance between actin and ABD. Our results provide valuable insights into the molecular details of this critical cellular phenomenon and further contribute to an understanding of cytoskeletal dynamics in health and disease.     

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.     

Highly sampled tetranucleotide and tetraloop motifs enable evaluation of common RNA force fields

Recent modifications and improvements to standard nucleic acid force fields have attempted to fix problems and issues that have been observed as longer timescale simulations have become routine. Although previous work has shown the ability to fold the UUCG stem–loop structure, until now no group has attempted to quantify the performance of current force fields using highly converged structural populations of the tetraloop conformational ensemble. In this study, we report the use of multiple independent sets of multidimensional replica exchange molecular dynamics (M-REMD) simulations with different initial conditions to generate well-converged conformational ensembles for the tetranucleotides r(GACC) and r(CCCC), as well as the larger UUCG tetraloop motif. By generating what is to our knowledge the most complete RNA structure ensembles reported to date for these systems, we remove the coupling between force field errors and errors due to incomplete sampling, providing a comprehensive comparison between current top-performing MD force fields for RNA. Of the RNA force fields tested in this study, none demonstrate the ability to correctly identify the most thermodynamically stable structure for all three systems. We discuss the deficiencies present in each potential function and suggest areas where improvements can be made. The results imply that although “short” (nsec-μsec timescale) simulations may stay close to their respective experimental structures and may well reproduce experimental observables, inevitably the current force fields will populate alternative incorrect structures that are more stable than those observed via experiment.     

The chaperone domain BRICHOS prevents CNS toxicity of amyloid-β peptide in Drosophila melanogaster

Aggregation of the amyloid-β peptide (Aβ) into toxic oligomers and amyloid fibrils is linked to the development of Alzheimer’s disease (AD). Mutations of the BRICHOS chaperone domain are associated with amyloid disease and recent in vitro data show that BRICHOS efficiently delays Aβ42 oligomerization and fibril formation. We have generated transgenic Drosophila melanogaster flies that express the Aβ42 peptide and the BRICHOS domain in the central nervous system (CNS). Co-expression of Aβ42 and BRICHOS resulted in delayed Aβ42 aggregation and dramatic improvements of both lifespan and locomotor function compared with flies expressing Aβ42 alone. Moreover, BRICHOS increased the ratio of soluble:insoluble Aβ42 and bound to deposits of Aβ42 in the fly brain. Our results show that the BRICHOS domain efficiently reduces the neurotoxic effects of Aβ42, although significant Aβ42 aggregation is taking place. We propose that BRICHOS-based approaches should be explored with an aim towards the future prevention and treatment of AD.KEY WORDS: Amyloid, Alzheimer’s disease, Protein misfolding, Chaperone     

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.     

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.