Attempting Well-Tempered Funnel Metadynamics Simulations for the Evaluation of the Binding Kinetics of Methionine Aminopeptidase-II Inhibitors.

Understanding the unbinding kinetics of protein-ligand complexes is considered a significant approach for the design of ligands with desired specificity and safety. In recent years, enhanced sampling methods have emerged as effective tools for studying the unbinding kinetics of protein-ligand complexes at the atomistic level. MetAP-II is a target for the treatment of cancer for which not a single effective drug is available yet. The identification of the dissociation rate of ligands from the complexes often serves as a better predictor for in vivo efficacy than the ligands' binding affinity. Here, funnel-based restraint well-tempered metadynamics simulations were applied to predict the residence time of two ligands bound to MetAP-II, along with the ligand association and dissociation mechanism involving the identification of the binding hotspot during ligand egress. The ligand-egressing route revealed by metadynamics simulations also correlated with the identified pathways from the CAVER analysis and by the enhanced sampling simulation using PLUMED. Ligand 1 formed a strong H-bond interaction with GLU364 estimating a higher residence time of 28.22 ± 5.29 ns in contrast to ligand 2 with a residence time of 19.05 ± 3.58 ns, which easily dissociated from the binding pocket of MetAP-II. The results obtained from the simulations were consistent to reveal ligand 1 being superior to ligand 2; however, the experimental data related to residence time were close for both ligands, and no kinetic data were available for ligand 2. The current study could be considered the first attempt to apply an enhanced sampling method for the evaluation of the binding kinetics and thermodynamics of two different classes of ligands to a binuclear metalloprotein.

Iron coordination to pyochelin siderophore influences dynamics of FptA receptor from Pseudomonas aeruginosa: a molecular dynamics simulation study.

FptA is a TonB-dependent transporter that permits the high affinity binding and transport of Fe(III)-pyochelin complex across the outer membrane of Pseudomonas aeruginosa. Molecular dynamics simulations were employed to FptA receptor and its complexes with pyochelin, and co-crystallized Fe(III)-pyochelin-ethanediol and Fe(III)-pyochelin-water embedded in dilauroyl phosphatidyl choline bilayer for the evaluation of their structural and dynamical properties. The evaluation of properties of the receptor bound to pyochelin molecule and Fe(III)-pyochelin complexes helped to figure out the iron coordination effect on the receptor properties. Moreover, comparison of these four simulation systems revealed further information on the dynamical changes occurred in extracellular loops, in particular loop-7 corresponding to the missing amino acid residues including the close-by loop-8 that was largely affected by the metal coordination to pyochelin. The binding of iron to pyochelin molecule affected the overall structure of the receptor therefore, evaluation fo the gyration radii and hydrogen bonding were evaluated as well as analysis of the pore size were also carried out to understand the effect of metal coordination on the dynamics of the helices which form a kind of translocation channel to transport the siderophore across the FptA protein into the periplasmic space. The properties of each component of the molecular systems were therefore observed to be perturbed by the incorporation of iron to the pyochelin molecule thus demonstrating that the bacteria use its receptor to abstract and transport iron from extracellular environment for its survival and that was made possible to understand at the molecular level through successful implementation of molecular dynamics simulations.

Solvation of cholesterol in different solvents: a molecular dynamics simulation study.

