Structure and dynamics of Trichoderma harzianum Cel7B suggest molecular architecture adaptations required for a wide spectrum of activities on plant cell wall polysaccharides.

Cellulases from glycoside hydrolase family 7 (GH7) play crucial roles in plant lignocellulose deconstruction by fungi, but structural information available for GH7 fungal endoglucanases is limited when compared to the number of known sequences in the family. Here, we report the X-ray structure of the glycosylated catalytic domain (CD) of Trichoderma harzianum endoglucanase, ThCel7B, solved and refined at 2.9 Šresolution. Additionally, our extensive molecular dynamics simulations of this enzyme in complex with a variety of oligosaccharides provide a better understanding of its promiscuous hydrolytic activities on plant cell wall polysaccharides. The simulations demonstrate the importance of the hydrogen bond between substrate O2 hydroxyl in the subsite -1 and a side chain of catalytic Glu196 which renders ThCel7B capable to catalytically cleave cello and xylooligosaccharides, but not mannooligosaccharides. Moreover, detailed structural analyses and MD simulations revealed an additional binding pocket, suitable for accommodation of oligosaccharide decorations and/or substrates with mixed glycoside bonds that abuts onto the binding cleft close to subsite +2.

Molecular dynamics simulations of the dielectric properties of fructose aqueous solutions.

The static dielectric permittivity and dielectric relaxation properties of fructose aqueous solutions of different concentrations ranging from 1.0 to 4.0 mol l(-1) are investigated by means of molecular dynamics simulations. The contributions from intra- and interspecies molecular correlations were computed individually for both the static and frequency-dependent dielectric properties, and the results were compared with the available experimental data. Simulation results in the time- and frequency-domains were analyzed and indicate that the presence of fructose has little effect on the position of the fast, high-frequency (>500 cm(-1)) components of the dielectric response spectrum. The low-frequency (<0.1 cm(-1)) components, however, are markedly influenced by sugar concentration. Our analysis indicates that fructose-fructose and fructose-water interactions strongly affect the rotational-diffusion regime of molecular motions in the solutions. Increasing fructose concentration not only enhances sugar-sugar and sugar-water low frequency contributions to the dielectric loss spectrum but also slows down the reorientational dynamics of water molecules. These results are consistent with previous computer simulations carried out for other disaccharide aqueous solutions.

Anomalous dynamics of hydration water in carbohydrate solutions.

Molecular dynamics simulations have been carried out to investigate the dynamics of fructose aqueous solutions up to 70 wt % concentration. We find that the hydrogen (H)-bonded network of fructose molecules extends with increasing sugar content and forms a structurally heterogeneous system around and above 45 wt % concentration, characterized as a percolated-like solute domain permeated by patchy regions of solvent. The presence of such aggregates in concentrated solutions promotes the slowing down of water translational, reorientational, and H-bonding dynamics, typical of many biomolecular environments. Analysis of the effects of the topological and energetic disorder of the sugar aggregates on vicinal water dynamics, similar to that recently carried out for the hydration layer of proteins by Pizzitutti et al. (J. Phys. Chem. B 2007, 111, 7584), reveals many similarities between the dynamical anomaly of the hydration layers of both systems. Like a protein surface, topological and energetic disorders of the sugar aggregates both contribute to the translational diffusion anomaly. However, unlike in the vicinity of a protein surface, the rotational relaxation is also hindered by the topological disorder created by the intertwined, percolating sugar clusters in concentrated solutions.

Ligand dissociation from estrogen receptor is mediated by receptor dimerization: evidence from molecular dynamics simulations.

