Study of the interaction of GFG tripeptide with cesium perfluorooctanoate micelles by means of NMR spectroscopy and MD simulations

The interaction of glycyl-phenylalanyl-glycine (GFG) with bilayers formed by cesium perfluorooctanoate (CsPFO) in water was investigated in the isotropic phase by means of 1H NMR and molecular dynamics (MD) simulations. Details on the preferential location of the different residues of GFG were obtained from selective variations of chemical shift with peptide concentration and of line width in the presence of the paramagnetic ion Mn2+. The analysis of 1H NMR spectra recorded at different concentrations and temperatures allowed the association constant and the enthalpy change upon binding to be evaluated. MD simulations highlighted the hydrogen bonds formed between the different GFG functional groups and the micelle. Both NMR and MD gave indications of high affinity of GFG with the micelle, with the N-terminal residue anchoring on the surface via hydrogen bonds with the micelle COO(-) groups.     

Coarse-grained molecular dynamics simulations of the sphere to rod transition in surfactant micelles

Surfactant molecules self-assemble in aqueous solutions to form various micellar structures such as spheres, rods, or lamellae. Although phase transitions in surfactant solutions have been studied experimentally, their molecular mechanisms are still not well understood. In this work, we show that molecular dynamics (MD) simulations using the coarse-grained (CG) MARTINI force field and explicit CG solvent, validated against atomistic MD studies, can accurately represent micellar assemblies of cetyltrimethylammonium chloride (CTAC). The effect of salt on micellar structures is studied for aromatic anionic salts, e.g., sodium salicylate, and simple inorganic salts, e.g., sodium chloride. Above a threshold concentration, sodium salicylate induces a sphere to rod transition in the micelle. CG MD simulations are shown to capture the dynamics of this shape transition and support a mechanism based on the reduction in the micelle-water interfacial tension induced by the adsorption of the amphiphilic salicylate ions. At the threshold salt concentration, the interface is nearly saturated with adsorbed salicylate ions. Predictions of the effect of salt on the micelle structure in different CG solvent models, namely, single-site standard water and three-site polarizable water, show qualitative agreement. This suggests that phase transitions in aqueous micelle solutions could be investigated by using standard CG water models which allow for 3 orders of magnitude reduction in the computational time as compared to that required for atomistic MD simulations.     

Interpreting protein structural dynamics from NMR chemical shifts

In this investigation, semiempirical NMR chemical shift prediction methods are used to evaluate the dynamically averaged values of backbone chemical shifts obtained from unbiased molecular dynamics (MD) simulations of proteins. MD-averaged chemical shift predictions generally improve agreement with experimental values when compared to predictions made from static X-ray structures. Improved chemical shift predictions result from population-weighted sampling of multiple conformational states and from sampling smaller fluctuations within conformational basins. Improved chemical shift predictions also result from discrete changes to conformations observed in X-ray structures, which may result from crystal contacts, and are not always reflective of conformational dynamics in solution. Chemical shifts are sensitive reporters of fluctuations in backbone and side chain torsional angles, and averaged (1)H chemical shifts are particularly sensitive reporters of fluctuations in aromatic ring positions and geometries of hydrogen bonds. In addition, poor predictions of MD-averaged chemical shifts can identify spurious conformations and motions observed in MD simulations that may result from force field deficiencies or insufficient sampling and can also suggest subsets of conformational space that are more consistent with experimental data. These results suggest that the analysis of dynamically averaged NMR chemical shifts from MD simulations can serve as a powerful approach for characterizing protein motions in atomistic detail.     

Comparison of properties of Aib-rich peptides in crystal and solution: a molecular dynamics study.

