DFT and MD study of adsorption sensitivity of aluminium phosphide nanotube towards some air pollutant gas molecules

To investigate the adsorption behaviour of CS₂, CO₂, SO₂, H₂Se and H₂S gas molecules on the external surface of (6, 0) single-walled aluminium phosphide nanotube (AlPNT), the density functional theory (DFT) calculations at the B3LYP level of theory are performed. The partial densities of states (PDOS) for the SO₂ molecule, the S and O atoms of SO₂ molecule before and after adsorption on the surface of AlPNT have been plotted. The vibrational frequencies and physical properties such as chemical potential, chemical hardness, dipole moment and chemical electrophilicity of all studied complexes have been systematically investigated. The electron density and the Laplacian of the electron density for bond critical points have been examined by the AIM theory. Also the molecular dynamics (MD) simulations of two complexes with the minimum and maximum negative interaction energies that is: AlPNT/CO₂ and AlPNT/SO₂ complexes, respectively, have been considered.     

Theoretical study of the surface properties of poly(dimethylsiloxane) and poly(tetrafluoroethylene)

Molecular dynamics (MD) simulations of poly(dimethylsiloxane) (PDMS) and poly(tetrafluoroethylene) (PTFE) were carried out to determine their surface properties and energies. This study helps to gain better insight into the molecular modeling of PDMS and PTFE, in particular how different approaches affect calculations of surface energy. Current experimental and theoretical data were used to further understand the surface properties of PDMS and PTFE as well as to validate and verify results obtained from the combination of density functional theory (DFT) calculations (including periodic boundary conditions) and MD simulations. Detailed analysis of the structure and electronic properties (by calculation of the projected density of states) of the bulk and surface models of PDMS and PTFE was performed. The sensitivity of the surface energy calculation of these two polymers to the chemistry and model preparation was indicated. The balance between the molecular density, weight (which also reflects bond orientation in the surface region), bond flexibility, and intramolecular interactions including bond stretching was revealed to govern the results obtained. In modeling, the structural organization of polymer near a given surface (types and number of end groups and broken bonds due to application of different cut offs of the periodic structure) also significantly affects the final results. Besides the structural differences, certain simulation parameters, such the DFT functionals and simulation boxes utilized, play an important role in determining surface energy. The models used here were shown to be sufficient due to their good agreement with experimental and other theoretical data related to surface properties and surface energies.     

Viral structural transition mechanisms revealed by multiscale molecular dynamics/order parameter extrapolation simulation

On the basis of an all‐atom multiscale analysis theory of nanosystem dynamics, a multiscale molecular dynamics/order parameter extrapolation (MD/OPX) approach has recently been developed. It accelerates MD for long‐time simulation of large bionanosystems and addresses rapid atomistic fluctuations and slowly varying coherent dynamics simultaneously. In this study, MD/OPX is optimized and implemented to simulate viral capsid structural transitions. Specifically, 200 ns MD/OPX simulation of the swollen state of cowpea chlorotic mottle virus capsid reveals that it undergoes significant energy‐driven shrinkage in vacuum, which is a symmetry‐breaking process involving local initiation and front propagation. © 2009 Wiley Periodicals, Inc. Biopolymers 93: 61–73, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Modeling and simulation of chemomechanics at the cell-matrix interface

Chemomechanical characteristics of the extracellular materials with which cells interact can have a profound impact on cell adhesion and migration. To understand and modulate such complex multiscale processes, a detailed understanding of the feedback between a cell and the adjacent microenvironment is crucial. Here, we use computational modeling and simulation to examine the cell-matrix interaction at both the molecular and continuum lengthscales. Using steered molecular dynamics, we consider how extracellular matrix (ECM) stiffness and extracellular pH influence the interaction between cell surface adhesion receptors and extracellular matrix ligands, and we predict potential consequences for focal adhesion formation and dissolution. Using continuum level finite element simulations and analytical methods to model cell-induced ECM deformation as a function of ECM stiffness and thickness, we consider the implications toward design of synthetic substrata for cell biology experiments that intend to decouple chemical and mechanical cues.Key words: cell adhesion, focal adhesion, steered molecular dynamics, finite element, chemomechanics, multiscale modeling, elasticity theory     

Multiscale modeling of cellular actin filaments: from atomistic molecular to coarse-grained dynamics

