Finding a suitable system scale to optimize program performance on software DSM systems

Recently, software distributed shared memory systems have successfully provided an easy user interface to parallel user applications on distributed systems. In order to prompt program performance, most of DSM systems usually were greedy to utilize all of available processors in a computer network to execute user programs. However, using more processors to execute programs cannot necessarily guarantee to obtain better program performance. The overhead of paralleling programs is increased by the addition in the number of processors used for program execution. If the performance gain from program parallel cannot compensate for the overhead, increasing the number of execution processors will result in performance degradation and resource waste. In this paper, we proposed a mechanism to dynamically find a suitable system scale to optimize performance for DSM applications according to run-time information. The experimental results show that the proposed mechanism can precisely predict the processor number that will result in the best performance and then effectively optimize the performance of the test applications by adapting system scale according to the predicted result. Yi-Chang Zhuang received his B.S., M.S. and Ph.D. degrees in electrical engineering from National Cheng Kung University in 1995, 1997, and 2004. He is currently working as an engineer at Industrial Technology Research Institute in Taiwan. His research interests include object-based storage, file systems, distributed systems, and grid computing. Jyh-Biau Chang is currently an assistant professor at the Information Management Department of Leader University in Taiwan. He received his B.S., M.S. and Ph.D. degrees from Electrical Engineering Department of National Cheng Kung University in 1994, 1996, and 2005. His research interest is focused on cluster and grid computing, parallel and distributed system, and operating system. Tyng-Yeu Liang is currently an assistant professor who teaches and studies at Department of Electrical Engineering, National Kaohsiung University of Applied Sciences in Taiwan. He received his B.S., M.S. and Ph.D. degrees from National Cheng Kung University in 1992, 1994, and 2000. His study is interested in cluster and grid computing, image processing and multimedia. Ce-Kuen Shieh currently is a professor at the Electrical Engineering Department of National Cheng Kung University in Taiwan. He is also the chief of computation center at National Cheng Kung University. He received his Ph.D. degree from the Department of Electrical Engineering of National Cheng Kung University in 1988. He was the chairman of the Electrical Engineering Department of National Cheng Kung University from 2002 to 2005. His research interest is focused on computer network, and parallel and distributed system. Laurence T. Yang is a professor at the Department of Computer Science, St. Francis Xavier University, Canada. His research includes high performance computing and networking, embedded systems, ubiquitous/pervasive computing and intelligence, and autonomic and trusted computing.     

A Parallel Algorithm for Nonequilibrium Molecular Dynamics Simulation of Shear Flow on Distributed Memory Machines

Abstract

An algorithm is described which allows Nonequilibrium Molecular Dynamics (NEMD) simulations of a fluid undergoing planar Couette flow (shear flow) to be carried out on a distributed memory parallel processor using a (spatial) domain decomposition technique. Unlike previous algorithms, this algorithm uses a co-moving, or Lagrangian, simulation box. Also, the shape of the simulation box changes throughout the course of the simulation. The algorithm, which can be used for two or three dimensional systems, has been tested on a Fujitsu AP1000 Parallel computer with 128 processors.     

Conformational ensemble of an intrinsically flexible loop in mitochondrial import protein Tim21 studied by modeling and molecular dynamics simulations

BackgroundTim21, a subunit of a highly dynamic translocase of the inner mitochondrial membrane (TIM23) complex, translocates proteins by interacting with subunits in the translocase of the outer membrane (TOM) complex and Tim23 channel in the TIM23 complex. A loop segment in Tim21, which is in close proximity of the binding site of Tim23, has different conformations in X-ray, NMR and new crystal contact-free space (CCFS) structures. MD simulations can provide information on the structure and dynamics of the loop in solution.MethodsThe conformational ensemble of the loop was characterized using loop modeling and molecular dynamics (MD) simulations.ResultsMD simulations confirmed mobility of the loop. Multidimensional scaling and clustering were used to characterize the dynamic conformational ensemble of the loop. Free energy landscape showed that the CCFS crystal structure occupied a low energy region as compared to the conventional X-ray crystal structure. Analysis of crystal packing indicates that the CCFS provides larger conformational space for the motions of the loop.ConclusionsOur work reported the conformational ensemble of the loop in solution, which is in agreement with the structure obtained from CCFS approach. The combination of the experimental techniques and computational methods is beneficial for studying highly flexible regions of proteins.General significanceComputational methods, such as loop modeling and MD simulations, have proved to be useful for studying conformational flexibility of proteins. These methods in integration with experimental techniques such as CCFS has the potential to transform the studies on flexible regions of proteins.     

