Mechanical behaviour of bio-inspired brittle-matrix nanocomposites under different strain rates using molecular dynamics

We study the generic mechanical behaviour of ceramic–ceramic nanocomposites inspired from biological materials. The nanocomposite models considered in our study are the regularly and stairwise staggered arrangements of stiff brittle platelets embedded in compliant brittle matrix. Molecular dynamics simulations are carried out to investigate the effect of strain rate on these nanocomposites. The variation in stress–strain behaviour and mechanical properties are analysed. The evolution of deformation processes is also investigated. Our results show the existence of different strain rate regimes separated by critical strain rate. Deformation mechanisms such as matrix cracking, crack bridging, interfacial debonding and hence platelet pullout are observed at lower strain rates. Amorphous deformation and direct debonding without matrix cracking are observed at higher strain rates.     

Molecular dynamics (MD) simulations and binding free energy calculation studies between inhibitors and type II dehydroquinase (DHQ2)

Type II dehydroquinase (DHQ2) is the third enzyme of the shikimic acid pathway, and it has been the effective target for tuberculosis (TB). So far, developing multiple potent inhibitors of the DHQ2 of Mycobacterium tuberculosis (DHQ2-Mt) has been considered to be the new therapy to TB. Molecular dynamics simulations followed by molecular mechanics-generalised Born surface area were carried out to calculate the free binding energy and to determine the affinity ability of the four chosen inhibitor molecules, L1, L2, L3 and L4. Energy decomposition analyses show the electrostatic interaction and van der Waals interaction of the ligands to every residue of the DHQ2-Mt. The results suggest that some important residues have different interactions with the four ligands, such as Arg19 and Tyr24. These interactions may have an effect on the ligand binding affinity. The binding affinity of monosubstituted inhibitor is higher than that of disubstituted inhibitor, due to some important interactions with the DHQ2-Mt residues. These computational works will be significant to the theoretical research in the future.     

Diffusional interaction behavior of NSAIDs in lipid bilayer membrane using molecular dynamics (MD) simulation: Aspirin and Ibuprofen

In this research, for the first time, molecular dynamics (MD) method was used to simulate aspirin and ibuprofen at various concentrations and in neutral and charged states. Effects of the concentration (dosage), charge state, and existence of an integral protein in the membrane on the diffusion rate of drug molecules into lipid bilayer membrane were investigated on 11 systems, for which the parameters indicating diffusion rate and those affecting the rate were evaluated. Considering the diffusion rate, a suitable score was assigned to each system, based on which, analysis of variance (ANOVA) was performed. By calculating the effect size of the indicative parameters and total scores, an optimum system with the highest diffusion rate was determined. Consequently, diffusion rate controlling parameters were obtained: the drug–water hydrogen bond in protein-free systems and protein–drug hydrogen bond in the systems containing protein.     

Deoxyribonuclease I Inhibitory Properties,Molecular Docking and Molecular Dynamics Simulations of 1-(Pyrrolidin-2-yl)propan-2-one Derivatives

Deoxyribonuclease I (DNase I) inhibitory properties of two 1-(pyrrolidin-2-yl)propan-2-one derivatives were examined in vitro. Determined IC₅₀ values of 1-[1-(4-methoxyphenyl)pyrrolidin-2-yl]propan-2-one ( 1 ) (192.13±16.95 μM) and 1-[1-(3,4,5-trimethoxyphenyl)pyrrolidin-2-yl]propan-2-one ( 2 ) (132.62±9.92 μM) exceed IC₅₀ value of crystal violet, used as a positive control, 1.89- and 2.73-times, respectively. Compounds are predicted to be nontoxic and to have favorable pharmacokinetic profiles, with high gastrointestinal absorption and blood-brain barrier permeability. Molecular docking and molecular dynamics simulations suggest that interactions with Glu 39, Glu 78, Arg 111, Pro 137, Asp 251 and His 252 are an important factor for inhibitors affinity toward the DNase I. Determined inhibitory properties along with predicted ADMET profiles and observed interactions would be beneficial for the discovery of new active 1-(pyrrolidin-2-yl)propan-2-one-based inhibitors of DNase I.     

