How do trehalose, maltose, and sucrose influence some structural and dynamical properties of lysozyme? Insight from molecular dynamics simulations

The influence of three well-known disaccharides, namely, trehalose, maltose, and sucrose, on some structural and dynamical properties of lysozyme has been investigated by means of molecular dynamics computer simulations in the 37-60 wt % concentration range. The effects of sugars on the protein conformation are found to be relatively weak, in agreement with the preferential hydration of lysozyme. Conversely, sugars seem to increase significantly the relaxation times of the protein. These effects are shown to be correlated to the fractional solvent accessibilities of lysozyme residues and further support the slaving of protein dynamics. Moreover, a significant increase in the relaxation times of lysozyme, sugars, and water molecules is observed within the studied concentration range and may result from the percolation of the hydrogen-bond network of sugar molecules. This percolation appears to be of primary importance to explain the influence of sugars on the dynamical properties of lysozyme and water.     

How sensitive are nanosecond molecular dynamics simulations of proteins to changes in the force field?

The sensitivity of molecular dynamics simulations to variations in the force field has been examined in relation to a set of 36 structures corresponding to 31 proteins simulated by using different versions of the GROMOS force field. The three parameter sets used (43a1, 53a5, and 53a6) differ significantly in regard to the nonbonded parameters for polar functional groups and their ability to reproduce the correct solvation and partitioning behavior of small molecular analogues of the amino acid side chains. Despite the differences in the force field parameters no major differences could be detected in a wide range of structural properties such as the root-mean-square deviation from the experimental structure, radii of gyration, solvent accessible surface, secondary structure, or hydrogen bond propensities on a 5 to 10 ns time scale. The small differences that were observed correlated primarily with the presence of charged residues as opposed to residues that differed most between the parameter sets. The work highlights the variation that can be observed in nanosecond simulations of protein systems and implications of this for force field validation, as well as for the analysis of protein simulations in general.     

How does the urea dynamics differ from water dynamics inside the reverse micelle?

In this study, the urea dynamics inside AOT reverse micelle (RM) has been monitored without intervention of water using time-resolved fluorescence techniques from the picosecond to nanosecond time regime. It has been observed that urea dynamics inside the reverse micelle is severely retarded compared to water RM due to the formation of highly networked urea cluster inside the RM. Time-resolved fluorescence anisotropy study also confirms the existence of a confined environment around the dye at higher concentrations of urea inside the reverse micelle. The dynamics of urea-water mixtures inside AOT reverse micelle has also been monitored with increasing urea concentration to get insight about the effect of urea on the overall solvation dynamics feature. It has been observed that with the increase in urea concentration, the overall dynamics becomes slower, and it infers the presence of few water or urea molecules, those strongly associated with surrounding urea and (or) water by hydrogen bonds.     

How does huperzine A enter and leave the binding gorge of acetylcholinesterase? Steered molecular dynamics simulations

The entering and leaving processes of Huperzine A (HupA) binding with the long active-site gorge of Torpedo californica acetylcholinesterase (TcAChE) have been investigated by using steered molecular dynamics simulations. The analysis of the force required along the pathway shows that it is easier for HupA to bind to the active site of AChE than to disassociate from it, which for the first time interprets at the atomic level the previous experimental result that unbinding process of HupA is much slower than its binding process to AChE. The direct hydrogen bonds, water bridges, and hydrophobic interactions were analyzed during two steered molecular dynamics (SMD) simulations. Break of the direct hydrogen bond needs a great pulling force. The steric hindrance of bottleneck might be the most important factor to produce the maximal rupture force for HupA to leave the binding site but it has a little effect on the binding process of HupA with AChE. Residue Asp72 forms a lot of water bridges with HupA leaving and entering the AChE binding gorge, acting as a clamp to take out HupA from or put HupA into the active site. The flip of the peptide bond between Gly117 and Gly118 has been detected during both the conventional MD and SMD simulations. The simulation results indicate that this flip phenomenon could be an intrinsic property of AChE and the Gly117-Gly118 peptide bond in both HupA bound and unbound AChE structures tends to adopt the native enzyme structure. At last, in a vacuum the rupture force is increased up to 1500 pN while in water solution the greatest rupture force is about 800 pN, which means water molecules in the binding gorge act as lubricant to facilitate HupA entering or leaving the binding gorge.     

