Adenosine triphosphate compartmentation in the rat heart: a 31P spin-lattice relaxation study

Longitudinal relaxation times (T1) of phosphorus compounds in the perfused rat heart and erythrocytes were measured using the 31P Driven-Equilibrium Single-Pulse Observation of T1 relaxation (DESPOT) method at 33 degrees C. Both creatine phosphate in the heart and the three phosphate groups of adenosine triphosphate (ATP) in erythrocytes showed single-exponential relaxation. The three phosphate groups of ATP in the heart, however, had two T1 components. The T1 values of the short and the long T1 components of the beta-phosphate of ATP were ca. 0.4 and 14 s, respectively. The fraction with the long T1 represented ca. 30% of the total ATP content. These results suggested that there were two major pools of intracellular ATP in the rat heart which could be determined by 31P NMR spectroscopy.     

Standards for molecular dynamics modelling and simulation of relaxation

An attempt is made to formulate a set of requirements for simulation and modelling of relaxation in dense media. Each requirement is illustrated by examples of numerical simulation of particles with different types of interaction given by soft-sphere, Lennard–Jones, embedded atom method or Coulomb potential. The approaches developed are expected to be universal for some classes of relaxation processes in liquids, fluids, crystals and plasmas.     

Molecular dynamics simulation of packing and mobility of lipids in bilayer membranes

A model of DSPC lipid membrane in gel and liquid-crystalline states has been developed. The parameters have been determined that enable one to calculate the molecular dynamics of lipid bilayers in the full-atromic approximation. The parameters of packing and mobility of lipid molecules for the liquid crystalline state of the bilayer have been calculated. The values agree well with experimental data. Based on the model of the liquid crystalline state of the membrane, a system in the gel-like state has been constructed. The model of the gel-like state reproduces well the packing of lipids in real bilayers, whereas the mobility of molecules in the gel-like state was found to be overestimated.     

Cholesterol dynamics in membranes of raft composition: a molecular point of view from 2H and 31P solid-state NMR

Lipidic membrane systems that have been reported to be composed of sphingomyelin (SM)-cholesterol (Chol) microdomains or "rafts" by Dietrich et al. [palmitoyloleoyl-phosphatidylcholine(POPC)/SM/Chol, 1/1/1; Dietrich, C., Bagatolli, L. A., Volovyk, Z. N., Thompson, N. L., Levi, M., Jacobson, K., and Gratton, E. (2001) Biophys. J. 80, 1417-1428] and by Schroeder et al. [SCRL: Liver-PC/Liver-phosphatidylethanolamine/SM/Cerebrosides/Chol, 1/1/1/1/2; Schroeder, R., London, E., and Brown, D. (1994) Proc. Natl. Acad. Sci. U.S.A. 91, 12130-12134] were investigated under the form of fully hydrated liposomes by the noninvasive solid-state (31)P and (2)H NMR method. Liposomes of binary lipid composition POPC/Chol and SM/Chol were also studied as boundary/control systems. All systems are found to be in the liquid-ordered phase (Lo) at physiological temperatures. Use of deuterium-labeled cholesterol afforded finding both the position of the sterol motional axis and its molecular order parameter. The axis of anisotropic rotation of cholesterol is such that the molecule is, on average, quasiperpendicular to the membrane plane, in all of the four systems investigated. Cholesterol order parameters greater than 0.8 are observed, indicating that the sterol is in a very motionally restricted environment in the temperature range 0-60 degrees C. The binary mixtures present "boundary" situations with the lowest values for POPC/Chol and the highest for SM/Chol. The SCRL raft mixture has the same ordering as the SM/Chol, i.e., the highest order parameter values over the temperature range. It demonstrates that in the SCRL mixture cholesterol dynamics is as in the binary system SM/Chol, therefore, suggesting that it might be depleted from the rest of the membrane to form complexes as if it were alone with SM. On the other hand, the mixture POPC/SM/Chol exhibits an intermediate ordering situation between those of SM/Chol and POPC/Chol. This strongly suggests that cholesterol could be in fast exchange, at the NMR time scale (milli- to microseconds), between two or more membrane regions of different dynamics and questions the statement of "rigid domains" made of SM and cholesterol in the model "raft" system POPC/SM/Chol.     