To the best of our knowledge, molecular dynamics simulations of an isolated cholesterol immersed in four different solvents of varying polarity, such as water, methanol, dimethyl sulfoxide and benzene, were reported for the first time to gain insights into the structural and dynamical properties. The study was mainly focused on the evaluation of solvation of cholesterol with respect to its hydrophilic and hydrophobic structural components in the form of respective functional groups interacting with the solvents. Structural evaluations suggested that both hydrophilic and hydrophobic groups of cholesterol were interacting with the solvents, in particular methanol and dimethyl sulfoxide, which presented both types of interactions that are polar and non-polar. On the other hand, the highly polar water and non-polar benzene demonstrated extreme solvation behavior, since water was involved only in hydrogen bonding to the solute hydroxyl group and non-polar benzene formed strong van der Waals interactions only. Furthermore, the hydrophobic effect of cholesterol was also analyzed mainly in polar solvents, as the effect was more pronounced in the polar environment thereby preventing the solvent mobility in the solvation layer(s). The dynamical properties in terms of lateral diffusion and hydrogen bond dynamics as well as free energies of solvation also corroborated the findings based on the structural data and the hydrophobic character of cholesterol was later quantified by the computation of the averaged solvent accessible surface area. The polarity effect of the solvents on the aggregation property of cholesterol was further investigated, which is of big concern from the clinical point of view due to its major role in cardiovascular ailments. It was another major finding of the present study that aggregation was shown to be facilitated by highly polar solvents like water.

Highly selective enrichment of phosphopeptides using poly(dibenzo-18-crown-6) as a solid-phase extraction material.

A poly(dibenzo-18-crown-6) was used as a new solid-phase extraction material for the selective enrichment of phosphopeptides. Isolation of phosphopeptides was achieved based on specific ionic interactions between poly(dibenzo-18-crown-6) and the phosphate group of phosphopeptides. Thus, a method was developed and optimized, including loading, washing and elution steps, for the selective enrichment of phosphopeptides. To assess this potential, tryptic digest of three proteins (α- casein, ß-casein and ovalbumin) was applied on poly(dibenzo-18-crown-6). The nonspecific products were removed by centrifugation and washing. The spectrometric analysis was performed using matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. Highly selective enrichment of both mono- and multiphosphorylated peptides was achieved using poly(dibenzo-18-crown-6) as solid-phase extraction material with minimum interference from nonspecific compounds. Furthermore, evaluation of the efficiency of the poly(dibenzo-18-crown-6) was performed by applying the digest of egg white. Finally, quantum mechanical calculations were performed to calculate the binding energies to predict the affinity between poly(dibenzo-18-crown-6) and various ligands. The newly identified solid-phase extraction material was found to be a highly efficient tool for phosphopeptide recovery from tryptic digest of proteins.

Hydration facilitates oxygenation of hemocyanin: perspectives from molecular dynamics simulations.

Molecular dynamics simulations were applied to deoxy- and oxy-hemocyanins using newly developed force field parameters for the dicopper site to evaluate their structural and dynamical properties. Data obtained from the simulations provided information of the oxygenation effect on the active site and overall topology of the protein that was analyzed by root-mean-square deviations, b-factors, and dicopper coordination geometries. Domain I of the protein was found to demonstrate higher flexibility with respect to domain II because of the interfacial rotation between domain I and II that was further endorsed by computing correlative domain movements for both forms of the protein. The oxygenation effect on the overall structure of the protein or polypeptide subunit was further explored via gyration radii evaluated for the metal-binding domain and for the whole subunit. The evaluation of hydration dynamics was carried out to understand the water mediated role of amino acid residues of the solvent tunnel facilitating the entry of oxygen molecule to the dicopper site of hemocyanin.

Hydration of iron-porphyrins: ab initio quantum mechanical charge field molecular dynamics simulation study.

The ab initio quantum mechanical charge field molecular dynamics (QMCF-MD) simulation approach was successfully applied to Fe2+-P and Fe3+-P in water to evaluate their structural, dynamical and energetic properties. Based on the structural data, it was found that Fe2+-P accommodates one water molecule in the first coordination sphere of the Fe2+ ion including the four nitrogen atoms of the porphyrin system coordinating with central metal species. On the other hand, two water molecules were coordinated to Fe3+-P, thus forming a hexa-coordinated species. Comparison of dynamical properties such as the vibrational power spectrum and ligand mean residence times to other metal-free porphyrin systems demonstrate the ions' influence on the hydration structure, enabling a characterisation of the strong interaction of the ions which greatly reduces the hydrogen bonding potential of the complex. The association of water molecules with the metal ions in both solutes was quantified by computing the free energy of binding obtained via the potential of mean force. This further confirmed the strong association of water to the metal ions which was conversely weak as inferred from the energetic data for the Fe2+-P system.