Car-Parinello molecular dynamics simulations of thionitroxide and S-nitrosothiol in the gas phase, methanol, and water—A theoretical study of S-nitrosylation

A dilemma about whether thionitroxide radical (RSNHO) or S-nitrosothiol (RSNO) is observed in protein S-nitrosylation has arisen recently. To illustrate the effect of chemical environment on these structures, this paper presents quantum mechanical molecular dynamics of thionitroxide, and cis-and trans-S-nitrosothiols in the gas phase, methanol, and water. By using Car-Parrinello molecular dynamics (CPMD), we have observed that there is free rotation about the S-N bond at 300 K in thionitroxide, but no such rotation is observed for S-nitrosothiol. The C-S-N-O torsion angle distribution in thionitroxide is s-ignificantly dependent upon the surrounding environment, leading to either gauche-, cis-, or trans-conformation. In the case of S-nitrosothiol the C-S-N-O plane is twisted slightly by 5°-15° in the cis-isomer, while the periplanar structure is well-retained in the trans-isomer. The calculated results are in agreement with the X-ray crystallographic data of small molecular RSNO species. Interestingly, for both compounds, the CPMD simulations show that solvation can cause a decrease in the S-N bond length. Moreover, the oxygen atom of thionitroxide is found to be a good hydrogen-bond acceptor, forming an oxyanion-hole-like hydrogen bonding network.     

What is the form of muscimol from fly agaric mushroom (Amanita muscaria) in water? An insight from NMR experiment supported by molecular modeling

The biologically active alkaloid muscimol is present in fly agaric mushroom (Amanita muscaria), and its structure and action is related to human neurotransmitter _γ-aminobutyric acid (GABA). The current study reports on determination of muscimol form present in water solution using multinuclear ¹H and ¹³C nuclear magnetic resonance (NMR) experiments supported by density functional theory molecular modeling. The structures of three forms of free muscimol molecule both in the gas phase and in the presence of water solvent, modeled by polarized continuous model, and nuclear magnetic isotropic shieldings, the corresponding chemical shifts, and indirect spin–spin coupling constants were calculated. Several J-couplings observed in proton and carbon NMR spectra, not available before, are reported. The obtained experimental spectra, supported by theoretical calculations, favor the zwitterion form of muscimol in water. This structure differs from NH isomer, previously determined in dimethyl sulfoxide (DMSO) solution. In addition, positions of signals C3 and C5 are reversed in both solvents.     

Can the dodecahedral water cluster naturally form in methane aqueous solutions? A molecular dynamics study on the hydrate nucleation mechanisms

By performing a large scale of molecular dynamics simulations, we analyze 60 x 10(6) hydration shells of methane to examine whether the dodecahedral water cluster (DWC) can naturally form in methane aqueous solutions--a fundamental question relevant to the nucleation mechanisms of methane hydrate. The analyzing method is based on identifying the incomplete cages (ICs) from the hydration shells and quantifying their cagelike degrees (zetaC=0-1). Here, the zetaC is calculated according to the H-bond topological network of IC and reflects how the IC resembles the complete polyhedral cage. In this study, we obtain the zetaC distributions of ICs in methane solutions and find the occurrence probabilities of ICs reduce with zetaC very rapidly. The ICs with zetaC>or=0.65 are studied, which can be regarded as the acceptable cagelike structures in appearance. Both increasing the methane concentration and lowering the temperature can increase their occurrence probabilities through slowing down the water molecules. Their shapes, cage-maker numbers, and average radii are also discussed. About 13-14 of these ICs are face saturated, meaning that every edges are shared by two faces. The face-saturated ICs have the potential to act as precursors of hydrate nucleus because they can prevent the encaged methane from directly contacting other dissolved methane when an event of methane aggregation occurs. The complete cages, i.e., the ICs with zetaC=1, form only in the solutions with high methane concentration, and their occurrence probabilities are about 10(-6). Most of their shapes are different from the known hydrate cages, but we indeed observe a standard 5(12)6(2) hydrate cage. We do not find the expected DWC, and its occurrence probability is estimated to be far less than 10(-7). Additionally, the IC analysis proposed in this work is also very useful in other studies not only on the formation, dissociation, and structural transition of hydrates but also on the hydrophobic hydration of apolar solutes.     

