肌球蛋白Ⅵ定向运动的蒙特卡罗模拟

根据myosin Ⅵ分子马达头部核苷的变化建立一个四态循环,然后利用蒙特卡罗方法模拟分子马达的动力学性质,讨论了myosin Ⅵ速度、滞留时间与溶液中ATP浓度、ADP浓度和负载力F的关系。模拟得到：myosin Ⅵ会在轨道微丝上进行梯跳运动并且每次在微丝上的滞留时间并不相同;另外,myosin Ⅵ的速度会随着ATP浓度的增加而增大,随着ADP浓度和负载力的增加而减小;负载力对myosin Ⅵ在胞内执行输运还是连接功能起一定的调节作用。     

Monte carlo simulation of 3-D buffered Ca(2+) diffusion in neuroendocrine cells

Buffered Ca(2+) diffusion in the cytosol of neuroendocrine cells is a plausible explanation for the slowness and latency in the secretion of hormones. We have developed a Monte Carlo simulation to treat the problem of 3-D diffusion and kinetic reactions of ions and buffers. The 3-D diffusion is modeled as a random walk process that follows the path of each ion and buffer molecule, combined locally with a stochastic treatment of the first-order kinetic reactions involved. Such modeling is able to predict [Ca(2+)] and buffer concentration time courses regardless of how low the calcium influx is, and it is therefore a convenient method for dealing with physiological calcium currents and concentrations. We study the effects of the diffusional and kinetic parameters of the model on the concentration time courses as well as on the local equilibrium of buffers with calcium. An in-mobile and fast endogenous buffer as described by, Biophys. J. 72:674-690) was able to reach local equilibrium with calcium; however, the exogenous buffers considered are displaced drastically from equilibrium at the start of the calcium pulse, particularly below the pores. The versatility of the method also allows the effect of different arrangements of calcium channels on submembrane gradients to be studied, including random distribution of calcium channels and channel clusters. The simulation shows how the particular distribution of channels or clusters can be of relevance for secretion in the case where the distribution of release granules is correlated with the channels or clusters.     

Monte carlo simulation of an ecosystem: A matching between two levels of observation

A Montes Carlo computer simulation is designed in order to described the evolution of an ecosystem under imposed environmental conditions of energy. Some preliminary results of a simple situation show a suggestive matching between the detailed random behaviors of individuals and a global deterministic model. This fact may provide a new insight into the understanding and control of biomass production processes and represents a new tool for the evaluation of the reliability of the classical treatments.     

Combined monte carlo and molecular dynamics simulation of fully hydrated dioleyl and palmitoyl-oleyl phosphatidylcholine lipid bilayers

We have applied a new equilibration procedure for the atomic level simulation of a hydrated lipid bilayer to hydrated bilayers of dioleyl-phosphatidylcholine (DOPC) and palmitoyl-oleyl phosphatidylcholine (POPC). The procedure consists of alternating molecular dynamics trajectory calculations in a constant surface tension and temperature ensemble with configurational bias Monte Carlo moves to different regions of the configuration space of the bilayer in a constant volume and temperature ensemble. The procedure is applied to bilayers of 128 molecules of POPC with 4628 water molecules, and 128 molecules of DOPC with 4825 water molecules. Progress toward equilibration is almost three times as fast in central processing unit (CPU) time compared with a purely molecular dynamics (MD) simulation. Equilibration is complete, as judged by the lack of energy drift in 200-ps runs of continuous MD. After the equilibrium state was reached, as determined by agreement between the simulation volume per lipid molecule with experiment, continuous MD was run in an ensemble in which the lateral area was restrained to fluctuate about a mean value and a pressure of 1 atm applied normal to the bilayer surface. Three separate continuous MD runs, 200 ps in duration each, separated by 10,000 CBMC steps, were carried out for each system. Properties of the systems were calculated and averaged over the three separate runs. Results of the simulations are presented and compared with experimental data and with other recent simulations of POPC and DOPC. Analysis of the hydration environment in the headgroups supports a mechanism by which unsaturation contributes to reduced transition temperatures. In this view, the relatively horizontal orientation of the unsaturated bond increases the area per lipid, resulting in increased water penetration between the headgroups. As a result the headgroup-headgroup interactions are attenuated and shielded, and this contributes to the lowered transition temperature.     