In this article, we present a computational multiscale model for the characterization of subcellular proteins. The model is encoded inside a simulation tool that builds coarse-grained (CG) force fields from atomistic simulations. Equilibrium molecular dynamics simulations on an all-atom model of the actin filament are performed. Then, using the statistical distribution of the distances between pairs of selected groups of atoms at the output of the MD simulations, the force field is parameterized using the Boltzmann inversion approach. This CG force field is further used to characterize the dynamics of the protein via Brownian dynamics simulations. This combination of methods into a single computational tool flow enables the simulation of actin filaments with length up to 400 nm, extending the time and length scales compared to state-of-the-art approaches. Moreover, the proposed multiscale modeling approach allows to investigate the relationship between atomistic structure and changes on the overall dynamics and mechanics of the filament and can be easily (i) extended to the characterization of other subcellular structures and (ii) used to investigate the cellular effects of molecular alterations due to pathological conditions.     

Assessment of the chitosan-functionalized graphene oxide as a carrier for loading thioguanine,an antitumor drug and effect of urea on adsorption process: Combination of DFT computational and molecular dynamics simulation studies

In this study, the interaction thioguanine (TG) anticancer drug with the functionalized graphene oxide (GO) nanosheet surface is theoretically studied in both gas phase and separately in physiological media using the density functional theory (DFT) calculations. DFT calculations indicated the adsorption and solvation energies are negative for f-GONS/TG complexes which propose the adsorption process of TG molecule onto the f-GONS surface is possible from the energetic viewpoint. QTAIM calculations confirm the nature of partially covalent-partially electrostatic between drug and nanosheet. These results are sorely relevant that an approach for loading of TG molecule is the chemical modification of GO using covalent functionalization which can serve as a nanocarrier to load drug molecules. Moreover, to understand the effect of urea on the nature of the interaction between TG and f-GONS, molecular dynamics (MD) simulation was employed. The results indicated that in the presence of urea the adsorption process gets affected and leads to instability of system, while the affinity of the TG for adsorption onto GO surface is increased in pure water.

Communicated by Ramaswamy H. Sarma     

Surface structures of oligoglycines: A molecular dynamics simulation

An atomistic molecular dynamics (MD) simulation of the adsorption of biantennary oligoglycine [H-Gly₄-NH(CH₂)₅]₂ onto a graphite and mica surface is described. The structure of the resultant adsorption layers is analyzed. The secondary structure motifs of peptide blocks are studied by the STRIDE (structural identification) and DSSP (dictionary of the secondary structure of proteins) methods. The results of the study confirm the possibility of forming a polyglycine-II (PGII) structure in the monolayers of biantennary oligoglycines (BAOG) on a graphite surface previously supposed from the data of atomic force microscopy.     

Multiscale modeling and mechanics of filamentous actin cytoskeleton

The adaptive structure and functional changes of the actin cytoskeleton are induced by its mechanical behavior at various temporal and spatial scales. In particular, the mechanical behaviors at different scales play important roles in the mechanical functions of various cells, and these multiscale phenomena require clarification. To establish a milestone toward achieving multiscale modeling and simulation, this paper reviews mathematical analyses and simulation methods applied to the mechanics of the filamentous actin cytoskeleton. The actin cytoskeleton demonstrates characteristic behaviors at every temporal and spatial scale, and mathematical models and simulation methods can be applied to each level of actin cytoskeletal structure ranging from the molecular to the network level. This paper considers studies on mathematical models and simulation methods based on the molecular dynamics, coarse-graining, and continuum dynamics approaches. Every temporal and spatial scale of actin cytoskeletal structure is considered, and it is expected that discrete and continuum dynamics ranging from functional expression at the molecular level to macroscopic functional expression at the whole cell level will be developed and applied to multiscale modeling and simulation.     

Multiscale modeling of lipids and lipid bilayers

A multiscale modeling approach is applied for simulations of lipids and lipid assemblies on mesoscale. First, molecular dynamics simulation of initially disordered system of lipid molecules in water within all-atomic model was carried out. On the next stage, structural data obtained from the molecular dynamics (MD) simulation were used to build a coarse-grained (ten sites) lipid model, with effective interaction potentials computed by the inverse Monte Carlo method. Finally, several simulations of the coarse-grained model on longer length- and time-scale were performed, both within Monte Carlo and molecular dynamics simulations: a periodical sample of lipid molecules ordered in bilayer, a free sheet of such bilayer without periodic boundary conditions, formation of vesicle from a plain membrane, process of self-assembly of lipids randomly dispersed in volume. It was shown that the coarse-grained model, developed exclusively from all-atomic simulation data, reproduces well all the basic features of lipids in water solution.     

