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.
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.
Project description:Quantitative bias analysis can be used to empirically assess how far study estimates are from the truth (i.e., an estimate that is free of bias). These methods can be used to explore the potential impact of confounding bias, selection bias (collider stratification bias), and information bias. Quantitative bias analysis includes methods that can be used to check the robustness of study findings to multiple types of bias and methods that use simulation studies to generate data and understand the hypothetical impact of specific types of bias in a simulated data set. In this article, we review 2 strategies for quantitative bias analysis: 1) traditional probabilistic quantitative bias analysis and 2) quantitative bias analysis with generated data. An important difference between the 2 strategies relates to the type of data (real vs. generated data) used in the analysis. Monte Carlo simulations are used in both approaches, but the simulation process is used for different purposes in each. For both approaches, we outline and describe the steps required to carry out the quantitative bias analysis and also present a bias-analysis tutorial demonstrating how both approaches can be applied in the context of an analysis for selection bias. Our goal is to highlight the utility of quantitative bias analysis for practicing epidemiologists and increase the use of these methods in the epidemiologic 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.
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.
Project description:Human Immunodeficiency Virus 1 (HIV-1) evades adaptive immunity by means of its extremely high mutation rate, which allows the HIV envelope glycoprotein to continuously escape from the action of antibodies. However, some broadly neutralizing antibodies (bNAbs) targeting specific viral regions show the ability to block the infectivity of a large number of viral variants. The discovery of these antibodies opens new avenues in anti-HIV therapy; however, they are still suboptimal tools as their amplitude of action ranges between 50% and 90% of viral variants. In this context, being able to discriminate between sensitive and resistant strains to an antibody would be of great interest for the design of optimal clinical antibody treatments and to engineer potent bNAbs for clinical use. Here, we describe a hierarchical procedure to predict the antibody neutralization efficacy of multiple viral isolates to three well-known anti-CD4bs bNAbs: VRC01, NIH45-46 and 3BNC117. Our method consists of simulating the three-dimensional binding process between the gp120 and the antibody by using Protein Energy Landscape Exploration (PELE), a Monte Carlo stochastic approach. Our results clearly indicate that the binding profiles of sensitive and resistant strains to a bNAb behave differently, showing the latter's weaker binding profiles, that can be exploited for predicting antibody neutralization efficacy in hypermutated HIV-1 strains.
Project description:Thermochemical heat-storage applications, based on the reversible endo-/exothermic hydration reaction of salts, are intensively investigated to search for compact heat-storage devices. To achieve a truly valuable storage system, progressively complex salts are investigated. For these salts, the equilibrium temperature and pressure conditions are not always easy to predict. However, these conditions are crucial for the design of thermochemical heat-storage systems. A biased grand-canonical Monte Carlo (GCMC) tool is developed, enabling the study of equilibrium conditions at the molecular level. The GCMC algorithm is combined with reactive force field molecular dynamics (ReaxFF), which allows bond formation within the simulation. The Weeks-Chandler-Andersen (WCA) potential is used to scan multiple trial positions for the GCMC algorithm at a small cost. The most promising trial positions can be selected for recomputation with the more expensive ReaxFF. The developed WCA-ReaxFF-GCMC tool was used to study the hydration of MgCl2·nH2O. The simulation results show a good agreement with experimental and thermodynamic equilibriums for multiple hydration levels. The hydration shows that water, present at the surface of crystalline salt, deforms the surface layers and promotes further hydration of these deformed layers. Additionally, the WCA-ReaxFF-GCMC algorithm can be used to study other, non-TCM-related, reactive sorption processes.
Project description:Many-body entanglement unveils additional aspects of quantum matter and offers insights into strongly correlated physics. While ground-state entanglement has received much attention in the past decade, the study of mixed-state quantum entanglement using negativity in interacting fermionic systems remains largely unexplored. We demonstrate that the partially transposed density matrix of interacting fermions, similar to their reduced density matrix, can be expressed as a weighted sum of Gaussian states describing free fermions, enabling the calculation of rank-n Rényi negativity within the determinant quantum Monte Carlo framework. We calculate the rank-two Rényi negativity for the half-filled Hubbard model and the spinless t-V model. Our calculation reveals that the area law coefficient of the Rényi negativity for the spinless t-V model has a logarithmic finite-size scaling at the finite-temperature transition point. Our work contributes to the calculation of entanglement and sets the stage for future studies on quantum entanglement in various fermionic many-body mixed states.
Project description:A strategy is presented to implement Gaussian process potentials in molecular simulations through parallel programming. Attention is focused on the three-body nonadditive energy, though all algorithms extend straightforwardly to the additive energy. The method to distribute pairs and triplets between processes is general to all potentials. Results are presented for a simulation box of argon, including full box and atom displacement calculations, which are relevant to Monte Carlo simulation. Data on speed-up are presented for up to 120 processes across four nodes. A 4-fold speed-up is observed over five processes, extending to 20-fold over 40 processes and 30-fold over 120 processes.
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.