In order to study the differences of the structural properties of Aib-rich peptides in solution and in the crystalline state, molecular dynamics (MD) simulations of the Aib-containing peptide II (pBrBz-(Aib)5-Leu-(Aib)2-OMe) were performed in the crystalline state, starting from two different conformers obtained experimentally by X-ray diffraction. The structural properties as derived from X-ray crystallography (e.g., torsional angles and hydrogen bonds) are well-reproduced in both constant-volume and constant-pressure simulations, although the force-field parameters used result in a too-high density of the crystals. Through comparison with the results from previous MD and nuclear magnetic resonance (NMR) studies of the very similar peptide I (Z-(Aib)s-Leu-(Aib)2-OMe) in dimethylsulfoxide (DMSO) solution, it is found that, in the crystal simulation, the conformational distribution of peptide II is much narrower than that in the solution simulation of peptide. I. This leads to a significant difference in 3 [symbol: see text] (HN, HC alpha) coupling constant values, in agreement with experimental data, whereas the NOE intensities or proton-proton distance bounds appear insensitive to the difference in conformational distribution. For small peptides the differences between their conformational distribution in the crystalline form and in solution may be much larger than for proteins, a fact which should be kept in mind when interpreting molecular properties in the solution state by using X-ray crystallographic data.     

Comparative analysis of the conformational profile of substance P using simulated annealing and molecular dynamics

The present study describes an extensive conformational search of substance P using two different computational methods. On the one hand, the peptide was studied using the iterative simulated annealing, and on the other, molecular dynamics simulations at 300 and 400 K. With the former method, the peptide was studied in vacuo with a dielectric constant of 80, whereas using the latter study the peptide was studied in a box of TIP3P water molecules. Analysis of the results obtained using both methodologies was carried out using an in-house methodology using a cluster analysis method based on information theory. Comparison of the two sampling methodologies and the different environment used in the calculations is also analyzed. Finally, the conformational motifs that are characteristic of substance P in a hydrophilic environment are presented and compared with the experimental results available in the literature.     

Molecular dynamics simulation and thermodynamic modeling of the self-assembly of the triterpenoids asiatic acid and madecassic acid in aqueous solution

The self-assembly behavior of the triterpenoids asiatic acid (AA) and madecassic acid (MA), both widely studied bioactive phytochemicals that are similar in structure to bile salts, were investigated in aqueous solution through atomistic-level molecular dynamics (MD) simulation. AA and MA molecules initially distributed randomly in solution were observed to aggregate into micelles during 75 ns of MD simulation. A "hydrophobic contact criterion" was developed to identify micellar aggregates from the computer simulation results. From the computer simulation data, the aggregation number of AA and MA micelles, the monomer concentration, the principal moments of the micelle radius of gyration tensor, the one-dimensional growth exhibited by AA and MA micelles as the aggregation number increases, the level of internal ordering within AA and MA micelles (quantified using two different orientational order parameters), the local environment of atoms within AA and MA in the micellar environment, and the total, hydrophilic, and hydrophobic solvent accessible surface areas of the AA and MA micelles were each evaluated. The MD simulations conducted provide insights into the self-assembly behavior of structurally complex, nontraditional surfactants in aqueous solution. Motivated by the high computational cost required to obtain an accurate estimate of the critical micelle concentrations (CMCs) of AA and MA from evaluation of the average monomer concentration present in the AA and MA simulation cells, a modified computer simulation/molecular-thermodynamic model (referred to as the MCS-MT model) was formulated to quantify the free-energy change associated with optimal AA and MA micelle formation in order to predict the CMCs of AA and MA. The predicted CMC of AA was found to be 59 microM, compared with the experimentally measured CMC of 17 microM, and the predicted CMC of MA was found to be 96 microM, compared with the experimentally measured CMC of 62 microM. The AA and MA CMCs predicted using the MCS-MT model are much more accurate than the CMCs inferred from the monomer concentrations of AA and MA present in the simulation cells after micelle self-assembly (2390 microM and 11,300 microM, respectively). The theoretical modeling results obtained for AA and MA indicate that, by combining computer simulation inputs with molecular-thermodynamic models of surfactant self-assembly, reasonably accurate estimates of surfactant CMCs can be obtained with a fraction of the computational expense that would be required by using computer simulations alone.     

Molecular Dynamics of Pectin Extension

