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.     

Density functional theory and molecular dynamics investigations on substituted banana-shaped compounds

Density functional theory (DFT) calculations and molecular dynamics (MD) simulations on the atomic level were performed on three different substituted banana-shaped compounds derived from 1,3-phenylene bis[4-(4-n-hexyloxyphenyliminomethyl)benzoate] (P-6-O-PIMB). The DFT studies were carried out on the isolated molecules, and in the MD simulations clusters were treated with up to 64 monomers. The effect of polar substituents, such as chlorine and the nitro group, on the central 1,3-phenylene unit of banana-shaped compounds was investigated. In particular, flexibility, polarity, electrostatic potential (ESP) group charge distributions, B-factors, bending angles and molecular lengths were considered. The MD results were analysed by trajectories of significant torsion angles as well as order parameters such as radial atom pair distribution functions g(r), orientational correlation functions g(o), diffusion coefficients (D) and root mean square deviations (RMSD) values. The g(r) and g(o) values show that a certain long range order is generated by the introduction of a NO₂ group in the 2-position of the central 1,3-phenylene ring. In contrast, the chlorination at the 4 and 6 positions of the central 1,3-phenylene unit decreases the long range order tendency by its perturbation effect on the conformations in such molecules. Moreover, g(r) and g(o) values, as well as diffusion coefficients, show that in the NO₂ substituted compound the formation of microphase areas is preferred. Finally, the aggregation effect in such compounds was studied in a systematic way by a comparison of the conformational properties of the isolated molecules and the monomers in the clusters. Figure Molecular dynamics (MD) simulations on the aggregation behaviour of substituted banana-shaped compounds Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Acetylene hydrogenation on anatase TiO2(101) supported Pd4 cluster: oxygen deficiency effect

Acetylene hydrogenation on both the perfect and oxygen defective anatase TiO(2)(101) surfaces supported Pd(4) cluster has been studied using density functional theory calculations with a Hubbard U correction (DFT + U). The adsorbed Pd(4) cluster on the perfect surface prefers to form a tetrahedral structure, while it likely moves to the oxygen defective site to form a distorted tetrahedral structure by removing a bridging oxygen atom. For the defective surface, it exhibits a stronger ability to capture Pd(4) cluster as charge transfer is significantly performed due to the oxygen deficiency. Moreover, it is found that the oxygen defective surface shows higher activity for acetylene hydrogenation, and the possible reason may lie in the weaker adsorption strength between the Pd cluster and the adsorbed molecules on the defective surface as compared to the case on the perfect surface.     

DROIDS 1.20: A GUI-Based Pipeline for GPU-Accelerated Comparative Protein Dynamics

Traditional informatics in comparative genomics work only with static representations of biomolecules (i.e., sequence and structure), thereby ignoring the molecular dynamics (MD) of proteins that define function in the cell. A comparative approach applied to MD would connect this very short timescale process, defined in femtoseconds, to one of the longest in the universe: molecular evolution measured in millions of years. Here, we leverage advances in graphics-processing-unit-accelerated MD simulation software to develop a comparative method of MD analysis and visualization that can be applied to any two homologous Protein Data Bank structures. Our open-source pipeline, DROIDS (Detecting Relative Outlier Impacts in Dynamic Simulations), works in conjunction with existing molecular modeling software to convert any Linux gaming personal computer into a “comparative computational microscope” for observing the biophysical effects of mutations and other chemical changes in proteins. DROIDS implements structural alignment and Benjamini-Hochberg-corrected Kolmogorov-Smirnov statistics to compare nanosecond-scale atom bond fluctuations on the protein backbone, color mapping the significant differences identified in protein MD with single-amino-acid resolution. DROIDS is simple to use, incorporating graphical user interface control for Amber16 MD simulations, cpptraj analysis, and the final statistical and visual representations in R graphics and UCSF Chimera. We demonstrate that DROIDS can be utilized to visually investigate molecular evolution and disease-related functional changes in MD due to genetic mutation and epigenetic modification. DROIDS can also be used to potentially investigate binding interactions of pharmaceuticals, toxins, or other biomolecules in a functional evolutionary context as well.     

