Local and long range potentials for heparin-protein systems for coarse-grained simulations

Heparin belongs to glycosaminoglycans (GAGs), a class of periodic linear anionic polysaccharides, which are functionally important components of the extracellular matrix owing to their interactions with various protein targets. Heparin is known to be involved in many cell signaling processes, while the experimental data available for heparin are significantly more abundant than for other GAGs. At the same time, the length and conformational flexibility of the heparin represent major challenges for its theoretical analysis. Coarse-grained (CG) approaches, which enable us to extend the size- and time-scale by orders of magnitude owing to reduction of system representation, appear, therefore, to be useful in simulating these systems. In this work, by using umbrella-sampling molecular dynamics simulations, we derived and parameterized the CG backbone-local potentials of heparin chains and the orientational potentials for the interactions of heparin with amino acid side chains to be further included in the physics-based Unified Coarse-Grained Model of biological macromolecules. With these potentials, simulations of extracellular matrix processes where both heparin and multiple proteins participate will be possible.     

Assessment of chemical-crosslink-assisted protein structure modeling in CASP13

With the advance of experimental procedures obtaining chemical crosslinking information is becoming a fast and routine practice. Information on crosslinks can greatly enhance the accuracy of protein structure modeling. Here, we review the current state of the art in modeling protein structures with the assistance of experimentally determined chemical crosslinks within the framework of the 13th meeting of Critical Assessment of Structure Prediction approaches. This largest-to-date blind assessment reveals benefits of using data assistance in difficult to model protein structure prediction cases. However, in a broader context, it also suggests that with the unprecedented advance in accuracy to predict contacts in recent years, experimental crosslinks will be useful only if their specificity and accuracy further improved and they are better integrated into computational workflows.     

Fluorescence decay time distribution analysis of cyclic enkephalin analogues. Influence of the solvents and configuration of amino acids in position 2 and 3 on changes in conformation

The lifetime distribution calculations were applied to study the influence of configuration of amino-acid residues in positions 2 and 3 on changes in conformation of the peptide chain of cyclic analogues of enkephalins containing a fluorescence energy donor and acceptor in different solvents. In all the solvents studied the lifetime distributions were bimodal. This testified to the presence of two families of conformations. In this paper the relationship between the population of each conformation and configuration of the residues in position 2 and 3, and the solvent used is discussed.     

Molecular simulation study of cooperativity in hydrophobic association: clusters of four hydrophobic particles

The multibody contribution to the potential of mean force (PMF) of hydrophobic association of four methane molecules in water was investigated by means of umbrella-sampling molecular dynamics. Two systems were considered: (i). a trigonal pyramid with three methane molecules at contact distance forming a fixed base, the fourth molecule being placed on the top with variable distance from the base; and (ii). a regular uniformly expanding tetrahedron. Methane-methane distances as far as 12.5 A, i.e. beyond the second solvent-separated minimum of the PMF, were considered to address the baseline problem. In contrast to the small effect in the three-body case studied previously (Protein Sci 9 (2000) 1235), the multibody contribution was found to amount to approximately 0.2 kcal/mol per methane-methane pair, or approximately 25% of the depth of the contact minimum in the PMF. The main effect of the multibody contribution to the PMF is a reduction of the height of the barrier between the contact and solvent separated minima and a narrowing of the region of its maximum, while the region of the contact minimum is affected only weakly. The reduction of the barrier is due to four-body contributions. The cooperative contributions to the PMF agree very well with those computed from the molecular surface of the systems under consideration, which further supports earlier observations that the molecular surface can be used with good accuracy to describe the energetics of hydrophobic association.     

Computer Modeling of the Solution Conformation of Cyclic Enkephalins

The probable conformations of two cyclic enkephalin analogs, DNS-cyclo[D-Dab-Gly-Trp-Leu] (I) and DNS-cyclo[D-Dab-Gly-Trp-D-Leu] (II) (DNS = dansyl), were determined by combining the results of NOE, vicinal coupling constant and fluorescence energy transfer measurements with theoretical calculations. The common feature of the conformations for both peptides is the presence of a -turn at residues 2 and 3.     

Computer Modeling of the Solution Conformation of Cyclic Enkephalins

Summary The probable conformations of two cyclic enkephalin analogs, DNS-cyclo[d-Dab-Gly-Trp-Leu] (I) and DNS-cyclo[d-Dab-Gly-Trp-d-Leu] (II) (DNS=dansyl), were determined by combining the results of NOE, vicinal coupling constant and fluorescence energy transfer measurements with theoretical calculations. The common feature of the conformations for both peptides is the presence of a β-turn at residues 2 and 3.     