Summary: A pectin 10mer under constant pulling speed and constant force was studied using the atomistic simulations. Molecular dynamics (MD) with the Amber99 and Amber-Glycam04 forcefields were performed. The main result of the present Amber-based MD simulations is that the two plateaux of the experimental force- extension dependence for pectin can be explained by transitions between three conformational states of pectin monomer ring (first from a chair (⁴C₁) to boat conformation and second from boat to an inverted chair (¹C₄) conformation). A multi-state dynamical model of single biopolymer extension under external force was elaborated and applied to extension of polymers with three-state monomers relevant to pectin.     

Molecular dynamics studies of cation aggregation in the room temperature ionic liquid [C10mim][Br] in aqueous solution

The structure of an aqueous 1-n-decyl-3-methylimidazolium bromide solution and its vapor-liquid interface has been studied using molecular dynamics (MD) simulations. Starting from an isotropic solution, spontaneous self-assembly of cations into small micellar aggregates has been observed. The decyl chains are buried inside the micelle to avoid unfavorable interactions with water, leaving the polar headgroups exposed to water. The cation aggregation numbers, ranging from 15 to 24 compare favorably with experimental estimates. Results are presented for the organization of solvent around the cations. The structure of the aggregates as determined from the present MD simulations does not support the staircase model proposed on the basis of nuclear magnetic resonance studies on similar aqueous ionic-liquid solutions. The distribution of ions in bulk solutions and at an air/water interface is also discussed.     

Protein folding in a reverse micelle environment: the role of confinement and dehydration

Characterization of the molecular interactions that stabilize the folded state of proteins including hydrogen bond formation, solvation, molecular crowding, and interaction with membrane environments is a fundamental goal of theoretical biophysics. Inspired by recent experimental studies by Gai and co-workers, we have used molecular dynamics simulations to explore the structure and dynamics of the alanine-rich AKA(2) peptide in bulk solution and in a reverse micelle environment. The simulated structure of the reverse micelle shows substantial deviations from a spherical geometry. The AKA(2) peptide is observed to (1) remain in a helical conformation within a spherically constrained reverse micelle and (2) partially unfold when simulated in an unconstrained reverse micelle environment, in agreement with experiment. While aqueous solvation is found to stabilize the N- and C-termini random coil portions of the peptide, the helical core region is stabilized by significant interaction between the nonpolar surface of the helix and the aliphatic chains of the AOT surfactant. The results suggest an important role for nonpolar peptide-surfactant and peptide-lipid interactions in stabilizing helical geometries of peptides in reverse micelle environments.     

Motion of a disordered polypeptide chain as studied by paramagnetic relaxation enhancements, 15N relaxation, and molecular dynamics simulations: how fast is segmental diffusion in denatured ubiquitin?

Molecular dynamics (MD) simulations have been widely used to analyze dynamic conformational equilibria of folded proteins, especially in relation to NMR observables. However, this approach found little use in the studies of disordered proteins, where the sampling of vast conformational space presents a serious problem. In this paper, we demonstrate that the latest advances in computation technology make it possible to overcome this limitation. The experimentally validated (calibrated) MD models allow for new insights into structure/dynamics of disordered proteins. As a test system, we have chosen denatured ubiquitin in solution with 8 M urea at pH 2. High-temperature MD simulations in implicit solvent have been carried out for the wild-type ubiquitin as well as MTSL-tagged Q2C, D32C, and R74C mutants. To recalibrate the MD data (500 K) in relation to the experimental conditions (278 K, 8 M urea), the time axes of the MD trajectories were rescaled. The scaling factor was adjusted such as to maximize the agreement between the simulated and experimental (15)N relaxation rates. The resulting effective length of the trajectories, 311 μs, ensures good convergence properties of the MD model. The constructed MD model was validated against the array of experimental data, including additional (15)N relaxation parameters, multiple sets of paramagnetic relaxation enhancements (PREs), and the radius of gyration. In each case, a near-quantitative agreement has been obtained, suggesting that the model is successful. Of note, the MD-based approach rigorously predicts the quantities that are inherently dynamic, i.e., dependent on the motional correlation times. This cannot be accomplished, other than in empirical fashion, on the basis of static structural models (conformational ensembles). The MD model was further used to investigate the relative translational motion of the MTSL label and the individual H(N) atoms. The derived segmental diffusion coefficients proved to be nearly uniform along the peptide chain, averaging to D = 0.49-0.55 × 10(-6) cm(2)/s. This result was verified by direct analysis of the experimental PRE data using the recently proposed Ullman-Podkorytov model. In this model, MTSL and H(N) moieties are treated as two tethered spheres undergoing mutual diffusion in a harmonic potential. The fitting of the experimental data involving D as a single adjustable parameter leads to D = 0.45 × 10(-6) cm(2)/s, in good agreement with the MD-based analyses. This result can be compared with the range of estimates obtained from the resonance energy transfer experiments, D = 0.2-6.0 × 10(-6) cm(2)/s.     