Molecular dynamics and atomic charge calculations in the study of heparin conformation in aqueous solution

HF/6-31G** and molecular dynamics (MD) simulations were used to evaluate the performance of different atomic charge basis sets (i.e., Mulliken, Lowdin, and Electrostatic Potential Derived Charges--ESP) in heparin simulations. HF/3-21 G calculations were also used to study the NMR conformation of the IdoA residue. The results thus obtained indicated that ESP and Lowdin charges gave the better results in heparin simulations, followed by Mulliken charges, and that the minimum-energy conformation of IdoA can be different from that observed by NMR spectroscopy by less than 1 Angstrom. However, it was found that this small conformational modification is capable of inducing a change of almost 200 kJ/mol in the interactions of heparin with the surrounding environment, which is a meaningful amount of energy in the context of ligand-receptor interactions. This information can be potentially of great relevance in the design of heparin-derived antithrombotic compounds.     

First Principles Modeling of Dissociative Adsorption at Crystal Surfaces: Hydrogen on Pt(111)

Multiscale simulation has the potential of becoming the new modeling paradigm in chemical sciences. An important class of multiscale models involves the mapping of a finer scale model into an approximate surface that is used by a coarser scale model. As a specific example of this class we present the case of the adsorption dynamics of diatomic molecules on single crystal catalyst surfaces. The prototype system studied is the dissociative adsorption of H₂ on Pt(111). The finer scale model consists of density functional theory (DFT) periodic slab calculations that provide a small dataset for training an atomistic scale potential energy surface. The coarser scale model uses a semi-classical molecular dynamics (MD) algorithm to obtain the sticking coefficient as a function of the incident energy. Comparison to experimental data and published simulation work is presented. Finally, major challenges in multiscale modeling of chemical reactivity in coupled DFT/MD simulations are discussed, specifically the need for a systematic method of assessing the accuracy of the coarse graining process.     

In silico mutations and molecular dynamics studies on a winged bean chymotrypsin inhibitor protein

Winged bean chymotrypsin inhibitor (WCI) has an intruding residue Asn14 that plays a crucial role in stabilizing the reactive site loop conformation. This residue is found to be conserved in the Kunitz (STI) family of serine protease inhibitors. To understand the contribution of this scaffolding residue on the stability of the reactive site loop, it was mutated in silico to Gly, Ala, Ser, Thr, Leu and Val and molecular dynamics (MD) simulations were carried out on the mutants. The results of MD simulations reveal the conformational variability and range of motions possible for the reactive site loop of different mutants. The N-terminus side of the scissile bond, which is close to a beta-barrel, is conformationally less variable, while the C-terminus side, which is relatively far from any such secondary structural element, is more variable and needs stability through hydrogen-bonding interactions. The simulated structures of WCI and the mutants were docked in the peptide-binding groove of the cognate enzyme chymotrypsin and the ability to form standard hydrogen-bonding interactions at P3, P1 and P2' residues were compared. The results of the MD simulations coupled with docking studies indicate that hydrophobic residues like Leu and Val at the 14th position are disruptive for the integrity of the reactive site loop, whereas a residue like Thr, which can stabilize the C-terminus side of the scissile bond, can be predicted at this position. However, the size and charge of the Asn residue made it most suitable for the best maintenance of the integrity of the reactive site loop, explaining its conserved nature in the family.     

Parameterization of aromatic azido groups: application as photoaffinity probes in molecular dynamics studies

