Coarse-grained molecular simulations of large biomolecules

Recently, we saw a dramatic increase in the number of researches that rely on coarse-grained (CG) simulations for large biomolecules. Here, first, we briefly describe recently developed and used CG models for proteins and nucleic acids. Balance between structure-based and physico-chemical terms is a key issue. We also discuss the multiscale algorithms used to derive CG parameters. Next, we comment on the dynamics used in CG simulations with an emphasis on the importance of hydrodynamic interactions. We then discuss the pros and cons of CG simulations. Finally, we overview recent exciting applications of CG simulations. Publicly available tools and software for CG simulations are also summarized.     

High speed atomic force microscopy of biomolecules by image tracking.

An image-tracking procedure for atomic force microscopy is proposed and tested, which allows repeated imaging of the same area without suffering from lateral drift. The drift correction procedure is based on on-line cross-correlation of succeeding images. Using the image-tracking procedure allows zooming in on a small scan area over a long period and thus increases the frame rate inversely proportional to the scan area. Application of the procedure is demonstrated for diffusion of 5.4-kb DNA plasmids. With a scan area of 500 * 500 nm(2), a single plasmid can be imaged for more than 30 min at 4 s per frame, with a drift less than 10 nm. The high temporal resolution allows detailed analysis of the diffusion of DNA molecules. A diffusion coefficient of 30 nm(2)/s is found for most DNA molecules, though many molecules are temporally pinned to the mica surface, restricting diffusion.     

High efficiency Hall effect micro-biosensor platform for detection of magnetically labeled biomolecules

Detection of magnetically labeled biomolecules using micro-Hall biosensors is a promising method for monitoring biomolecular recognition processes. The measurement efficiency of standard systems is limited by the time taken for magnetic beads to reach the sensing area of the Hall devices. Here, micro-current lines were integrated with Hall effect structures to manipulate the position of magnetic beads via field gradients generated by localized currents flowing in the current lines. Beads were accumulated onto the sensor surface within seconds of passing currents through the current lines. Real-time detection of magnetic beads using current lines integrated with Hall biosensors was achieved. These results are promising in establishing Hall biosensor platforms as efficient and inexpensive means of monitoring biomolecular reactions for medical applications.     

Effect of artificial periodicity in simulations of biomolecules under Ewald boundary conditions: a continuum electrostatics study.

Ewald and related methods are nowadays routinely used in explicit-solvent simulations of biomolecules, although they impose an artificial periodicity in systems which are inherently non-periodic. The consequences of this approximation should be assessed, since they may crucially affect the reliability of computer simulations under Ewald boundary conditions. In the present study we use a method based on continuum electrostatics to investigate the nature and magnitude of possible periodicity-induced artifacts on the potentials of mean force for conformational equilibria in biomolecules. Three model systems and pathways are considered: polyalanine oligopeptides (unfolding), a DNA tetranucleotide (separation of the strands), and the protein Sac7d (conformations from a molecular dynamics simulation). Artificial periodicity may significantly affect these conformational equilibria, in each case stabilizing the most compact conformation of the biomolecule. Three factors enhance periodicity-induced artifacts: (i) a solvent of low dielectric permittivity; (ii) a solute size which is non-negligible compared to the size of the unit cell; and (iii) a non-neutral solute. Neither the neutrality of the solute nor the absence of charge pairs at distances exceeding half the edge of the unit cell do guarantee the absence of artifacts.     

High sensitivity detection of molecular recognition using magnetically labelled biomolecules and magnetoresistive sensors

Small magnetoresistive spin valve sensors (2 x 6 microm(2)) were used to detect the binding of single streptavidin functionalized 2 microm magnetic microspheres to a biotinylated sensor surface. The sensor signals, using 8 mA sense current, were in the order of 150-400 microV for a single microsphere depending on sensor sensitivity and the thickness of the passivation layer over the sensor surface. Sensor saturation signals were 1-2 mV representing an estimated 6-20 microspheres, with a noise level of approximately 10 microV. The detection of biomolecular recognition for the streptavidin-biotin model was shown using both single and differential sensor architectures. The signal data compares favourably with previously reported signals for high numbers of magnetic microspheres detected using larger multilayered giant magnetoresistance sensors. A wide range of applications is foreseen for this system in the development of biochips, high sensitivity biosensors and the detection of single molecules and single molecule interactions.     

High temperature unfolding simulations of the TRPZ1 peptide

We report high temperature molecular dynamics simulations of the unfolding of the TRPZ1 peptide using an explicit model for the solvent. The system has been simulated for a total of 6 μs with 100-ns minimal continuous stretches of trajectory. The populated states along the simulations are identified by monitoring multiple observables, probing both the structure and the flexibility of the conformations. Several unfolding and refolding transition pathways are sampled and analyzed. The unfolding process of the peptide occurs in two steps because of the accumulation of a metastable on-pathway intermediate state stabilized by two native backbone hydrogen bonds assisted by nonnative hydrophobic interactions between the tryptophan side chains. Analysis of the un/folding kinetics and classical commitment probability calculations on the conformations extracted from the transition pathways show that the rate-limiting step for unfolding is the disruption of the ordered native hydrophobic packing (Trp-zip motif) leading from the native to the intermediate state. But, the speed of the folding process is mainly determined by the transition from the completely unfolded state to the intermediate and specifically by the closure of the hairpin loop driven by formation of two native backbone hydrogen bonds and hydrophobic contacts between tryptophan residues. The temperature dependence of the unfolding time provides an estimate of the unfolding activation enthalpy that is in agreement with experiments. The unfolding time extrapolated to room temperature is in agreement with the experimental data as well, thus providing a further validation to the analysis reported here.     