Molecular dynamics simulations reveal the disruption mechanism of a 2,4-thiazolidinedione derivative C30 against tau hexapeptide (PHF6) oligomer

Derivatives of 2,4-thiazolidinedione have been reported to inhibit the aggregation of tau protein, in which compound 30 (C30) not only inhibit 80% of paired helical filament 6 (PHF6) aggregation, but also inhibit K18 and full-length tau aggregation. However, its inhibitory mechanism is unclear. In this study, to investigate the effect of C30 on tau protein, all-atom molecular dynamics simulation was performed on the PHF6 oligomer with and without C30. The results show that C30 can cause significant conformational changes in the PHF6 oligomer. The nematic order parameter P2 and secondary structure analyses show that C30 destroys the ordered structure of PHF6 oligomer, reduces the content of β-sheet structure, and transforms β-sheet into random coil structure. By clustering analysis, it was found that C30 has four possible binding sites on the PFH6 oligomer, and the binding ability order is S1 > S2 > S4 > S3. Following a more in-depth analyses of each site, it was determined that the S1 site is the most possible binding site mainly located between layers of L1 and L3. The hydrophobic interaction is the driving force for the binding of C30 to PHF6 oligomer. In addition, L1P4_Y310, L1P5_Y310, L3P1_V309, and L3P2_V309 are key residues for C30 binding to oligomer. Moreover, π-π interaction formed by L1P4_Y310 and L1P5_Y310 with C30 and the hydrogen bonding interaction formed by C30 with L3P3_Q307 are beneficial to the combination of C30 and oligomer. The fully understanding disrupt the mechanism of 2,4-thiazolidinedione derivative on PHF6 oligomer and the identification of binding sites will help design and discover new AD inhibitors in the future.     

Molecular dynamics simulations on the conformational transitions from the GA98 (GA88) to GB98 (GB88) proteins

We performed conventional and targeted molecular dynamics simulations to address the dynamic transition mechanisms of the conformational transitions from the G_A98 protein with only 1 mutation of Leu45Tyr to G_B98 and from the G_A88 protein with 7 mutations of Gly24Ala, Ile25Thr, Ile30Phe, Ile33Tyr, Leu45Tyr, Ile49Thr, and Leu50Lys to G_B88. The results show that the conformational transition mechanism from the mutated 3α G_A98 (G_A88) state to the α+4β G_B98 (G_B88) state via several intermediate conformations involves the bending of loops at the N and C termini firstly, the unfolding of αA and αC, then the traversing of αB, and the formation of the 4β layer with the conversion of the hydrophobic core. The bending of loops at the N and C termini and the formation of the crucial transition conformation with the full unfolded structure are key factors in their transition processes. The communication of the interaction network, the bending directions of loops, and the traversing site of αB in the transition of G_A98 to G_B98 are markedly different from those in G_A88 to G_B88 because of the different mutated residues. The analysis of the correlations and the calculated mass center distances between some segments further supported their conformational transition mechanisms. These results could help people to better understand the Paracelsus challenge. Copyright © 2016 John Wiley & Sons, Ltd.     

Partitioning of 2,6-Bis(1H-Benzimidazol-2-yl)pyridine fluorophore into a phospholipid bilayer: complementary use of fluorescence quenching studies and molecular dynamics simulations