Common functionally important motions of the nucleotide‐binding domain of Hsp70

The 70 kDa heat shock proteins (Hsp70) are a family of molecular chaperones involved in protein folding, aggregate prevention, and protein disaggregation. They consist of the substrate‐binding domain (SBD) that binds client substrates, and the nucleotide‐binding domain (NBD), whose cycles of nucleotide hydrolysis and exchange underpin the activity of the chaperone. To characterize the structure–function relationships that link the binding state of the NBD to its conformational behavior, we analyzed the dynamics of the NBD of the Hsp70 chaperone from Bos taurus (PDB 3C7N:B) by all‐atom canonical molecular dynamics simulations. It was found that essential motions within the NBD fall into three major classes: the mutual class, reflecting tendencies common to all binding states, and the ADP‐ and ATP‐unique classes, which reflect conformational trends that are unique to either the ADP‐ or ATP‐bound states, respectively. “Mutual” class motions generally describe “in‐plane” and/or “out‐of‐plane” (scissor‐like) rotation of the subdomains within the NBD. This result is consistent with experimental nuclear magnetic resonance data on the NBD. The “unique” class motions target specific regions on the NBD, usually surface loops or sites involved in nucleotide binding and are, therefore, expected to be involved in allostery and signal transmission. For all classes, and especially for those of the “unique” type, regions of enhanced mobility can be identified; these are termed “hot spots,” and their locations generally parallel those found by NMR spectroscopy. The presence of magnesium and potassium cations in the nucleotide‐binding pocket was also found to influence the dynamics of the NBD significantly. Proteins 2015; 83:282–299. © 2014 Wiley Periodicals, Inc.     

Acidic-basic properties of three alanine-based peptides containing acidic and basic side chains: comparison between theory and experiment

The purpose of this work was to evaluate the effect of the nature of the ionizable end groups, and the solvent, on their acid-base properties in alanine-based peptides. Hence, the acid-base properties of three alanine-based peptides: Ac-KK-(A)(7)-KK-NH(2) (KAK), Ac-OO-(A)(7)-DD-NH(2) (OAD), Ac-KK-(A)(7)-EE-NH(2) (KAE), where A, D, E, K, and O denote alanine, aspartic acid, glutamic acid, lysine, and ornithine, respectively, were determined in water and in methanol by potentiometry. With the availability of these data, the ability of two theoretical methods to simulate pH-metric titration of those peptides was assessed: (i) the electrostatically driven Monte Carlo method with the ECEPP/3 force field and the Poisson-Boltzmann approach to compute solvation energy (EDMC/PB/pH), and (ii) the molecular dynamics method with the AMBER force field and the Generalized Born model (MD/GB/pH). For OAD and KAE, pK(a1) and pK(a2) correspond to the acidic side chains. For all three compounds in both solvents, the pK(a1) value is remarkably lower than the pK(a) of a compound modeling the respective isolated side chain, which can be explained by the influence of the electrostatic field from positively charged ornithine or lysine side chains. The experimental titration curves are reproduced well by the MD/GB/pH approach, the agreement being better if restraints derived from NMR measurements are incorporated in the conformational search. Poorer agreement is achieved by the EDMC/PB/pH method.     

Further evidence for the absence of polyproline II stretch in the XAO peptide

It has been suggested that the alanine-based peptide with sequence Ac-XX-[A](7)-OO-NH(2), termed XAO where X denotes diaminobutyric acid and O denotes ornithine, exists in a predominantly polyproline-helix (P(II)) conformation in aqueous solution. In our recent work, we demonstrated that this "polyproline conformation" should be regarded as a set of local conformational states rather than as the overall conformation of the molecule. In this work, we present further evidence to support this statement. Differential scanning calorimetry measurements showed only a very small peak in the heat capacity of an aqueous solution of XAO at 57 degrees C, whereas the suggested transition to the P(II) structure should occur at approximately 30 degrees C. We also demonstrate that the temperature dependence of the (3)J(HNHalpha) coupling constants of the alanine residues can be explained qualitatively in terms of Boltzmann averaging over all local conformational states; therefore, this temperature dependence proves that a conformational transition does not occur. Canonical MD simulations with the solvent represented by the generalized Born model, and with time-averaged NMR-derived restraints, demonstrate the presence of an ensemble of structures with a substantial amount of local P(II) conformational states but not with an overall P(II) conformation.