The accuracy of molecular dynamics (MD) simulations is limited by the availability of parameters for the molecular system of interest. In most force fields, parameters of common chemical groups are already present. With the development of novel small organic molecules as probes to study biological systems, more chemical groups require parameterization. An azide group is often used in studies of biological systems but computational studies are still impeded by the lack of parameters. In this paper, we present a set of molecular mechanics (MM) parameters for aromatic and aliphatic azido groups, and their application in MD simulations of a photoaffinity probe currently used in our laboratory for mapping binding modes available in the active site of histone deacetylases. The parameters were developed for the generalized Amber force field (GAFF) using density functional theory (DFT) calculations at B3LYP 6-311G(d) level. The parameters were validated by geometry optimization and MD simulations.     

Force field development on pigments of photosystem 2 reaction centre

We developed new parameters for molecular dynamics (MD) simulations, namely partial atomic charges, equilibrium bond-lengths, angles, dihedrals, atom types, and force constants of chlorophyll a (Chl) and plastoquinone (PQ), and both reduced and neutral form of primary acceptor (PHO) molecule. These parameters are essential for MD simulations that can interpret various structure functional relationships during primary processes of charge separation and stabilization in photosystem 2 reaction centres.     

Searching the conformational complexity and binding properties of HDAC6 through docking and molecular dynamic simulations

Histone deacetylases (HDACs) are a family of proteins involved in the deacetylation of histones and other non-histones substrates. HDAC6 belongs to class II and shares similar biological functions with others of its class. Nevertheless, its three-dimensional structure that involves the catalytic site remains unknown for exploring the ligand recognition properties. Therefore, in this contribution, homology modeling, 100-ns-long Molecular Dynamics (MD) simulation and docking calculations were combined to explore the conformational complexity and binding properties of the catalytic domain 2 from HDAC6 (DD2-HDAC6), for which activity and affinity toward five different ligands have been reported. Clustering analysis allowed identifying the most populated conformers present during the MD simulation, which were used as starting models to perform docking calculations with five DD2-HDAC6 inhibitors: Cay10603 (CAY), Rocilinostat (RCT), Tubastatin A (TBA), Tubacin (TBC), and Nexturastat (NXT), and then were also submitted to 100-ns-long MD simulations. Docking calculations revealed that the five inhibitors bind at the DD2-HDAC6 binding site with the lowest binding free energy, the same binding mode is maintained along the 100-ns-long MD simulations. Overall, our results provide structural information about the molecular flexibility of apo and holo DD2-HDAC6 states as well as insight of the map of interactions between DD2-HDAC6 and five well-known DD2-HDAC6 inhibitors allowing structural details to guide the drug design. Finally, we highlight the importance of combining different theoretical approaches to provide suitable structural models for structure-based drug design.     

Comparative molecular dynamics simulations of HIV-1 integrase and the T66I/M154I mutant: binding modes and drug resistance to a diketo acid inhibitor

HIV-1 IN is an essential enzyme for viral replication and an interesting target for the design of new pharmaceuticals for use in multidrug therapy of AIDS. L-731,988 is one of the most active molecules of the class of beta-diketo acids. Individual and combined mutations of HIV-1 IN at residues T66, S153, and M154 confer important degrees of resistance to one or more inhibitors belonging to this class. In an effort to understand the molecular mechanism of the resistance of T66I/M154I IN to the inhibitor L-731,988 and its specific binding modes, we have carried out docking studies, explicit solvent MD simulations, and binding free energy calculations. The inhibitor was docked against different protein conformations chosen from prior MD trajectories, resulting in 2 major orientations within the active site. MD simulations have been carried out for the T66I/M154I DM IN, DM IN in complex with L-731,988 in 2 different orientations, and 1QS4 IN in complex with L-731,988. The results of these simulations show a similar dynamical behavior between T66I/M154I IN alone and in complex with L-731,988, while significant differences are observed in the mobility of the IN catalytic loop (residues 138-149). Water molecules bridging the inhibitor to residues from the active site have been identified, and residue Gln62 has been found to play an important role in the interactions between the inhibitor and the protein. This work provides information about the binding modes of L-731,988, as well as insight into the mechanism of inhibitor-resistance in HIV-1 integrase.     

Theoretical investigations on Azotobacter vinelandii ferredoxin I: effects of electron transfer on protein dynamics