Structural insight into the role of Gln293Met mutation on the Peloruside A/Laulimalide association with αβ-tubulin from molecular dynamics simulations,binding free energy calculations and weak interactions analysis

Peloruside A (PLA) and Laulimalide (LAU) are novel microtubule-stabilizing agents with promising properties against different cancer types. These ligands share a non-taxoid binding site at the outer surface of β-tubulin and promote microtubule stabilization by bridging two adjacent αβ-tubulin dimers from parallel protofilaments. Recent site-directed mutagenesis experiments confirmed the existence of a unique β-tubulin site mutation (Gln293Met) that specifically increased the activity of PLA and caused resistance to LAU, without affecting the stability of microtubules in the absence of the ligands. In this work, fully atomistic molecular dynamics simulations were carried out to examine the PLA and LAU association with native and mutated αβ-tubulin in the search for structural and energetic evidence to explain the role of Gln293Met mutation on determining the activity of these ligands. Our results revealed that Gln293Met mutation induced the loss of relevant LAU–tubulin contacts but exerted negligible changes in the interaction networks responsible for PLA–tubulin association. Binding free energy calculations (MM/GBSA and MM/PBSA), and weak interaction analysis (aNCI) predicted an increased affinity for PLA, and a weakened association for LAU after mutation, thus suggesting that Gln293Met mutation exerts its action by a modulation of drug–tubulin interactions. These results are valuable to increase understanding about PLA and LAU activity and to assist the future design of novel agents targeting the PLA/LAU binding pocket.     

Flexibility of a biotinylated ligand in artificial metalloenzymes based on streptavidin—an insight from molecular dynamics simulations with classical and ab initio force fields

In the field of enzymatic catalysis, creating activity from a non catalytic scaffold is a daunting task. Introduction of a catalytically active moiety within a protein scaffold offers an attractive means for the creation of artificial metalloenzymes. With this goal in mind, introduction of a biotinylated d⁶-piano-stool complex within streptavidin (SAV) affords enantioselective artificial transfer-hydrogenases for the reduction of prochiral ketones. Based on an X-ray crystal structure of a highly selective hybrid catalyst, displaying significant disorder around the biotinylated catalyst [η⁶-(p-cymene)Ru(Biot-p-L)Cl], we report on molecular dynamics simulations to shed light on the protein–cofactor interactions and contacts. The results of these simulations with classical force field indicate that the SAV-biotin and SAV-catalyst complexes are more stable than ligand-free SAV. The point mutations introduced did not affect significantly the overall behavior of SAV and, unexpectedly, the P64G substitution did not provide additional flexibility to the protein scaffold. The metal-cofactor proved to be conformationally flexible, and the S112K or P64G mutants proved to enhance this effect in the most pronounced way. The network of intermolecular hydrogen bonds is efficient at stabilizing the position of biotin, but much less at fixing the conformation of an extended biotinylated ligand. This leads to a relative conformational freedom of the metal-cofactor, and a poorly localized catalytic metal moiety. MD calculations with ab initio potential function suggest that the hydrogen bonds alone are not sufficient factors for full stabilization of the biotin. The hydrophobic biotin-binding pocket (and generally protein scaffold) maintains the hydrogen bonds between biotin and protein.     

What can you learn from a molecular probe? New insights on the behavior of C343 in homogeneous solutions and AOT reverse micelles

