Abiological catalysis by myoglobin mutant with a genetically incorporated unnatural amino acid.

To inculcate biocatalytic activity in the oxygen-storage protein myoglobin (Mb), a genetically engineered myoglobin mutant H64DOPA (DOPA = L-3,4-dihydroxyphenylalanine) has been created. Incorporation of unnatural amino acids has already demonstrated their ability to accomplish many non-natural functions in proteins efficiently. Herein, the presence of redox-active DOPA residue in the active site of mutant Mb presumably stabilizes the compound I in the catalytic oxidation process by participating in an additional hydrogen bonding (H-bonding) as compared to the WT Mb. Specifically, a general acid-base catalytic pathway was achieved due to the availability of the hydroxyl moieties of DOPA. The reduction potential values of WT (E° = -260 mV) and mutant Mb (E° = -300 mV), w.r.t. Ag/AgCl reference electrode, in the presence of hydrogen peroxide, indicated an additional H-bonding in the mutant protein, which is responsible for the peroxidase activity of the mutant Mb. We observed that in the presence of 5 mM H2O2, H64DOPA Mb oxidizes thioanisole and benzaldehyde with a 10 and 54 folds higher rate, respectively, as opposed to WT Mb. Based on spectroscopic, kinetic, and electrochemical studies, we deduce that DOPA residue, when present within the distal pocket of mutant Mb, alone serves the role of His/Arg-pair of peroxidases.

Getting Drugs through Small Pores: Exploiting the Porins Pathway in Pseudomonas aeruginosa.

Understanding molecular properties of outer membrane channels of Gram-negative bacteria is of fundamental significance as they are the entry point of polar antibiotics into bacteria. Outer membrane proteomics revealed OccK8 (OprE) to be among the five most expressed substrate specific channels of the clinically important Pseudomonas aeruginosa. The high-resolution X-ray structure and electrophysiology highlighted a very narrow pore. However, experimental in vitro methods showed the transport of natural amino acids and antibiotics, among them ceftazidime. We used molecular dynamics simulations to reveal the importance of the physicochemical properties of ceftazidime in modulating the translocation through OccK8, proposing a structure-function relationship. As in general porins, the internal electric field favors the translocation of polar molecules by gainful energy compensation in the central constriction region. Importantly, the comparatively narrow OccK8 pore can undergo a substrate-induced expansion to accommodate relatively large-sized substrates.

Rationalizing the permeation of polar antibiotics into Gram-negative bacteria.

The increasing level of antibiotic resistance in Gram-negative bacteria, together with the lack of new potential drug scaffolds in the pipeline, make the problem of infectious diseases a global challenge for modern medicine. The main reason that Gram-negative bacteria are particularly challenging is the presence of an outer cell-protecting membrane, which is not present in Gram-positive species. Such an asymmetric bilayer is a highly effective barrier for polar molecules. Several protein systems are expressed in the outer membrane to control the internal concentration of both nutrients and noxious species, in particular: (i) water-filled channels that modulate the permeation of polar molecules and ions according to concentration gradients, and (ii) efflux pumps to actively expel toxic compounds. Thus, besides expressing specific enzymes for drugs degradation, Gram-negative bacteria can also resist by modulating the influx and efflux of antibiotics, keeping the internal concentration low. However, there are no direct and robust experimental methods capable of measuring the permeability of small molecules, thus severely limiting our knowledge of the molecular mechanisms that ultimately control the permeation of antibiotics through the outer membrane. This is the innovation gap to be filled for Gram-negative bacteria. This review is focused on the permeation of small molecules through porins, considered the main path for the entry of polar antibiotics into Gram-negative bacteria. A fundamental understanding of how these proteins are able to filter small molecules is a prerequisite to design/optimize antibacterials with improved permeation. The level of sophistication of modern molecular modeling algorithms and the advances in new computer hardware has made the simulation of such complex processes possible at the molecular level. In this work we aim to share our experience and perspectives in the context of a multidisciplinary extended collaboration within the IMI-Translocation consortium. The synergistic combination of structural data, in vitro assays and computer simulations has proven to give new insights towards the identification and description of physico-chemical properties modulating permeation. Once similar general rules are identified, we believe that the use of virtual screening techniques will be very helpful in searching for new molecular scaffolds with enhanced permeation, and that molecular modeling will be of fundamental assistance to the optimization stage.

Molecular Properties of Astaxanthin in Water/Ethanol Solutions from Computer Simulations.

Astaxanthin (AXT) is a reference model of xanthophyll carotenoids, which is used in medicine and food industry, and has potential applications in nanotechnology. Because of its importance, there is a great interest in understanding its molecular properties and aggregation mechanism in water and mixed solvents. In this paper, we report a novel model of AXT for molecular dynamics simulation. The model is used to estimate different properties of the molecule in pure solutions and in water/ethanol mixtures. The calculated diffusion coefficients of AXT in pure water and ethanol are (3.22 ± 0.01) × 10(-6) cm(2) s(-1) and (2.7 ± 0.4) × 10(-6) cm(2) s(-1), respectively. Our simulations also show that the content of water plays a clear effect on the morphology of the AXT aggregation in water/ethanol mixture. In up to 75% (v/v) water concentration, a loosely connected network of dimers and trimers and two-dimensional array structures are observed. At higher water concentrations, AXT molecules form more compact three-dimensional structures that are preferentially solvated by the ethanol molecules. The ethanol preferential binding and the formation of a well connected hydrogen bonding network on these AXT clusters, suggest that such preferential solvation can play an important role in controlling the aggregate structure.