The Fip35 WW Domain Folds with Structural and Mechanistic Heterogeneity in Molecular Dynamics Simulations

We describe molecular dynamics simulations resulting in the folding the Fip35 Hpin1 WW domain. The simulations were run on a distributed set of graphics processors, which are capable of providing up to two orders of magnitude faster computation than conventional processors. Using the Folding@home distributed computing system, we generated thousands of independent trajectories in an implicit solvent model, totaling over 2.73 ms of simulations. A small number of these trajectories folded; the folding proceeded along several distinct routes and the system folded into two distinct three-stranded β-sheet conformations, showing that the folding mechanism of this system is distinctly heterogeneous.     

Brownian and Essential Dynamics Studies of the HIV-1 Integrase Catalytic Domain

Abstract

The three-dimensional structure of the active site region of the enzyme HIV-1 integrase is not unambiguously known. This region includes a flexible peptide loop that cannot be well resolved in crystallographic determinations. Here we present two different computional approaches with different levels of resolution and on different time-scales to understand this flexibility and to analyze the dynamics of this part of the protein. We have used molecular dynamics simulations with an atomic model to simulate the region in a realistic and reliable way for 1 ns. It is found that parts of the loop wind up after 300 ps to extend an existing helix. This indicates that the helix is longer than in the earlier crystal structures that were used as basis for this study. Very recent crystal data confirms this finding, underlining the predictive value of accurate MD simulations. Essential dynamics analysis of the MD trajectory yields an anharmonic motion of this loop. We have supplemented the MD data with a much lower resolution Brownian dynamics simulation of 600 ns length. It provides ideas about the slow-motion dynamics of the loop. It is found that the loop explores a conformational space much larger than in the MD trajectory, leading to a “gating”-like motion with respect to the active site.     

Molecular Dynamics Simulations of Phospholipid Bilayers

Abstract

Molecular dynamics (MD) simulations at 37°C have been performed on three phospholipid bilayer systems composed of the lipids DLPE, DOPE, and DOPC. The model used included 24 explicit lipid molecules and explicit waters of solvation in the polar head group regions, together with constant-pressure periodic boundary conditions in three dimensions. Using this model, a MD simulation samples part of an infinite planar lipid bilayer. The lipid dynamics and packing behavior were characterized. Furthermore, using the results of the simulations, a number of diverse properties including bilayer structural parameters, hydrocarbon chain order parameters, dihedral conformations, electron density profile, hydration per lipid, and water distribution along the bilayer normal were calculated. Many of these properties are available for the three lipid systems chosen, making them well suited for evaluating the model and protocols used in these simulations by direct comparisons with experimental data. The calculated MD behavior, chain disorder, and lipid packing parameter, i.e. the ratio of the effective areas of hydrocarbon tails and head group per lipid (a_t/a_h), correctly predict the aggregation preferences of the three lipids observed experimentally at 37°C, namely: a gel bilayer for DLPE, a hexagonal tube for DOPE, and a liquid crystalline bilayer for DOPC. In addition, the model and conditions used in the MD simulations led to good agreement of the calculated properties of the bilayers with available experimental results, demonstrating the reliability of the simulations. The effects of the cis unsaturation in the hydrocarbon chains of DOPE and DOPC, compared to the fully saturated one in DLPE, as well as the effects of the different polar head groups of PC and PE with the same unsaturated chains on the lipid packing and bilayer structure have been investigated. The results of these studies indicate the ability of MD methods to provide molecular-level insights into the structure and dynamics of lipid assemblies.     

Combined use of periodic reaction field and coarse-grained molecular dynamics simulations. I. phospholipid monolayer systems

Abstract

A periodic reaction field based on a linear-combination-based isotropic periodic sum (LIPS) method was applied for coarse-grained molecular dynamics simulations of zwitterionic lipid systems. In phospholipid monolayer systems with various number of lipid molecules, the density profile, lipid orientation and surface tension were mainly calculated using the periodic reaction field and Ewald sum. The results from the periodic reaction field were almost equal to that from the Ewald sum. It is concluded that the periodic reaction field method has a great possibility to provide a high accuracy in determining coarse-grained zwitterionic lipid systems.     

Crystal contact-free conformation of an intrinsically flexible loop in protein crystal: Tim21 as the case study

