A molecular dynamics study of the 41-56 beta-hairpin from B1 domain of protein G.

The structural and dynamical behavior of the 41-56 beta-hairpin from the protein G B1 domain (GB1) has been studied at different temperatures using molecular dynamics (MD) simulations in an aqueous environment. The purpose of these simulations is to establish the stability of this hairpin in view of its possible role as a nucleation site for protein folding. The conformation of the peptide in the crystallographic structure of the protein GB1 (native conformation) was lost in all simulations. The new equilibrium conformations are stable for several nanoseconds at 300K (>10 ns), 350 K (>6.5 ns), and even at 450 K (up to 2.5 ns). The new structures have very similar hairpin-like conformations with properties in agreement with available experimental nuclear Overhauser effect (NOE) data. The stability of the structure in the hydrophobic core region during the simulations is consistent with the experimental data and provides further evidence for the role played by hydrophobic interactions in hairpin structures. Essential dynamics analysis shows that the dynamics of the peptide at different temperatures spans basically the same essential subspace. The main equilibrium motions in this subspace involve large fluctuations of the residues in the turn and ends regions. Of the six interchain hydrogen bonds, the inner four remain stable during the simulations. The space spanned by the first two eigenvectors, as sampled at 450 K, includes almost all of the 47 different hairpin structures found in the database. Finally, analysis of the hydration of the 300 K average conformations shows that the hydration sites observed in the native conformation are still well hydrated in the equilibrium MD ensemble.     

Role of Dynamics in the Autoinhibition and Activation of the Hyperpolarization-activated Cyclic Nucleotide-modulated (HCN) Ion Channels

The hyperpolarization-activated cyclic nucleotide-modulated (HCN) ion channels control rhythmicity in neurons and cardiomyocytes. Cyclic AMP allosterically modulates HCN through the cAMP-dependent formation of a tetrameric gating ring spanning the intracellular region (IR) of HCN, to which cAMP binds. Although the apo versus holo conformational changes of the cAMP-binding domain (CBD) have been previously mapped, only limited information is currently available on the HCN IR dynamics, which have been hypothesized to play a critical role in the cAMP-dependent gating of HCN. Here, using molecular dynamics simulations validated and complemented by experimental NMR and CD data, we comparatively analyze HCN IR dynamics in the four states of the thermodynamic cycle arising from the coupling between cAMP binding and tetramerization equilibria. This extensive set of molecular dynamics trajectories captures the active-to-inactive transition that had remained elusive for other CBDs, and it provides unprecedented insight on the role of IR dynamics in HCN autoinhibition and its release by cAMP. Specifically, the IR tetramerization domain becomes more flexible in the monomeric states, removing steric clashes that the apo-CDB structure would otherwise impose. Furthermore, the simulations reveal that the active/inactive structural transition for the apo-monomeric CBD occurs through a manifold of pathways that are more divergent than previously anticipated. Upon cAMP binding, these pathways become disallowed, pre-confining the CBD conformational ensemble to a tetramer-compatible state. This conformational confinement primes the IR for tetramerization and thus provides a model of how cAMP controls HCN channel gating.     

Experiment and theory highlight role of native state topology in SH3 folding.

We use a combination of experiments, computer simulations and simple model calculations to characterize, first, the folding transition state ensemble of the src SH3 domain, and second, the features of the protein that determine its folding mechanism. Kinetic analysis of mutations at 52 of the 57 residues in the src SH3 domain revealed that the transition state ensemble is even more polarized than suspected earlier: no single alanine substitution in the N-terminal 15 residues or the C-terminal 9 residues has more than a two-fold effect on the folding rate, while such substitutions at 15 sites in the central three-stranded beta-sheet cause significant decreases in the folding rate. Molecular dynamics (MD) unfolding simulations and ab initio folding simulations on the src SH3 domain exhibit a hierarchy of folding similar to that observed in the experiments. The similarity in folding mechanism of different SH3 domains and the similar hierarchy of structure formation observed in the experiments and the simulations can be largely accounted for by a simple native state topology-based model of protein folding energy landscapes.     

