Effects of water content on the tetrahedral intermediate of chymotrypsin - trifluoromethylketone in polar and non-polar media: observations from molecular dynamics simulation

The work uses MD simulation to study effects of five water contents (3 %, 10 %, 20 %, 50 %, 100 %?v/v) on the tetrahedral intermediate of chymotrypsin - trifluoromethyl ketone in polar acetonitrile and non-polar hexane media. The water content induced changes in the structure of the intermediate, solvent distribution and H-bonding are analyzed in the two organic media. Our results show that the changes in overall structure of the protein almost display a clear correlation with the water content in hexane media while to some extent U-shaped/bell-shaped dependence on the water content is observed in acetonitrile media with a minimum/maximum at 10–20 % water content. In contrast, the water content change in the two organic solvents does not play an observable role in the stability of catalytic hydrogen-bond network, which still exhibits high stability in all hydration levels, different from observations on the free enzyme system [Zhu L, Yang W, Meng YY, Xiao X, Guo Y, Pu X, Li M (2012) J Phys Chem B 116(10):3292–3304]. In low hydration levels, most water molecules mainly distribute near the protein surface and an increase in the water content could not fully exclude the organic solvent from the protein surface. However, the acetonitrile solvent displays a stronger ability to strip off water molecules from the protein than the hexane. In a summary, the difference in the calculated properties between the two organic solvents is almost significant in low water content (<10 %) and become to be small with increasing water content. In addition, some structural properties at 10?~?20 %?v/v hydration zone, to large extent, approach to those in aqueous solution.

Theoretical study of the solvatochromism of a donor-acceptor bithiophene

The solvation and the solvatochromic behavior of the 5-(methylthio)-5′-nitro-2,2′-bithiophene 1 in diethyl ether, dichloromethane, acetonitrile, methanol and formamide was theoretically investigated with an iterative molecular and quantum mechanics (QM/MM) approach. Calculated longest-wavelength solvatochromic absorption band of 1, obtained as averages of statistically uncorrelated configurations, including the solute and explicit solvent molecules of the first and second solvation layer, were in excellent agreement with the experimental results.

Effect of non aqueous solvent on structural stability of α-amylase: A cost-effective prospective for protein stabilization

The aim of the present study is to assess the effect of non-aqueous organic solvent on structural stability, molecular integrity and structure of α-amylase. The activity and thermal stability of the enzyme was measured before and after treatment with non polar solvent (i.e. hexane). The activity was found to be marginally affected and thermal stability was found to be significantly increased after treatment with hexane. The enzyme was found to be more resistant to thermal inactivation in hexane compared to in an aqueous buffer. The fluorescence measurement indicated a blue shift of 3 nm in the emission maximum (λ_max) probably due to a minor change in the polarity of aromatic amino acid residues after treatment with a non-aqueous solvent. Assessment of thermal denaturation profile, 1-anilino-8-naphthalene-sulfonate (ANS) binding and acrylamide quenching of the enzyme suggested an increase in the molecular integrity and overall stability of the enzyme after treatment with hexane. However, these entire molecular events were not accompanied by any major change in the secondary structure. Our findings suggest that treatment of proteins or enzymes in non-aqueous solvents could be an attractive and cost-effective strategy to improve their structural stability without compromising their biological functions.     

Molecular dynamics analysis of a series of 22 potential farnesyltransferase substrates containing a CaaX-motif

Protein farnesyltransferase (FTase) is an important target in many research fields, more markedly so in cancer investigation since several proteins known to be involved in human cancer development are thought to serve as substrates for FTase and to require farnesylation for proper biological activity. Several FTase inhibitors (FTIs) have advanced into clinical testing. Nevertheless, despite the progress in the field several functional and mechanistic doubts on the FTase catalytic activity have persisted. This work provides some crucial information on this important enzyme by describing the application of molecular dynamics simulations using specifically designed molecular mechanical parameters for a variety of 22 CaaX peptides known to work as natural substrates or inhibitors for this enzyme. The study involves a comparative analysis of several important molecular aspects, at the mechanistic level, of the behavior of substrates and inhibitors at the dynamic level, including the behavior of the enzyme and peptides, as well as their interaction, together with the effect of the solvent. Properties evaluated include the radial distribution function of the water molecules around the catalytically important zinc metal atom and cysteine sulfur of CaaX, the conformations of the substrate and inhibitor and the corresponding RMSF values, critical hydrogen bonds, and several catalytically relevant distances. These results are discussed in light of recent experimental and computational evidence that provides new insights into the activity of this enzyme.

