Sucrose in aqueous solution revisited, Part 2: adaptively biased molecular dynamics simulations and computational analysis of NMR relaxation

We report 100 ns molecular dynamics simulations, at various temperatures, of sucrose in water (with concentrations of sucrose ranging from 0.02 to 4M), and in a 7:3 water‐DMSO mixture. Convergence of the resulting conformational ensembles was checked using adaptive‐biased simulations along the glycosidic Φ and ψ torsion angles. NMR relaxation parameters, including longitudinal (R₁) and transverse (R₂) relaxation rates, nuclear Overhauser enhancements (NOE), and generalized order parameter (S²) were computed from the resulting time‐correlation functions. The amplitude and time scales of molecular motions change with temperature and concentration in ways that track closely with experimental results, and are consistent with a model in which sucrose conformational fluctuations are limited (with 80–90% of the conformations having ??ψ values within 20° of an average conformation), but with some important differences in conformation between pure water and DMSO‐water mixtures. © 2011 Wiley Periodicals, Inc. Biopolymers 97: 289–302, 2012.     

How does it really move? Recent progress in the investigation of protein nanosecond dynamics by NMR and simulation

Nuclear magnetic resonance (NMR) spin relaxation experiments currently probe molecular motions on timescales from picoseconds to nanoseconds. The detailed interpretation of these motions in atomic detail benefits from complementarity with the results from molecular dynamics (MD) simulations. In this mini-review, we describe the recent developments in experimental techniques to study the backbone dynamics from ¹⁵N relaxation and side-chain dynamics from ¹³C relaxation, discuss the different analysis approaches from model-free to dynamics detectors, and highlight the many ways that NMR relaxation experiments and MD simulations can be used together to improve the interpretation and gain insights into protein dynamics.     

FF12MC: A revised AMBER forcefield and new protein simulation protocol

Specialized to simulate proteins in molecular dynamics (MD) simulations with explicit solvation, FF12MC is a combination of a new protein simulation protocol employing uniformly reduced atomic masses by tenfold and a revised AMBER forcefield FF99 with (i) shortened C? H bonds, (ii) removal of torsions involving a nonperipheral sp³ atom, and (iii) reduced 1–4 interaction scaling factors of torsions ? and ψ. This article reports that in multiple, distinct, independent, unrestricted, unbiased, isobaric–isothermal, and classical MD simulations FF12MC can (i) simulate the experimentally observed flipping between left‐ and right‐handed configurations for C14–C38 of BPTI in solution, (ii) autonomously fold chignolin, CLN025, and Trp‐cage with folding times that agree with the experimental values, (iii) simulate subsequent unfolding and refolding of these miniproteins, and (iv) achieve a robust Z score of 1.33 for refining protein models TMR01, TMR04, and TMR07. By comparison, the latest general‐purpose AMBER forcefield FF14SB locks the C14–C38 bond to the right‐handed configuration in solution under the same protein simulation conditions. Statistical survival analysis shows that FF12MC folds chignolin and CLN025 in isobaric–isothermal MD simulations 2–4 times faster than FF14SB under the same protein simulation conditions. These results suggest that FF12MC may be used for protein simulations to study kinetics and thermodynamics of miniprotein folding as well as protein structure and dynamics. Proteins 2016; 84:1490–1516. © 2016 The Authors Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.     

Motional timescale predictions by molecular dynamics simulations: Case study using proline and hydroxyproline sidechain dynamics