Molecular basis of substrate translocation through the outer membrane channel OprD of Pseudomonas aeruginosa.

The objective of this study is to identify the structural features governing the transport of molecules through nanometric channel proteins at a molecular level. Our focus is to come up with a precise understanding of the structure and dynamics of the outer membrane porin OprD of the Gram-negative bacterium Pseudomonas aeruginosa by studying the translocation of natural amino acid residues/substrates through it. We used in silico electrophysiology and metadynamics simulation techniques as they can provide precise information on the molecule/channel interactions at the atomic scale that allows testing quantitative structure-function relationships. We performed our simulations on the whole OprD protein, with all loops modelled and without any constraints to keep the channel open. Dynamics of both internal and external loops and the polar nature of the eyelet region play important roles in modulating the translocation of molecules through OprD by creating two alternative paths for translocation. All positive residues take the main path upon binding in the negative pocket just above the constriction region. The same factor is unfavourable for negative substrates and hence they have a relatively high barrier for translocation. Differently, neutral substrates do not show any specificity and they can follow the two alternative paths.

Theoretical study of binding and permeation of ether-based polymers through interfaces.

We present a molecular dynamics simulation study on the interactions of poly(ethylene oxide) (PEO), poly(propylene oxide) (PPO), and their ABA-type block copolymer, poloxamers, at water/n-heptane and 1,2-dimyristoyl-sn-glycero-3-phospatidycholine (DMPC) lipid bilayer/water interfaces. The partition coefficients in water/1-octanol of the linear polyethers up to three monomers were calculated. The partition coefficients evidenced a higher hydrophobicity of the PPO in comparison to PEO. At the water/n-heptane interface, the polymers tend to adopt elongated conformations in agreement with similar experimental ellipsometry studies of different poloxamers. In the case of the poloxamers at the n-heptane/water interface, the stronger preference of the PPO block for the hydrophobic phase resulted in bottle-brush-type polymer conformations. At lipid bilayer/water interface, the PEO polymers, as expected from their hydrophilic nature, are weakly adsorbed on the surface of the lipid bilayer and locate in the water phase close to the headgroups. The free energy barriers of permeation calculated for short polymer chains suggest a thermodynamics propensity for the water phase that increase with the chain length. The lower affinity of PEO for the hydrophobic interior of the lipid bilayer resulted in the spontaneous expulsion within the simulation time. On the contrary, PPO chains and poloxamers have a longer residence time inside the bilayer, and they tend to concentrate in the tail region of the bilayer near the polar headgroups. In addition, polymers with PPO unit length comparable to the thickness of the hydrophobic region of the bilayer tend to span across the bilayer.

Interaction of Curcumin with PEO-PPO-PEO block copolymers: a molecular dynamics study.

Curcumin, a naturally occurring drug molecule, has been extensively investigated for its various potential usages in medicine. Its water insolubility and high metabolism rate require the use of drug delivery systems to make it effective in the human body. Among various types of nanocarriers, block copolymer based ones are the most effective. These polymers are broadly used as drug-delivery systems, but the nature of this process is poorly understood. In this paper, we propose a molecular dynamics simulation study of the interaction of Curcumin with block copolymer based on polyethylene oxide (PEO) and polypropylene oxide (PPO). The study has been conducted considering the smallest PEO and PPO oligomers and multiple chains of the block copolymer Pluronic P85. Our study shows that the more hydrophobic 1,2-dimethoxypropane (DMP) molecules and PPO block preferentially coat the Curcumin molecule. In the case of the Pluronic P85, simulation shows formation of a drug-polymer aggregate within 50 ns. This process leaves exposed the PEO part of the polymers, resulting in better solvation and stability of the drug in water.

Understanding the interaction of block copolymers with DMPC lipid bilayer using coarse-grained molecular dynamics simulations.

In this paper, we present a computational model of the adsorption and percolation mechanism of poloxamers (poly(ethylene oxide) (PEO) and poly(propylene oxide) (PPO) triblock copolymers) across a 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) lipid bilayer. A coarse-grained model was used to cope with the long time scale of the percolation process. The simulations have provided details of the interaction mechanism of Pluronics with lipid bilayer. In particular, the results have shown that polymer chains containing a PPO block with a length comparable to the DMPC bilayer thickness, such as P85, tends to percolate across the lipid bilayer. On the contrary, Pluronics with a shorter PPO chain, such as L64 and F38, insert partially into the membrane with the PPO block part while the PEO blocks remain in water on one side of the lipid bilayer. The percolation of the polymers into the lipid tail groups reduces the membrane thickness and increases the area per lipid. These effects are more evident for P85 than L64 or F38. Our findings are qualitatively in good agreement with published small-angle X-ray scattering experiments that have evidenced a thinning effect of Pluronics on the lipid bilayer as well as the role of the length of the PPO block on the permeation process of the polymer through the lipid bilayer. Our theoretical results complement the experimental data with a detailed structural and dynamic model of poloxamers at the interface and inside the lipid bilayer.

