Prediction of Mechanical Properties of Polymers With Various Force Fields

The effect of force field type on the predicted elastic properties of a polyimide is examined using a multiscale modeling technique. Molecular Dynamics simulations are used to predict the atomic structure and elastic properties of the polymer by subjecting a representative volume element of the material to bulk and shear finite deformations. The elastic properties of the polyimide are determined using three force fields: AMBER, OPLS-AA, and MM3. The predicted values of Young s modulus and shear modulus of the polyimide are compared with experimental values. The results indicate that the mechanical properties of the polyimide predicted with the OPLS-AA force field most closely matched those from experiment. The results also indicate that while the complexity of the force field does not have a significant effect on the accuracy of predicted properties, small differences in the force constants and the functional form of individual terms in the force fields determine the accuracy of the force field in predicting the elastic properties of the polyimide

On Conduction in a Bacterial Sodium Channel

Voltage-gated Na+-channels are transmembrane proteins that are responsible for the fast depolarizing phase of the action potential in nerve and muscular cells. Selective permeability of Na+ over Ca2+ or K+ ions is essential for the biological function of Na+-channels. After the emergence of the first high-resolution structure of a Na+-channel, an anionic coordination site was proposed to confer Na+ selectivity through partial dehydration of Na+ via its direct interaction with conserved glutamate side chains. By combining molecular dynamics simulations and free-energy calculations, a low-energy permeation pathway for Na+ ion translocation through the selectivity filter of the recently determined crystal structure of a prokaryotic sodium channel from Arcobacter butzleri is characterised. The picture that emerges is that of a pore preferentially occupied by two ions, which can switch between different configurations by crossing low free-energy barriers. In contrast to K+-channels, the movements of the ions appear to be weakly coupled in Na+-channels. When the free-energy maps for Na+ and K+ ions are compared, a selective site is characterised in the narrowest region of the filter, where a hydrated Na+ ion, and not a hydrated K+ ion, is energetically stable

Free energy of binding of coiled-coil complexes with different electrostatic environments: the influence of force field polarisation and capping

Coiled-coils are well known protein–protein interaction motifs, with the leucine zipper region of activator protein-1 (AP-1) consisting of the c-Jun and c-Fos proteins being a typical example. Molecular dynamics (MD) simulations using the MM/GBSA method have been used to predict the free energy of interaction of these proteins. The influence of force field polarisation and capping on the predicted free energy of binding of complexes with different electrostatic environments (net charge) were investigated. Although both force field polarisation and peptide capping are important for the prediction of the absolute free energy of binding, peptide capping has the largest influence on the predicted free energy of binding. Polarisable simulations appear better suited to determine structural properties of the complexes of these proteins while non-polarisable simulations seem to give better predictions of the associated free energies of bindin

Insights into the Complex Formed by Matrix Metalloproteinase-2 and Alloxan Inhibitors: Molecular Dynamics Simulations and Free Energy Calculations

Matrix metalloproteinases (MMP) are well-known biological targets implicated in tumour progression, homeostatic regulation, innate immunity, impaired delivery of pro-apoptotic ligands, and the release and cleavage of cell-surface receptors. Hence, the development of potent and selective inhibitors targeting these enzymes continues to be eagerly sought. In this paper, a number of alloxan-based compounds, initially conceived to bias other therapeutically relevant enzymes, were rationally modified and successfully repurposed to inhibit MMP-2 (also named gelatinase A) in the nanomolar range. Importantly, the alloxan core makes its debut as zinc binding group since it ensures a stable tetrahedral coordination of the catalytic zinc ion in concert with the three histidines of the HExxHxxGxxH metzincin signature motif, further stabilized by a hydrogen bond with the glutamate residue belonging to the same motif. The molecular decoration of the alloxan core with a biphenyl privileged structure allowed to sample the deep S1′ specificity pocket of MMP-2 and to relate the high affinity towards this enzyme with the chance of forming a hydrogen bond network with the backbone of Leu116 and Asn147 and the side chains of Tyr144, Thr145 and Arg149 at the bottom of the pocket. The effect of even slight structural changes in determining the interaction at the S1′ subsite of MMP-2 as well as the nature and strength of the binding is elucidated via molecular dynamics simulations and free energy calculations. Among the herein presented compounds, the highest affinity (pIC50 = 7.06) is found for BAM, a compound exhibiting also selectivity (>20) towards MMP-2, as compared to MMP-9, the other member of the gelatinases

