scholarly journals Structural and dynamical anomalies of soft particles interacting through harmonic repulsions

2015 ◽  
Vol 17 (8) ◽  
pp. 5734-5742 ◽  
Author(s):  
Wenze Ouyang ◽  
Bin Sun ◽  
Zhiwei Sun ◽  
Shenghua Xu

Molecular dynamics (MD) simulations are carried out to investigate the structural and dynamical anomalies in the core-softened fluid with harmonic repulsions.

Polymers ◽  
2021 ◽  
Vol 13 (4) ◽  
pp. 552
Author(s):  
Andrey O. Kurbatov ◽  
Nikolay K. Balabaev ◽  
Mikhail A. Mazo ◽  
Elena Yu. Kramarenko

We studied the conformational behavior of silicon-containing dendrimers during their adsorption onto a flat impenetrable surface by molecular dynamics (MD) simulations. Four homologous series of dendrimers from the 4th up to the 7th generations were modeled, namely, two types of carbosilane dendrimers differing by the functionality of the core Si atom and two types of siloxane dendrimers with different lengths of the spacers. Comparative analysis of the fractions of adsorbed atoms belonging to various structural layers within dendrimers as well as density profiles allowed us to elucidate not only some general trends but also the effects determined by dendrimer specificity. In particular, it was found that in contrast to the carbosilane dendrimers interacting with the adsorbing surface mainly by their peripheral layers, the siloxane dendrimers with the longer –O–Si(CH3)2–O spacers expose atoms from their interior to the surface spreading out on it. These findings are important for the design of functional materials on the basis of silicon-containing dendrimers.


2019 ◽  
Vol 3 (4) ◽  
pp. 24
Author(s):  
Marie-Anne Hervé du Penhoat ◽  
Anis Hamila ◽  
Marie-Pierre Gaigeot ◽  
Rodolphe Vuilleumier ◽  
Kentaro Fujii ◽  
...  

It has been suggested that core ionization in DNA atoms could induce complex, irreparable damage. Synchrotron soft X-rays have been used to probe the damage induced by such events in thin films of DNA components. In a complementary approach, we investigate the fragmentation dynamics following a carbon or oxygen K-shell ionization in 2-deoxy-D-ribose (DR), a major component in the DNA chain. Core ionization of the sugars hydration layer is also studied. To that aim, we use state-of-the-art ab initio Density Functional Theory-based Molecular Dynamics (MD) simulations. The ultrafast dissociation dynamics of the core ionized molecule, prior Auger decay, is modeled for about 10 fs. We show that the core-ionization of oxygen atoms within DR or its hydration layer may induce proton transfers towards nearby molecules, before Auger decay. In a second step, we model an Auger effect occurring either at the beginning or at the end of the core–hole dynamics. Two electrons are removed from the deepest valence molecular orbitals localized on the initially core-ionized oxygen atom (O*), and this electronic state is propagated by means of Ehrenfest MD. We show an ultrafast dissociation of the DR2+ molecule C-O* bonds, which, in most cases, seems independent of the time at which Auger decay occurs.


2000 ◽  
Vol 653 ◽  
Author(s):  
Celeste Sagui ◽  
Thoma Darden

AbstractFixed and induced point dipoles have been implemented in the Ewald and Particle-Mesh Ewald (PME) formalisms. During molecular dynamics (MD) the induced dipoles can be propagated along with the atomic positions either by interation to self-consistency at each time step, or by a Car-Parrinello (CP) technique using an extended Lagrangian formalism. The use of PME for electrostatics of fixed charges and induced dipoles together with a CP treatment of dipole propagation in MD simulations leads to a cost overhead of only 33% above that of MD simulations using standard PME with fixed charges, allowing the study of polarizability in largemacromolecular systems.


2020 ◽  
Author(s):  
Matías R. Machado ◽  
Sergio Pantano