Molecular dynamics simulation study of solvent effects on conformation and dynamics of polyethylene oxide and polypropylene oxide chains in water and in common organic solvents.

In this paper, the conformation and dynamics properties of polyethylene oxide (PEO) and polypropylene oxide (PPO) polymer chains at 298 K have been studied in the melt and at infinite dilution condition in water, methanol, chloroform, carbon tetrachloride, and n-heptane using molecular dynamics simulations. The calculated density of PEO melt with chain lengths of n = 2, 3, 4, 5 and, for PPO, n = 7 are in good agreement with the available experimental data. The conformational properties of PEO and PPO show an increasing gauche preference for the O-C-C-O dihedral in the following order water>methanol>chloroform>carbon tetrachloride = n-heptane. On the contrary, the preference for trans conformation has a maximum in carbon tetrachloride and n-heptane followed in the order by chloroform, methanol, and water. The PEO conformational preferences are in qualitative agreement with results of NMR studies. PEO chains formed different types of hydrogen bonds with polar solvent molecules. In particular, the occurrence of bifurcated hydrogen bonding in chloroform was also observed. Radii of gyration of PEO chains of length larger than n = 9 monomers showed a good agreement with light scattering data in water and in methanol. For the shorter chains the observed deviations are probably due to the enhanced hydrophobic effects caused by the terminal methyl groups. For PEO the fitting of end-to-end distance distributions with the semi-flexible chain model at 298 K provided persistence lengths of 0.375 and 0.387 nm in water and methanol, respectively. Finally, the radius of gyration of Pluronic P85 turned out to be 2.25 ± 0.4 nm at 293 K in water in agreement with experimental data.

Diffusion of 1,2-dimethoxyethane and 1,2-dimethoxypropane through phosphatidycholine bilayers: a molecular dynamics study.

In this paper, a theoretical study of 1,2-dimethoxyethane (DME) and 1,2-dimethoxypropane (DMP) at water/n-heptane and 1,2-dimyristoyl-sn-glycero-3-phospatidycholine (DMPC) lipid bilayer/water interfaces using the umbrella sampling method is reported. Recently proposed GROMOS96/OPLS compatible models for DME and DMP have been used for the simulation studies. The percolation free energy barrier of one DME and DMP molecule from water to n-heptane phase calculated using the umbrella sampling method turned out to be equal to ~18.5 kJ/mol and ~6 kJ/mol, respectively. In the case of the DMPC lipid bilayer, overall free energy barriers of ~20 kJ/mol and ∼12 kJ/mol were obtained for DME and DMP, respectively. The spontaneous diffusion of DME and DMP in the lipid bilayer has also been investigated using unconstrained molecular dynamics simulations at the water/DMPC interface and inside the lipid bilayer. As expected from the estimated percolation barriers, simulation results show that DME, contrary to DMP, spontaneously diffuse into the aqueous solution from the lipid interior. In addition, simulations with multiple DME or DMP molecules at the interface show spontaneous diffusion within 50 ns inside the DMPC layer only for DMP.

Structure and dynamics of 1,2-dimethoxyethane and 1,2-dimethoxypropane in aqueous and non-aqueous solutions: a molecular dynamics study.

Herein, we report a comparative modelling study of 1,2-dimethoxyethane (DME) and 1,2-dimethoxypropane (DMP) at 298 K and 318 K in the liquid state, water mixtures, and at infinite dilution condition in water, methanol, carbon tetrachloride, and n-heptane. Both DME and DMP are united-atom models compatible with GROMOS∕OPLS force fields. Calculated thermodynamic and structural properties of the pure DME and DMP liquids resulted in excellent agreement with the experimental data. In aqueous solutions, densities, diffusion coefficients, and concentration dependent conformers of DME, were in agreement with experimental data. The calculated free energy of solvation (ΔG(hyd)) at 298 K is equal to -22.1 ± 0.8 kJ mol(-1) in good agreement with the experimental value of 20.2 kJ mol(-1). In addition, the free energy of solvation of DME in non-aqueous solvents follows the trend methanol ≈ water < carbon tetrachloride < n-heptane, consistently with the dielectric constant of the solvents. On contrary, the presence of an extra methyl group on chiral carbon makes DMP less soluble than DME in water (ΔG(hyd) = -16.0 ± 1.1 kJ mol(-1)) but more soluble in non-polar solvents as n-heptane. Finally, for the DMP the chiral discrimination of the two enantiomers was calculated as solvation free energy difference of one DMP isomer in the solution of the other. The obtained value of ΔΔG(RS) = -3.7 ± 1.4 kJ mol(-1) indicates a net chiral discrimination of the two enantiomers.