We propose a new approach for force field optimizations which aims at reproducing dynamics characteristics using biomolecular MD simulations, in addition to improved prediction of motionally averaged structural properties available from experiment. As the source of experimental data for dynamics fittings, we use ¹³C NMR spin‐lattice relaxation times T₁ of backbone and sidechain carbons, which allow to determine correlation times of both overall molecular and intramolecular motions. For structural fittings, we use motionally averaged experimental values of NMR J couplings. The proline residue and its derivative 4‐hydroxyproline with relatively simple cyclic structure and sidechain dynamics were chosen for the assessment of the new approach in this work. Initially, grid search and simplexed MD simulations identified large number of parameter sets which fit equally well experimental J couplings. Using the Arrhenius‐type relationship between the force constant and the correlation time, the available MD data for a series of parameter sets were analyzed to predict the value of the force constant that best reproduces experimental timescale of the sidechain dynamics. Verification of the new force‐field (termed as AMBER99SB‐ILDNP) against NMR J couplings and correlation times showed consistent and significant improvements compared to the original force field in reproducing both structural and dynamics properties. The results suggest that matching experimental timescales of motions together with motionally averaged characteristics is the valid approach for force field parameter optimization. Such a comprehensive approach is not restricted to cyclic residues and can be extended to other amino acid residues, as well as to the backbone. Proteins 2014; 82:195–215. © 2013 Wiley Periodicals, Inc.     

Dual effects of familial Alzheimer's disease mutations (D7H,D7N,and H6R) on amyloid β peptide: Correlation dynamics and zinc binding

Although the N‐terminal region of Amyloid β (Aβ) peptides plays dual roles as metal‐coordinating sites and conformational modulator, few studies have been performed to explore the effects of mutations at this region on the overall conformational ensemble of Aβ and the binding propensity of metal ions. In this work, we focus on how three familial Alzheimer's disease mutations (D7H, D7N, and H6R) alter the structural characteristics and thermodynamic stabilities of Aβ42 using molecular dynamics simulations. We observe that each mutation displays increased β‐sheet structures in both N and C termini. In particular, both the N terminus and central hydrophobic region of D7H can form stable β‐hairpin structures with its C terminus. The conserved turn structure at Val²⁴–Lys²⁸ in all peptides and Zn²⁺‐bound Aβ42 is confirmed as the common structural motif to nucleate folding of Aβ. Each mutant can significantly increase the solvation free energy and thus enhance the aggregation of Aβ monomers. The correlation dynamics between Aβ(1–16) and Aβ(17–42) fragments are elucidated by linking the domain motions with the corresponding structured conformations. We characterize the different populations of correlated domain motions for each mutant from a more macroscopic perspective, and unexpectedly find that Zn²⁺‐bound Aβ42 ensemble shares the same populations as Aβ42, indicating that the binding of Zn²⁺ to Aβ follows the conformational selection mechanism, and thus is independent of domain motions, even though the structures of Aβ have been modified at a residue level. Proteins 2014; 82:3286–3297. © 2014 Wiley Periodicals, Inc.     

Constant pH molecular dynamics of proteins in explicit solvent with proton tautomerism

pH is a ubiquitous regulator of biological activity, including protein‐folding, protein‐protein interactions, and enzymatic activity. Existing constant pH molecular dynamics (CPHMD) models that were developed to address questions related to the pH‐dependent properties of proteins are largely based on implicit solvent models. However, implicit solvent models are known to underestimate the desolvation energy of buried charged residues, increasing the error associated with predictions that involve internal ionizable residue that are important in processes like hydrogen transport and electron transfer. Furthermore, discrete water and ions cannot be modeled in implicit solvent, which are important in systems like membrane proteins and ion channels. We report on an explicit solvent constant pH molecular dynamics framework based on multi‐site λ‐dynamics (CPHMD^MSλD). In the CPHMD^MSλD framework, we performed seamless alchemical transitions between protonation and tautomeric states using multi‐site λ‐dynamics, and designed novel biasing potentials to ensure that the physical end‐states are predominantly sampled. We show that explicit solvent CPHMD^MSλD simulations model realistic pH‐dependent properties of proteins such as the Hen‐Egg White Lysozyme (HEWL), binding domain of 2‐oxoglutarate dehydrogenase (BBL) and N‐terminal domain of ribosomal protein L9 (NTL9), and the pK_a predictions are in excellent agreement with experimental values, with a RMSE ranging from 0.72 to 0.84 pK_a units. With the recent development of the explicit solvent CPHMD^MSλD framework for nucleic acids, accurate modeling of pH‐dependent properties of both major class of biomolecules—proteins and nucleic acids is now possible. Proteins 2014; 82:1319–1331. © 2013 Wiley Periodicals, Inc.     

Bacterial collagen‐binding domain targets undertwisted regions of collagen