Cholesterol in model membranes. A molecular dynamics simulation.

Molecular dynamics simulations of a model membrane with inserted cholesterol molecules have been performed to study the perturbing influence of cholesterol. In the fluid phase of a lipid bilayer at 13 mol% concentration of cholesterol, local ordering of the hydrocarbon chains is induced. This perturbation decays with the distance from the cholesterol, and the effect extends 1.25 nm. It can be monitored in several ways, e.g., by an order parameter corresponding to deuterium nuclear magnetic resonance quadrupolar splittings, by the fraction of gauche bonds, or by the local bilayer thickness. At constant surface density, the local ordering is accompanied by disordering of the bulk phase, and, consequently, the net ordering effect is small. After compressing the system laterally in accordance with experimentally known surface areas, the bulk order parameters agree with those of a pure system, and the average order parameters are in accordance with experimental data. The necessity for this lateral compression is supported by calculated lateral pressures. At lower cholesterol concentration (3%), no direct perturbing effect is observed. A smaller lateral pressure than in a pure system indicates that the system with cholesterol is expected to have a smaller surface area, which would result in an increase of the order parameters, thus accounting for the experimental observations. The lack of spatial variation is, however, puzzling and may indicate a cooperative ordering effect.     

Charge pairing of headgroups in phosphatidylcholine membranes: A molecular dynamics simulation study

Molecular dynamics simulation of the hydrated dimyristoylphosphatidylcholine (DMPC) bilayer membrane in the liquid-crystalline phase was carried out for 5 ns to study the interaction among DMPC headgroups in the membrane/water interface region. The phosphatidylcholine headgroup contains a positively charged choline group and negatively charged phosphate and carbonyl groups, although it is a neutral molecule as a whole. Our previous study (Pasenkiewicz-Gierula, M., Y. Takaoka, H. Miyagawa, K. Kitamura, and A. Kusumi. 1997. J. Phys. Chem. 101:3677-3691) showed the formation of water cross-bridges between negatively charged groups in which a water molecule is simultaneously hydrogen bonded to two DMPC molecules. Water bridges link 76% of DMPC molecules in the membrane. In the present study we show that relatively stable charge associations (charge pairs) are formed between the positively and negatively charged groups of two DMPC molecules. Charge pairs link 93% of DMPC molecules in the membrane. Water bridges and charge pairs together form an extended network of interactions among DMPC headgroups linking 98% of all membrane phospholipids. The average lifetimes of DMPC-DMPC associations via charge pairs, water bridges and both, are at least 730, 1400, and over 1500 ps, respectively. However, these associations are dynamic states and they break and re-form several times during their lifetime.     

A 31P-NMR spin-lattice relaxation and 31P[1H] nuclear Overhauser effect study of sonicated small unilamellar phosphatidylcholine vesicles.

The motional properties of the inner and outer monolayer headgroups of egg phosphatidylcholine (PC) small unilamellar vesicles (SUV) were investigated by 31P-NMR temperature-dependent spin-lattice relaxation time constant (T1) and 31P[1H] nuclear Overhauser effect (NOE) analyses. Three different aspects of the dynamics of PC headgroups were investigated using the T1 analysis. First, differences in the dynamics of the headgroup region of both surfaces of the SUV were measured after application of a chemical shift reagent, PrCl3, to either the extra- or intravesicular volumes. Second, the ability of the T1 experiment to resolve the different motional states was evaluated in the absence of shift reagent. Third, comparison between correlation times obtained from a resonance frequency dependent 31P[1H] NOE analysis allowed a determination of the applicability of a simplified motional model to describe phosphorus dipolar relaxation. Temperature-dependent 31P-NMR T1 values obtained for the individual monolayers at 81.0 and 162.0 MHz were modelled assuming that phosphorus undergoes both a dipolar and an anisotropic chemical shielding relaxation mechanism, each being described by the same correlation time, tau. At 162.0 MHz, the position of the T1 minimum for the inner monolayer was 9 degrees higher than that of the outer region, indicating a higher level of motional restriction for the inner leaflet, in agreement with 31P[1H] NOE measurements. The 162.0 MHz T1 profile of the combined SUV monolayers exhibited a smooth minimum located at the midpoint of the monolayer minima positions, effectively masking the presence of the individual surfaces. 31P[1H] NOE results obtained at 32.3, 81.0 and 162.0 MHz did not agree with those predicted from a simple dipolar relaxation model. These results suggest a T1-temperature method can neither discriminate two or more closely related motional time scales in a heterogeneous environment (such as incorporation of protein into lipid bilayers) nor allow accurate determination of the correlation time at the position of the minimum when the dipolar relaxation rate makes a significant contribution to the overall rate.     