A conformational dynamics study of alpha-l-Rhap-(1-->2)[alpha-l-Rhap-(1-->3)]-alpha-l-Rhap-OMe in solution by NMR experiments and molecular simulations

The conformational preference of alpha-l-Rhap-(1-->2)[alpha-l-Rhap-(1-->3)]-alpha-l-Rhap-OMe in solution has been studied by NMR spectroscopy using one-dimensional (1)H,(1)H T-ROESY experiments and measurement of trans-glycosidic (3)J(C,H) coupling constants. Molecular dynamics (MD) simulations with a CHARMM22 type of force field modified for carbohydrates were performed with water as the explicit solvent. The homonuclear cross-relaxation rates, interpreted as effective proton-proton distances, were compared to those obtained from simulation. Via a Karplus torsional relationship, (3)J(C,H) values were calculated from simulation and compared to experimental data. Good agreement was observed between experimental data and the MD simulation, except for one inter-residue T-ROE between protons in the terminal sugar residues. The results show that the trisaccharide exhibits substantial conformational flexibility, in particular along the psi glycosidic torsion angles. Notably, for these torsions, a high degree of correlation (77%) was observed in the MD simulation revealing either psi(2)(+) psi(3)(+) or psi(2)(-)psi(3)(-) states. The simulations also showed that non-exoanomeric conformations were present at the phi torsion angles, but to a limited extent, with the phi(3) state populated to a larger extent than the phi(2) state. Further NMR analysis of the trisaccharide by translational diffusion measurements and (13)C T(1) relaxation experiments quantified global reorientation using an anisotropic model together with interpretation of the internal dynamics via the "model-free" approach. Fitting of the dynamically averaged states to experimental data showed that the psi(2)(+)psi(3)(+) state is present to approximately 49%, psi(2)(-) psi(3)(-) to approximately 39%, and phi(3) (non-exo) to approximately 12%. Finally, using a dynamic and population-averaged model, (1)H,(1)H T-ROE buildup curves were calculated using a full relaxation matrix approach and were found to be in excellent agreement with experimental data, in particular for the above inter-residue proton-proton interaction between the terminal residues.     

Molecular dynamics simulations of angiotensin II in aqueous and dimethyl sulfoxide environments

Angiotensin II (Ang II) is an octapeptidic hormone, which plays an important role in the mechanisms of blood pressure control. In this work, extensive molecular dynamics (MD) simulations have been carried out on this peptide, both in aqueous and in dimethyl sulfoxide (DMSO) environments. Experimentally proposed models for the structure of angiotensin II in both environments are not consensual and the results obtained have provided some further insight about the structural properties of this hormone. In these simulations, the N-terminus of Ang II in the aqueous environment has been associated with a considerable larger flexibility than the correspondent C-terminus, but this was not found in the case of the DMSO environment. This is consistent with the assumption that the biological activity of Ang II is associated with its C-terminal residues embedded in a hydrophobic environment of its AT1 receptor. Other features detected in DMSO environment were an H(His6 imidazole)-O(Phe8 carboxylate) hydrogen bond and a salt-bridge structure involving the Asp1 and Arg2 side chains. An additional important conformational feature is the spatial proximity between Tyr4 and His6 in both water and DMSO environments. This molecular feature may trigger the interest for the synthetic chemists to apply rational design for the synthesis of novel AT1 antagonists.     

Morphology and mechanical properties of surfactant aggregates at water-silica interfaces: molecular dynamics simulations