Structural, energetic, and dynamical studies of Azotobacter vinelandii ferredoxin I are presented for native and mutant forms. The protein contains two iron-sulfur clusters, one of which ([3Fe-4S]) is believed to play a central role in the electron-coupled proton transfer. Different charge sets for the [3Fe-4S] cluster in its reduced and oxidized state are calculated with broken symmetry ab initio density functional theory methods and used in molecular dynamics (MD) simulations. The validity of the ab initio calculations is assessed by comparing partially optimized structures of the [3Fe-4S] clusters with x-ray structures. Possible proton transfer pathways between the protein and the iron-sulfur cluster are examined by both MD simulations and ab initio calculations. The MD simulations identify three main-chain hydrogen atoms--HN(13), HN(14), and HN(16)--that are within H-bonding distance of the [3Fe-4S] cluster throughout the MD simulations. They could thus play a role in the proton transfer from the protein to the iron-sulfur cluster. By contrast, the HD2(15) atom of the Asp-15 is seldom close enough to the [3Fe-4S] cluster to transfer a proton. Poisson-Boltzmann calculations indicate that there is a low, but nonzero probability, that Asp-15 is protonated at pH 7; this is a requirement for it to serve as a proton donor. Ab initio calculations with a fragment model for the protein find similar behavior for the transfer of a proton from the OH of the protonated side chain and the main-chain NH of Asp-15. The existence of a stable salt bridge between Asp-15 and Lys-84 in the D15E mutant, versus its absence in the wild-type, has been suggested as the cause of the difference in the rate of proton transfer. Extensive MD simulations were done to test this idea; the results do not support the proposal. The present findings, together with the available data, serve as the basis for an alternative proposal for the mechanism of the coupled electron-proton transfer reaction in ferredoxin I.     

DFT: a dynamic study of the interaction of ethanol and methanol with platinum

The interaction between alcohol molecules and platinum (Pt) was studied using molecular dynamics (MD; Born-Oppenheimer method). Alcohol molecules like ethanol and methanol present a similar molecular structure, with a methyl group (CH₃) at one end and a fragment of hydroxyl (OH) at the other. This fact generates two orientations that are considered in the interaction with Pt. The MD calculation results for these two orientations indicate a preferential orientation due to energy interactions. A plausible reaction mechanism that takes into account the interaction between Pt and alcohol is presented. The charge transference obtained from the Pt–alcohol interaction was also analyzed. The energy for the two orientations was calculated by indicating the preferential orientation. The methyl and hydroxyl groups are involved in heterolytic breakage of hydrogen bonds, joined to a carbon atom in the former and to an oxygen atom in the latter; however, the methyl group reaction seems to be the most important.     

Thermal decomposition behaviour of polyethylene in oxygen-free and low oxygen content circumstances by reactive molecular dynamic simulation

In the present study, the thermal decomposition characteristics of polyethylene (PE) in oxygen-free and low oxygen content circumstances were examined by molecular dynamic (MD) simulations at atomic scale using reactive force field (ReaxFF). Temporal evolutions of species were captured reasonably during the processes of thermal decomposition. The effects of oxygen content, temperature and heating rate were also analysed. In addition, the kinetic properties were predicted with reliable parameters. The results show good agreements with the available ones, which illustrate that the species with two carbon atoms are the vast majority of final products. Higher oxygen content and temperature promote the generation of small molecules with carbon atom number less than or equal to 10. In the presence of oxygen, greater activation energy span and reaction order are calculated with lower adjust R², which indicates complex reactions according to kinetic analysis. The initial decomposition temperature of PE is proportional to the heating rate owing to heat transfer lag.     

A comparison of four different conformations adopted by human telomeric G‐quadruplex using computer simulations