Head group and chain behavior in biological membranes: a molecular dynamics computer simulation.

A computer-modeled hydrated bilayer model of the lipid 2,3-dimyristoyl-D-glycero-1-phosphorylcholine in the L alpha phase was built. Particular care was taken in building the starting structure with the inclusion of structural detail reported in experiments on the L alpha phase. Molecular dynamics simulations using the molecular dynamics and energy refinement program AMBER 3.1 force field with an optimized parameters for liquid simulation parameter set were run to study the motions and conformations of the lipid molecules and characterize the behavior and structure of the head groups and the hydrocarbon lipid chains. Although the head groups were observed to show great flexibility, certain head-group torsion combinations appeared favored. The observed tilt of the lipid chains is discussed and is consistent with previous experimental findings. Motion of the lipid chains is shown to be correlated with those chains immediately surrounding, but correlation with chains more distant varies with time.     

Long-term molecular dynamics simulation of copper azurin: structure, dynamics and functionality

A long-term molecular dynamics simulation (1.1 ns), at 300 K, of fully hydrated azurin has been performed to put into relationship the protein dynamics to functional properties with particular attention to those structural elements involved in the electron transfer process. A detailed analysis of the root mean square deviations and fluctuations and of the intraprotein H-bonding pattern has allowed us to demonstrate that a rigid arrangement of the beta-stranded protein skeleton is maintained during the simulation run, while a large mobility is registered in the solvent-exposed connecting regions (turns) and in the alpha-helix. Moreover, the structural elements, likely involved in the electron transfer path, show a stable H-bonding arrangement and low fluctuations. Analysis of the dynamical cross-correlation map has revealed the existence of correlated motions among residues connected by hydrogen bonds and of correlated and anti-correlated motions between regions which are supposed to be involved in the functional process, namely the hydrophobic patch and the regions close to the copper reaction center. The results are briefly discussed also in connection to the current through-bond tunneling model for the electron transfer process. Finally, a comparison with the structural and the dynamical behaviour of plastocyanin, whose structure and functional role are very similar to those of azurin, has been performed.     

NMR relaxation processes of 31P in macromolecules

³¹P-Nmr relaxation parameters (spin-lattice relaxation time, linewidth, and nuclear Overhauser effect) were obtained at three different frequencies for poly(U) and a well-defined (145 ± 3 base-pair) fragment of DNA in solution. Data sets for the two samples were analyzed by theories which included relaxation by the mechanisms of ³¹P chemical shift anisotropy as well as by ¹H-³¹P dipole–dipole interaction. Neither data set could be satisfactorily described by a single correlation time. A model of a rigid rotor most nearly fits the data for the DNA molecule. Parameters obtained from the least-square fit indicate (1) that the DNA undergoes anisotropic reorientation with a correlation time τ₀ = 6.5 × 10^?7 sec for the end-to-end motion, (2) the ratio of diffusion constants D_∥/D_⊥ is 91, and (3) that the linewidth is due to chemical shift dispersion to the extent of 0.5 ppm. Some deviations of the calculated from the observed values suggested that significant torsional and bending motions may also take place for this DNA. Another model which contains isotropic motion but with a broad distribution of correlation times was required to fit the data for poly(U). A log ? χ² distribution function of correlation times [Scheafer, J. (1973) Macromolecules 6 , 881–888] described well the motion of poly(U) with the average correlation time τ = 3.3 × 10^?9 sec and a distribution parameter p = 14.     