Successful use of fluorescence sensing in elucidating the biophysical properties of lipid membranes requires knowledge of the distribution and location of an emitting molecule in the bilayer. We report here that 2,6-bis(1H-benzimidazol-2-yl)pyridine (BBP), which is almost non-fluorescent in aqueous solutions, reveals a strong emission enhancement in a hydrophobic environment of a phospholipid bilayer, making it interesting for fluorescence probing of water content in a lipid membrane. Comparing the fluorescence behavior of BBP in a wide variety of solvents with those in phospholipid vesicles, we suggest that the hydrogen bonding interactions between a BBP fluorophore and water molecules play a crucial role in the observed “light switch effect”. Therefore, the loss of water-induced fluorescence quenching inside a membrane are thought to be due to deep penetration of BBP into the hydrophobic, water-free region of a bilayer. Characterized by strong quenching by transition metal ions in solution, BBP also demonstrated significant shielding from the action of the quencher in the presence of phospholipid vesicles. We used the increase in fluorescence intensity, measured upon titration of probe molecules with lipid vesicles, to estimate the partition constant and the Gibbs free energy (ΔG) of transfer of BBP from aqueous buffer into a membrane. Partitioning BBP revealed strongly favorable ΔG, which depends only slightly on the lipid composition of a bilayer, varying in a range from − 6.5 to − 7.0 kcal/mol. To elucidate the binding interactions of the probe with a membrane on the molecular level, a distribution and favorable location of BBP in a POPC bilayer were modeled via atomistic molecular dynamics (MD) simulations using two different approaches: (i) free, diffusion-driven partitioning of the probe molecules into a bilayer and (ii) constrained umbrella sampling of a penetration profile of the dye molecule across a bilayer. Both of these MD approaches agreed with regard to the preferred location of a BBP fluorophore within the interfacial region of a bilayer, located between the hydrocarbon acyl tails and the initial portion of the lipid headgroups. MD simulations also revealed restricted permeability of water molecules into this region of a POPC bilayer, determining the strong fluorescence enhancement observed experimentally for the membrane-partitioned form of BBP.     

Molecular dynamics simulations of a cyclic-beta-(1-->2) glucan containing an alpha-(1-->6) linkage as a 'molecular alleviator' for the macrocyclic conformational strain

The conformational preferences of a cyclic osmoregulated periplasmic glucan of Ralstonia solanacearum (OPGR), which is composed of 13 glucose units and linked entirely via beta-(1-->2) linkages excluding one alpha-(1-->6) linkage, were characterized by molecular dynamics simulations. Of the three force fields modified for carbohydrates that were applied to select a suitable one for the cyclic glucan, the carbohydrate solution force field (CSFF) was found to most accurately simulate the cyclic molecule. To determine the conformational characteristics of OPGR, we investigated the glycosidic dihedral angle distribution, fluctuation, and the potential energy of the glucan and constructed hypothetical cyclic (CYS13) and linear (LINEAR) glucans. All beta-(1-->2)-glycosidic linkages of OPGR adopted stable conformations, and the dihedral angles fluctuated in this energy region with some flexibility. However, despite the inherent flexibility of the alpha-(1-->6) linkage, the dihedral angles have no transition and are more rigid than that in a linear glucan. CYS13, which consists of only beta-(1-->2) linkages, is somewhat less flexible than other glycans, and one of its linkages adopts a higher energy conformation. In addition, the root-mean-square fluctuation of this linkage is lower than that of other linkages. Furthermore, the potential energy of glucans increases in the order of LINEAR, OPGR, and CYS13. These results provide evidence of the existence of conformational constraints in the cyclic glucan. The alpha-(1-->6)-glycosidic linkage can relieve this constraint more efficiently than the beta-(1-->2) linkage. The conformation of OPGR can reconcile the tendency for individual glycosidic bonds to adopt energetically favorable conformations with the requirement for closure of the macrocyclic ring by losing the inherent flexibility of the alpha-(1-->6)-glycosidic linkage.     

Molecular dynamic simulations of cisplatin- and oxaliplatin-d(GG) intrastand cross-links reveal differences in their conformational dynamics