<p> Despite the relevance of properly setting ionic concentrations in Molecular Dynamics (MD) simulations, methods or practical rules to set ionic strength are scarce and rarely documented. Based on a recently proposed thermodynamics method we provide an accurate rule of thumb to define the electrolytic content in simulation boxes. Extending the use of good practices in setting up MD systems is promptly needed to ensure reproducibility and consistency in molecular simulations.</p>


2019 ◽  
Vol 16 (3) ◽  
pp. 291-300
Author(s):  
Saumya K. Patel ◽  
Mohd Athar ◽  
Prakash C. Jha ◽  
Vijay M. Khedkar ◽  
Yogesh Jasrai ◽  
...  

Background: Combined in-silico and in-vitro approaches were adopted to investigate the antiplasmodial activity of Catharanthus roseus and Tylophora indica plant extracts as well as their isolated components (vinblastine, vincristine and tylophorine). </P><P> Methods: We employed molecular docking to prioritize phytochemicals from a library of 26 compounds against Plasmodium falciparum multidrug-resistance protein 1 (PfMDR1). Furthermore, Molecular Dynamics (MD) simulations were performed for a duration of 10 ns to estimate the dynamical structural integrity of ligand-receptor complexes. </P><P> Results: The retrieved bioactive compounds viz. tylophorine, vinblastin and vincristine were found to exhibit significant interacting behaviour; as validated by in-vitro studies on chloroquine sensitive (3D7) as well as chloroquine resistant (RKL9) strain. Moreover, they also displayed stable trajectory (RMSD, RMSF) and molecular properties with consistent interaction profile in molecular dynamics simulations. </P><P> Conclusion: We anticipate that the retrieved phytochemicals can serve as the potential hits and presented findings would be helpful for the designing of malarial therapeutics.


Polymers ◽  
2020 ◽  
Vol 13 (1) ◽  
pp. 99
Author(s):  
Cristian Privat ◽  
Sergio Madurga ◽  
Francesc Mas ◽  
Jaime Rubio-Martínez

Solvent pH is an important property that defines the protonation state of the amino acids and, therefore, modulates the interactions and the conformational space of the biochemical systems. Generally, this thermodynamic variable is poorly considered in Molecular Dynamics (MD) simulations. Fortunately, this lack has been overcome by means of the Constant pH Molecular Dynamics (CPHMD) methods in the recent decades. Several studies have reported promising results from these approaches that include pH in simulations but focus on the prediction of the effective pKa of the amino acids. In this work, we want to shed some light on the CPHMD method and its implementation in the AMBER suitcase from a conformational point of view. To achieve this goal, we performed CPHMD and conventional MD (CMD) simulations of six protonatable amino acids in a blocked tripeptide structure to compare the conformational sampling and energy distributions of both methods. The results reveal strengths and weaknesses of the CPHMD method in the implementation of AMBER18 version. The change of the protonation state according to the chemical environment is presumably an improvement in the accuracy of the simulations. However, the simulations of the deprotonated forms are not consistent, which is related to an inaccurate assignment of the partial charges of the backbone atoms in the CPHMD residues. Therefore, we recommend the CPHMD methods of AMBER program but pointing out the need to compare structural properties with experimental data to bring reliability to the conformational sampling of the simulations.


Molecules ◽  
2021 ◽  
Vol 26 (6) ◽  
pp. 1711
Author(s):  
Mohamed Ahmed Khaireh ◽  
Marie Angot ◽  
Clara Cilindre ◽  
Gérard Liger-Belair ◽  
David A. Bonhommeau