Structural, stability, dynamic and binding properties of the ALS-causing T46I mutant of the hVAPB MSP domain as revealed by NMR and MD simulations

T46I is the second mutation on the hVAPB MSP domain which was recently identified from non-Brazilian kindred to cause a familial amyotrophic lateral sclerosis (ALS). Here using CD, NMR and molecular dynamics (MD) simulations, we characterized the structure, stability, dynamics and binding capacity of the T46I-MSP domain. The results reveal: 1) unlike P56S which we previously showed to completely eliminate the native MSP structure, T46I leads to no significant disruption of the native secondary and tertiary structures, as evidenced from its far-UV CD spectrum, as well as Cα and Cβ NMR chemical shifts. 2) Nevertheless, T46I does result in a reduced thermodynamic stability and loss of the cooperative urea-unfolding transition. As such, the T46I-MSP domain is more prone to aggregation than WT at high protein concentrations and temperatures in vitro, which may become more severe in the crowded cellular environments. 3) T46I only causes a 3-fold affinity reduction to the Nir2 peptide, but a significant elimination of its binding to EphA4. 4) EphA4 and Nir2 peptide appear to have overlapped binding interfaces on the MSP domain, which strongly implies that two signaling networks may have a functional interplay in vivo. 5) As explored by both H/D exchange and MD simulations, the MSP domain is very dynamic, with most loop residues and many residues on secondary structures highly fluctuated or/and exposed to bulk solvent. Although T46I does not alter overall dynamics, it does trigger increased dynamics of several local regions of the MSP domain which are implicated in binding to EphA4 and Nir2 peptide. Our study provides the structural and dynamic understanding of the T46I-causing ALS; and strongly highlights the possibility that the interplay of two signaling networks mediated by the FFAT-containing proteins and Eph receptors may play a key role in ALS pathogenesis.     

Structural Characterization of N-WASP Domain V Using MD Simulations with NMR and SAXS Data

Because of their large conformational heterogeneity, structural characterization of intrinsically disordered proteins (IDPs) is very challenging using classical experimental methods alone. In this study, we use NMR and small-angle x-ray scattering (SAXS) data with multiple molecular dynamics (MD) simulations to describe the conformational ensemble of the fully disordered verprolin homology domain of the neural Aldrich syndrome protein involved in the regulation of actin polymerization. First, we studied several back-calculation software of SAXS scattering intensity and optimized the adjustable parameters to accurately calculate the SAXS intensity from an atomic structure. We also identified the most appropriate force fields for MD simulations of this IDP. Then, we analyzed four conformational ensembles of neural Aldrich syndrome protein verprolin homology domain, two generated with the program flexible-meccano with or without NMR-derived information as input and two others generated by MD simulations with two different force fields. These four conformational ensembles were compared to available NMR and SAXS data for validation. We found that MD simulations with the AMBER-03w force field and the TIP4P/2005s water model are able to correctly describe the conformational ensemble of this 67-residue IDP at both local and global level.     

Molecular dynamics simulations of the Bcl-2 protein to predict the structure of its unordered flexible loop domain