Syntheses of 4,6-dihydroxypyrimidine diones, their urease inhibition, in vitro, in silico, and kinetic studies.

A library of 4,6-dihydroxypyrimidine diones (1-35) were synthesized and evaluated for their urease inhibitory activity. Structure-activity relationships, and mechanism of inhibition were also studied. All compounds were found to be active with IC50 values between 22.6±1.14-117.4±0.73µM, in comparison to standard, thiourea (IC50=21.2±1.3µM). Kinetics studies on the most active compounds 2-7, 16, 17, 28, and 33 were performed to investigate their modes of inhibition, and dissociation constants Ki. Compounds 2, 3, 7, 16, 28, and 33 were found to be mixed-type of inhibitors with Ki values in the range of 7.91±0.024-13.03±0.013µM, whereas, compounds 4-6, and 17 were found to be non-competitive inhibitors with Ki values in the range of 9.28±0.019-13.05±0.023µM. In silico study was also performed, and a good correlation was observed between experimental and docking studies. This study is continuation of our previously reported urease inhibitory activity of pyrimidine diones, representing potential leads for further research as possible treatment of diseases caused by ureolytic bacteria.

Synthesis and in vitro α-chymotrypsin inhibitory activity of 6-chlorobenzimidazole derivatives.

A library of benzimidazole derivatives 1-20 were synthesized, and studied for their α-chymotrypsin (α-CT) inhibitory activity in vitro. Kinetics and molecular docking studies were performed to identify the type of inhibition. Compound 1 was found to be a good inhibitor of α-chymotrypsin enzyme (IC50=14.8±0.1µM, Ki=16.4µM), when compared with standard chymostatin (IC50=5.7±0.13µM). Compounds 2-8, 15, 17, and 18 showed significant inhibitory activities. All the inhibitors were found to be competitive inhibitors, except compound 17, which was a mixed type inhibitor. The substituents (R) in para and ortho positions of phenyl ring B, apparently played a key role in the inhibitory potential of the series. Compounds 1-20 were also studied for their cytotoxicity profile by using 3T3 mouse fibroblast cells and compounds 3, 5, 6, 8, 12-14, 16, 17, 19, and 20 were found to be cytotoxic. Molecular docking was performed on the most active members of the series in comparison to the standard compound, chymostatin, to identify the most likely binding modes. The compounds reported here can serve as templates for further studies for new inhibitors of α-chymotrypsin and other chymotrypsin-like serine proteases enzymes.

2-Arylquinazolin-4(3H)-ones: A new class of α-glucosidase inhibitors.

Twenty-five derivatives of 2-arylquinazolin-4(3H)-ones (1-25) were evaluated for their yeast (Saccharomyces cerevisiae) α-glucosidase inhibitory activities. All synthetic compounds, except 1 and 6, were found to be several hundred fold more active (IC50 values in the range of 0.3±0.01-117.9±1.76µM), than the standard drug, acarbose (IC50=840±1.73µM). The enzyme kinetic studies on the most active compounds 12, 4, 19, and 13 were performed for the determination of their modes of inhibition and dissociation constants Ki. Study of the modes of inhibition of compounds 12, and 4 were also performed using molecular modeling techniques. In brief, current study identifies a novel class of α-glucosidase inhibitors which can be further studied for the treatment of hyperglycemia and obesity.

Hydration of the cyanide ion: an ab initio quantum mechanical charge field molecular dynamics study.