Clostridium histolyticum collagenase causes extensive degradation of collagen in connective tissue that results in gas gangrene. The C‐terminal collagen‐binding domain (CBD) of these enzymes is the minimal segment required to bind to a collagen fibril. CBD binds unidirectionally to the undertwisted C‐terminus of triple helical collagen. Here, we examine whether CBD could also target undertwisted regions even in the middle of the triple helix. Collageneous peptides with an additional undertwisted region were synthesized by introducing a Gly → Ala substitution [(POG)_xPOA(POG)_y]₃, where x + y = 9 and x > 3). ¹H–¹⁵N heteronuclear single quantum coherence nuclear magnetic resonance (HSQC NMR) titration studies with ¹⁵N‐labeled CBD demonstrated that the minicollagen binds to a 10 Å wide 25 Å long cleft. Six collagenous peptides each labeled with a nitroxide radical were then titrated with ¹⁵N‐labeled CBD. CBD binds to either the Gly → Ala substitution site or to the C‐terminus of each minicollagen. Small‐angle X‐ray scattering measurements revealed that CBD prefers to bind the Gly → Ala site to the C‐terminus. The HSQC NMR spectra of ¹⁵N‐labeled minicollagen and minicollagen with undertwisted regions were unaffected by the titration of unlabeled CBD. The results imply that CBD binds to the undertwisted region of the minicollagen but does not actively unwind the triple helix.     

Sop‐GPU: Accelerating biomolecular simulations in the centisecond timescale using graphics processors

Theoretical exploration of fundamental biological processes involving the forced unraveling of multimeric proteins, the sliding motion in protein fibers and the mechanical deformation of biomolecular assemblies under physiological force loads is challenging even for distributed computing systems. Using a C_α‐based coarse‐grained self organized polymer (SOP) model, we implemented the Langevin simulations of proteins on graphics processing units (SOP‐GPU program). We assessed the computational performance of an end‐to‐end application of the program, where all the steps of the algorithm are running on a GPU, by profiling the simulation time and memory usage for a number of test systems. The ～90‐fold computational speedup on a GPU, compared with an optimized central processing unit program, enabled us to follow the dynamics in the centisecond timescale, and to obtain the force‐extension profiles using experimental pulling speeds (v_f = 1–10 μm/s) employed in atomic force microscopy and in optical tweezers‐based dynamic force spectroscopy. We found that the mechanical molecular response critically depends on the conditions of force application and that the kinetics and pathways for unfolding change drastically even upon a modest 10‐fold increase in v_f. This implies that, to resolve accurately the free energy landscape and to relate the results of single‐molecule experiments in vitro and in silico, molecular simulations should be carried out under the experimentally relevant force loads. This can be accomplished in reasonable wall‐clock time for biomolecules of size as large as 10⁵ residues using the SOP‐GPU package. Proteins 2010; © 2010 Wiley‐Liss, Inc.     

Mode‐coupling Smoluchowski dynamics of a double‐stranded DNA oligomer

The local dynamics of a double‐stranded DNA d(TpCpGpCpG)₂ is obtained to second order in the mode‐coupling expansion of the Smoluchowski diffusion theory. The time correlation functions of bond variables are derived and the ¹³C‐nmr spin–lattice relaxation times T₁ of different ¹³C along the chains are calculated and compared to experimental data from the literature at three frequencies. The DNA is considered as a fluctuating three‐dimensional structure undergoing rotational diffusion. The fluctuations are evaluated using molecular dynamics simulations, with the ensemble averages approximated by time averages along a trajectory of length 1 ns. Any technique for sampling the configurational space can be used as an alternative. For a fluctuating three‐dimensional (3D) structure using the three first‐order vector modes of lower rates, higher order basis sets of second‐rank tensor are built to give the required mode coupling dynamics. Second‐ and even first‐order theories are found to be in close agreement with the experimental results, especially at high frequency, where the differences in T₁ for ¹³C in the base pairs, sugar, and backbone are well described. These atomistic calculations are of general application for studying, on a molecular basis, the local dynamics of fluctuating 3D structures such as double‐helix DNA fragments, proteins, and protein–DNA complexes. © 1999 John Wiley & Sons, Inc. Biopoly 50: 613–629, 1999     