B-cell lymphoma (Bcl-2) protein is an anti-apoptotic member of the Bcl-2 family. It is functionally demarcated into four Bcl-2 homology (BH) domains: BH1, BH2, BH3, BH4, one flexible loop domain (FLD), a transmembrane domain (TM), and an X domain. Bcl-2’s BH domains have clearly been elucidated from a structural perspective, whereas the conformation of FLD has not yet been predicted, despite its important role in regulating apoptosis through its interactions with JNK-1, PKC, PP2A phosphatase, caspase 3, MAP kinase, ubiquitin, PS1, and FKBP38. Many important residues that regulate Bcl-2 anti-apoptotic activity are present in this domain, for example Asp34, Thr56, Thr69, Ser70, Thr74, and Ser87. The structural elucidation of the FLD would likely help in attempts to accurately predict the effect of mutating these residues on the overall structure of the protein and the interactions of other proteins in this domain. Therefore, we have generated an increased quality model of the Bcl-2 protein including the FLD through modeling. Further, molecular dynamics (MD) simulations were used for FLD optimization, to predict the flexibility, and to determine the stability of the folded FLD. In addition, essential dynamics (ED) was used to predict the collective motions and the essential subspace relevant to Bcl-2 protein function. The predicted average structure and ensemble of MD-simulated structures were submitted to the Protein Model Database (PMDB), and the Bcl-2 structures obtained exhibited enhanced quality. This study should help to elucidate the structural basis for Bcl-2 anti-apoptotic activity regulation through its binding to other proteins via the FLD.     

Conformational interconversion in compstatin probed with molecular dynamics simulations

Compstatin is a 13-residue cyclic peptide that has the potential to become a therapeutic agent against unregulated complement activation. In our effort to understand the structural and dynamic characteristics of compstatin that form the basis for rational and combinatorial optimization of structure and activity, we performed 1-ns molecular dynamics (MD) simulations. We used as input in the MD simulations the ensemble of 21 lowest energy NMR structures, the average minimized structure, and a global optimization structure. At the end of the MD simulations we identified five conformations, with populations ranging between 9% and 44%. These conformations are as follows: 1) coil with alphaR-alphaR beta-turn, as was the conformation of the initial ensemble of NMR structures; 2) beta-hairpin with epsilon-alphaR beta-turn; 3) beta-hairpin with alphaR-alphaR beta-turn; 4) beta-hairpin with alphaR-beta beta-turn; and 5) alpha-helical. Conformational switch was possible with small amplitude backbone motions of the order of 0.1-0.4 A and free energy barrier crossing of 2-11 kcal/mol. All of the 21 MD structures corresponding to the NMR ensemble possessed a beta-turn, with 14 structures retaining the alphaR-alphaR beta-turn type, but the average minimized structure and the global optimization structures were converted to alpha-helical conformations. Overall, the MD simulations have aided to gain insight into the conformational space sampled by compstatin and have provided a measure of conformational interconversion. The calculated conformers will be useful as structural and possibly dynamic templates for optimization in the design of compstatin using structure-activity relations (SAR) or dynamics-activity relations (DAR).     

Dimerization of the p53 oligomerization domain: identification of a folding nucleus by molecular dynamics simulations

Dimerization of the p53 oligomerization domain involves coupled folding and binding of monomers. To examine the dimerization, we have performed molecular dynamics (MD) simulations of dimer folding from the rate-limiting transition state ensemble (TSE). Among 799 putative transition state structures that were selected from a large ensemble of high-temperature unfolding trajectories, 129 were identified as members of the TSE via calculation of a 50% transmission coefficient from at least 20 room-temperature simulations. This study is the first to examine the refolding of a protein dimer using MD simulations in explicit water, revealing a folding nucleus for dimerization. Our atomistic simulations are consistent with experiment and offer insight that was previously unobtainable.     

Thermal stability of single‐domain antibodies estimated by molecular dynamics simulations

Single‐domain antibodies (sdAbs) function like regular antibodies, however, consist of only one domain. Because of their low molecular weight, sdAbs have advantages with respect to production and delivery to their targets and for applications such as antibody drugs and biosensors. Thus, sdAbs with high thermal stability are required. In this work, we chose seven sdAbs, which have a wide range of melting temperature (T_m) values and known structures. We applied molecular dynamics (MD) simulations to estimate their relative stability and compared them with the experimental data. High‐temperature MD simulations at 400 K and 500 K were executed with simulations at 300 K as a control. The fraction of native atomic contacts, Q, measured for the 400 K simulations showed a fairly good correlation with the T_m values. Interestingly, when the residues were classified by their hydrophobicity and size, the Q values of hydrophilic residues exhibited an even better correlation, suggesting that stabilization is correlated with favorable interactions of hydrophilic residues. Measuring the Q value on a per‐residue level enabled us to identify residues that contribute significantly to the instability and thus demonstrating how our analysis can be used in a mutant case study.     