A simple method for the simulation of steady-state diffusion through membranes: pressure-tuned,boundary-driven molecular dynamics

We present a novel molecular dynamics-based simulation technique for investigating gas transport through membranes. In our simulations, the main control parameters are the partial pressure for the components on the input side of the membrane and the total pressure on the output side. The essential point of our scheme is that this pressure control should be realised by adjusting the particle numbers in the input and output side control cells indirectly. Although this perturbation is applied sufficiently far from the membrane, the bulk-phase properties are well controlled in a simulation cell of common size. Numerical results are given for silicalite-1 membrane with permeating CH₄, CO₂, H₂ and N₂ gases as well as with binary mixtures of CO₂ with the other three components. To describe interactions between particles, we used the simple shifted and cut Lennard–Jones potential with parameters available in the literature. It is expected that the proposed technique can be applied to several other types of membranes and transported fluids in order to support the development of a deeper understanding of separation processes.     

Boundary lipids of the nicotinic acetylcholine receptor: Spontaneous partitioning via coarse-grained molecular dynamics simulation

Reconstituted nicotinic acetylcholine receptors (nAChRs) exhibit significant gain-of-function upon addition of cholesterol to reconstitution mixtures, and cholesterol affects the organization of nAChRs within domain-forming membranes, but whether nAChR partitions to cholesterol-rich liquid-ordered (“raft” or l_o) domains or cholesterol-poor liquid-disordered (l_do) domains is unknown. We use coarse-grained molecular dynamics simulations to observe spontaneous interactions of cholesterol, saturated lipids, and polyunsaturated (PUFA) lipids with nAChRs. In binary Dipalmitoylphosphatidylcholine:Cholesterol (DPPC:CHOL) mixtures, both CHOL and DPPC acyl chains were observed spontaneously entering deep “non-annular” cavities in the nAChR TMD, particularly at the subunit interface and the β subunit center, facilitated by the low amino acid density in the cryo-EM structure of nAChR in a native membrane. Cholesterol was highly enriched in the annulus around the TMD, but this effect extended over (at most) 5–10 Å. In domain-forming ternary mixtures containing PUFAs, the presence of a single receptor did not significantly affect the likelihood of domain formation. nAChR partitioned to any cholesterol-poor l_do domain that was present, regardless of whether the l_do or l_o domain lipids had PC or PE headgroups. Enrichment of PUFAs among boundary lipids was positively correlated with their propensity for demixing from cholesterol-rich phases. Long n-3 chains (tested here with Docosahexaenoic Acid, DHA) were highly enriched in annular and non-annular embedded sites, partially displacing cholesterol and completely displacing DPPC, and occupying sites even deeper within the bundle. Shorter n-6 chains were far less effective at displacing cholesterol from non-annular sites.     

Anisotropic 2H-nuclear magnetic resonance spin-lattice relaxation in cerebroside- and phospholipid-cholesterol bilayer membranes.

The axially symmetric powder pattern 2H-nuclear magnetic resonance (NMR) lineshapes observed in the liquid crystalline phase of pure lipid or lipid/cholesterol bilayers are essentially invariant to temperature, or, equivalently, to variations in the correlation times characterizing C-2H bond reorientations. In either of these melted phases, where correlation times for C-2H bond motions are shorter than 10(-7) s, information on the molecular dynamics of the saturated hydrocarbon chain would be difficult to obtain using lineshape analyses alone, and one must resort to other methods, such as the measurement of 2H spin-lattice relaxation rates, in order to obtain dynamic information. In pure lipid bilayers, the full power of the spin-lattice relaxation technique has yet to be realized, since an important piece of information, namely the orientation dependence of the 2H spin-lattice relaxation rates is usually lost due to orientational averaging of T1 by rapid lateral diffusion. Under more favorable circumstances, such as those encountered in the lipid/cholesterol mixtures of this study, the effects of orientational averaging by lateral diffusion are nullified, due to either a marked reduction (by at least an order of magnitude) in the diffusion rate, or a marked increase in the radii of curvature of the liposomes. In either case, the angular dependence of 2H spin-lattice relaxation is accessible to experimental study, and can be used to test models of molecular dynamics in these systems. Simulations of the partially recovered lineshapes indicate that the observed T1 anisotropies are consistent with large amplitude molecular reorientation of the C-2H bond among a finite number of sites.(ABSTRACT TRUNCATED AT 250 WORDS)     