The behavior of C343, a common molecular probe utilized in solvation dynamics experiments, was studied in homogeneous media and in aqueous and nonaqueous reverse micelles (RMs). In homogeneous media, the Kamlet and Taft solvatochromic comparison method quantified solute-solvent interactions from the absorption and emission bands showing that the solvatochromic behavior of the dye depends not only on the polarity of the medium but also on the hydrogen-bonding properties of the solvent. Specifically, in the ground state the molecule displays a bathochromic shift with the polarity polarizability (pi) and the H-bond acceptor (beta) ability of the solvents and a hypsochromic shift with the hydrogen donor ability (alpha) of the media. The carboxylic acid group causes C343 to display greater sensitivity to the beta than to the pi polarity parameter; this sensitivity increases in the excited state, while the dependence on alpha vanishes. This demonstrates that C343 forms a stable H-bond complex with solvents with high H-bond acceptor ability (high beta) and low H-bond donor character (low alpha). Spectroscopy in nonpolar solvents reveals J-aggregate formation. With information from the Kamlet-Taft analysis, C343 was used to explore RMs composed of water or polar solvents/sodium 1,4-bis-2-ethylhexylsulfosuccinate (AOT)/isooctane using absorption, emission, and time-resolved spectroscopies. Sequestered polar solvents included ethylene glycol (EG), formamide (FA), N,N-dimethylformamide (DMF), and N,N-dimethylacetamide (DMA). Dissolved in the AOT RM systems at low concentration, C343 exists as a monomer, and when introduced to the RM samples in its protonated form, C343 remains protonated driving it to reside in the interface rather than the water pool. The solvathochromic behavior of the dye depends the specific polar solvent encapsulated in the RMs, revealing different types of interactions between the solvents and the surfactant. EG and water H-bond with the AOT sulfonate group destroying their bulk H-bonded structures. While water remains well segregated from the nonpolar regions, EG appears to penetrate into the oil side of the interface. In aqueous AOT RMs, C343 interacts with neither the sulfonate group nor the water, perhaps because of intramolecular H-bonding in the dye. DMF and DMA interact primarily through dipole-dipole forces, and the strong interactions with AOT sodium counterions destroy their bulk structure. FA also interacts with the Na+ counterions but retains its H-bond network present in bulk solvent. Surprisingly, FA appears to be the only polar solvent other than water forming a "polar-solvent pool" with macroscopic properties similar to the bulk.     

Modeling the molecular dynamics of cytochrome C in aqueous and water–methanol environment

One-microsecond molecular dynamics of horse heart cytochrome C was modeled in aqueous and water–methanol environment. It was shown that the coordination bond between Met-80 sulfur and heme iron is broken in water– methanol solution.     

How does solvent molecular size affect the microscopic structure in polymer solutions?

Monte Carlo simulation has been used to investigate the effects of linear solvent molecular size on polymer chain conformation in solutions. Increasing the solvent molecular size leads to shrinkage of the polymer chains and increase of the critical overlap concentrations. The root-mean-square radius of gyration of polymer chains (R(g)) is less sensitive to the variation of polymer concentration in solutions of larger solvent molecules. In addition, the dependency of R(g) on polymer concentration under normal solvent conditions and solvent molecular size is in good agreement with scaling laws. When the solvent molecular size approaches the ideal end-to-end distance of the polymer chain, an extra aggregation of polymer chains occurs, and the solvent becomes the so-called medium-sized solvent. When the size of solvent molecules is smaller than the medium size, the polymer chains are swollen or partially swollen. However, when the size of solvent molecules is larger than the medium size, the polymer coils shrink and segregate, enwrapped by the large solvent molecules.     

Structural changes in the molecular sheets along (hk0) planes derived from cellulose Iβ by molecular dynamics simulations