Interactions of calcium with the external surfaces of fullerenes and endofullerenes doped with radioactive sodium iodide

We report first-principles calculations carried out to analyze the adsorption of calcium on the outer surface of the fullerene C₆₀, yielding [C₆₀?+?mCa]. Geometric optimization (GO) and molecular dynamics (MD) simulation were performed using the plane-wave pseudopotential method within the framework of density functional theory (DFT) and time-dependent DFT (TD-DFT) to investigate the configurations, the associated energies in the ground state, and the stabilities of fullerenes and endofullerenes doped with radioactive sodium iodide when they interact with calcium atoms on the outer fullerene surface (i.e., [nNa¹³¹I@C₆₀?+?mCa]). The reason for investigating these calcium-functionalized (endo)fullerene systems was to gauge their potential stability when used as vectors to deliver radioiodine to cancerous tissue in the human body. In the simulations, we found that the geometric limit on the number of calcium atoms that can be physisorbed on the outer surface of an empty fullerene while maintaining its structural stability is 28 calcium atoms, which also takes into account the proportional expansion of the fullerene as the number of absorbed calcium atoms increases. However, the stability of a fullerene system during calcium adsorption also strongly depends on whether any atoms or molecules are being encapsulated by the fullerene, as these encapsulated atoms/molecules can also interact with the fullerene and influence its stability. A Mulliken electronegativity analysis revealed that, when atoms inside and/or outside the fullerene donate charge (electrons) to the fullerene, the fullerene expands. The excess charge on the carbon atoms of the fullerene weakens some of the carbon–carbon bonds, potentially causing them to break, in which case the fullerene loses its ability to encapsulate molecules and releases them.

Open image in new window

Graphical Abstract DFT simulation of a endo fullerene doped with radioactive sodium iodide interacting with 28 calcium atoms in a geometric arrangement

     

Prediction of chlortetracycline adsorption on the Fe3O4 nanoparticle using molecular dynamics simulation

Abstract

Molecular dynamics (MD) simulation was applied to investigate the adsorption mechanism of chlortetracycline (CTC) antibiotic molecule as the aqueous pollutant on the Fe₃O₄ nanoparticle (NP). Two different NP sizes with a diameter of about 1.4?nm and 3.5?nm were selected. Initially, the stability of both NPs in water was investigated by calculating radial distribution function curves of NP atoms. Simulation results confirmed the stable crystallographic structures of both NPs. However, small NP induce greater structural stabilization. Then, CTC molecules were adsorbed on NPs surface in various pollutant concentrations. Electrostatic and hydrogen bond were the major types of interactions between CTC molecules and the adsorbent surface. CTC molecules formed a complex with NP surface from their amine side chains; while they were parallel to each other in their aromatic rings and π-π bond between two CTC molecules was formed. Diffusion rate of CTC molecules could predict the adsorption mechanism. At lower concentration of CTC, CTC molecules tend to adsorb on the NP surface. At these concentrations, the diffusion rate of CTC was high. By increasing the CTC concentration, the pollutant agglomeration was enhanced which decreased the diffusion rate. At this time, the surface of NP was saturated. In addition, the results of isotherm curves showed that CTC adsorption on small NPs could be defined with both Langmuir and Freundlich isotherm models, while Freundlich isotherm model was more appropriate for larger NPs. In conclusion, observations confirmed that MD simulation could successfully predict the behavior of CTC adsorption on the Fe₃O₄ NP surface.

Communicated by Ramaswamy H. Sarma     

Theoretical evaluation of flotation performance of carboxyl hydroxamic acids with different number of polar groups on the surfaces of diaspore (010) and kaolinite (001)