The diffusion of carbon dioxide (CO2) and ethanol (EtOH) is a fundamental transport process behind the formation and growth of CO2 bubbles in sparkling beverages and the release of organoleptic compounds at the liquid free surface. In the present study, CO2 and EtOH diffusion coefficients are computed from molecular dynamics (MD) simulations and compared with experimental values derived from the Stokes-Einstein (SE) relation on the basis of viscometry experiments and hydrodynamic radii deduced from former nuclear magnetic resonance (NMR) measurements. These diffusion coefficients steadily increase with temperature and decrease as the concentration of ethanol rises. The agreement between theory and experiment is suitable for CO2. Theoretical EtOH diffusion coefficients tend to overestimate slightly experimental values, although the agreement can be improved by changing the hydrodynamic radius used to evaluate experimental diffusion coefficients. This apparent disagreement should not rely on limitations of the MD simulations nor on the approximations made to evaluate theoretical diffusion coefficients. Improvement of the molecular models, as well as additional NMR measurements on sparkling beverages at several temperatures and ethanol concentrations, would help solve this issue.


CrystEngComm ◽  
2021 ◽  
Author(s):  
Andrey Sarikov ◽  
Anna Marzegalli ◽  
Luca Barbisan ◽  
Massimo Zimbone ◽  
Corrado Bongiorno ◽  
...  

In this work, annihilation mechanism of stacking faults (SFs) in epitaxial 3C-SiC layers grown on Si(001) substrates is studied by molecular dynamics (MD) simulations. The evolution of SFs located in...


Sensors ◽  
2021 ◽  
Vol 21 (8) ◽  
pp. 2621
Author(s):  
Seunghwa Yang

Here, we systematically interrogate the effects of grafting single-walled (SWNT) and multi-walled carbon nanotubes (MWNT) to polymer matrices by using molecular dynamics (MD) simulations. We specifically investigate key material properties that include interfacial load transfer, alteration of nanotube properties, and dispersion of nanotubes in the polymer matrix. Simulations are conducted on a periodic unit cell model of the nanocomposite with a straight carbon nanotube and an amorphous polyethylene terephthalate (PET) matrix. For each type of nanotube, either 0%, 1.55%, or 3.1% of the carbon atoms in the outermost nanotubes are covalently grafted onto the carbon atoms of the PET matrix. Stress-strain curves and the elastic moduli of nanotubes and nanocomposites are determined based on the density of covalent grafting. Covalent grafting promotes two rivalling effects with respect to altering nanotube properties, and improvements in interfacial load transfer in the nanocomposites are clearly observed. The enhanced interface enables external loads applied to the nanocomposites to be efficiently transferred to the grafted nanotubes. Covalent functionalization of the nanotube surface with PET molecules can alter the solubility of nanotubes and improve dispersibility. Finally, we discuss the current limitations and challenges in using molecular modelling strategies to accurately predict properties on the nanotube and polymers systems studied here.


2021 ◽  
pp. 1-12
Author(s):  
Haiyan Li ◽  
Zanxia Cao ◽  
Guodong Hu ◽  
Liling Zhao ◽  
Chunling Wang ◽  
...  

BACKGROUND: The ribose-binding protein (RBP) from Escherichia coli is one of the representative structures of periplasmic binding proteins. Binding of ribose at the cleft between two domains causes a conformational change corresponding to a closure of two domains around the ligand. The RBP has been crystallized in the open and closed conformations. OBJECTIVE: With the complex trajectory as a control, our goal was to study the conformation changes induced by the detachment of the ligand, and the results have been revealed from two computational tools, MD simulations and elastic network models. METHODS: Molecular dynamics (MD) simulations were performed to study the conformation changes of RBP starting from the open-apo, closed-holo and closed-apo conformations. RESULTS: The evolution of the domain opening angle θ clearly indicates large structural changes. The simulations indicate that the closed states in the absence of ribose are inclined to transition to the open states and that ribose-free RBP exists in a wide range of conformations. The first three dominant principal motions derived from the closed-apo trajectories, consisting of rotating, bending and twisting motions, account for the major rearrangement of the domains from the closed to the open conformation. CONCLUSIONS: The motions showed a strong one-to-one correspondence with the slowest modes from our previous study of RBP with the anisotropic network model (ANM). The results obtained for RBP contribute to the generalization of robustness for protein domain motion studies using either the ANM or PCA for trajectories obtained from MD.


Sign in / Sign up

Export Citation Format

Share Document