Effect of Warmup Protocol and Sampling Time on Convergence of Molecular Dynamics Simulations of a DNA Dodecamer Using AMBER 4.1 and Particle-Mesh Ewald Method

Abstract

This report describes one 3000 ps and two 1500 ps molecular dynamic simulations on a TATA box containing dodecamer DNA duplex in a periodic box of TIP3P water molecules, using the AMBER 4.1 implementation of the particle-mesh Ewald method. We compare the effect of warmup protocol and simulation time length on the root-mean square deviation (RMSD) parameter. For the longer simulation, the RMSD computed for the 500–1000 ps time interval is representative of longer time intervals, including 500–3000 ps. The various warmup protocols do not appear to have a significant effect on the simulation results. Based on the present results, DNA sequence-dependent differences in RMSD, or related properties, should exceed two standard deviations before being attributed to non-simulation factors, such as warmup protocol and sampling time effects; we recommend a minimum criterion of at least a three standard deviation difference with a sampling period of at least 500–1000 ps. In addition, while end effects appear negligible there is a consistent dependence of RMSD on DNA helix length.     

A role for Ca2+ stores in kainate receptor-dependent synaptic facilitation and LTP at mossy fiber synapses in the hippocampus

Compared with NMDA receptor-dependent LTP, much less is known about the mechanism of induction of NMDA receptor-independent LTP; the most extensively studied form of which is mossy fiber LTP in the hippocampus. In the present study we show that Ca2+-induced Ca2+ release from intracellular stores is involved in the induction of mossy fiber LTP. This release also contributes to the kainate receptor-dependent component of the pronounced synaptic facilitation that occurs during high-frequency stimulation. We also present evidence that the trigger for this Ca2+ release is Ca2+ permeation through kainate receptors. However, these novel synaptic mechanisms can be bypassed when the Ca2+ concentration is raised (from 2 to 4 mM), via a compensatory involvement of L-type Ca2+ channels. These findings suggest that presynaptic kainate receptors at mossy fiber synapses can initiate a cascade involving Ca2+ release from intracellular stores that is important in both short-term and long-term plasticity.     

Lumbar muscle electromyographic dynamic topography during flexion-extension

The objective of this study is to introduce dynamic topography of surface electromyography (SEMG) to visualize lumbar muscle myoelectric activity and provides a new view to analyze muscle activity in vivo. A total of 20 healthy male subjects and 15 males LBP were enrolled. An electrode-array was applied to the lumbar region to collect SEMG. The root mean square (RMS) value was calculated for each channel, and then a 160×120 matrix was constructed using a linear cubic spline interpolation of each scan to create a 2-D color topographic image. Along a definite interval of action, a series of RMS topography matrices was concatenated as a function of position and time, to form a dynamic topographical video of lumbar muscle activity. Relative area (RA), relative width (RW), relative height (RH) and Width-to-Height Ratio (W/H) were chosen as the four quantitative parameters in measuring topographic features. Normal RMS dynamic topography was found to have a consistent, symmetric pattern with a high intensity area in the paraspinal area. LBP patients had a different RMS dynamic topography, with an asymmetric, broad, or disorganized distribution. Quantitative SEMG features were found significantly different between normal control and LBP. After physiotherapy rehabilitation, the dynamic topography images of LBP tended towards the normal pattern.There are obvious differences in lumbar muscle coordination between healthy subjects and LBP patients. The dynamic topography allows the continuous visualization of the distribution of surface EMG signals and the coordination of muscular contractions.     

Synthesis and structure-activity relationships of novel dipeptides and reduced dipeptides as ligands for melanocortin subtype-4 receptor

A series of benzylic piperazines (e.g., 4 and 5) attached to an 'address element', the dipeptide H-D-Tic-D-p-Cl-Phe-OH, 3 has been identified as ligands for the melanocortin subtype-4 receptor (MC4R). We describe herein the structure-activity relationship (SAR) studies on the N-terminal residue of the 'address element'. Several novel dipeptides and reduced dipeptides with high MC4R binding affinities and selectivity emerged from this SAR study.     

Privileged structure based ligands for melanocortin-4 receptors--aliphatic piperazine derivatives

Aliphatic carbocyclic replacement of the benzyl group of compound 1 yielded compounds with high affinity for the melanocortin-4 receptor (MC4R). Compounds with a cyclohexyl group showed a consistent high affinity, while different polar groups with less basicity were good replacements for the original diethyl amines. Substitution of the polar group found in these privileged structures with an aliphatic moiety produced compounds with high affinity for MC4R.     

Using Molecular Dynamics Simulations on Crambin to Evaluate the Suitability of Different Continuum Dielectric and Hydrogen Atom Models for Protein Simulations

Abstract

Molecular dynamics simulations of enzymes with enough explicit waters of solvation to realistically account for solute-solvent interactions can burden the computational resources required to perform the simulation by more than two orders of magnitude. Since enzyme simulations even with an implicit solvation model can be imposing for a supercomputer, it is important to assess the suitability of different continuum dielectric models for protein simulations. A series of 100-picosecond molecular dynamics simulations were performed on the X-ray crystal structure of the protein crambin to examine how well computed structures, obtained using seven continuum dielectric and two hydrogen atom models, agreed with the X-ray structure. The best level of agreement between computed and experimental structures was obtained using a constant dielectric of 2 and the all-hydrogen model. Continuum dielectric models of 1,1^*r, and 2^*r also led to computed structures in reasonably good agreement with the X-ray structure. In all cases, the all-hydrogen model gave better agreement than the united atom model, although, in one case, the difference was not significant. Dielectric models of 4, 80, and 4^*r with either hydrogen model yielded significantly poorer fits. It is especially noteworthy that the observed trends did not semiquantitatively converge until about 50 picoseconds into the simulations, suggesting that validation studies for protein calculations based on energy minimizations or short simulations should be viewed with caution.