Spectral densities of nitrogen nuclei in Escherichia coli ribonuclease HI obtained by 15N NMR relaxation and molecular dynamics

Summary Spectral densities of the ¹⁵N amide in Escherichia coli ribonuclease HI, obtained from NMR relaxation experiments, were compared with those calculated using a molecular dynamics (MD) simulation. All calculations and comparisons assumed that the auto-correlation function describing the internal motions of the molecule was independent of the auto-correlation function associated with overall rotational diffusion. Comparisons were limited to those residues for which the auto-correlation function of internal motions rapidly relaxed and reached a steady state within 205 ps. The results show the importance of frequency components as well as amplitudes of internal motions in order to obtain a meaningful comparison of MD simulations with NMR data.     

Global fold of human cannabinoid type 2 receptor probed by solid‐state 13C‐, 15N‐MAS NMR and molecular dynamics simulations

The global fold of human cannabinoid type 2 (CB₂) receptor in the agonist‐bound active state in lipid bilayers was investigated by solid‐state ¹³C‐ and ¹⁵N magic‐angle spinning (MAS) NMR, in combination with chemical‐shift prediction from a structural model of the receptor obtained by microsecond‐long molecular dynamics (MD) simulations. Uniformly ¹³C‐ and ¹⁵N‐labeled CB₂ receptor was expressed in milligram quantities by bacterial fermentation, purified, and functionally reconstituted into liposomes. ¹³C MAS NMR spectra were recorded without sensitivity enhancement for direct comparison of C_α, C_β, and C?O bands of superimposed resonances with predictions from protein structures generated by MD. The experimental NMR spectra matched the calculated spectra reasonably well indicating agreement of the global fold of the protein between experiment and simulations. In particular, the ¹³C chemical shift distribution of C_α resonances was shown to be very sensitive to both the primary amino acid sequence and the secondary structure of CB₂. Thus the shape of the C_α band can be used as an indicator of CB₂ global fold. The prediction from MD simulations indicated that upon receptor activation a rather limited number of amino acid residues, mainly located in the extracellular Loop 2 and the second half of intracellular Loop 3, change their chemical shifts significantly (≥1.5 ppm for carbons and ≥5.0 ppm for nitrogens). Simulated two‐dimensional ¹³C_α(i)? ¹³C?O(i) and ¹³C?O(i)? ¹⁵NH(i + 1) dipolar‐interaction correlation spectra provide guidance for selective amino acid labeling and signal assignment schemes to study the molecular mechanism of activation of CB₂ by solid‐state MAS NMR. Proteins 2014; 82:452–465. © 2013 Wiley Periodicals, Inc.     

Conformational analysis of potent sweet taste ligands by NMR and computer simulations

We report the conformational analysis by ¹H‐nmr and computer simulations of five potent sweet molecules, N‐(3,3‐dimethylbutyl)‐L ‐aspartyl‐S‐(α‐methyl)phenylalanine methylester (1; 5000 times more potent than sucrose), L ‐aspartyl‐D ‐valine (S)‐α‐methoxycarbonylmethylbenzylamide (2; 1400 times more potent than sucrose), L ‐aspartyl‐D ‐valine α‐phenylcyclopentylamide (3; 1200 times more potent than sucrose), L ‐aspartyl‐D ‐α‐aminobutyric acid (S)‐α‐cyclohexylpropylamide (4; 2300 times more potent than sucrose), and L ‐aspartyl‐D ‐valine (R)‐α‐methylthiomethylbenzylamide (5; 3000 times more potent than sucrose). The “L‐shaped” structure, which we believe to be responsible for sweet taste, is accessible to all five sweet compounds in solution. This structure is characterized by a zwitterionic ring formed by the A‐H and B containing moieties located in the +y axis and by the hydrophobic group X pointing into the +x axis. Other accessible conformations of these flexible molecules are extended conformations with the A‐H and B containing moieties in the +y axis and the hydrophobic group X pointing in the −y axis and reversed L‐shaped structures with the hydrophobic group X projecting along the −x axis. The remarkable potency of the N‐alkylated compound 1 supports our recent hypothesis that a second hydrophobic binding domain in addition to interactions arising from the L‐shaped structure leads to an enhancement of sweetness potency. © 1999 John Wiley & Sons, Inc. Biopoly 49: 525–539, 1999     