Vibrational circular dichroism and IR spectral analysis as a test of theoretical conformational modeling for a cyclic hexapeptide

A model cyclohexapeptide, cyclo-(Phe-(D)Pro-Gly-Arg-Gly-Asp) was synthesized and its IR and VCD spectra were used as a test of density functional theory (DFT) level predictions of spectral intensities for a peptide with a nonrepeating but partially constricted conformation. Peptide structure and flexibility was estimated by molecular dynamics (MD) simulations and the spectra were simulated using full quantum mechanical (QM) approaches for the complete peptide and for simplified models with truncated side chains. After simulated annealing, the backbone conformation of the ring structure is relatively stable, consisting of a normal beta-turn and a tight loop (no H-bond) which does not vary over short trajectories. Only in quite long MD runs at high temperatures do other conformations appear. MD simulations were carried out for the cyclic peptide in water and in TFE, which match experimental solvents, as well as with and without protonation of the Asp carboxyl group. DFT spectral simulations were made using the annealed structure and were extended to include basis set variation, to determine an optimal computational approach, and solvent simulation with a polarized continuum model (PCM). Stepwise full DFT simulation of spectra was done for various sequences with the same backbone geometry but based on (1) solely Gly residues, (2) Ala substitution except Gly and Pro, and (3) complete sequences with side chains. Additionally, a selection of structures was used to compute IR and VCD spectra with the optimal method to determine structural variation effects. The side chains, especially the Asp-COOH and Arg-NH(2) transitions, had an impact on the computed amide frequencies, IR intensities and VCD pattern. Since experimentally these groups would have little chirality, due to conformational variation, they do not impact the observed VCD spectra. Correcting for frequency shifts, the Ala model for the cyclopeptide gives the clearest representation of the amide VCD. The experimental sign pattern for the amide I' band in D(2)O and also the sharper, more intense amide I VCD band in TFE was seen to some degree in one conformer with Type II' turns, but the data favor a mix of structures.     

Structure of hen egg-white lysozyme solvated in TFE/water: a molecular dynamics simulation study based on NMR data

Various experimental studies of hen egg white lysozyme (HEWL) in water and TFE/water clearly indicate structural differences between the native state and TFE state of HEWL, e.g. the helical content of the protein in the TFE state is much higher than in the native state. However, the available detailed NMR studies were not sufficient to determine fully a structure of HEWL in the TFE state. Different molecular dynamics (MD) simulations, i.e. at room temperature, at increased temperature and using proton–proton distance restraints derived from NMR NOE data, have been used to generate configurational ensembles corresponding to the TFE state of HEWL. The configurational ensemble obtained at room temperature using atom-atom distance restraints measured for HEWL in TFE/water solution satisfies the experimental data and has the lowest protein energy. In this ensemble residues 50–58, which are part of the β-sheet in native HEWL, adopt fluctuating α-helical secondary structure.     

Molecular dynamics simulations of peptides and proteins with a continuum electrostatic model based on screened Coulomb potentials