Recognition of 5-Hydroxymethylcytosine by the Uhrf1 SRA Domain

Recent discovery of 5-hydroxymethylcytosine (5hmC) in genomic DNA raises the question how this sixth base is recognized by cellular proteins. In contrast to the methyl-CpG binding domain (MBD) of MeCP2, we found that the SRA domain of Uhrf1, an essential factor in DNA maintenance methylation, binds 5hmC and 5-methylcytosine containing substrates with similar affinity. Based on the co-crystal structure, we performed molecular dynamics simulations of the SRA:DNA complex with the flipped cytosine base carrying either of these epigenetic modifications. Our data indicate that the SRA binding pocket can accommodate 5hmC and stabilizes the flipped base by hydrogen bond formation with the hydroxyl group

Molecular dynamics simulation studies of the interactions between ionic liquids and amino acids in aqueous solution

Although the understanding of the influence of ionic liquids (ILs) on the solubility behavior of biomolecules in aqueous solutions is relevant for the design and optimization of novel biotechnological processes, the underlying molecular-level mechanisms are not yet consensual or clearly elucidated. In order to contribute to the understanding of the molecular interactions established between amino acids and ILs in aqueous media, classical molecular dynamics (MD) simulations were performed for aqueous solutions of five amino acids with different structural characteristics (glycine, alanine, valine, isoleucine, and glutamic acid) in the presence of 1-butyl-3-methylimidazolium bis(trifluoromethyl)sulfonyl imide. The results from MD simulations enable to relate the properties of the amino acids, namely their hydrophobicity, to the type and strength of their interactions with ILs in aqueous solutions and provide an explanation for the direction and magnitude of the solubility phenomena observed in [IL + amino acid + water] systems by a mechanism governed by a balance between competitive interactions of the IL cation, IL anion, and water with the amino acids

Perturbation of the Dimer Interface of Triosephosphate Isomerase and its Effect on Trypanosoma cruzi

Most of the enzymes of parasites have their counterpart in the host. Throughout evolution, the three-dimensional architecture of enzymes and their catalytic sites are highly conserved. Thus, identifying molecules that act exclusively on the active sites of the enzymes from parasites is a difficult task. However, it is documented that the majority of enzymes consist of various subunits, and that conservation in the interface of the subunits is lower than in the catalytic site. Indeed, we found that there are significant differences in the interface between the two subunits of triosephosphate isomerase from Homo sapiens and Trypanosoma cruzi (TcTIM), which causes Chagas disease in the American continent. In the search for agents that specifically inhibit TcTIM, we found that 2,2′-dithioaniline (DTDA) is far more effective in inactivating TcTIM than the human enzyme, and that its detrimental effect is due to perturbation of the dimer interface. Remarkably, DTDA prevented the growth of Escherichia coli cells that had TcTIM instead of their own TIM and killed T. cruzi epimastigotes in culture. Thus, this study highlights a new approach base of targeting molecular interfaces of dimers

Direct knock-on of desolvated ions governs strict ion selectivity in K+ channels

The seeming contradiction that K+ channels conduct K+ ions at maximal throughput rates while not permeating slightly smaller Na+ ions has perplexed scientists for decades. Although numerous models have addressed selective permeation in K+ channels, the combination of conduction efficiency and ion selectivity has not yet been linked through a unified functional model. Here, we investigate the mechanism of ion selectivity through atomistic simulations totalling more than 400 μs in length, which include over 7,000 permeation events. Together with free-energy calculations, our simulations show that both rapid permeation of K+ and ion selectivity are ultimately based on a single principle: the direct knock-on of completely desolvated ions in the channels' selectivity filter. Herein, the strong interactions between multiple 'naked' ions in the four filter binding sites give rise to a natural exclusion of any competing ions. Our results are in excellent agreement with experimental selectivity data, measured ion interaction energies and recent two-dimensional infrared spectra of filter ion configurations