Domain motions of hyaluronan lyase underlying processive hyaluronan translocation

Hyaluronan lyase (Hyal) is a surface enzyme occurring in many bacterial organisms including members of Streptococcus species. Streptococcal Hyal primarily degrades hyaluronan‐substrate (HA) of the extracellular matrix. This degradation appears to facilitate the spread of this bacterium throughout host tissues. Unlike purely endolytic degradation of its other substrates, unsulfated chondroitin or some chondroitin sulfates, the degradation of HA by Hyal proceeds by processive exolytic cleavage of one disaccharide at a time following an initial endolytic cut. Molecular dynamics (MD) studies of Hyal from Streptococcus pneumoniae are presented that address the enzyme's molecular mechanism of action and the role of domain motions for processive functionality. The analysis of extensive sub‐microsecond MD simulations of this enzyme action on HA‐substrates of different lengths and the connection between the domain dynamics of Hyal and the translocation of the HA‐substrate reveals that opening/closing and twisting domain motions of the Hyal are intimately linked to processive HA degradation. Enforced simulations confirmed this finding as the domain motions in SpnHyal were found to be induced by enforced substrate translocation. These results establish the dynamic interplay between Hyal flexibility and substrate translocation and provide insight into the processive mechanism of Hyal. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Tunable white light emission from single‐phased Li2SrSiO4:Dy3+ phosphors by co‐doping with Eu3+

A series of single‐phase full‐color emitting Li₂Sr_1−x−ySiO₄:xDy³⁺,yEu³⁺ phosphors were synthesized by solid‐state reaction and characterized by X‐ray diffraction and photoluminescence analyses. The samples showed emission peaks at 488 nm (blue), 572 nm (yellow), 592 nm (orange) and 617 nm (red) under 393 nm excitation. The photoluminescence excitation spectra, comprising the Eu–O charge transfer band and 4f–4f transition bands of Dy³⁺ and Eu³⁺, range from 200 to 500 nm. The Commission Internationale de I'Eclairage chromaticity coordinates for Li₂Sr_0.98−xSiO₄:0.02Dy³⁺,xEu³⁺ phosphors were simulated. By manipulating Eu³⁺ and Dy³⁺ concentrations, the color points of Li₂Sr_1−x−ySiO₄:xDy³⁺,yEu³⁺ were tuned from the greenish‐white region to white light and eventually to reddish‐white region, demonstrating that a tunable white light can be obtained by Li₂Sr_1−x−ySiO₄:xDy³⁺,yEu³⁺ phosphors. Li₂Sr_0.98−xSiO₄:0.02Dy³⁺, xEu³⁺ can serve as a white‐light‐emitting phosphor for phosphor‐converted light‐emitting diode. Copyright © 2014 John Wiley & Sons, Ltd.     

Molecular dynamics simulations reveal that apo‐HisJ can sample a closed conformation

The Escherichia coli histidine binding protein HisJ is a type II periplasmic binding protein (PBP) that preferentially binds histidine and interacts with its cytoplasmic membrane ABC transporter, HisQMP₂, to initiate histidine transport. HisJ is a bilobal protein where the larger Domain 1 is connected to the smaller Domain 2 via two linking strands. Type II PBPs are thought to undergo “Venus flytrap” movements where the protein is able to reversibly transition from an open to a closed conformation. To explore the accessibility of the closed conformation to the apo state of the protein, we performed a set of all‐atom molecular dynamics simulations of HisJ starting from four different conformations: apo‐open, apo‐closed, apo‐semiopen, and holo‐closed. The simulations reveal that the closed conformation is less dynamic than the open one. HisJ experienced closing motions and explored semiopen conformations that reverted to closed forms resembling those found in the holo‐closed state. Essential dynamics analysis of the simulations identified domain closing/opening and twisting as main motions. The formation of specific inter‐hinge strand and interdomain polar interactions contributed to the adoption of the closed apo‐conformations although they are up to 2.5‐fold less prevalent compared with the holo‐closed simulations. The overall sampling of the closed form by apo‐HisJ provides a rationale for the binding of unliganded PBPs with their cytoplasmic membrane ABC transporters. Proteins 2014; 82:386–398. © 2013 Wiley Periodicals, Inc.     