A continuum electrostatics approach for molecular dynamics (MD) simulations of macromolecules is presented and analyzed for its performance on a peptide and a globular protein. The approach incorporates the screened Coulomb potential (SCP) continuum model of electrostatics, which was reported earlier. The model was validated in a broad set of tests some of which were based on Monte Carlo simulations that included single amino acids, peptides, and proteins. The implementation for large-scale MD simulations presented in this article is based on a pairwise potential that makes the electrostatic model suitable for fast analytical calculation of forces. To assess the suitability of the approach, a preliminary validation is conducted, which consists of (i) a 3-ns MD simulation of the immunoglobulin-binding domain of streptococcal protein G, a 56-residue globular protein and (ii) a 3-ns simulation of Dynorphin, a biological peptide of 17 amino acids. In both cases, the results are compared with those obtained from MD simulations using explicit water (EW) molecules in an all-atom representation. The initial structure of Dynorphin was assumed to be an alpha-helix between residues 1 and 9 as suggested from NMR measurements in micelles. The results obtained in the MD simulations show that the helical structure collapses early in the simulation, a behavior observed in the EW simulation and consistent with spectroscopic data that suggest that the peptide may adopt mainly an extended conformation in water. The dynamics of protein G calculated with the SCP implicit solvent model (SCP-ISM) reveals a stable structure that conserves all the elements of secondary structure throughout the entire simulation time. The average structures calculated from the trajectories with the implicit and explicit solvent models had a cRMSD of 1.1 A, whereas each average structure had a cRMSD of about 0.8A with respect to the X-ray structure. The main conformational differences of the average structures with respect to the crystal structure occur in the loop involving residues 8-14. Despite the overall similarity of the simulated dynamics with EW and SCP models, fluctuations of side-chains are larger when the implicit solvent is used, especially in solvent exposed side-chains. The MD simulation of Dynorphin was extended to 40 ns to study its behavior in an aqueous environment. This long simulation showed that the peptide has a tendency to form an alpha-helical structure in water, but the stabilization free energy is too weak, resulting in frequent interconversions between random and helical conformations during the simulation time. The results reported here suggest that the SCP implicit solvent model is adequate to describe electrostatic effects in MD simulation of both peptides and proteins using the same set of parameters. It is suggested that the present approach could form the basis for the development of a reliable and general continuum approach for use in molecular biology, and directions are outlined for attaining this long-term goal.     

Bending of the calmodulin central helix: a theoretical study.

The crystal structure of calcium-calmodulin (CaM) reveals a protein with a typical dumbbell structure. Various spectroscopic studies have suggested that the central linker region of CaM, which is alpha-helical in the crystal structure, is flexible in solution. In particular, NMR studies have indicated the presence of a flexible backbone between residues Lys 77 and Asp 80. This flexibility is related directly to the function of the protein because it enables the N- and C-terminal domains of the protein to move toward each other and bind to the CaM-binding domain of a target protein. We have investigated the flexibility of the CaM central helix by a variety of computational techniques: molecular dynamics (MD) simulations, normal mode analysis (NMA), and essential dynamics (ED) analysis. Our MD results reproduce the experimentally determined location of the bend in a simulation of only the CaM central helix, indicating that the bending point is an intrinsic property of the alpha-helix, for which the remainder of the protein is not important. Interestingly, the modes found by the ED analysis of the MD trajectory are very similar to the lowest frequency modes from the NM analysis and to modes found by an ED analysis of different structures in a set of NMR structures. Electrostatic interactions involving residues Arg 74 and Asp 80 seem to be important for these bending motions and unfolding, which is in line with pH-dependent NMR and CD studies.     

Molecular Dynamics Ensemble Refinement of Intrinsically Disordered Peptides According to Deconvoluted Spectra from Circular Dichroism