Estrogen Receptor (ER) is an important target for pharmaceutical design. Like other ligand-dependent transcription factors, hormone binding regulates ER transcriptional activity. Nevertheless, the mechanisms by which ligands enter and leave ERs and other nuclear receptors remain poorly understood. Here, we report results of locally enhanced sampling molecular dynamics simulations to identify dissociation pathways of two ER ligands [the natural hormone 17beta-estradiol (E(2)) and the selective ER modulator raloxifene (RAL)] from the human ERalpha ligand-binding domain in monomeric and dimeric forms. E(2) dissociation occurs via three different pathways in ER monomers. One resembles the mousetrap mechanism (Path I), involving repositioning of helix 12 (H12), others involve the separation of H8 and H11 (Path II), and a variant of this pathway at the bottom of the ligand-binding domain (Path II'). RAL leaves the receptor through Path I and a Path I variant in which the ligand leaves the receptor through the loop region between H11 and H12 (Path I'). Remarkably, ER dimerization strongly suppresses Paths II and II' for E(2) dissociation and modifies RAL escape routes. We propose that differences in ligand release pathways detected in the simulations for ER monomers and dimers provide an explanation for previously observed effects of ER quaternary state on ligand dissociation rates and suggest that dimerization may play an important, and hitherto unexpected, role in regulation of ligand dissociation rates throughout the nuclear receptor family.

Carbohydrate clustering in aqueous solutions and the dynamics of confined water.

We use molecular dynamics simulations to investigate structure and dynamics of fructose aqueous solutions in the 1-5 M concentration range at ambient conditions. We analyze hydration structures, H-bond statistics, and size distribution of H-bonded carbohydrate clusters as functions of concentration. We find that the local tetrahedral order of water is reasonably well-preserved and that the solute tends to appear as scattered "isolated" molecules at low concentrations and as H-bonded clusters for less diluted solutions. The sugar cluster size distribution exhibits a sharp transition to a percolated cluster between 3.5 and 3.8 M. The percolated cluster forms an intertwined network of H-bonded saccharides that imprisons water. For the dynamics, we find good agreement between simulation and available experimental results for the self-diffusion coefficients. Water librational dynamics is little affected by sugar concentration, whereas reorientational relaxation is described by a concentration-independent bulk-like component attributed to noninterfacial water molecules and a slower component (strongly concentration dependent) that arises from interfacial solvent molecules and, hence, depends on the dynamics of the cluster structure itself. Analysis of H-bonding survival probability functions indicates that the formation of carbohydrate clusters upon increasing concentration enhances the H-bond relaxation time and slows down the entire system dynamics. We find that multiexponential or stretched-exponential fits alone cannot describe the H-bond survival probabilities for the entire postlibrational time span of our data (0.1-100 ps), as opposed to a combined stretched-plus-biexponential function, which provides excellent fits. Our results suggest that water dynamics in concentrated fructose solutions resembles in many ways that of protein hydration water.

Development of new Cd2+ and Pb2+ Lennard-Jones parameters for liquid simulations.

We present new Lennard-Jones parameters for Cd2+ and Pb2+ ion-water interactions and describe a general methodology to obtain these parameters for any ion. Our strategy is based on the adjustment of ion parameters to reproduce simultaneously experimental absolute hydration free energy and structural properties, namely, g(r) and coordination numbers, obtained from X-ray liquid scattering and quantum mechanical/molecular mechanical (QM/MM) calculations. The validation of the obtained parameters is made by the calculation of dynamical properties and comparing them with experimental values and theoretical results from the literature. The transferability of parameters is checked by the calculation of thermodynamic, structural, and dynamical properties cited above with four different water models. The results obtained for Cd2+ and Pb2+ show an overall agreement with reference values. The absolute hydration free energy calculated with the TIP3P, SPC/E, SPC, and TIP4P water models presents, respectively, percent differences of 3.8, 3.0, 4.3, and 7.2% for lead(II) and 9.8, 8.4, 10.2, and 14.1% for cadmium(II) when compared with experimental values. Ion-water mean distance and coordination numbers for the first coordination shell are in good agreement with experimental and QM/MM results for both ions. Cd2+ shows a lesser diffusion coefficient compared to that of Pb2+ despite its smaller atomic radius, indicating a more persistent first coordination shell for the cadmium(II) ion, a result confirmed with calculations of the mean residence time of water molecules in the first coordination shell.

Molecular dynamics simulations reveal multiple pathways of ligand dissociation from thyroid hormone receptors.