NMR investigations of structural and dynamics features of natively unstructured drug peptide – salmon calcitonin: implication to rational design of potent sCT analogs

Backbone dynamics and conformational properties of drug peptide salmon calcitonin have been studied in aqueous solution using nuclear magnetic resonance (NMR). Although salmon calcitonin (sCT) is largely unfolded in solution (as has been reported in several circular dichroism studies), the secondary H^α chemical shifts and three bond H^N–H^α coupling constants indicated that most of the residues of the peptide are populating the α‐helical region of the Ramachandran (?, ψ) map. Further, the peptide in solution has been found to exhibit multiple conformational states exchanging slowly on the NMR timescale (10²–10³ s^?1), inferred by the multiple chemical shift assignments in the region Leu4–Leu12 and around Pro23 (for residues Gln20–Tyr22 and Arg24). Possibly, these slowly exchanging multiple conformational states might inhibit symmetric self‐association of the peptide and, in part, may account for its reduced aggregation propensity compared with human calcitonin (which lacks this property). The ¹⁵N NMR‐relaxation data revealed (i) the presence of slow (microsecond‐to‐millisecond) timescale dynamics in the N‐terminal region (Cys1–Ser5) and core residues His17 and Asn26 and (ii) the presence of high frequency (nanosecond‐to‐picosecond) motions in the C‐terminal arm. Put together, the various results suggested that (i) the flexible C‐terminal of sCT (from Thr25–Thr31) is involved in identification of specific target receptors, (ii) whereas the N‐terminal of sCT (from Cys1–Gln20) in solution – exhibiting significant amount of conformational plasticity and strong bias towards biologically active α‐helical structure – facilitates favorable conformational adaptations while interacting with the intermembrane domains of these target receptors. Thus, we believe that the structural and dynamics features of sCT presented here will be useful guiding attributes for the rational design of biologically active sCT analogs. Copyright © 2012 European Peptide Society and John Wiley & Sons, Ltd.     

The recognition of membrane‐bound PtdIns3P by PX domains

Phox‐homology (PX) domains target proteins to the organelles of the secretary and endocytic systems by binding to phosphatidylinositol phospholipids (PIPs). Among all the structures of PX domains that have been solved, only three have been solved in a complex with the main physiological ligand: PtdIns3P. In this work, molecular dynamic simulations have been used to explore the structure and dynamics of the p40^phox–PX domain and the SNX17–PX domain and their interaction with membrane‐bound PtdIns3P. In the simulations, both PX domains associated spontaneously with the membrane‐bound PtdIns3P and formed stable complexes. The interaction between the p40^phox–PX domain and PtdIns3P in the membrane was found to be similar to the crystal structure of the p40^phox–PX–PtdIns3P complex that is available. The interaction between the SNX17–PX domain and PtdIns3P was similar to that observed in the p40^phox–PX–PtdIns3P complex; however, some residues adopted different orientations. The simulations also showed that nonspecific interactions between the β1–β2 loop and the membrane play an important role in the interaction of membrane bound PtdIns3P and different PX domains. The behaviour of unbound PtdIns3P within a 2‐oleoyl‐1‐palmitoyl‐sn‐glycero‐3‐phosphocholine (POPC) membrane environment was also examined and compared to the available experimental data and simulation studies of related molecules. Proteins 2014; 82:2332–2342. © 2014 Wiley Periodicals, Inc.     

A First‐Principles and Experimental Investigation of Nickel Solubility into the P2 NaxCoO2 Sodium‐Ion Cathode