We have developed a computational method of atomistically refining the structural ensemble of intrinsically disordered peptides (IDPs) facilitated by experimental measurements using circular dichroism spectroscopy (CD). A major challenge surrounding this approach stems from the deconvolution of experimental CD spectra into secondary structure features of the IDP ensemble. Currently available algorithms for CD deconvolution were designed to analyze the spectra of proteins with stable secondary structures. Herein, our work aims to minimize any bias from the peptide deconvolution analysis by implementing a non-negative linear least-squares fitting algorithm in conjunction with a CD reference data set that contains soluble and denatured proteins (SDP48). The non-negative linear least-squares method yields the best results for deconvolution of proteins with higher disordered content than currently available methods, according to a validation analysis of a set of protein spectra with Protein Data Bank entries. We subsequently used this analysis to deconvolute our experimental CD data to refine our computational model of the peptide secondary structure ensemble produced by all-atom molecular dynamics simulations with implicit solvent. We applied this approach to determine the ensemble structures of a set of short IDPs, that mimic the calmodulin binding domain of calcium/calmodulin-dependent protein kinase II and its 1-amino-acid and 3-amino-acid mutants. Our study offers a, to our knowledge, novel way to solve the ensemble secondary structures of IDPs in solution, which is important to advance the understanding of their roles in regulating signaling pathways through the formation of complexes with multiple partners.     

All-atom structure ensembles of islet amyloid polypeptides determined by enhanced sampling and experiment data restraints

Exploring the accurate structure ensembles are critical to understand the functions of intrinsically disordered proteins (IDPs). As a well-known IDP, islet amyloid polypeptide (IAPP) plays important roles in the development of human type II diabetes (T2D). The toxicity of human IAPP (hIAPP) is induced by the amyloidosis of the peptide, however, its aggregation mechanism remains ambiguous. The characterization of structure ensemble of hIAPP, as well as the differences between hIAPP and its non-amyloidogenic homologous such as rat IAPP (rIAPP), would greatly help to illuminate the amyloidosis mechanism of IAPP. In this study, the atomic structure ensembles of hIAPP and rIAPP were characterized by all-atom molecular dynamics (MD) simulations combined with enhanced sampling technology and experiment data restraints. The obtained structure ensembles were firstly compared with those determined by the conventional MD (cMD) and enhanced sampling without experiment data restraints. The results showed that the enhanced sampling and experiment data restraints would improve the simulation accuracy. The transient N-terminal α-helix structures were adopted by the sub-states of both hIAPP and rIAPP, however, the C-terminal helical structures were only present on hIAPP. The hydrophobic residues in the amyloid-core region of hIAPP are exposed to the solvent. The structure ensemble differences between hIAPP and rIAPP revealed in this work provide potential explain to the amyloidogenic mechanism and would be helpful for the design of drugs to combat T2D.     

Ensemble MD simulations restrained via crystallographic data: Accurate structure leads to accurate dynamics

Currently, the best existing molecular dynamics (MD) force fields cannot accurately reproduce the global free‐energy minimum which realizes the experimental protein structure. As a result, long MD trajectories tend to drift away from the starting coordinates (e.g., crystallographic structures). To address this problem, we have devised a new simulation strategy aimed at protein crystals. An MD simulation of protein crystal is essentially an ensemble simulation involving multiple protein molecules in a crystal unit cell (or a block of unit cells). To ensure that average protein coordinates remain correct during the simulation, we introduced crystallography‐based restraints into the MD protocol. Because these restraints are aimed at the ensemble‐average structure, they have only minimal impact on conformational dynamics of the individual protein molecules. So long as the average structure remains reasonable, the proteins move in a native‐like fashion as dictated by the original force field. To validate this approach, we have used the data from solid‐state NMR spectroscopy, which is the orthogonal experimental technique uniquely sensitive to protein local dynamics. The new method has been tested on the well‐established model protein, ubiquitin. The ensemble‐restrained MD simulations produced lower crystallographic R factors than conventional simulations; they also led to more accurate predictions for crystallographic temperature factors, solid‐state chemical shifts, and backbone order parameters. The predictions for ¹⁵N relaxation rates are at least as accurate as those obtained from conventional simulations. Taken together, these results suggest that the presented trajectories may be among the most realistic protein MD simulations ever reported. In this context, the ensemble restraints based on high‐resolution crystallographic data can be viewed as protein‐specific empirical corrections to the standard force fields.