BackgroundIn protein crystals, flexible loops are frequently deformed by crystal contacts, whereas in solution, the large motions result in the poor convergence of such flexible loops in NMR structure determinations. We need an experimental technique to characterize the structural and dynamic properties of intrinsically flexible loops of protein molecules.MethodsWe designed an intended crystal contact-free space (CCFS) in protein crystals, and arranged the flexible loop of interest in the CCFS. The yeast Tim 21 protein was chosen as the model protein, because one of the loops (loop 2) is distorted by crystal contacts in the conventional crystal.ResultsYeast Tim21 was fused to the MBP protein by a rigid α-helical linker. The space created between the two proteins was used as the CCFS. The linker length provides adjustable freedom to arrange loop 2 in the CCFS. We re-determined the NMR structure of yeast Tim21, and conducted MD simulations for comparison. Multidimensional scaling was used to visualize the conformational similarity of loop 2. We found that the crystal contact-free conformation of loop 2 is located close to the center of the ensembles of the loop 2 conformations in the NMR and MD structures.ConclusionsLoop 2 of yeast Tim21 in the CCFS adopts a representative, dominant conformation in solution.General significanceNo single powerful technique is available for the characterization of flexible structures in protein molecules. NMR analyses and MD simulations provide useful, but incomplete information. CCFS crystallography offers a third route to this goal.     

The Stabilizing Effects of O-glycosylation on the Secondary Structural Integrity of the Designed α-loop-α motif by Molecular Dynamics Simulations

Abstract

In this study, various 400 ps molecular dynamics simulations were conducted to determine the stabilizing effect of O-glycosylation on the secondary structural integrity of the design α-loop-α motif, which has the optimal loop length of 7 Gly residues (denoted as N-A₁₆G₇A₁₆-C). In general, O-glycosylation stabilizes the structural integrity of the model peptide regardless of the length and position of glycosylation sites because it decreases the opportunity for water molecules to compete for the intramolecular hydrogen bonds. The designed peptide exhibits the highest helicity when residues 11 and 31 are replaced with Ser residues followed by O-linked with 3 galactose residues, representing the “face-to-face” glycosylation near the loop. In this case, the loop exhibits an extended conformation and several new hydrogen bonds are observed between the main chain of the loop and the galactose residues, resulting in decreasing the fluctuation and increasing the stability of the entire peptide. When the glycosylation are made close to the loop, the secondary structural integrity of the α-loop-α motif increases with the number of galactose residues. In addition, “face- to-face” glycosylation increases the structural integrity of this motif to a greater extent than “back-to-back” glycosylation. However, when the glycosylation are created away from the loop and near the N- and C-termini, no general rule is found for the stabilizing effect.     

Predicting doxorubicin drug delivery by single-walled carbon nanotube through cell membrane in the absence and presence of nicotine molecules: a molecular dynamics simulation study

Abstract

In order to study the interaction of the anticancer agent Doxorubicin with the single-walled carbon nanotubes with different diameters as drug delivery systems, the molecular dynamics (MD) simulations have been used. Also, for design and development of intracellular Doxorubicin drug delivery systems, a series of steered MD simulations are applied to explore the possibility of encapsulated Doxorubicin–carbon nanotube penetration through a lipid bilayer in presence and absence of Nicotine molecules at different pulling rates. Our simulation results showed that in spite of the adsorption of drug molecules on the outer sidewall of the nanotubes, the spontaneous localization of one Doxorubicin molecule into the cavity of the nanovectors with larger diameters is observed. It is found that the presence of Nicotine molecules in extracellular medium increases the required force for pulling nanotube-encapsulated drug as well as the required time for penetration process, especially at higher velocity. Also, the entering process of the Nicotine molecules into the carbon nanotube causes that the encapsulated drug molecule is fully released in the hydrophobic phase of the lipid bilayer.

Communicated by Ramaswamy H. Sarma     

Molecular dynamics simulation of proteins under high pressure: Structure,function and thermodynamics

BackgroundMolecular dynamics (MD) simulation is well-recognized as a powerful tool to investigate protein structure, function, and thermodynamics. MD simulation is also used to investigate high pressure effects on proteins. For conducting better MD simulation under high pressure, the main issues to be addressed are: (i) protein force fields and water models were originally developed to reproduce experimental properties obtained at ambient pressure; and (ii) the timescale to observe the pressure effect is often much longer than that of conventional MD simulations.Scope of reviewFirst, we describe recent developments in MD simulation methodologies for studying the high-pressure structure and dynamics of protein molecules. These developments include force fields for proteins and water molecules, and enhanced simulation techniques. Then, we summarize recent studies of MD simulations of proteins in water under high pressure.Major conclusionsRecent MD simulations of proteins in solution under pressure have reproduced various phenomena identified by experiments using high pressure, such as hydration, water penetration, conformational change, helix stabilization, and molecular stiffening.General significanceMD simulations demonstrate differences in the properties of proteins and water molecules between ambient and high-pressure conditions. Comparing the results obtained by MD calculations with those obtained experimentally could reveal the mechanism by which biological molecular machines work well in collaboration with water molecules.     