The adsorption behaviors of three carboxyl hydroxamic acids on diaspore (010) and kaolinite (001) have been studied by density functional theory (DFT) and molecular dynamics (MD) method. The results indicated that carboxyl hydroxamic acids could adsorb on diaspore surface by ionic bonds and hydrogen bonds, and adsorb on kaolinite surface by hydrogen bonds. The models of carboxyl hydroxamic acids adsorbed on diaspore and kaolinite surfaces are proposed.

An explorative study on Staphylococcus aureus MurE inhibitor: induced fit docking,binding free energy calculation,and molecular dynamics

Staphylococcus aureus MurE enzyme catalyzes the addition of l-lysine as third residue of the peptidoglycan peptide moiety. Due to the high substrate specificity and its ubiquitous nature among bacteria, MurE enzyme is considered as one of the potential target for the development of new therapeutic agents. In the present work, induced fit docking (IFD), binding free energy calculation, and molecular dynamics (MD) simulation were carried out to elucidate the inhibition potential of 2-thioxothiazolidin-4-one based inhibitor 1 against S. aureus MurE enzyme. The inhibitor 1 formed majority of hydrogen bonds with the central domain residues Asn151, Thr152, Ser180, Arg187, and Lys219. Binding free-energy calculation by MM-GBSA approach showed that van der Waals (ΔG_vdW, ?57.30?kcal/mol) and electrostatic solvation (ΔG_solv, ?36.86?kcal/mol) energy terms are major contributors for the inhibitor binding. Further, 30-ns MD simulation was performed to validate the stability of ligand–protein complex and also to get structural insight into mode of binding. Based on the IFD and MD simulation results, we designed four new compounds D1–D4 with promising binding affinity for the S. aureus MurE enzyme. The designed compounds were subjected to the extra-precision docking and binding free energy was calculated for complexes. Further, a 30-ns MD simulation was performed for D1/4C13 complex.     

Performance of the DLPNO-CCSD and recent DFT methods in the calculation of isotropic and dipolar contributions to 14N hyperfine coupling constants of nitroxide radicals

Adsorption mechanism and competitive adsorption of NH₃ and As₂O₃ molecules on CuO (111) surface were investigated via density functional theory (DFT) simulations. The adsorption configuration, adsorption energy, electronic gains/losses, and projected density of states were thoroughly discussed. Results showed that As₂O₃ molecule was more likely to adsorb on CuO surface than NH₃ molecule. For the adsorption of As₂O₃ on CuO surface, the “O_suf-O_suf bridge” site was the active site with the chemical adsorption energy of −1.39 eV. In the process of NH₃ adsorption on CuO surface, N-Cu bond was formed in chemical adsorption, which was mainly attributed to the charge transfer from NH₃ molecule to substrate surface. However, no obvious structural change occurred in physisorption. Thermodynamic analysis indicated that physisorption was the primary form for the adsorption of NH₃ and As₂O₃ molecules on CuO surface under the actual SCR condition.

     

Developing polarized protein-specific charges for protein dynamics: MD free energy calculation of pKa shifts for Asp26/Asp20 in thioredoxin

Ab initio quantum mechanical calculation of protein in solution is carried out to generate polarized protein-specific charge(s) (PPC) for molecular dynamics (MD) stimulation of protein. The quantum calculation of protein is made possible by developing a fragment-based quantum chemistry approach in combination with the implicit continuum solvent model. The computed electron density of protein is utilized to derive PPCs that represent the polarized electrostatic state of protein near the native structure. These PPCs are atom-centered like those in the standard force fields and are thus computationally attractive for molecular dynamics simulation of protein. Extensive MD simulations have been carried out to investigate the effect of electronic polarization on the structure and dynamics of thioredoxin. Our study shows that the dynamics of thioredoxin is stabilized by electronic polarization through explicit comparison between MD results using PPC and AMBER charges. In particular, MD free-energy calculation using PPCs accurately reproduced the experimental value of pK_a shift for ionizable residue Asp²⁶ buried inside thioredoxin, whereas previous calculations using standard force fields overestimated pK_a shift by twice as much. Accurate prediction of pK_a shifts by rigorous MD free energy simulation for ionizable residues buried inside protein has been a significant challenge in computational biology for decades. This study presented strong evidence that electronic polarization of protein plays an important role in protein dynamics.     

A force field for naproxen