The telomeric G‐quadruplexes for their unique structural features are considered as potential anticancer drug targets. These, however, exhibit structural polymorphism as different topology types for the intra‐molecular G‐quadruplexes from human telomeric G‐rich sequences have been reported based on NMR spectroscopy and X‐ray crystallography. These techniques provide detailed atomic‐level information about the molecule but relative conformational stability of the different topologies remains unsolved. Therefore, to understand the conformational preference, we have carried out quantum chemical calculations on G‐quartets; used all‐atom molecular dynamics (MD) simulations and steered molecular dynamics (SMD) simulations to characterize the four human telomeric G‐quadruplex topologies based on its G‐tetrad core‐types, viz., parallel, anti‐parallel, mixed‐(3 + 1)‐form1 and mixed‐(3 + 1)‐form2. We have also studied a non‐telomeric sequence along with these telomeric forms giving a comparison between the two G‐rich forms. The structural properties such as base pairing, stacking geometry and backbone conformations have been analyzed. The quantum calculations indicate that presence of a sodium ion inside the G‐tetrad plane or two potassium ions on both sides of the plane give it an overall planarity which is much needed for good stacking to form a helix. MD simulations indicate that capping of the G‐tetrad core by the TTA loops keep the terminal guanine bases away from water. The SMD simulations along with equilibrium MD studies indicate that the parallel and non‐telomeric forms are comparatively less stable. We could come to the conclusion that the anti‐parallel form and also the mixed‐(3 + 1)‐form1 topology are most likely to represent the major conformation., 2016. © 2015 Wiley Periodicals, Inc. Biopolymers 105: 83–99, 2016     

Thermal and pressure-induced martensitic phase transformations in a Ni–Al alloy modelled by Sutton–Chen embedded atom method

The Ni–Al alloys which exhibit the thermoelastic martensitic phase transformations in the composition range from 60 to 65 atomic percentage (at.%) of Ni are widely used in the high technology applications. In this study, both thermal and pressure-induced phase transformations in Ni-37.5 at.%Al alloy model were investigated by a molecular dynamics (MD) method. Physical interactions between atoms in the alloy system were modelled using the Sutton–Chen version of the embedded atom method based on many-body interactions. The potential parameters for cross interactions between Ni and Al atoms were estimated by optimising the results obtained from the MD simulations, taking into account the experimental data including the crystal lattice properties of the model alloy in high temperature phase.     

Deep analysis of N-cadherin/ADH-1 interaction: a computational survey

Due to the considerable role of N-cadherin in cancer metastasis, tumor growth, and progression, inhibition of this protein has been highly regarded in recent years. Although ADH-1 has been known as an appropriate inhibitor of N-cadherin in clinical trials, its chemical nature and binding mode with N-cadherin have not been precisely specified yet. Accordingly, in this study, quantum mechanics calculations were used to investigate the chemical nature of ADH-1. These calculations clarify the molecular properties of ADH-1 and determine its reactive sites. Based on the results, the oxygen atoms are suitable for electrophilic reactivity, while the hydrogen atoms that are connected to nitrogen atoms are the favorite sites for nucleophilic reactivity. The higher electronegativity of the oxygen atoms makes them the most reactive portions in this molecule. Molecular docking and molecular dynamics (MD) simulation have also been applied to specify the binding mode of ADH-1 with N-cadherin and determine the important residues of N-cadherin involving in the interaction with ADH-1. Moreover, the verified model by MD simulation has been studied to extract the free energy value and find driving forces. These calculations and molecular electrostatic potential map of ADH-1 indicated that hydrophobic and electrostatic interactions are almost equally involved in the implantation of ADH-1 in the N-cadherin binding site. The presented results not only enable a closer examination of N-cadherin in complex with ADH-1 molecule, but also are very beneficial in designing new inhibitors for N-cadherin and can help to save time and cost in this field.     

Understanding the conformational flexibility and electrostatic properties of curcumin in the active site of rhAChE via molecular docking,molecular dynamics,and charge density analysis