Coarse-grained simulations of the salt dependence of the radius of gyration of polyelectrolytes as models for biomolecules in aqueous solution

The salt dependent radius of gyration of a polyelectrolyte in aqueous solution is calculated in an environment where the polyelectrolyte is surrounded by a permeable membrane that exchanges only solvent particles with the bulk. We obtain additionally the scaling exponent of the gyration radius as a function of the polymerization degree, and find that the polyelectrolyte retains a stretched conformation during the condensation and re-expansion process, indicating that these effects are of an electrostatic nature. The solvent quality is also shown to affect the polyelectrolyte conformation, especially for the poor solvent case. These results are obtained using a hybridized Monte Carlo technique with the coarse-grained, dissipative particle dynamics method with fluctuating number of solvent particles. The full range of the electrostatic interactions is included in the simulations, using the Ewald sum method, and the counterions and solvent molecules are included explicitly. In the complex systems mentioned above, the electrostatic interactions and the solvent quality play a key role in understanding phenomena that do not occur in uncharged systems. Our results are compared and validated with the behavior of some biomolecules under similar environments.     

Comparative analyses of biomolecules

Systematic examination of the gene encoding CYP2B6, a human cytochrome P450, has characterized genetic polymorphisms that might account for its variability in expression and function between individuals.     

High efficient sedimentary basin simulations on hybrid CPU-GPU clusters

The key to achieving high performance on a GPU-enhanced cluster is efficient exploitation of each GPU’s powerful computing capability. Moreover, rationally balancing the workload between CPUs and GPUs can release additional computing power, which arises from the CPUs. In this paper, we extend our earlier work on using a hybrid CPU-GPU cluster for real-world sedimentary basin simulation, by further improving the involved CUDA implementations. A thorough analysis of the achieved new performance is also carried out. By using 1024 GPUs and 12288 CPU cores together, our best CPU-GPU hybrid implementation is able to achieve a double-precision performance of 72.8 TFlops, in connection with simulations on a huge 131072×131072 mesh.     

High pressure modulates amyloid formation

A common mechanism of conformational changes and pathological aggregation of proteins associated with amyloid diseases remains to be proven. High pressure is emerging as a new strategy for studying aspects of amyloid formation. Pressure provides a convenient means to populate and characterize partially folded states, which are thought to have a key role in assembly processes of proteins into amyloid fibrils. High pressure can also be used to dissociate aggregates and amyloid fibrils or on the opposite to generate such species.     

High pressure studies on transthyretin

High hydrostatic pressure (HHP) is a powerful tool to study protein folding and the dynamics and structure of folding intermediates. Aggregates and amyloids, derived from partially folding intermediates at the junction between productive and off-pathway folding, have been studied as well, which promises better understanding of the protein misfolding diseases. Here is summarized the recent data we have collected with transthyretin under pressure.     

Pressure-temperature phase diagrams of biomolecules

The pressure-temperature phase diagram of various biomolecules is reviewed. Special attention is focused on the elliptic phase diagram of proteins. The phenomenological thermodynamic theory describing this diagram explains the heat, cold and pressure denaturations in a unified picture. The limitations and possible developments of this theory are discussed as well. It is pointed out that a more complex diagram can be obtained when the intermolecular interactions are also taken into account. In this case metastable states appear on the pressure-temperature (p-T) diagram due to intermolecular interactions. Pressure-temperature phase diagrams of other biopolymers are also discussed. While the p-T diagrams of helix-coil transition of nucleic acids and of gel-liquid crystal transition of lipid bilayers are non-elliptical, those of gelatinization of starch and of phase separation of some synthetic polymers show an elliptic profile, similar to that of proteins. Finally, the p-T diagram of bacterial inactivation is shown to be elliptic. From the point of view of basic science, this fact shows that the key factor of inactivation should be the protein type, and from the viewpoint of practical applications, it serves as the theoretical basis of pressure treatment of biosystems.     

Single molecule nanomanipulation of biomolecules

The development of nanomanipulation techniques has given investigators the ability to manipulate single biomolecules and to record mechanical events of biomolecules at the single molecule level. The techniques were developed to elucidate the mechanism of molecular motors. We can directly monitor the unitary process of the mechanical work and the energy conversion processes by combining these techniques with the single molecule imaging techniques. Our results strongly suggest that the sliding movement of the actomyosin motor is driven by Brownian movement. Other groups have reported data that are more consistent with the lever arm model. These methods and imaging techniques enable us to monitor the behavior of biomolecules at work and will be applied to other molecular machines.