A united-atom potential model for naproxen suitable for molecular dynamics (MD) simulation has been developed. The charge distribution is approximated by point charges obtained from ab initio calculations using the CHELPG method. Also the intramolecular interactions such as bond and angle vibration, and the torsion potential are obtained from ab initio calculations. The dispersive interaction contribution is taken from the literature. By MD simulation using a naproxen film in slap geometry, the temperature dependence of the density, surface tension and self-diffusion coefficient as well as the melting temperature for the developed potential model are obtained.     

Molecular insight into adsorption affinities of Carmustine drug on boron and nitrogen doped functionalized single-walled carbon nanotubes using density functional theory including dispersion correction calculations and molecular dynamics simulation

Abstract

We report a quantum mechanics calculation and molecular dynamics simulation study of Carmustine drug (BNU) adsorption on the surface of nitrogen (N) and boron (B) doped-functionalized single-walled carbon nanotubes. The stability of the optimized complexes is determined on the basis of relative adsorption energy (ΔE_ads). The ΔE_ads results claim that drug molecule tends to adsorb on the nitrogen and boron doped functionalized tubes with the energy values in the range of ?61.177 to ?95.806?kJ/mol. Based on the obtained results, it is observed that N-doping compared with B-doping has improved more effectively drug absorption on the surface of functionalized nanotube. The results of Atoms in Molecule calculations indicate that drug adsorbs molecularly via hydrogen bonds interactions on the surface doped-functionalized carbon nanotubes. Moreover, molecular dynamics simulation is performed to investigate the dynamics behavior of the drug molecules on the nitrogen-doped functionalized carbon nanotube (f-NNT) and functionalized carbon nanotube (f-CNT). The higher average calculated electrostatic and van der Waals energies as well as higher number of intermolecular hydrogen bonds in BNU-f-NNT in comparison with BNU-f-CNT model suggest the more effectual interaction between drug molecules and nitrogen-doped functionalized carbon nanotube.

Communicated by Ramaswamy H. Sarma     

Vibrational spectra of an RDX film over an aluminum substrate from molecular dynamics simulations and density functional theory

We report calculated vibrational spectra in the range of 0–3,500 cm^?1 of RDX (hexahydro-1,3,5-trinitro-1,3,5-triazine) molecules adsorbed on a model aluminum surface. A molecular film was modeled using two approaches: (1) density functional theory (DFT) was used to optimize a single RDX molecule interacting with its periodic images, and (2) a group of nine molecules extracted from the crystal structure was deposited on the surface and interacted with its periodic images via molecular dynamics (MD) simulations. In both cases, the molecule was initialized in the AAA conformer geometry having the three nitro groups in axial positions, and kept that conformation in the DFT examination, but some molecules were found to change to the AAE conformer (two nitro groups in axial and one in equatorial position) in the MD analysis. The vibrational spectra obtained from both methods are similar to each other, except in the regions where collective RDX intermolecular interactions (captured by MD simulations) are important, and compare fairly well with experimental findings.

Binding Mode Analyses and Pharmacophore Model Development for Stilbene Derivatives as a Novel and Competitive Class of α-Glucosidase Inhibitors

Stilbene urea derivatives as a novel and competitive class of non-glycosidic α-glucosidase inhibitors are effective for the treatment of type II diabetes and obesity. The main purposes of our molecular modeling study are to explore the most suitable binding poses of stilbene derivatives with analyzing the binding affinity differences and finally to develop a pharmacophore model which would represents critical features responsible for α-glucosidase inhibitory activity. Three-dimensional structure of S. cerevisiae α-glucosidase was built by homology modeling method and the structure was used for the molecular docking study to find out the initial binding mode of compound 12, which is the most highly active one. The initial structure was subjected to molecular dynamics (MD) simulations for protein structure adjustment at compound 12-bound state. Based on the adjusted conformation, the more reasonable binding modes of the stilbene urea derivatives were obtained from molecular docking and MD simulations. The binding mode of the derivatives was validated by correlation analysis between experimental K_i value and interaction energy. Our results revealed that the binding modes of the potent inhibitors were engaged with important hydrogen bond, hydrophobic, and π-interactions. With the validated compound 12-bound structure obtained from combining approach of docking and MD simulation, a proper four featured pharmacophore model was generated. It was also validated by comparison of fit values with the K_i values. Thus, these results will be helpful for understanding the relationship between binding mode and bioactivity and for designing better inhibitors from stilbene derivatives.