Parameterization and Application of an Implicit Solvent Model for Macromolecules

Abstract

As the field of theoretical biophysics begins to recognize systems of longer timescales and larger magnitude, rapid approaches for investigating these systems are required. One promising simplification of the typical system of a solute surrounded by water is the use of implicit solvation models. The generalized Born implicit solvent offers a rapid approach for computing the electrostatic effects of bulk solvent without the explicit representation of water molecules. This report describes the parameterization of a generalized Born (GB) model for protein and nucleic acid structures. As a demonstration of the usefulness of this approach, the GB model is applied toward the discrimination of misfolded and properly folded protein structures. This study attempts to illustrate the potential of the GB model for molecular dynamics simulations over longer timescales as well as for screening large structural databases.     

Transcriptional Bursting Diversifies the Behaviour of a Toggle Switch: Hybrid Simulation of Stochastic Gene Expression

Hybrid models for gene expression combine stochastic and deterministic representations of the underlying biophysical mechanisms. According to one of the simplest hybrid formalisms, protein molecules are produced in randomly occurring bursts of a randomly distributed size while they are degraded deterministically. Here, we use this particular formalism to study two key regulatory motifs—the autoregulation loop and the toggle switch. The distribution of burst times is determined and used as a basis for the development of exact simulation algorithms for gene expression dynamics. For the autoregulation loop, the simulations are compared to an analytic solution of a master equation. Simulations of the toggle switch reveal a number of qualitatively distinct scenarios with implications for the modelling of cell-fate selection.     

The Implementation of the Iterative Diagonalization Scheme and ab initio Molecular Dynamics Simulation with the LAPW Method

Abstract

A new ab initio molecular dynamics method based on the full-potential linearized-augmented-plane-wave (LAPW) basis set has been implemented. The LAPW basis set has been successfully employed for systems containing localized electrons such as first row atoms and transition metals. In our implementation of the LAPW-MD scheme, iterative residual minimization algorithm is used to solve the electronic states problem. The atoms are moved according to forces derived from the Hellman–Feynman theorem and incomplete basis set correction terms. The performance of the program is further enhanced by parallelization. We will discuss technical details of the program implementation and present results obtained from this code to the equilibrium structures and vibrational properties of simple diatomic molecules.     

Large Scale Molecular Dynamics on Parallel Computers using the Link-cell Algorithm

Parallel computers offer a more cost-effective route to high performance computing than traditional single processor machines. Software for such machines is still in its infancy and they are often much more difficult to program than sequential machines. In addition many of the algorithms which are successful with sequential and vector processors are no longer appropriate. Both the force calculation and integration steps of molecular dynamics are parallel in nature and for that reason we have developed a parallel algorithm based on the link cell technique. This method is particularly efficient when the range of intermolecular potential is much smaller than the dimensions of the simulation box. The details of the algorithm are presented for systems of atoms in two and three dimensions using a number of decompositions into sub-units. The algorithm has been tested on an Intel iPSC/2 and a Cray X-MP/416 and the results are presented for simulations of up to 2 · 10⁶ atoms.     

Fluorescence microscopy for visualizing single-molecule protein dynamics

BackgroundSingle-molecule fluorescence imaging (smFI) has evolved into a valuable method used in biophysical and biochemical studies as it can observe the real-time behavior of individual protein molecules, enabling understanding of their detailed dynamic features. smFI is also closely related to other state-of-the-art microscopic methods, optics, and nanomaterials in that smFI and these technologies have developed synergistically.Scope of reviewThis paper provides an overview of the recently developed single-molecule fluorescence microscopy methods, focusing on critical techniques employed in higher-precision measurements in vitro and fluorescent nanodiamond, an emerging promising fluorophore that will improve single-molecule fluorescence microscopy.Major conclusionssmFI will continue to improve regarding the photostability of fluorophores and will develop via combination with other techniques based on nanofabrication, single-molecule manipulation, and so on.General significanceQuantitative, high-resolution single-molecule studies will help establish an understanding of protein dynamics and complex biomolecular systems.     

Investigation of molecular dissolution mechanism of ketoprofen binary and ternary solid dispersions by molecular dynamics simulations

Abstract