Dilute and concentrated surfactant systems at the solid-liquid interface are examined using classical molecular dynamics simulations. Particular emphasis is placed on understanding how surfactants aggregate and form the micellar structure, how micelles change shape at high concentrations in aqueous media and in the presence of hydrophilic surfaces, and at what force this micellar structure breaks apart during indentation of micelle-covered surfaces with a proximal probe microscope tip. The specific system of interest is C12TAB (n-dodecyltrimethylammonium bromide) surfactant in an aqueous medium that is modeled with empirical potentials. The simulations predict that the micelle structure in water is compact and either spherical or elliptical in shape. In the presence of a hydrophilic surface of silica, the structure evolves into a flat elliptical shape, in agreement with experimental findings. The simulated indentation of the micelle/silica system causes the micelle to break apart at an indentation force of about 1 nN and form a surfactant monolayer. The predicted force curve is in excellent agreement with experimental measurements.     

DIRECT‐ID: An automated method to identify and quantify conformational variations—application to β2‐adrenergic GPCR

The conformational dynamics of a macromolecule can be modulated by a number of factors, including changes in environment, ligand binding, and interactions with other macromolecules, among others. We present a method that quantifies the differences in macromolecular conformational dynamics and automatically extracts the structural features responsible for these changes. Given a set of molecular dynamics (MD) simulations of a macromolecule, the norms of the differences in covariance matrices are calculated for each pair of trajectories. A matrix of these norms thus quantifies the differences in conformational dynamics across the set of simulations. For each pair of trajectories, covariance difference matrices are parsed to extract structural elements that undergo changes in conformational properties. As a demonstration of its applicability to biomacromolecular systems, the method, referred to as DIRECT‐ID, was used to identify relevant ligand‐modulated structural variations in the β₂‐adrenergic (β₂AR) G‐protein coupled receptor. Micro‐second MD simulations of the β₂AR in an explicit lipid bilayer were run in the apo state and complexed with the ligands: BI‐167107 (agonist), epinephrine (agonist), salbutamol (long‐acting partial agonist), or carazolol (inverse agonist). Each ligand modulated the conformational dynamics of β₂AR differently and DIRECT‐ID analysis of the inverse‐agonist vs. agonist‐modulated β₂AR identified residues known through previous studies to selectively propagate deactivation/activation information, along with some previously unidentified ligand‐specific microswitches across the GPCR. This study demonstrates the utility of DIRECT‐ID to rapidly extract functionally relevant conformational dynamics information from extended MD simulations of large and complex macromolecular systems. © 2015 Wiley Periodicals, Inc.     

Fluorescence quenching induced by conformational fluctuations in unsolvated polypeptides

Time-resolved measurements were conducted to relate the fluorescence lifetimes of dye-derivatized polypeptides to local conformational dynamics in trapped, unsolvated peptide ions. This research was performed to better understand the intramolecular interactions leading to the observed increase of fluorescence quenching with temperature and, in particular, how this quenching is related to conformational fluctuations. Dye-derivatized polyproline ions, Dye-[Pro] n -Arg (+)-Trp, are formed by electrospray ionization and trapped in a variable-temperature quadrupole ion trap where they are exposed to a pulsed laser which excites fluorescence. Lifetime data exhibit fluorescence quenching as a result of an interaction between the dye and tryptophan (Trp) side chain. This result is consistent with solution measurements performed for comparison. The lifetime temperature dependence is closely fit over the range 150-463 K by an Arrhenius model of the ensemble averaged quenching rate, k q. Model fits of the measured lifetimes yield a frequency prefactor of approximately 10 (11) s (-1) for k q characteristic of collective motions of the side chains identified in molecular dynamics (MD) simulations. The data fits also yield activation barriers of approximately 0.3 eV, which are comparable to intramolecular electrostatic interactions calculated between the unshielded charge on the Arg residue and the dye. As a result, the quenching rate appears to be determined by the rate of conformational fluctuations and not by the rate of a specific quenching mechanism. The peptide sequence of Dye-Trp-[Pro] n -Arg (+) was also studied and identified a dependence of the quenching rate on the electrostatic field in the vicinity of the dye, Trp pair. Molecular dynamics simulations were performed over the range of experimental measurements to study trajectories relevant to the quenching interaction. The MD simulations indicate that as the temperature is increased, conformational fluctuations in the presence of strong electrostatic fields of the charged Arg (+) residue can result in both (a) an increased number of dye and Trp separations <8 A and (b) increased exothermicity for electron transfer reactions between the dye and Trp. Consequently, the MD simulations are consistent with increased fluorescence quenching with temperature resulting from the occurrence of conformers having specific positions of the dye, Trp, and Arg (+). As a result, the fluorescence lifetime provides a local probe of conformational fluctuations averaged over the ion ensemble.     