Acetylcholinesterase (AChE) is an important enzyme responsible for Alzheimer’s disease, as per report, keto-enol form of curcumin inhibits this enzyme. The present study aims to understand the binding mechanism of keto-enol curcumin with the recombinant human Acetylcholinesterase (rhAChE) from its conformational flexibility, intermolecular interactions, charge density distribution, and the electrostatic properties at the active site of rhAChE. To accomplish this, a molecular docking analysis of curcumin with the rhAChE was performed, which gives the structure and conformation of curcumin in the active site of rhAChE. Further, the charge density distribution and the electrostatic properties of curcumin molecule (lifted from the active site of rhAChE) were determined from the high level density functional theory (DFT) calculations coupled with the charge density analysis. On the other hand, the curcumin molecule was optimized (gas phase) using DFT method and further, the structure and charge density analysis were also carried out. On comparing the conformation, charge density distribution and the electrostatic potential of the active site form of curcumin with the corresponding gas phase form reveals that the above said properties are significantly altered when curcumin is present in the active site of rhAChE. The conformational stability and the interaction of curcumin in the active site are also studied using molecular dynamics simulation, which shows a large variation in the conformational geometry of curcumin as well as the intermolecular interactions.     

A comparative study of trypsin specificity based on QM/MM molecular dynamics simulation and QM/MM GBSA calculation

Hydrogen bonding and polar interactions play a key role in identification of protein-inhibitor binding specificity. Quantum mechanics/molecular mechanics molecular dynamics (QM/MM MD) simulations combined with DFT and semi-empirical Hamiltonian (AM1d, RM1, PM3, and PM6) methods were performed to study the hydrogen bonding and polar interactions of two inhibitors BEN and BEN1 with trypsin. The results show that the accuracy of treating the hydrogen bonding and polar interactions using QM/MM MD simulation of PM6 can reach the one obtained by the DFT QM/MM MD simulation. Quantum mechanics/molecular mechanics generalized Born surface area (QM/MM-GBSA) method was applied to calculate binding affinities of inhibitors to trypsin and the results suggest that the accuracy of binding affinity prediction can be significantly affected by the accurate treatment of the hydrogen bonding and polar interactions. In addition, the calculated results also reveal the binding specificity of trypsin: (1) the amidinium groups of two inhibitors generate favorable salt bridge interaction with Asp189 and form hydrogen bonding interactions with Ser190 and Gly214, (2) the phenyl of inhibitors can produce favorable van der Waals interactions with the residues His58, Cys191, Gln192, Trp211, Gly212, and Cys215. This systematic and comparative study can provide guidance for the choice of QM/MM MD methods and the designs of new potent inhibitors targeting trypsin.     

Binding free energies for nicotine analogs inhibiting cytochrome P450 2A6 by a combined use of molecular dynamics simulations and QM/MM-PBSA calculations

Molecular dynamics (MD) simulations and hybrid quantum mechanical/molecular mechanical (QM/MM) calculations have been performed to explore the dynamic behaviors of cytochrome P450 2A6 (CYP2A6) binding with nicotine analogs (that are typical inhibitors) and to calculate their binding free energies in combination with Poisson–Boltzmann surface area (PBSA) calculations. The combined MD simulations and QM/MM-PBSA calculations reveal that the most important structural parameters affecting the CYP2A6-inhibitor binding affinity are two crucial internuclear distances, that is, the distance between the heme iron atom of CYP2A6 and the coordinating atom of the inhibitor, and the hydrogen-bonding distance between the N297 side chain of CYP2A6 and the pyridine nitrogen of the inhibitor. The combined MD simulations and QM/MM-PBSA calculations have led to dynamic CYP2A6-inhibitor binding structures that are consistent with the observed dynamic behaviors and structural features of CYP2A6-inhibitor binding, and led to the binding free energies that are in good agreement with the experimentally-derived binding free energies. The agreement between the calculated binding free energies and the experimentally-derived binding free energies suggests that the combined MD and QM/MM-PBSA approach may be used as a valuable tool to accurately predict the CYP2A6-inhibitor binding affinities in future computational design of new, potent and selective CYP2A6 inhibitors.