Monte carlo simulations of protein-like heteropolymers.

Properties of a simple model of polypeptide chains were studied by the means of the Monte Carlo method. The chains were built on the (310) hybrid lattice. The residues interacted with long-range potential. There were two kinds of residues: hydrophobic and hydrophilic forming a typical helical pattern -HHPPHPP-. Short range potential was used to prefer helical conformations of the chain. It was found that at low temperatures the model chain formes dense and partially ordered structures (non-unique). The presence of the local potential led to an increase of helicity. The effect of the interplay between the two potentials was studied. After the collapse of the chain further annealing caused rearrangement of helical structures. Dynamic properties of the chain at low temperature depended strongly on the local chain ordering.     

Monte carlo study of the percolation in two-dimensional polymer systems

The structure of a two-dimensional film formed by adsorbed polymer chains was studied by means of Monte Carlo simulations. The polymer chains were represented by linear sequences of lattice beads and positions of these beads were restricted to vertices of a two-dimensional square lattice. Two different Monte Carlo methods were employed to determine the properties of the model system. The first was the random sequential adsorption (RSA) and the second one was based on Monte Carlo simulations with a Verdier-Stockmayer sampling algorithm. The methodology concerning the determination of the percolation thresholds for an infinite chain system was discussed. The influence of the chain length on both thresholds was presented and discussed. It was shown that the RSA method gave considerably lower thresholds for longer chains. This behavior can be explained by a different pool of chain conformations used in the calculations in both methods under consideration.

Monte carlo simulation of bacterial disinfection: Nonlinear and time‐explicit intensity of transition

Bacteria being disinfected in fluid media are discrete entities and mesoscopic in size; moreover, they are incessantly as well as irregularly in motion and in collision among themselves or with the surrounding solid surfaces. As such, it is highly likely that some of the attributes of the bacterial population, for example, their number concentration, will fluctuate randomly. This is especially the case at the tail‐end of disinfection when the population of bacteria is sparse. It might be effectual, therefore, to explore the resultant random fluctuations via a stochastic paradigm. Proposed herein is a Markovian stochastic model for the rate of bacterial disinfection, whose intensity of transition takes into account the contact time of the bacteria with the disinfecting agent to eliminate any given percentage of the bacteria in terms of a nonlinear function of time. The model's master equation has been simulated by resorting to the Monte Carlo method to circumvent the undue complexities in solving it analytically or numerically via conventional numerical techniques. For illustration, the mean, the variance (standard deviation), and the coefficient of variation of the number concentration of bacteria during disinfection have been estimated through Monte Carlo simulation. The results of simulation compare favorably with the available experimental data as well as with those computed from the corresponding deterministic model. © 2010 American Institute of Chemical Engineers Biotechnol. Prog., 2010     

Melting of the solvent structure around a RNA duplex: a molecular dynamics simulation study

From three 2.4-ns molecular dynamics simulations of the r(CpG)(12) duplex conducted at 5, 25 and 37 degrees C, a strong temperature dependence of the dynamics of the water molecules and ions located in the first nucleic acid coordination shell is observed. At 5 degrees C, the highest residence times of bound water molecules exceed 1 ns while, at 37 degrees C, they decrease to 0.5 ns in agreement with available NMR data. Similar temperature dependencies are observed for the potassium ions bound to the duplex. In this temperature range, the structure of the RNA helix remains essentially unchanged. Thus, the observed alterations correspond to a 'premelting' of the solvent structure around the duplex. It is proposed that, before the nucleic acid structure melts, the entropy of the solvent increases to a point where it is no longer compensated by the enthalpic contribution of solute-solute and solute-solvent interactions. At this stage, the weakest structural elements start to melt. In other terms, the experimentally observed melting processes are preceded by a melting of the more labile solvent structure.     

Monte Carlo simulation of myoglobin primary to helical structure

A Monte Carlo simulation is presented of the formation of the individual helices of myoglobin from their primary to their helical structures. A simplified model in which each amino acid residue is replaced by a single interaction center is used. The small helices formed are in good agreement with experiment, while the larger helices are moderately well reproduced.     