The research aimed to investigate the molecular dissolution mechanism of both binary and ternary solid dispersions by molecular dynamics simulations. The simulation results indicated that the drug molecules were much easier to be released from surfactant-containing ternary systems than from binary ones. Moreover, sodium dodecyl sulfate as an additive in ternary systems had better effects than Tween 60. The simulation results were in well agreement with the experimental results. This research presented a reasonable explanation of molecular dissolution mechanism for both binary and ternary solid dispersions, which may benefit the future development of solid dispersion formulations.     

Isobaric and Isothermal Molecular Dynamics Simulations of Diatomic Systems

Abstract

Isobaric molecular dynamics simulations were carried out for diatomic systems using different algorithms available in the literature. Two-centered Lennard-Jones potentials with and without quadrupolar interactions were used. Thermodynamic properties obtained from the isobaric algorithms compared very well with those of an equivalent simulation in the microcanonical ensemble; however, some differences were observed when similar comparisons were carried out for dynamic properties. More specifically, the constant pressure constraint affects the translational dynamics of the system because of the non-negligible differences between the momenta and the instantaneous velocities of the molecules.

Furthermore, the following studies were carried out using isobaric MD simulations: 1. Low temperature spontaneous FCC-orthorhombic (and vice versa) transition of a diatomic system with quadrupolar interactions as a function of the molecular bond length. 2. Effect of quadrupolar interaction on isobaric melting of a model diatomic system. 3. Effect of pressure on melting properties of a model diatomic system with quadrupolar interactions.     

Structural and dynamical correlations in PfHGXPRT oligomers: A molecular dynamics simulation study

PfHGXPRT is a key enzyme involved in purine nucleotide salvage pathway of the malarial parasite, Plasmodium falciparum. Atomistic molecular dynamics simulations have been performed on two types of PfHGXPRT dimers (D1 and D3) and its tetramer in their apo and ligand-bound states. A significant event in the catalytic cycle is the dynamics of a gate that provides access for the ligand molecules to the reaction center. The gate is formed by loops II and IV, the former being the most flexible. Large amplitude conformational changes have been observed in active site loop II. Upon complete occupancy of the active site, loop II gets stabilized due to specific interactions between its residues and the ligand molecules. Remote loop, X, is seen to be less fluxional in the D3 dimer than in D1 which is rationalized as due to the greater number of inter-subunit contacts in the former. The presence of ligand molecules in subunits of the tetramer further reduces the flexibility of loop X epitomizing a communication between this region and the active sites in the tetramer. These observations are in accordance with the outcomes of several experimental investigations. Participation of loop X in the oligomerization process has also been discerned. Between the two types of dimers in solution, D1 tetramerizes readily and thus would not be present as free dimers. We conjecture an equilibrium to exist between D3 and the tetramer in solution; upon binding of the ligand molecules to the D3 dimer, this equilibrium shifts toward the tetramer.     

Dynamics and Stability of GCAA Tetraloops with 2-Aminopurine and Purine Substitutions

Abstract

The contributions of various interactions in the GGCGCAAGCC hairpin containing a GCAA tetraloop were studied by computer simulations using the substitutions of functional groups. The guanosine (G) in the first tetraloop position or in the C-G closing base pair was replaced by 2-aminopurine (AP), and the individual tetraloop's adenosines (A) were replaced by purine (PUR). These substitutions eliminated particular hydrogen bonds thought to stabilize the GCAA tetraloop. For each substitution, molecular dynamics (MD) simulations were carried out in an aqueous solution with sodium counterions, using the CHARMM27 force field. The MD simulations showed that the substitutions in the first (G→AP) and the third (A→PUR) position of the GCAA tetraloop did not significantly influence the conformation of the hairpin. A long-lived bridging water molecule observed in the GCAA loop was present in both modified loops. The substitutions made in the last loop position (A→PUR) or in the C-G base pair closing the tetraloop (G→AP) to some extent influenced the loop structure and dynamics. These loops did not display the long- lived bridging water molecules. When the second A in the GCAA loop was replaced by PUR, the first A in the loop was observed in the anti or in the syn orientation about the gly- cosyl bond. The G to AP substitution in C-G base pair led to a change of their arrangement from the Watson-Crick to wobble. The MD simulations of the hairpin with C-AP wobble closing base pair showed increased conformational dynamics of the hairpin. The changes of hairpin formation free energy associated with the substitutions of individual bases were calculated by the free energy perturbation method. Our theoretical estimates suggest a larger destabilization for the G to AP substitutions in GCAA loop than for the substitutions of individual A's by PUR, which is in accordance with experimental tendency. The calculations predicted a similar free energy change for G to AP substitutions in the GCAA tetraloop and in the C-G closing base pair.