This paper presents an ab initio quantum mechanical charge field molecular dynamics simulation study of the cyanide anion (CN(-)) in aqueous solution where hydrogen bond formation plays a dominant role in the hydration process. Preferential orientation of water hydrogens compared to oxygen atoms was quantified in terms of radial, angular as well as coordination number distributions. All structural results indicate that the water hydrogens are attracted towards CN(-) atoms, thus contributing to the formation of the hydration layer. Moreover, a clear picture of the local arrangement of water molecules around the ellipsoidal CN(-) ion is provided via angular-radial distribution and spatial distribution functions. Apart from the structural analysis, the evaluation of water dynamics in terms of ligand mean residence times and H-bond correlation functions indicates the weak structure making capacity of the CN(-) ion. The similar values of H-bond lifetimes obtained for the NHwat and CHwat bonds indicate an isokinetic behaviour of these H-bonds, since there is a very small difference in the magnitude of the lifetimes. On the other hand, the H-bond lifetimes between water molecules of the hydration shell, and between solute and solvent evidence the slightly stable hydration of the CN(-). Overall, the H-bonding dominates in the hydration process of the cyanide anion enabling it to become soluble in the aqueous environment associated to chemical and biological processes.

Dynamics of ligand exchange mechanism at Cu(II) in water: an ab initio quantum mechanical charge field molecular dynamics study with extended quantum mechanical region.

Ab initio quantum mechanical charge field molecular dynamics (QMCF-MD) were successfully applied to Cu(II) embedded in water to elucidate structure and to understand dynamics of ligand exchange mechanism. From the simulation studies, it was found that using an extended large quantum mechanical region including two shells of hydration is required for a better description of the dynamics of exchanging water molecules. The structural features characterized by radial distribution function, angular distribution function and other analytical parameters were consistent with experimental data. The major outcome of this study was the dynamics of exchange mechanism and reactions in the first hydration shell that could not be studied so far. The dynamical data such as mean residence time of the first shell water molecules and other relevant data from the simulations are close to the results determined experimentally. Another major characteristic of hydrated Cu(II) is the Jahn-Teller distortion which was also successfully reproduced, leading to the final conclusion that the dominating aqua complex is a 6-coordinated species. The ab initio QMCF-MD formalism proved again its capabilities of unraveling even ambiguous properties of hydrated species that are far difficult to explore by any conventional quantum mechanics/molecular mechanics (QM/MM) approach or experiment.

Structure and dynamics of 5-lipoxygenase and its complexes - a molecular dynamics simulation study.

Molecular dynamics simulations were applied to human 5-LOX to obtain detailed information on its structure and dynamics with and without ligands. The dynamical properties evaluated based on root mean square deviations, root mean square fluctuations and secondary structure prediction helped decipher the contrast dynamic behavior of the systems pointing toward the ligand binding effect. The ligand binding to the protein also perturbed other properties of the protein such as the central bending of the protein and water coordination to the metal ion. The central bending in the protein was reported to be very significant that was associated with the allosteric modulation in the lipoxygenases; therefore, on a similar line, the central bending was evaluated in terms of hinge angle analysis which showed substantial bending between the C-terminal and the N-terminal domain via the linker residues which connects the two domains. On the other hand, the suspected water coordination to the metal ion in the protein was ruled out by computing the iron-water radial distribution function which showed that the water molecule was not found to be in the vicinity of the metal ion. Finally, the binding free energy was estimated for Zileuton and CAPE1 inhibitors bound to 5-LOX via the thermodynamic integration approach which showed that CAPE1 had a strong binding potential for the active site of the protein compared to Zileuton, and the free energy data correlated well with their IC50 values corresponding to the high inhibition potential of CAPE1 compared to Zileuton.

Interaction of gentamicin and gentamicin-AOT with poly-(lactide-co-glycolate) in a drug delivery system - density functional theory calculations and molecular dynamics simulation.