Combined Monte Carlo and molecular dynamics simulation of hydrated lipid-cholesterol lipid bilayers at low cholesterol concentration

We have applied a hybrid equilibration and sampling procedure for the atomic level simulation of a hydrated lipid bilayer to systems consisting of dipalmitoyl phosphatidylcholine (DPPC) and cholesterol, and palmitoyl-oleyl phosphatidylcholine (POPC) at low (approximately 6%) cholesterol concentration. The procedure is applied to bilayers of 94 molecules of DPPC, 6 molecules of cholesterol, and 3205 water molecules, and to bilayers of 120 molecules of POPC, 8 molecules of cholesterol, and 4268 water molecules, at a temperature of 325 K. After equilibration, three separate 400-ps continuous molecular dynamics runs, separated by 10,000 configurational bias Monte Carlo steps, were carried out for each system. Properties of the systems were calculated and averaged over the three separate runs. Results of the simulations are presented and compared with experimental data and with other recent simulations of DPPC and cholesterol, and of pure DPPC, and pure POPC. Certain properties of the bilayers are indistinguishable from cholesterol-free bilayers, including lateral diffusion and electron density. Other properties, most notably the order parameter profile, show the effect of cholesterol even at low concentrations.     

Molecular simulation of the effects of alcohols on peptide structure

The effects of alcohols on local protein structure have been simulated using computational approaches and model peptides. Molecular simulations were carried out on a 7-residue peptide created in both an extended conformation and an alpha-helix to explore alcohol-induced changes in peptide structure. It was assumed that alcohols hydrogen bond at peptide carbonyl groups with an optimum geometry and compete with water molecules at these site. Energy minimization of the peptide/alcohol assemblies revealed that alcohols induced a twist in the peptide backbone as a function of (1) the methylene chain length, (2) the hydrogen-bond geometry, (3) halogenation of the molecule, (4) concentration, and (5) the dielectric constant. The rank ordering of the potencies of the alcohols was hexafluoroisopropanol > trifluoroethanol approximately pentanol > butanol > ethanol > methanol. Helix destabilization by cosolvent was measured by examining the hydrogen-bond lengths in peptide structures that resulted from a combination of energy minimization and molecular dynamics simulations. Destabilization was also found to be dependent upon the chemical nature of the alcohol and the hydrogen-bond geometry. The data suggest that alcohols at low concentrations affect protein structure mainly through a combination of hydrogen-bonding and hydrophobic interactions that are influenced by the properties of the solvent.     

Monte carlo simulations of tau proteins: effect of phosphorylation

We perform Monte Carlo simulations of tau proteins bound to a cylinder that mimics a microtubule (MT), and then study them in solution. Tau protein binds to a highly anionic MT surface to stabilize the cylindrical structure of MT. The negatively charged tail domain floats away from the anionic MT surface while positively charged tau segments localize near the MT surface. Monte Carlo simulations demonstrate that, in 3RS tau isoform (which has three imperfect repeats (R) short (S) isoform), amino acids are more condensed near a highly charged interface compared to 4RL isoform (which has four imperfect repeats (R) long (L) isoform). In 4RL isoform, amino acids in tail domain stay mostly apart from the MT surface. In the bulk solution, dephosphorylated taus are separated due to Coulomb repulsion between similarly charged isoforms. Moderate phosphorylation of 3RS isoform decreases average intermolecular distance between dephosphorylated and phosphorylated taus and lead to their overlap. Further phosphorylation does not change noticeably the intermolecular distances.     

Monte carlo simulation-based algorithms for analysis of shotgun proteomic data

Two new statistical models based on Monte Carlo Simulation (MCS) have been developed to score peptide matches in shotgun proteomic data and incorporated in a database search program, MassMatrix (www.massmatrix.net). The first model evaluates peptide matches based on the total abundance of matched peaks in the experimental spectra. The second model evaluates amino acid residue tags within MS/MS spectra. The two models provide complementary scores for peptide matches that result in higher confidence in peptide identification when significant scores are returned from both models. The MCS-based models use a variance reduction technique that improves estimation precision. Due to the high computational expense of MCS-based models, peptide matches were prefiltered by other statistical models before further evaluation by the MCS-based models. Receiver operating characteristic analysis of the data sets confirmed that MCS-based models improved the overall performance of the MassMatrix search software, especially for low-mass accuracy data sets.     