Theoretical evaluation of flotation performance of carboxyl hydroxamic acids with different number of polar groups on the surfaces of diaspore (010) and kaolinite (001)

The adsorption behaviors of three carboxyl hydroxamic acids on diaspore (010) and kaolinite (001) have been studied by density functional theory (DFT) and molecular dynamics (MD) method. The results indicated that carboxyl hydroxamic acids could adsorb on diaspore surface by ionic bonds and hydrogen bonds, and adsorb on kaolinite surface by hydrogen bonds. The models of carboxyl hydroxamic acids adsorbed on diaspore and kaolinite surfaces are proposed.

Determination of key receptor–ligand interactions of dopaminergic arylpiperazines and the dopamine D2 receptor homology model

Interest in structure-based G-protein-coupled receptor (GPCR) ligand discovery is huge, given that almost 30 % of all approved drugs belong to this category of active compounds. The GPCR family includes the dopamine receptor subtype D2 (D2DR), but unfortunately—as is true of most GPCRs—no experimental structures are available for these receptors. In this publication, we present the molecular model of D2DR based on the previously published crystal structure of the dopamine D3 receptor (D3DR). A molecular modeling study using homology modeling and docking simulation provided a rational explanation for the behavior of the arylpiperazine ligand. The observed binding modes and receptor–ligand interactions provided us with fresh clues about how to optimize selectivity for D2DR receptors.

A molecular dynamics study on sI hydrogen hydrate

A molecular dynamics simulation is carried out to explore the possibility of using sI clathrate hydrate as hydrogen storage material. Metastable hydrogen hydrate structures are generated using the LAMMPS software. Different binding energies and radial distribution functions provide important insights into the behavior of the various types of hydrogen and oxygen atoms present in the system. Clathrate hydrate cages become more stable in the presence of guest molecules like hydrogen.

Enzymes in non-aqueous solvents

For a variety of reasons including increased recognition of the large degree of association, by non-polar interaction, of enzymes with other cellular components such as membranes, enzymes are increasingly being investigated in mixed solvents less polar than water. Such solvents may be quite relevant because their polarity more nearly resembles the natural cellular microenvironment than does pure water. The single most important criterion in selecting a non-aqueous solvent is its compatibility with the maintenance of the enzyme's catalytic activity, which must be determined experimentally for each enzyme. Non-aqueous solvents have a variety of effects on enzymes: they may bind specifically, compete with substrate binding, dissociate multimers, shift an equilibrium between two enzyme conformations, alter the amount of helix, react with the enzyme, stabilize or destabilize the enzyme, and affect the rate of the catalytic reaction in several different ways. Typically, modest concentrations of hydroxylic solvents have little effect on rates, and may even enhance the rate significantly. Higher concentrations give lower rates, in a solvent-specific and enzyme-specific manner. Hydroxylic solvents may replace water as acceptor of a phosphoryl, glycosyl, or acyl group produced by a hydrolytic enzyme. Non-aqueous solvents also make it possible to run hydrolytic reactions in the reverse direction, forming a condensation product and water as a by-product. Non-aqueous solvents are being extensively used in cryoenzymology as antifreeze agents, in solubilizing and purifying enzymes, and to a lesser degree in two-phase systems in which the non-polar substrate is dissolved in the non-aqueous phase. Liquefied aqueous phenol is an extraordinary solvent for enzymes and other proteins. It is a powerful denaturant which rapidly and irreversibly extracts the enzyme into the phenol-rich phase of a phenol-water system. This property makes phenol useful for removing protein contaminants, and for detecting labile enzyme-substrate intermediates, by extracting substrate covalently bound to enzyme into the phenol-rich phase away from all other substrate, which generally remains in the aqueous phase.     

Investigation of the binding network of IGF-I on the cavity surface of IGFBP4