Gentamicin is used to treat brucellosis, an infectious disease caused by the Brucella species but the drug faces several issues such as low efficacy, instability, low solubility, and toxicity. It also has a very short half-life, therefore, requiring frequent dosing. Consequently, several other antibiotics are also being used for the treatment of brucellosis as a single dose as well as in combination with other antibiotics but none of these therapies are satisfactory. Nanoparticles in particular polymer-based ones utilizing polymers that are biodegradable and biocompatible for instance PLGA are a method of choice to overcome such drug delivery issues and enable potential targeted delivery. The current study focuses on the evaluation of the structural and dynamical properties of a drug-polymer system consisting of gentamicin drug and PLGA polymer nanoparticles in the water representing a targeted drug delivery system for the treatment of brucellosis. For this purpose, all-atom molecular dynamics simulations were carried out on the drug-polymer systems in the absence and presence of the surfactant bis(2-Ethylhexyl) sulfosuccinate (AOT) to determine the structural and dynamical properties as well as the effect of the surfactant on these properties. We also investigated systems in which the polymer constituents were in the form of monomeric units toward decoupling the primary interactions of the monomer units and polymer effects. The simulation results explain the nature of the interactions between the drug and the polymer as well as transport properties in terms of drug diffusion coefficients, which characterize the molecular behavior of gentamicin-polymer nanoparticles for use in brucellosis.

Deciphering Structural and Dynamical Properties of Hydrated Cobalt Porphyrins via Ab Initio Quantum Mechanical Charge Field Molecular Dynamics Simulation.

The present study successfully implemented the ab initio quantum mechanical charge field molecular dynamics (QMCF MD) formalism for the investigation of structural and dynamical properties of hydrated cobalt-porphyrin complexes. Considering the significance of cobalt ions in biological systems (for instance, vitamin B12), which reportedly incorporate cobalt ions in a d6, low spin, +3 state chelated in the corrin ring, an analog of porphyrin, the current study is focused on cobalt in the oxidation states +2 and +3 bound to the parent porphyrin lead structures embedded in an aqueous solution. These cobalt-porphyrin complexes were investigated in terms of their structural and dynamical properties at the quantum chemical level. The structural attributes of these hydrated complexes revealed the contrasting features of the water binding to these solutes, including a detailed evaluation of the associated dynamics. The study also yielded notable findings in regard to the respective electronic configurations vs coordination, which suggested that Co(II)-POR possesses a 5-fold square pyramidal coordination geometry in an aqueous solution containing the metal ion coordinating to four nitrogen atoms of the porphyrin ring and one axial water as the fifth ligand. On the other hand, high-spin Co(III)-POR was hypothesized to be more stable due to the smaller size-to-charge ratio of the cobalt ion, but the high-spin complex demonstrated unstable structural and dynamical behavior. However, the corresponding properties of the hydrated Co(III)LS-POR revealed a stable structure in an aqueous solution, thus suggesting the Co(III) ion to be in a low-spin state when bound to the porphyrin ring. Moreover, the structural and dynamical data were augmented by computing the free energy of water binding to the cobalt ions and the solvent-accessible surface area, which provide further information on thermochemical properties of the metal-water interaction and the hydrogen bonding potential of the porphyrin ring in these hydrated systems.

Genomic analysis of Chryseobacterium indologenes and conformational dynamics of the selected DD-peptidase.

Chrysobacterium indologenes is an emerging MDR pathogen that belongs to the family Flavobacteriaceae. The genome of the C. indologenes, isolated from the nephrotic patient, was sequenced through Illumina MiSeq. The pangenomics of available 56 C. indologenes strains using BPGA revealed an open pangenome (n=5553 CDS), core genome (2141), and accessory genome (2013). The CEG/DEG database identified 662 essential genes that drastically reduced to 68 genes after non-homology analyses towards human and gut microbiome. Further filtering the data for other drug target prioritizing parameters resulted in 32 putative targets. Keeping in view the crucial role played in cell wall biosynthesis, dacB was selected as the final target that encodes D-alanyl-d-alanine carboxypeptidase/endopeptidase (DD-peptidase). The 3D structure of dacB was modelled and rendered to docking analyses against two compound libraries of African plants (n=6842) and Tibetan medicines (n=52). The ADMET profiling exhibited the physicochemical properties of final compounds. The MD simulations showed the stability of inhibitor-DD-peptidase complex and interactions in terms of RMSD, RMSF, binding free energy calculation and H-bonding. We propose that the novel compounds Leptopene and ZINC95486338 from our findings might be potent DD-peptidase inhibitors that could aid in the development of new antibiotic-resistant therapy for the emerging MDR C. indologenes.