Cytosine ribose flexibility in DNA: a combined NMR 13C spin relaxation and molecular dynamics simulation study

Using (13)C spin relaxation NMR in combination with molecular dynamic (MD) simulations, we characterized internal motions within double-stranded DNA on the pico- to nano-second time scale. We found that the C-H vectors in all cytosine ribose moieties within the Dickerson-Drew dodecamer (5'-CGCGAATTCGCG-3') are subject to high amplitude motions, while the other nucleotides are essentially rigid. MD simulations showed that repuckering is a likely motional model for the cytosine ribose moiety. Repuckering occurs with a time constant of around 100 ps. Knowledge of DNA dynamics will contribute to our understanding of the recognition specificity of DNA-binding proteins such as cytosine methyltransferase.     

Coarse-grained molecular dynamics simulation of interactions between cyclic lipopeptide Bacillomycin D and cell membranes

According to experimental studies, Bacillomycin D has strong antimicrobial activities, but the antimicrobial mechanism is still unknown. In this paper, the interaction mechanisms between this cyclic lipopeptide and three different charged cell membranes are studied via Coarse-Grained Molecular Dynamics (CG MD) simulations. A specific CG model for the cyclic lipopeptide Bacillomycin D was developed. The insertion of cyclic lipopeptide Bacillomycin D into DOPC, DOPC/DPPA and DOPC/DOTAP cell membranes was investigated. The position distribution and stability of Bacillomycin D in the three different cell membranes were analysed and compared based on density profile calculations. Additionally, we focused on the Radial Distribution Function (RDF) curves between amino acid residues with negative charges or strong hydrophobic properties and the head groups of two different cell membranes. Based on changes in the curvature of the three membranes, the cyclic lipopeptide Bacillomycin D can cause localised surface protrusions in DOPC/DOTAP membranes, inward depressions in the surface of DOPC/DPPA membranes and inhibition deformation in the surface of DOPC membranes. This study will help to further understand the antibacterial mechanism of the cyclic lipopeptide Bacillomycin D and provide a basis for the development of new antibiotics.     

Molecular dynamics simulation and steered molecular dynamics simulation on irisin dimers

Irisin is found closely associated with promoting the browning of beige fat cells in white adipose tissue. The crystal structure reveals that irisin forms a continuous inter-subunit β-sheet dimer. Here, molecular dynamics (MD) simulation and steered molecular dynamics (SMD) simulation were performed to investigate the dissociation process and the intricate interactions between the two irisin monomers. In the process of MD, the interactions between the monomers were roughly analyzed through the average numbers of both hydrophobic contacts and H-bonds. Then, SMD was performed to investigate the accurate interaction energy between the monomers. By the analysis of dissociation energy, the van der Waals (vdW) force was identified as the major energy to maintain the dimer structure, which also verified the results of MD simulation. Meanwhile, 11 essential residues were discovered by the magnitude of rupture force during dissociation. Among them, residues Arg75, Glu79, Ile77, Ala88, and Trp90 were reported in a previous study using the method of mutagenesis and size exclusion chromatography, and several new important residues (Arg72, Leu74, Phe76, Gln78, Val80, and Asp91) were also identified. Interestingly, the new important residues that we discovered and the important residues that were reported are located in the opposite side of the β-sheet of the dimer.     

Oligoglycine surface structures: molecular dynamics simulation