Insulin-like growth factor-binding proteins (IGFBPs) control bioactivity and distribution of insulin-like growth factors (IGFs) through high-affinity complex of IGFBP and IGF. To get more insight into the binding interaction of IGF system, the site-directed mutagenesis and force-driving desorption methods were employed to study the interaction mechanism of IGFBP4 and IGF-I by molecular dynamics (MD) simulation. In IGF-I, residues Gly7 to Asp12 were found to be the hot spots and they mainly anchored on the N-domain of IGFBP4. The contact area, the shape and size of protein, the surroundings of the binding site, the hydrophobic and electrostatic interaction between the two proteins worked as a complex network to regulate the protein-protein interaction. It was also found that the unfolding of the helix was not inevitable in the mutant, and it could be regulated by careful selection of the substituted amino acid.

Modeling hydration mechanisms of enzymes in nonpolar and polar organic solvents

A comprehensive study of the hydration mechanism of an enzyme in nonaqueous media was done using molecular dynamics simulations in five organic solvents with different polarities, namely, hexane, 3-pentanone, diisopropyl ether, ethanol, and acetonitrile. In these solvents, the serine protease cutinase from Fusarium solani pisi was increasingly hydrated with 12 different hydration levels ranging from 5% to 100% (w/w) (weight of water/weight of protein). The ability of organic solvents to 'strip off' water from the enzyme surface was clearly dependent on the nature of the organic solvent. The rmsd of the enzyme from the crystal structure was shown to be lower at specific hydration levels, depending on the organic solvent used. It was also shown that organic solvents determine the structure and dynamics of water at the enzyme surface. Nonpolar solvents enhance the formation of large clusters of water that are tightly bound to the enzyme, whereas water in polar organic solvents is fragmented in small clusters loosely bound to the enzyme surface. Ions seem to play an important role in the stabilization of exposed charged residues, mainly at low hydration levels. A common feature is found for the preferential localization of water molecules at particular regions of the enzyme surface in all organic solvents: water seems to be localized at equivalent regions of the enzyme surface independently of the organic solvent employed.     

Cytotoxic effect and molecular docking of 4-ethoxycarbonylmethyl-1-(piperidin-4-ylcarbonyl)-thiosemicarbazide—a novel topoisomerase II inhibitor

The preliminary cytotoxic effect of 4-ethoxycarbonylmethyl-1-(piperidin-4-ylcarbonyl)-thiosemicarbazide hydrochloride (1)—a potent topoisomerase II inhibitor—was measured using a MTT assay. It was found that the compound decreased the number of viable cells in both estrogen receptor-positive MCF-7 and estrogen receptor-negative MDA-MB-231breast cancer cells, with IC₅₀ values of 146?±?2 and 132?±?2 μM, respectively. To clarify the molecular basis of the inhibitory action of 1, molecular docking studies were carried out. The results suggest that 1 targets the ATP binding pocket.

TDDFT prediction of UV–vis absorption and emission spectra of tocopherols in different media

We use the TDDFT/PBE0/6-31+G* method to determine the electronic absorption and emission energies, in different media, of the four forms of tocopherol, which differ by the number and the position of methyl groups on the chromanol. Geometries of the ground state S₀ and the first singlet excited state S₁ were optimized in the gas phase, and various solvents. The solvent effect is evaluated using an implicit solvation model (IEF-PCM). Our results are compared to the experimental ones obtained for the vitamin E content in several vegetable oils. For all forms of tocopherols, the HOMO–LUMO first vertical excitation is a π–π* transition. Gas phase and non-polar solvents (benzene and toluene) give higher absorption wavelengths than polar solvents (acetone, ethanol, methanol, DMSO, and water); this can be interpreted by a coplanarity between the O-H group and the chroman, allowing a better electronic resonance of the oxygen lone pairs and the aromatic ring, and therefore giving an important absorption wavelength, whereas the polar solvents give high emission wavelengths comparatively to gas phase and non-polar solvents. Fluorescence spectra permit the determination, the separation, and the identification of the four forms of tocopherols by a large difference in emission wavelength values.

Open image in new window

Graphical Abstract Scheme from process methodological to obtain the absorption and emission spectra for tocopherols

     

Kinetic behavior of phosphotriesterase and its non-covalent complexes with polyelectrolyte in systems with polar and non-polar organic solvents