Dynamic changes in the secondary structure of ECE-1 and XCE account for their different substrate specificities.

BACKGROUND: X-converting enzyme (XCE) involved in nervous control of respiration, is a member of the M13 family of zinc peptidases, for which no natural substrate has been identified yet. In contrast, it's well characterized homologue endothelin-converting enzyme-1 (ECE-1) showed broad substrate specificity and acts as endopeptidase as well as dipeptidase. To explore the structural differences between XCE and ECE-1, homology model of XCE was built using the complex structure of ECE-1 with phosphoramidon (pdb-id: 3DWB) as template. Phosphoramidon was docked into the binding site of XCE whereas phosphate oxygen of the inhibitor was used as water molecule to design the apo forms of both enzymes. Molecular dynamics simulation of both enzymes was performed to analyze the dynamic nature of their active site residues in the absence and presence of the inhibitor. RESULTS: Homology model of XCE explained the role of non-conserved residues of its S2' subsite. Molecular dynamics (MD) simulations identified the flexible transitions of F149/I150, N566/N571, W714/W719, and R145/R723 residues of ECE-1/XCE for the strong binding of the inhibitor. Secondary structure calculations using DSSP method reveals the folding of R145/R723 residue of ECE-1/XCE into ß-sheet structure while unfolding of the S2' subsite residues in aECE-1 and sustained compact folding of that of aXCE. The results evaluated are in good agreement with available experimental data, thus providing detailed molecular models which can explain the structural and specificities differences between both zinc peptidases. CONCLUSIONS: Secondary structure changes of both enzymes during the simulation time revealed the importance of ß-sheet structure of R145/R723 for its binding with the terminal carboxylate group of the inhibitor. Unfolding of the α-helix comprising the S2' subsite residues in aECE-1 correlate well with its endopeptidase activity while their compact folding in aXCE may account for the inactivity of the enzyme towards large C-terminal containing substrates.

Activation of TrkB receptors by NGFß mimetic peptide conjugated polymersome nanoparticles.

Activation of tyrosine kinase receptor B (TrkB), a neurotrophin receptor, has been shown to increase neuronal cell survival and promote regeneration. Stimulation of the TrkB receptor by neurotrophic growth factors has been identified as a possible therapeutic target for the treatment of neurodegenerative disorders. However, growth factor delivery is problematic because of a short half-life in vivo. We have conjugated hNgf-EE, a short peptide mimetic of NGFß to the surface of polymersome nanoparticles and shown that they are capable of activating the TrkB receptor in vitro in the SHSY-G7 cell line. We propose that polymersomes could act as a scaffold for the delivery of TrkB activating moieties and that the polymersome size and polyethylene glycol surface have been shown to increase in vivo retention time. These multifunctional nanoparticles have potential for the treatment of neurodegenerative disorders by TrkB activation. From the ClinicaL Editor: Tyrosine kinase receptor B activation has been shown to promote regeneration and survival of neurons. However, growth factor delivery to stimulate these receptors remains problematic. The authors demonstrate that a peptide mimetic of NGFß conjugated to the surface of polymersome nanoparticles is capable of activating the TrkB receptors. These nanoparticles may offer a novel treatment strategy for a variety of neurodegenerative disorders.