We studied the stability of molecular sheets with four cellotetraoses in an aqueous environment by molecular dynamics simulation to identify the molecular details of first structure as one of the possibilities in the course of crystallization of cellulose I. After simulation, the molecular sheets formed by van der Waals forces along the (11?0) and (110) crystal plane did not change their structures in an aqueous environment, whereas the other ones formed by hydrogen bonds along the (100) and (200) crystal plane changed into a van der Waals associated molecular sheet, similar to the former. These simulated molecular sheets formed by van der Waals forces were structurally stable in water because of their hydrophilic exterior and hydrophobic interior. Therefore, if the molecular sheet structures are formed in the real system, the sheets formed by van der Waals forces are probably the initial structure of crystallization. A close analysis indicated that these sheets could be classified into two groups in terms of the hydrogen bonding networks, camber angle, and main and side chain conformations. One group was the molecular sheets corresponding to the (110) after simulation. This sheet is probably rigid because intramolecular hydrogen bonds of the chains in the sheet are highly developed. The other group was the molecular sheets corresponding to (200), (100), and (11?0) crystal plane: the chains in these sheets seemed to be rather flexible due to their moderately developed intramolecular hydrogen bonds.     

Temperature dependent dielectric relaxation spectroscopy of amyl acetate-xylene solutions – An approach to molecular and co-operative dynamics

Dielectric relaxation parameters of AMA-Xylene solutions have been investigated thoroughly to discern the associativeness as well as heterogeneous behavior among the hetero molecules in the solutions such that dielectric constant (ε_o) and relaxation time (τ) of solutions has been determined. Complex permittivity spectra (CPS) for Amyl Acetate (AMA)-Xylene solutions were measured in 10 MHz–50 GHz frequency range using Time domain reflectometry (TDR) technique. It is exciting to explore the cooperative nature and association between AMA-Xylene. With co-operative domains (CDs), dynamics in AMA-AMA and AMA-Xylene molecules have been explicated while Kirkwood validates CDs with diverse exchanges through bonding. The experimental values of ε_o are compared with theoretical values obtained from the molecular parameters suggested by the Luzar. Thermodynamic parameters for AMA-Xylene show positive values of Enthalpy and Entropy for all concentrations, indicating that the system is endothermic and less ordered. Gibbs free energy decreases with increase in volume fraction of amyl acetate (V_AMA).     

How do SET-domain protein lysine methyltransferases achieve the methylation state specificity? Revisited by Ab initio QM/MM molecular dynamics simulations

A distinct protein lysine methyltransferase (PKMT) only transfers a certain number of methyl group(s) to its target lysine residue in spite of the fact that a lysine residue can be either mono-, di-, or tri-methylated. In order to elucidate how such a remarkable product specificity is achieved, we have carried out ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations on two SET-domain PKMTs: SET7/9 and Rubisco large subunit methyltransferase (LSMT). The results indicate that the methylation state specificity is mainly controlled by the methyl-transfer reaction step, and confirm that SET7/9 is a mono-methyltransferase while LSMT has both mono-and di-methylation activities. It is found that the binding of the methylated lysine substrate in the active site of SET7/ 9 opens up the cofactor AdoMet binding channel so that solvent water molecules get access to the active site. This disrupts the catalytic machinery of SET7/9 for the di-methylation reaction, which leads to a higher activation barrier, whereas for the LSMT, its active site is more spacious than that of SET7/9, so that the methylated lysine substrate can be accommodated without interfering with its catalytic power. These detailed insights take account of protein dynamics and are consistent with available experimental results as well as recent theoretical findings regarding the catalytic power of SET7/9.     

How delocalized is N,N,N',N'-tetraphenylphenylenediamine radical cation? An experimental and theoretical study on the electronic and molecular structure

The electronic and molecular structure of N,N,N',N'-tetraphenylphenylenediamine radical cation 1(+) is in focus of this study. Resonance Raman experiments showed that at least eight vibrational modes are strongly coupled to the optical charge resonance band which is seen in the NIR. With the help of a DFT-based vibrational analysis, these eight modes were assigned to symmetric vibrations. The contribution of these symmetric modes to the total vibrational reorganization energy is dominant. These findings are in agreement with the conclusions from a simple two-state two-mode Marcus-Hush analysis which yields a tiny electron-transfer barrier. The excellent agreement of the X-ray crystal structure analysis and the DFT computed molecular structure of 1(+) on one hand as well as the solvent and solid-state IR spectra and the DFT-calculated IR active vibrations on the other hand prove 1(+) adopts a symmetrical delocalized Robin-Day class III structure both in the solid state and in solution.