Soluble phosphotriesterase from E. coli DH5 together with E. coli DH5 cells with the phosphotriesterase activity were co-immobilized into poly(vinyl alcohol) (PVA) cryogel and studied in water/organic systems with polar and non-polar organic solvents. The phosphotriesterase activity was competitively inhibited by polar organic solvents. The inhibition constant correlated with the dielectric constant () of the solvent. The rate of the enzyme-catalyzed reaction in biphasic non-polar solvent/water systems was independent of water/organic ratio and the hydrophobicity of the solvent. Formation of the non-covalent complexes with polyelectrolytes was suggested to enhance the resistance of the phosphotriesterase towards inactivation by organic solvents in their homogeneous mixtures with water.     

Lipid Extraction from <Emphasis Type="Italic">Tetraselmis</Emphasis> sp. Microalgae for Biodiesel Production Using Hexane-based Solvent Mixtures

Lipid extraction is a critical step in the downstream processing of biodiesel production from microalgae. Solvent extraction using mixtures of non-polar and polar solvents is one of the most well-known processes for this purpose. Hexane is the most common solvent of choice for large-scale lipid extractions due to its technical and economic advantages, especially its high selectivity toward lipids and low cost. In this study, extractions using mixtures of hexane and polar solvents were evaluated for their performance in order to develop a more efficient method for large-scale lipid extraction from microalgae. The combination of hexane and methanol resulted in the highest fatty acid methyl ester (FAME) yield for lipids from Tetraselmis sp. The effects of extraction conditions, including proportions of methanol to hexane, ratios of total solvent volume to dry biomass, and extraction time, on extraction yields were evaluated to determine optimum conditions providing higher lipid and FAME yields. The optimal conditions were as follows: proportion of hexane to methanol of 1:1, ratio of total solvent volume to dry biomass of 10 mL/g, and extraction time of 120 min. Finally, the selected solvent mixture and optimal conditions were applied to larger scale extraction experiments with scale-up factors of 10, 50, and 100. FAME yields of large-scale extractions were almost completely consistent with increasing scale-up factors. The results of this study suggest that a hexane and methanol mixture is a promising solvent for large-scale lipid extraction from microalgae.     

Metal ions in sugar binding, sugar specificity and structural stability of Spatholobus parviflorus seed lectin

Spatholobus parviflorus seed lectin (SPL) is a heterotetrameric lectin, with two α and two β monomers. In the crystal structure of SPL α monomer, two residues at positions 240 and 241 are missing. This region was modeled based on the positional and sequence similarities. The role of metal ions in SPL structure was analyzed by 10 ns molecular dynamics simulation. MD simulations were performed in the presence and absence of metal ions to explain the loss of haemagglutinating property of the lectin due to demetallization. Demetallized structure was found to deviate drastically at the metal binding loop region. Affinity of different sugars like N-acetyl galactosamine (GalNAc), D-galactose and lactose towards the native and demetallized protein was calculated by molecular docking studies. It was found that the sugar binding site got severely distorted in demetallized lectin. Consequently, sugar binding ability of lectin might be decreasing in the demetallized condition. Isothermal titration calorimetric (ITC) analysis of the sugars in the presence of native and demetallized protein confirmed the in silico results. It was observed after molecular dynamics simulations, that significant structural deviations were not caused in the quaternary structure of demetallized lectin. It was confirmed that the structural changes modified the sugar binding ability, as well as sugar specificity of the present lectin. The role of metal ions in sugar binding is described based on the in silico studies and ITC analysis. A comprehensive analysis of the ITC data suggests that the sugar specificity of the metal bound lectin and the loss of sugar specificity due to metal chelation are not linear.

Retrospective molecular docking study of WY-25105 ligand to β-secretase and bias of the three-dimensional structure flexibility