Investigation of the structural and dynamical properties of human uncoupling protein 2 through molecular dynamics simulations.

Uncoupling protein 2 (UCP2) is an integral membrane protein that belongs to the family of mitochondrial anion carrier proteins. The absence of human UCP2 structure, lack of understanding of Cl- ion transport mechanism in the UCP2 and the associated biological functions motivated us to model the protein and investigate its structural and dynamical properties in a realistic mitochondrial lipid membrane system. The lipid-protein and protein-protein interactions were probed since they were found to be responsible for the conformational changes of the transmembrane (TM) helices which are involved in facilitating Cl- ion transport. Here, we employed multiscale molecular dynamics simulations including unbiased and biased MD for the investigation of the transport pathway in hUCP2 and interactions of the ion with TM helices within a membrane environment. We initially validated the hUCP2 model in the lipid membrane and then explored the transport pathway of Cl- ion and its interaction with positive residues of TM2 helix that have been reported to play a major role in the Cl- ion transport along with other TM helices of the protein. The simulation results suggest that the TM2 helix plays an important role in the formation of a stable ion channel due to the presence of arginine residues, in particular Arg88 which was found to be a key residue to maintain the channel pore through which the movement of Cl- ions occurs. Based on the results, it can be said that the study provides an atomic-level description of the Cl- ion transport mechanism in hUCP2 embedded in the mitochondrial lipid membrane.

Dynamics of metal binding and mutation in yybP-ykoY riboswitch of Lactococcus lactis.

Riboswitch is a regulatory segment of messenger RNA (mRNA), which by binding to various cellular metabolites regulates the activity of mRNA via modulating transcription, translation, alternative splicing, and stability of the mRNA. yybP-ykoY riboswitch of Lactococcus lactis, which is present upstream of the yoaB gene, functions as a Mn2+-specific genetic ON-switch, and modulates expression of proteins which are significant for Mn2+ homeostasis. The P1.1 switch helix of the aptamer domain of the riboswitch contains an intrinsic transcription terminator structure, which gets stabilized with Mn2+ binding and causes disruption of terminator structure and allows the continuation of transcription. The current research work involved the evaluation of structural and dynamical properties of the yybP-ykoY riboswitch of L. lactis in its Mn2+-free, Mn2+-bound (wild-type), and Mn2+-bound mutant (A41U) states by applying molecular dynamics simulations. Based on the simulations, the effects of Mn2+ absence and A41U mutation were evaluated on the structure and dynamics of the riboswitches followed by the computation of the free energy of metal binding in the wild-type and the mutant riboswitches. The simulation results provided insights into the properties of the riboswitch with the focus on the dynamics of the P1.1 switch helix, and the manganese binding site designated as MB site, as well as the relative stability of the wild-type and the mutant riboswitches, which helped to understand the structural and dynamical role of the metal ion involved in the function of Mn2+-sensing riboswitch.

Interaction of Phthalates with Lipid Bilayer Membranes.

Phthalates are esters of phthalic acid, widely used as additives in the manufacture of plastics. They are not covalently linked to polymer chains and can easily leach out, disperse in the environment, and get into contact with living organisms. Several short chain phthalates are classified as endocrine disruptors or hormonal active agents, and have also been reported to promote various kinds of cancer. However, the biological effects of longer chain analogues are less well known. Moreover, little is known on the permeation of phthalates and their metabolites through biological membranes and on their effects on the physical properties of membranes. Here we explore the interaction of a group of phthalates and their main metabolites with model biological membranes. We focus on three industrially relevant phthalates, with acyl chains of different sizes, and their monoester metabolites. We use molecular dynamics simulations to predict the distribution in model membranes, as well as permeabilities and effects on the structural, dynamic, and elastic properties of the membranes. We find that alterations of membrane properties are significant and only weakly affected by the size of acyl chains, suggesting that modifications of molecular size may not be sufficient to reduce the impact of this class of molecules on the environment and health.