Refinement of the thrombin-bound structure of a hirudin peptide by a restrained electrostatically driven Monte Carlo method.

Energy refinement of the structure of a linear peptide, hirudin56-65, bound to thrombin was carried out using a conformational search method in combination with restrained minimization. Five conformations originated from nmr data and distance geometry calculations having a similar global folding pattern but quite different backbone conformations were used as the starting structures. As a result of this approach, a series of low-energy conformations compatible with a set of upper and lower bounds of interproton distances determined from transferred nuclear Overhauser effects were found. A comparison among the lowest energy conformations of each run showed that the combination of energy refinement plus distance constraints led to a very well-defined structure for both the backbone and the side chains of the last 7 residues of the polypeptide. Furthermore, the low-energy conformations generated with this technique contain a segment of 3(10)-helix involving the last 5 residues at the COOH terminal end.     

A computer simulation study of the hydrated proton in a synthetic proton channel

Classical molecular dynamics simulations using the multistate empirical valence bond model for aqueous proton transport were performed to characterize the hydration structure of an excess proton inside a leucine-serine synthetic ion channel, LS2. For such a nonuniform pore size ion channel, it is found that the Zundel ion (H(5)O(2)(+)) solvation structure is generally more stable in narrow channel regions than in wider channel regions, which is in agreement with a recent study on idealized hydrophobic proton channels. However, considerable diversity in the relative stability of the Zundel to Eigen cation (H(9)O(4)(+)) was observed. Three of the five wide channel regions, one located at the channel's center and the other two located near the channel mouths, are found to show extraordinary preference for the Eigen solvation structure. This implies that proton hopping is inhibited in these regions and therefore suggests that these regions may behave as barriers in the proton conducting pathway inside the channel. The proton solvation is also greatly influenced by the local molecular environment of the protein. In particular, the polar side chains of the Ser residues, which are intimately involved in the solvation structure, can greatly influence proton solvation. However, no preference of the influence by the various Ser side chains was found; they can either promote or prevent the formation of certain solvation structures.     

Molecular dynamics simulation of trp-aporepressor in a solvent.

Molecular dynamics simulations of Escherichia coli trp-aporepressor were carried out in the absence and presence of explicit water molecules. The vacuum simulations resulted in significant deformation of the initial X-ray structure. A solvated simulation with a nonbonded cut-off radius of 9 A gave a better result, and the most satisfactory result was obtained when electrostatic interactions within a cut-off radius of 18 A were considered by a twin-range method. The trajectory from the last simulation was used to analyze the dynamical properties of the aporepressor. The root-mean-square fluctuations of the residues showed the rigidity of the central core and the flexibility of the DNA-binding sites, consistent with the X-ray temperature factors. The dynamical cross-correlation map indicated a significant negative correlation between the central core and the two DNA-binding sites, and thus reproduced the three-domain format (a central core and two DNA-binding heads) from a dynamical point of view. The core region showed weak, but many, intra- and inter-molecular correlations, while the helix-turn-helix DNA-binding motifs were free from correlations with other regions.     

Effects of solvent on the structure of the Alzheimer amyloid-beta(25-35) peptide

The free energy landscape for folding of the Alzheimer's amyloid-beta(25-35) peptide is explored using replica exchange molecular dynamics in both pure water and in HFIP/water cosolvent. This amphiphilic peptide is a natural by-product of the Alzheimer's amyloid-beta(1-40) peptide and retains the toxicity of its full-length counterpart as well as the ability to aggregate into beta-sheet-rich fibrils. Our simulations reveal that the peptide preferentially populates a helical structure in apolar organic solvent, while in pure water, the peptide adopts collapsed coil conformations and to a lesser extent beta-hairpin conformations. The beta-hairpin is characterized by a type II' beta-turn involving residues G29 and A30 and two short beta-strands involving residues N27, K28, I31, and I32. The hairpin is stabilized by backbone hydrogen-bonding interactions between residues K28 and I31; S26 and G33; and by side-chain-to-side-chain interactions between N27 and I32. Implications regarding the mechanism of aggregation of this peptide into fibrils and the role of the environment in modulating secondary structure are discussed.