Mercedes-Benz water molecules near hydrophobic wall: integral equation theories vs Monte Carlo simulations.
Ontology highlight
ABSTRACT:
SUBMITTER: Urbic T
PROVIDER: S-EPMC3206897 | biostudies-other | 2011 Oct
REPOSITORIES: biostudies-other
Similar Datasets
Project description:We introduce here a novel Monte Carlo simulation method for studying the interactions of hydrophobic peptides with lipid membranes. Each of the peptide's amino acids is represented as two interaction sites: one corresponding to the backbone alpha-carbon and the other to the side chain, with the membrane represented as a hydrophobic profile. Peptide conformations and locations in the membrane and changes in the membrane width are sampled using the Metropolis criterion, taking into account the underlying energetics. Using this method we investigate the interactions between the hydrophobic peptide M2delta and a model membrane. The simulations show that starting from an extended conformation in the aqueous phase, the peptide first adsorbs onto the membrane surface, while acquiring an ordered helical structure. This is followed by formation of a helical-hairpin and insertion into the membrane. The observed path is in agreement with contemporary understanding of peptide insertion into biological membranes. Two stable orientations of membrane-associated M2delta were obtained: transmembrane (TM) and surface, and the value of the water-to-membrane transfer free energy of each of them is in agreement with calculations and measurements on similar cases. M2delta is most stable in the TM orientation, where it assumes a helical conformation with a tilt of 14 degrees between the helix principal axis and the membrane normal. The peptide conformation agrees well with the experimental data; average root-mean-square deviations of 2.1 A compared to nuclear magnetic resonance structures obtained in detergent micelles and supported lipid bilayers. The average orientation of the peptide in the membrane in the most stable configurations reported here, and in particular the value of the tilt angle, are in excellent agreement with the ones calculated using the continuum-solvent model and the ones observed in the nuclear magnetic resonance studies. This suggests that the method may be used to predict the three-dimensional structure of TM peptides.
| S-EPMC1303652 | biostudies-literature
Project description:A Monte Carlo algorithm was developed to predict the most likely orientations of protein G B1, an immunoglobulin G (IgG) antibody-binding domain of protein G, adsorbed onto a hydrophobic surface. At each Monte Carlo step, the protein was rotated and translated as a rigid body. The assumption about rigidity was supported by quartz crystal microbalance with dissipation monitoring experiments, which indicated that protein G B1 adsorbed on a polystyrene surface with its native structure conserved and showed that its IgG antibody-binding activity was retained. The Monte Carlo simulations predicted that protein G B1 is likely adsorbed onto a hydrophobic surface in two different orientations, characterized as two mutually exclusive sets of amino acids contacting the surface. This was consistent with sum frequency generation (SFG) vibrational spectroscopy results. In fact, theoretical SFG spectra calculated from an equal combination of the two predicted orientations exhibited reasonable agreement with measured spectra of protein G B1 on polystyrene surfaces. Also, in explicit solvent molecular dynamics simulations, protein G B1 maintained its predicted orientation in three out of four runs. This work shows that using a Monte Carlo approach can provide an accurate estimate of a protein orientation on a hydrophobic surface, which complements experimental surface analysis techniques and provides an initial system to study the interaction between a protein and a surface in molecular dynamics simulations.
| S-EPMC5148762 | biostudies-literature
Project description:Monte Carlo simulations and integral equation theory were used to study the thermodynamics and structure of particles interacting through the smooth version of Stell-Hemmer interaction. We checked the possibility that a fluid with a core-softened potential reproduces anomalies of liquid water such as the density anomaly, the minimum in the isothermal compressibility as a function of temperature, and others. Critical points of the fluid were also determined. We showed that a potential with two characteristic distances is sufficient for the system to exhibit water-like behavior and anomalies, including the famous density maximum. We also showed that some versions of the integral equation theory completely fail to predict structure of such system, while others only predict it qualitatively.
| S-EPMC3829925 | biostudies-other
Project description:Using computer simulations, we systematically studied the influence of different design parameters of a spherical nanoparticle tethered with monovalent ligands on its efficiency of targeting planar cell surfaces containing mobile receptors. We investigate how the nanoparticle affinity can be affected by changing the binding energy, the percent of functionalization by ligands, tether length, grafting density, and nanoparticle core size. In general, using a longer tether length or increasing the number of tethered chains without increasing the number of ligands increases the conformational penalty for tether stretching/compression near the cell surface and leads to a decrease in targeting efficiency. At the same time, using longer tethers or a larger core size allows ligands to interact with receptors over a larger cell surface area, which can enhance the nanoparticle affinity toward the cell surface. We also discuss the selectivity of nanoparticle targeting of cells with a high receptor density. Based on the obtained results, we provide recommendations for improving the nanoparticle binding affinity and selectivity, which can guide future nanoparticle development for diagnostic and therapeutic purposes.
| S-EPMC2999362 | biostudies-literature
Project description:This review article is intended as a practical guide for newcomers to the field of kinetic Monte Carlo (KMC) simulations, and specifically to lattice KMC simulations as prevalently used for surface and interface applications. We will provide worked out examples using the kmos code, where we highlight the central approximations made in implementing a KMC model as well as possible pitfalls. This includes the mapping of the problem onto a lattice and the derivation of rate constant expressions for various elementary processes. Example KMC models will be presented within the application areas surface diffusion, crystal growth and heterogeneous catalysis, covering both transient and steady-state kinetics as well as the preparation of various initial states of the system. We highlight the sensitivity of KMC models to the elementary processes included, as well as to possible errors in the rate constants. For catalysis models in particular, a recurrent challenge is the occurrence of processes at very different timescales, e.g., fast diffusion processes and slow chemical reactions. We demonstrate how to overcome this timescale disparity problem using recently developed acceleration algorithms. Finally, we will discuss how to account for lateral interactions between the species adsorbed to the lattice, which can play an important role in all application areas covered here.
| S-EPMC6465329 | biostudies-literature
Project description:Integral equation approaches, based on the Ornstein-Zernike equation, provide a fast way to calculate phase diagrams and thermodynamic properties of systems as opposed to time-consuming and computationally expensive computer simulations. However, when employing integral equations it is necessary to introduce simplifications. The Ornstein-Zernike equation merely relates two unknown functions h(r) and c(r), and another relation (closer) between these two functions is needed. The later function cannot be obtained in a closed form and it is always in some approximations. Various approximations exist with each of its own advantages and disadvantages. In this work we extensively tested hyper-netted chain, Percus-Yevick, Kovalenko-Hirata, and Rogers-Young closure on an interaction model with core-softened potential. Convergence domain was established for each method. We calculated pair distribution functions, pressure, and excess energy. Results were compared with Monte Carlo simulation results and literature data from molecular dynamics simulations.
| S-EPMC3695973 | biostudies-other
Project description:The placebo response in epilepsy randomized clinical trials (RCTs) has recently been shown to largely reflect underlying natural variability in seizure frequency. Based on this observation, we sought to explore the parameter space of RCT design to optimize trial efficiency and cost.We used one of the world's largest patient reported seizure diary databases, SeizureTracker.com to derive virtual patients for simulated RCTs. We ran 1000 randomly generated simulated trials using bootstrapping (sampling with replacement) for each unique combination of trial parameters, sweeping a large set of parameters in durations of the baseline and test periods, number of patients, eligibility criteria, drug effect size, and patient dropout. We studied the resulting trial efficiency and cost.A total of 6,732,000 trials were simulated, drawing from 5097 patients in the database. We found that the strongest regression predictors of placebo response were durations of baseline and test periods. Drug effect size had a major impact on trial efficiency and cost. Dropout did not have a major impact on trial efficiency or cost. Eligibility requirements impacted trial efficiency to a limited extent. Cost was minimized while maintaining statistical integrity with very short RCT durations.This study suggests that RCT parameters can be improved over current practice to reduce costs while maintaining statistical power. In addition, use of a large-scale population dataset in a massively parallel computing analysis allows exploration of the wider parameter space of RCT design prior to running a trial, which could help accelerate drug discovery and approval.
| S-EPMC5553226 | biostudies-other
Project description:Biomacromolecule activity is usually related to its ability to keep a specific structure. However, in solution, many parameters (pH, ionic strength) and external compounds (polyelectrolytes, nanoparticles) can modify biomacromolecule structure as well as acid/base properties, thus resulting in a loss of activity and denaturation. In this paper, the impact of neutral and charged nanoparticles (NPs) is investigated by Monte Carlo simulations on polypeptide (PP) chains with primary structure based on bovine serum albumin. The influence of pH, salt valency, and NP surface charge density is systematically studied. It is found that the PP is extended at extreme pH, when no complex formation is observed, and folded at physiological pH. PP adsorption around oppositely-charged NPs strongly limits chain structural changes and modifies its acid/base properties. At physiological pH, the complex formation occurs only with positively-charged NPs. The presence of salts, in particular those with trivalent cations, introduces additional electrostatic interactions, resulting in a mitigation of the impact of negative NPs. Thus, the corona structure is less dense with locally-desorbed segments. On the contrary, very limited impact of salt cation valency is observed when NPs are positive, due to the absence of competitive effects between multivalent cations and NP.
| S-EPMC6432162 | biostudies-literature
Project description:A simple coarse-grained model of mucus structure and dynamics is proposed and evaluated. The model is based on simple cubic, face-centered lattice representation. Mucins are simulated as lattice chains in which each bead of the model chains represents a mucin domain, equivalent to its Kuhn segment. The remaining lattice sites are considered to be occupied by the solvent. Model mucins consist of three types of domains: polar (glycosylated central segments), hydrophobic, and cysteine-rich, located at the terminal part of the mucin chains. The sequence of these domains mimics the sequence of real mucins. Static and dynamic properties of the system were studied by means of Monte Carlo dynamics. It was shown that the model system undergoes sol-gel transition and that the interactions between hydrophobic domains are responsible for the transition and characteristic properties of the dynamic network in the gel phase. Cysteine-rich domains are essential for frictional properties of the system. Structural and dynamic properties of the model mucus observed in simulations are in qualitative agreement with known experimental facts and provide mechanistic explanation of complex properties of real mucus.
| S-EPMC2998598 | biostudies-literature
Project description:Predictive modeling of reaction equilibria presents one of the grand challenges in the field of molecular simulation. Difficulties in the study of such systems arise from the need (i) to accurately model both strong, short-ranged interactions leading to the formation of chemical bonds and weak interactions arising from the environment, and (ii) to sample the range of time scales involving frequent molecular collisions, slow diffusion, and infrequent reactive events. Here we present a novel reactive first-principles Monte Carlo (RxFPMC) approach that allows for investigation of reaction equilibria without the need to prespecify a set of chemical reactions and their ideal-gas equilibrium constants. We apply RxFPMC to investigate a nitrogen/oxygen mixture at T = 3000 K and p = 30 GPa, i.e., conditions that are present in atmospheric lightning strikes and explosions. The RxFPMC simulations show that the solvation environment leads to a significantly enhanced NO concentration that reaches a maximum when oxygen is present in slight excess. In addition, the RxFPMC simulations indicate the formation of NO2 and N2O in mole fractions approaching 1%, whereas N3 and O3 are not observed. The equilibrium distributions obtained from the RxFPMC simulations agree well with those from a thermochemical computer code parametrized to experimental data.
| S-EPMC4919768 | biostudies-literature