Nuclear receptor (NR) ligands occupy a pocket that lies within the core of the NR ligand-binding domain (LBD), and most NR LBDs lack obvious entry/exit routes upon the protein surface. Thus, significant NR conformational rearrangements must accompany ligand binding and release. The precise nature of these processes, however, remains poorly understood. Here, we utilize locally enhanced sampling (LES) molecular dynamics computer simulations to predict molecular motions of x-ray structures of thyroid hormone receptor (TR) LBDs and determine events that permit ligand escape. We find that the natural ligand 3,5,3'-triiodo-L-thyronine (T(3)) dissociates from the TRalpha1 LBD along three competing pathways generated through i), opening of helix (H) 12; ii), separation of H8 and H11 and the Omega-loop between H2 and H3; and iii), opening of H2 and H3, and the intervening beta-strand. Similar pathways are involved in dissociation of T(3) and the TRbeta-selective ligand GC24 from TRbeta; the TR agonist IH5 from the alpha- and beta-TR forms; and Triac from two natural human TRbeta mutants, A317T and A234T, but are detected with different frequencies in simulations performed with the different structures. Path I was previously suggested to represent a major pathway for NR ligand dissociation. We propose here that Paths II and III are also likely ligand escape routes for TRs and other NRs. We also propose that different escape paths are preferred in different situations, implying that it will be possible to design NR ligands that only associate stably with their cognate receptors in specific cellular contexts.

Optical Kerr effect in supercooled water.

We present molecular dynamics simulations of the optical Kerr effect in liquid and supercooled water and compare with recent time-resolved Kerr spectroscopy measurements [R. Torre, Nature (London) 428, 296 (2004)]. The short time features of the Kerr response, characterized by peaks near 15, 60, and 160 fs, are weakly temperature dependent. The long-time decay is well described by a stretched exponential with a nearly constant stretch parameter and relaxation times that follow a power law approximately (T-T(S))(-gamma), with T(S)=198.3 K and gamma=2.35. Our findings are discussed in the light of the spectroscopy data and previous simulation analyzes of the structural relaxation in supercooled water.

Solvation of coumarin 314 at water/air interfaces containing anionic surfactants. I. Low coverage.

Through the use of molecular dynamics techniques, we analyze equilibrium and dynamical aspects of the solvation of Coumarin 314 adsorbed at water/air interfaces in the presence of sodium dodecyl sulfate surfactant molecules. Three different coverages in the submonolayer regime were considered, 500, 250, and 100 A(2)/SDS molecule. The surfactant promotes two well-differentiated solvation environments, which can be clearly distinguished in terms of their structures for the largest surfactant coverage considered. The first one is characterized by the probe lying adjacent or exterior to two-dimensional spatial domains formed by clustered surfactant molecules. A second type of solvation environment is found in which the coumarin appears embedded within compact surfactant domains. Equilibrium and dynamical aspects of the interfacial orientation of the probe are investigated. Our results show a gradual transition from parallel to perpendicular dipolar alignment of the probe with respect to the interface as the concentration of surfactant rho(s) increases. The presence of the surfactant leads to an increase in the roughness and in the characteristic width of the water/air interface. These modifications are also manifested by the decorrelation times for the probe reorientational dynamics, which become progressively slower with rho(s) in both solvation states, although much more pronounced for the embedded ones. The dynamical characteristics of the solvation responses of the charged interfaces are also analyzed, and the implications of our findings to the interpretation of available experimental measurements are discussed.

A simulation study of the optical Kerr effect in liquid water.

A molecular dynamics simulation study is presented for the dynamics of the polarizability anisotropy of liquid water using the SPC/E model and a dipolar induction scheme that involves the intrinsic polarizability and first hyperpolarizability tensors obtained from ab initio quantum chemical calculations at the MP2/6-311++G(d,p) level. The time-correlation functions for the collective polarizability anisotropy, the optical Kerr effect response, and the frequency spectra are analyzed in terms of the intrinsic and induced polarizability contributions. At short times, the simulated Kerr nuclear response exhibits maxima near 15, 50 and 180 fs, followed by a diffusive tail which has been fitted by a bi-exponential with time constants ca. 0.4 and 2.5 ps. The short time features are in good agreement with available simulation and experimental results. The agreement with experiments is less satisfactory for the diffusive components. The main features of the frequency spectrum include a rotational-diffusion peak centered around 3 cm(-1), a collision-induced (hindered translations) band near 200 cm(-1), and a broad librational band at 450 cm(-1). The simulation results are in good agreement with experimental frequency spectra obtained from Kerr effect and related spectroscopies, but fail to reproduce the experimental band near 60 cm(-1).