The full-atomic molecular dynamics (MD) simulation of adsorption mode for diantennary oligoglycines [H-Gly4-NH(CH2)5]2 onto graphite and mica surface is described. The resulting structure of adsorption layers is analyzed. The peptide second structure motives have been studied by both STRIDE (structural identification) and DSSP (dictionary of secondary structure of proteins) methods. The obtained results confirm the possibility of polyglycine II (PGII) structure formation in diantennary oligoglycine (DAOG) monolayers deposited onto graphite surface, which was earlier estimated based on atomic-force microscopy measurements.     

Membrane electroporation: a molecular dynamics simulation

We present results of molecular dynamics simulations of lipid bilayers under a high transverse electrical field aimed at investigating their electroporation. Several systems are studied, namely 1), a bare bilayer, 2), a bilayer containing a peptide nanotube channel, and 3), a system with a peripheral DNA double strand. In all systems, the applied transmembrane electric fields (0.5 V.nm(-1) and 1.0 V.nm(-1)) induce an electroporation of the lipid bilayer manifested by the formation of water wires and water channels across the membrane. The internal structures of the peptide nanotube assembly and that of the DNA strand are hardly modified under field. For system 2, no perturbation of the membrane is witnessed at the vicinity of the channel, which indicates that the interactions of the peptide with the nearby lipids stabilize the bilayer. For system 3, the DNA strand migrates to the interior of the membrane only after electroporation. Interestingly enough, switching of the external transmembrane potential in cases 1 and 2 for few nanoseconds is enough to allow for complete resealing and reconstitution of the bilayer. We provide evidence that the electric field induces a significant lateral stress on the bilayer, manifested by surface tensions of magnitudes in the order of 1 mN.m(-1). This study is believed to capture the essence of several dynamical phenomena observed experimentally and provides a framework for further developments and for new applications.     

Statistical prediction and molecular dynamics simulation

We describe a statistical approach to the validation and improvement of molecular dynamics simulations of macromolecules. We emphasize the use of molecular dynamics simulations to calculate thermodynamic quantities that may be compared to experimental measurements, and the use of a common set of energetic parameters across multiple distinct molecules. We briefly review relevant results from the theory of stochastic processes and discuss the monitoring of convergence to equilibrium, the obtaining of confidence intervals for summary statistics corresponding to measured quantities, and an approach to validation and improvement of simulations based on out-of-sample prediction. We apply these methods to replica exchange molecular dynamics simulations of a set of eight helical peptides under the AMBER potential using implicit solvent. We evaluate the ability of these simulations to quantitatively reproduce experimental helicity measurements obtained by circular dichroism. In addition, we introduce notions of statistical predictive estimation for force-field parameter refinement. We perform a sensitivity analysis to identify key parameters of the potential, and introduce Bayesian updating of these parameters. We demonstrate the effect of parameter updating applied to the internal dielectric constant parameter on the out-of-sample prediction accuracy as measured by cross-validation.     

Conformational dynamics of RNA-peptide binding: a molecular dynamics simulation study

Molecular dynamics simulations of the binding of the heterochiral tripeptide KkN to the transactivation responsive (TAR) RNA of HIV-1 is presented, using an all-atom force field with explicit water. To obtain starting structures for the TAR-KkN complex, semirigid docking calculations were performed that employ an NMR structure of free TAR RNA. The molecular dynamics simulations show that the starting structures in which KkN binds to the major groove of TAR (as it is the case for the Tat-TAR complex of HIV-1) are unstable. On the other hand, the minor-groove starting structures are found to lead to several binding modes, which are stabilized by a complex interplay of stacking, hydrogen bonding, and electrostatic interactions. Although the ligand does not occupy the binding position of Tat protein, it is shown to hinder the interhelical motion of free TAR RNA. The latter is presumably necessary to achieve the conformational change of TAR RNA to bind Tat protein. Considering the time evolution of the trajectories, the binding process is found to be ligand-induced and cooperative. That is, the conformational rearrangement only occurs in the presence of the ligand and the concerted motion of the ligand and a large part of the RNA binding site is necessary to achieve the final low-energy binding state.