Evaluation of efficiency of reference potential spatial warping algorithm in conformational sampling of peptides

Molecular dynamics (MD) simulation using the reference potential spatial warping algorithm (REPSWA) is a promising method for computing the conformational equilibrium of a system with a rugged energy surface. Its effectiveness has previously been demonstrated using only a simplified model system of a hydrocarbon chain omitting nonbonded interaction terms from the potential energy function. To evaluate the applicability of REPSWA MD simulation to more realistic problems, we applied it to small peptide systems in an aqueous environment. Difficulties were encountered, however, forcing us to devise several modifications. We investigated their effectiveness in comparison to conventional constant-temperature and multicanonical MD simulations. We found that the sampling efficiency of the modified REPSWA MD after a careful optimization of its parameters was better than that of the constant-temperature MD and comparable to that of the multicanonical MD in several cases.     

Comprehensive Characterization of the Self-Folding Cavitand Dynamics

The conformational equilibria and guest exchange process of a resorcin[4]arene derived self-folding cavitand receptor have been characterized in detail by molecular dynamics simulations (MD) and ¹H EXSY NMR experiments. A multi-timescale strategy for exploring the fluxional behaviour of this system has been constructed, exploiting conventional MD and accelerated MD (aMD) techniques. The use of aMD allows the reconstruction of the folding/unfolding process of the receptor by sampling high-energy barrier processes unattainable by conventional MD simulations. We obtained MD trajectories sampling events occurring at different timescales from ns to s: 1) rearrangement of the directional hydrogen bond seam stabilizing the receptor, 2) folding/unfolding of the structure transiting partially open intermediates, and 3) guest departure from different folding stages. Most remarkably, reweighing of the biased aMD simulations provided kinetic barriers that are in very good agreement with those determined experimentally by ¹H NMR. These results constitute the first comprehensive characterization of the complex dynamic features of cavitand receptors. Our approach emerges as a valuable rational design tool for synthetic host-guest systems     

Molecular dynamics simulation of the A-DNA to B-DNA transition in aqueous RbCl solution

Unrestrained molecular dynamics (MD) simulations have been carried out to characterize the stability of DNA conformations and the dynamics of A-DNA→B-DNA conformational transitions in aqueous RbCl solutions. The PARM99 force field in the AMBER8 package was used to investigate the effect of RbCl concentration on the dynamics of the A→B conformational transition in the DNA duplex d(CGCGAATTCGCG)2 . Canonical Aand B-form DNA were assumed for the initial conformation and the final conformation had a length per complete turn that matched the canonical B-DNA. The DNA structure was monitored for 3.0 ns and the distances between the C5′ atoms were obtained from the simulations. It was found that all of the double stranded DNA strands of A-DNA converged to the structure of B-form DNA within 1.0 ns during the unrestrained MD simulations. In addition, increasing the RbCl concentration in aqueous solution hindered the A→B conformational transition and the transition in aqueous RbCl solution was faster than that in aqueous NaCl solution for the same electrolyte strength. The effects of the types and concentrations of counterions on the dynamics of the A→B conformational transition can be understood in terms of the variation in water activity and the number of accumulated counterions in the major grooves of A-DNA. The rubidium ion distributions around both fixed A-DNA and B-DNA were obtained using the restrained MD simulations to help explain the effect of RbCl concentration on the dynamics of the A→B conformational transition.