The difficulty in finding positive electrode materials for sodium‐ion (Na‐ion) batteries with a large specific energy has slowed down their commercialization. Layered transition metal (M) oxides Na_xMO₂ with a two‐layer oxygen stacking (P2, 0.6 ≤ x ≤ 0.75), are promising candidates. However, the high average metal oxidation state needed during synthesis means that P2 Na_xMO₂ cathodes often require the introduction of high‐valent cations (Mn⁴⁺, Ti⁴⁺, Sn⁵⁺, or Te⁶⁺), limiting the cathode's performance. Using a combination of first‐principles calculations and experiments, the feasibility of P2 cathodes containing only electrochemically active nickel and cobalt cations is investigated. It is found that P2 Na_xNi_yCo_1–yO₂ materials with x = 0.66, 0.75, and 0 ≤ y ≤ 0.33 are either thermodynamically stable or metastable yet close to the convex hull at typical P2 synthesis temperatures (≈1000 K). It is demonstrated that a novel P2 compound with y = 0.22 and both Ni^3+/4+ and Co^3+/4+ can be successfully synthesized. It is studied electrochemically and structurally, using in situ and ex situ X‐ray diffraction. It is demonstrated that the chemical space of P2 layered compounds is not fully explored yet and that ab initio phase diagrams allow the determination of new high‐specific energy positive electrodes to be targeted experimentally.     

All‐atom and coarse‐grained simulations of the forced unfolding pathways of the SNARE complex

The SNARE complex, consisting of three proteins (VAMP2, syntaxin, and SNAP‐25), is thought to drive membrane fusion by assembling into a four‐helix bundle through a zippering process. In support of the above zippering model, a recent single‐molecule optical tweezers experiment by Gao et al. revealed a sequential unzipping of SNARE along VAMP2 in the order of the linker domain → the C‐terminal domain → the N‐terminal domain. To offer detailed structural insights to this unzipping process, we have performed all‐atom and coarse‐grained steered molecular dynamics (sMD) simulations of the forced unfolding pathways of SNARE using different models and force fields. Our findings are summarized as follows: First, the sMD simulations based on either an all‐atom force field (with an implicit solvent model) or a coarse‐grained Go model were unable to capture the forced unfolding pathway of SNARE as observed by Gao et al., which may be attributed to insufficient simulation time and inaccurate force fields. Second, the sMD simulations based on a reparameterized coarse‐grained model (i.e., modified elastic network model) were able to predict a sequential unzipping of SNARE in good agreement with the findings by Gao et al. The key to this success is to reparameterize the intrahelix and interhelix nonbonded force constants against the pair‐wise residue–residue distance fluctuations collected from all‐atom MD simulations of SNARE. Therefore, our finding supports the importance of accurately describing the inherent dynamics/flexibility of SNARE (in the absence of force), in order to correctly simulate its unfolding behaviors under force. This study has established a useful computational framework for future studies of the zippering function of SNARE and its perturbations by point mutations with amino‐acid level of details, and more generally the forced unfolding pathways of other helix bundle proteins. Proteins 2014; 82:1376–1386. © 2014 Wiley Periodicals, Inc.     

Nuclear magnetic resonance‐based determination of dioxygen binding sites in protein cavities

The location and ligand accessibility of internal cavities in cysteine‐free wild‐type T4 lysozyme was investigated using O₂ gas‐pressure NMR spectroscopy and molecular dynamics (MD) simulation. Upon increasing the concentration of dissolved O₂ in solvent to 8.9 mM, O₂‐induced paramagnetic relaxation enhancements (PREs) to the backbone amide and side chain methyl protons were observed, specifically around two cavities in the C‐terminal domain. To determine the number of O₂ binding sites and their atomic coordinates from the 1/r⁶ distance dependence of the PREs, we established an analytical procedure using Akaike's Information Criterion, in combination with a grid‐search. Two O₂‐accessible sites were identified in internal cavities: One site was consistent with the xenon‐binding site in the protein in crystal, and the other site was established to be a novel ligand‐binding site. MD simulations performed at 10 and 100 mM O₂ revealed dioxygen ingress and egress as well as rotational and translational motions of O₂ in the cavities. It is therefore suggested that conformational fluctuations within the ground‐state ensemble transiently develop channels for O₂ association with the internal protein cavities.