Mismatch repair proteins, DNA damage-recognition proteins and translesion DNA polymerases discriminate between Pt-GG adducts containing cis-diammine ligands (formed by cisplatin (CP) and carboplatin) and trans-RR-diaminocyclohexane ligands (formed by oxaliplatin (OX)) and this discrimination is thought to be important in determining differences in the efficacy, toxicity and mutagenicity of these platinum anticancer agents. We have postulated that these proteins recognize differences in conformation and/or conformational dynamics of the DNA containing the adducts. We have previously determined the NMR solution structure of OX-DNA, CP-DNA and undamaged duplex DNA in the 5'-d(CCTCAGGCCTCC)-3' sequence context and have shown the existence of several conformational differences in the vicinity of the Pt-GG adduct. Here we have used molecular dynamics simulations to explore differences in the conformational dynamics between OX-DNA, CP-DNA and undamaged DNA in the same sequence context. Twenty-five 10 ns unrestrained fully solvated molecular dynamics simulations were performed starting from two different DNA conformations using AMBER v8.0. All 25 simulations reached equilibrium within 4 ns, were independent of the starting structure and were in close agreement with previous crystal and NMR structures. Our data show that the cis-diammine (CP) ligand preferentially forms hydrogen bonds on the 5' side of the Pt-GG adduct, while the trans-RR-diaminocyclohexane (OX) ligand preferentially forms hydrogen bonds on the 3' side of the adduct. In addition, our data show that these differences in hydrogen bond formation are strongly correlated with differences in conformational dynamics, specifically the fraction of time spent in different DNA conformations in the vicinity of the adduct, for CP- and OX-DNA adducts. We postulate that differential recognition of CP- and OX-GG adducts by mismatch repair proteins, DNA damage-recognition proteins and DNA polymerases may be due, in part, to differences in the fraction of time that the adducts spend in a conformation favorable for protein binding.     

Investigation of TiO2 anatase (1 0 1), (1 0 0) and (1 1 0) facets as immobilizer for a potential anticancer RNA aptamer: a classical molecular dynamics simulation

One of the most outstanding properties of TiO₂ nanosheets is their lack of harmful effects on the public health and environment, which makes them an appropriate agent for medical applications such as drug delivery. Interaction of an RNA aptamer with (1 0 1), (1 0 0) and (1 1 0) surfaces of TiO₂ anatase were investigated using the molecular dynamics simulation. The structural parameters including root-mean-square deviation and fluctuation, and the distance between the center-of-mass of RNA aptamer and the considered surfaces were discussed in detail. Besides, the effect of water between adsorbed aptamer and surface was investigated and analyzed by the help of dipole moment orientation, hydrogen bonds and density profile of these water molecules. Analysis of the structural parameters and interaction energies shows that the (1 1 0) surface is energetically more favorable for the adsorption of considered RNA aptamer than the (1 0 0) and (1 0 1) surfaces. Consequently, our results suggest a great potential of (1 1 0) surface of TiO₂ as an efficient candidate for drug delivery applications.     

Insights into the binding specificity of wild type and mutated wheat germ agglutinin towards Neu5Acα(2‐3)Gal: a study by in silico mutations and molecular dynamics simulations

Wheat germ agglutinin (WGA) is a plant lectin, which specifically recognizes the sugars NeuNAc and GlcNAc. Mutated WGA with enhanced binding specificity can be used as biomarkers for cancer. In silico mutations are performed at the active site of WGA to enhance the binding specificity towards sialylglycans, and molecular dynamics simulations of 20 ns are carried out for wild type and mutated WGAs (WGA1, WGA2, and WGA3) in complex with sialylgalactose to examine the change in binding specificity. MD simulations reveal the change in binding specificity of wild type and mutated WGAs towards sialylgalactose and bound conformational flexibility of sialylgalactose. The mutated polar amino acid residues Asn114 (S114N), Lys118 (G118K), and Arg118 (G118R) make direct and water mediated hydrogen bonds and hydrophobic interactions with sialylgalactose. An analysis of possible hydrogen bonds, hydrophobic interactions, total pair wise interaction energy between active site residues and sialylgalactose and MM‐PBSA free energy calculation reveals the plausible binding modes and the role of water in stabilizing different binding modes. An interesting observation is that the binding specificity of mutated WGAs (cyborg lectin) towards sialylgalactose is found to be higher in double point mutation (WGA3). One of the substituted residues Arg118 plays a crucial role in sugar binding. Based on the interactions and energy calculations, it is concluded that the order of binding specificity of WGAs towards sialylgalactose is WGA3 > WGA1 > WGA2 > WGA. On comparing with the wild type, double point mutated WGA (WGA3) exhibits increased specificity towards sialylgalactose, and thus, it can be effectively used in targeted drug delivery and as biological cell marker in cancer therapeutics. Copyright © 2014 John Wiley & Sons, Ltd.