β-Secretase (BACE) is a very promising target in the search for a treatment for Alzheimer’s disease using a protein–ligand inhibition approach. Given the many published X-ray structures of BACE protein, structure-based drug design has been used extensively to support new inhibitor discovery programs. Due to the high flexibility and large catalytic site of this protein, sampling of the huge conformational space of the binding site is the big challenge to overcome and is the main limitation of the most widely used docking programs. Incorrect treatment of these pitfalls can introduce bias into ligand docking and could affect the results. This is especially the case with the WY-25105 compound reported by the Wyeth Corporation as a BACE ligand that did not fit into any of the known crystal structures. In the present retrospective study, a set of available X-ray enzyme structures was selected and molecular dynamics simulations were conducted to generate more diverse representative BACE protein conformations. These conformations were then used for a docking study of the WY-25105 compound. The results confirmed the need to use an ensemble of structures in protein–ligand docking for identification of new binding modes in structure-based drug design of BACE inhibitors.

Conformational flexibility of the leucine binding protein examined by protein domain coarse-grained molecular dynamics

Periplasmic binding proteins are the initial receptors for the transport of various substrates over the inner membrane of gram-negative bacteria. The binding proteins are composed of two domains, and the substrate is entrapped between these domains. For several of the binding proteins it has been established that a closed-up conformation exists even without substrate present, suggesting a highly flexible apo-structure which would compete with the ligand-bound protein for the transporter interaction. For the leucine binding protein (LBP), structures of both open and closed conformations are known, but no closed-up structure without substrate has been reported. Here we present molecular dynamics simulations exploring the conformational flexibility of LBP. Coarse grained models based on the MARTINI force field are used to access the microsecond timescale. We show that a standard MARTINI model cannot maintain the structural stability of the protein whereas the ELNEDIN extension to MARTINI enables simulations showing a stable protein structure and nanosecond dynamics comparable to atomistic simulations, but does not allow the simulation of conformational flexibility. A modification to the MARTINI-ELNEDIN setup, referred to as domELNEDIN, is therefore presented. The domELNEDIN setup allows the protein domains to move independently and thus allows for the simulation of conformational changes. Microsecond domELNEDIN simulations starting from either the open or the closed conformations consistently show that also for LBP, the apo-structure is flexible and can exist in a closed form.

Effects of protein binding on a lipid bilayer containing local anesthetic articaine, and the potential of mean force calculation: a molecular dynamics simulation approach

Articaine, as a local anesthetic drug has been simulated in neutral and charged forms, and its interaction with the dimyristoylphosphatidylcholine (DMPC) lipid bilayer membrane is investigated by molecular dynamics simulation using GROMACS software. In order to obtain the optimum location of the drug molecules, as they penetrate into the membrane, umbrella sampling is applied and the free energy is calculated. The effect of protein binding to DMPC membrane on the process of drug diffusion through the membrane is considered. Five simulation systems are designed and by applying the potential of mean force, the molecular dynamics simulation on the system is performed. In light of the obtained results, the electrostatic potential, variation of lipid bilayer’s order parameter and the diffusion coefficient of drug are discussed.

Comparison of the structural characteristics of Cu2+-bound and unbound α-syn12 peptide obtained in simulations using different force fields

The effects of Cu²⁺ binding and the utilization of different force fields when modeling the structural characteristics of α-syn12 peptide were investigated. To this end, we performed extensive temperature replica exchange molecular dynamics (T-REMD) simulations on Cu²⁺-bound and unbound α-syn12 peptide using the GROMOS 43A1, OPLS-AA, and AMBER03 force fields. Each replica was run for 300 ns. The structural characteristics of α-syn12 peptide were studied based on backbone dihedral angle distributions, free-energy surfaces obtained with different reaction coordinates, favored conformations, the formation of different Turn structures, and the solvent exposure of the hydrophobic residues. The findings show that AMBER03 prefers to sample helical structures for the unbound α-syn12 peptide and does not sample any β-hairpin structure for the Cu²⁺-bound α-syn12 peptide. In contrast, the central structure of the major conformational clusters for the Cu²⁺-bound and unbound α-syn12 peptide according to simulations performed using the GROMOS 43A1 and OPLS-AA force fields is a β-hairpin with Turn_9-6. Cu²⁺ can also promote the formation of the β-hairpin and increase the solvent exposure of hydrophobic residues, which promotes the aggregation of α-syn12 peptide. This study can help us to understand the mechanisms through which Cu²⁺ participates in the fibrillation of α-syn12 peptide at the atomic level, which in turn represents a step towards elucidating the nosogenesis of Parkinson’s disease.