Search

The interaction between xylan and cellulose microfibrils is important for secondary cell wall properties in vascular plants; however, the molecular arrangement of xylan in the cell wall and the nature of the molecular bonding between the polysaccharides are unknown. In dicots, the xylan backbone of β-(1,4)-linked xylosyl residues is decorated by occasional glucuronic acid, and approximately one-half of the xylosyl residues are O-acetylated at C-2 or C-3. We recently proposed that the even, periodic spacing of GlcA residues in the major domain of dicot xylan might allow the xylan backbone to fold as a twofold helical screw to facilitate alignment along, and stable interaction with, cellulose fibrils; however, such an interaction might be adversely impacted by random acetylation of the xylan backbone. Here, we investigated the arrangement of acetyl residues in Arabidopsis xylan using mass spectrometry and NMR. Alternate xylosyl residues along the backbone are acetylated. Using molecular dynamics simulation, we found that a twofold helical screw conformation of xylan is stable in interactions with both hydrophilic and hydrophobic cellulose faces. Tight docking of xylan on the hydrophilic faces is feasible only for xylan decorated on alternate residues and folded as a twofold helical screw. The findings suggest an explanation for the importance of acetylation for xylan-cellulose interactions, and also have implications for our understanding of cell wall molecular architecture and properties, and biological degradation by pathogens and fungi. They will also impact strategies to improve lignocellulose processing for biorefining and bioenergy., (© 2014 The Authors. The Plant Journal published by Society for Experimental Biology and John Wiley & Sons Ltd.)

Aiming to contribute toward the characterization of new, biotechnologically relevant cellulolytic enzymes, we report here the first crystal structure of the catalytic core domain of Cel7A (cellobiohydrolase I) from the filamentous fungus Trichoderma harzianum IOC 3844. Our structural studies and molecular dynamics simulations show that the flexibility of Tyr260, in comparison with Tyr247 from the homologous Trichoderma reesei Cel7A, is enhanced as a result of the short side-chains of adjacent Val216 and Ala384 residues and creates an additional gap at the side face of the catalytic tunnel. T. harzianum cellobiohydrolase I also has a shortened loop at the entrance of the cellulose-binding tunnel, which has been described to interact with the substrate in T. reesei Cel7A. These structural features might explain why T. harzianum Cel7A displays higher k(cat) and K(m) values, and lower product inhibition on both glucoside and lactoside substrates, compared with T. reesei Cel7A., (© 2012 The Authors Journal compilation © 2012 FEBS.)

Fullerenes and cubane (C(8)H(8)) can be arranged to form heteromolecular crystals that exhibit interesting crystal phases. Experimental measurements indicate a rotor-stator phase for C(60)-cubane crystals in which the C(60) molecules rotate freely whereas cubane molecules are essentially static. A similar phase is found for C(70)-cubane crystals but, due to C(70)'s asymmetry, hindered rotations can be observed in specific crystal phases. Details of the rotational dynamics of the fullerenes in these heteromolecular crystals are difficult to be completely assessed by experiments. To this end, we have performed classical molecular dynamics simulations of C(70)-cubane crystals to investigate the behavior of C(70) fullerenes and cubanes in the face-centered cubic and body-centered tetragonal crystallographic phases. Our simulations show that, in the cubic phase, C(70) molecules are allowed to freely rotate whereas cubanes act as molecular bearings. In the tetragonal phase, the cubane molecules also remain practically fixed and the rotation of C(70) fullerenes becomes hindered. In this phase, C(70) molecules rotate around the fivefold axis, which in turn precesses about the c crystallographic direction of the unit cell. Details regarding the dynamics (e.g., energy barriers, reorientational relaxation processes, and phonon-libration coupling) of the C(70) molecules in both crystal phases are discussed. In general, our results agree with previous experimental findings for C(70)-cubane crystals.

Understanding the molecular basis of the binding modes of natural and synthetic ligands to nuclear receptors is fundamental to our comprehension of the activation mechanism of this important class of hormone regulated transcription factors and to the development of new ligands. Thyroid hormone receptors (TR) are particularly important targets for pharmaceuticals development because TRs are associated with the regulation of metabolic rates, body weight, and circulating levels of cholesterol and triglycerides in humans. While several high-affinity ligands are known, structural information is only partially available. In this work we obtain structural models of several TR-ligand complexes with unknown structure by docking high affinity ligands to the receptors' ligand binding domain with subsequent relaxation by molecular dynamics simulations. The binding modes of these ligands are discussed providing novel insights into the development of TR ligands. The experimental binding free energies are reasonably well-reproduced from the proposed models using a simple linear interaction energy free-energy calculation scheme.

Nuclear receptors are important targets for pharmaceuticals, but similarities between family members cause difficulties in obtaining highly selective compounds. Synthetic ligands that are selective for thyroid hormone (TH) receptor beta (TRbeta) vs. TRalpha reduce cholesterol and fat without effects on heart rate; thus, it is important to understand TRbeta-selective binding. Binding of 3 selective ligands (GC-1, KB141, and GC-24) is characterized at the atomic level; preferential binding depends on a nonconserved residue (Asn-331beta) in the TRbeta ligand-binding cavity (LBC), and GC-24 gains extra selectivity from insertion of a bulky side group into an extension of the LBC that only opens up with this ligand. Here we report that the natural TH 3,5,3'-triodothyroacetic acid (Triac) exhibits a previously unrecognized mechanism of TRbeta selectivity. TR x-ray structures reveal better fit of ligand with the TRalpha LBC. The TRbeta LBC, however, expands relative to TRalpha in the presence of Triac (549 A(3) vs. 461 A(3)), and molecular dynamics simulations reveal that water occupies the extra space. Increased solvation compensates for weaker interactions of ligand with TRbeta and permits greater flexibility of the Triac carboxylate group in TRbeta than in TRalpha. We propose that this effect results in lower entropic restraint and decreases free energy of interactions between Triac and TRbeta, explaining subtype-selective binding. Similar effects could potentially be exploited in nuclear receptor drug design.

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.

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.

Molecular dynamics simulations have been performed to investigate the interactions between chlorpromazine (CPZ) and Langmuir monolayers of the zwitterionic dipalmitoylphosphatidylcholine (DPPC) and the anionic dipalmitoylphosphatidylglycerol (DPPG). Simulations for a fixed surface density and different charge states - neutral and protonated CPZ - were able to capture important features of the CPZ-phospholipid monolayer interaction. Neutral CPZ is predominantly found in the hydrophobic tail region, whereas protonated CPZ is located at the lipid-water interface. Specific interactions (hydrogen bonds) between protonated CPZ and the lipid head groups were found for both zwitterionic and anionic monolayers. We computed lipid tail order parameters and investigated the effects of the drug upon tail ordering. We also computed electrostatic surface potentials and found qualitative good agreement with experimental results.

Molecular dynamics simulations have been performed to investigate the interactions between chlorpromazine (CPZ), a neuroleptic drug used in the treatment of psychiatric disorders, and dipalmitoylphosphatidylcholine (DPPC), a zwitterionic phospholipid, in Langmuir monolayers. The results from simulations carried out at different monolayer surface densities were able to capture important features of the CPZ-lipid interaction. We find that neutral (unprotonated) CPZ is preferentially located in the lipid tail region of the phospholipids, in little contact with the aqueous phase, and that the orientation of its rigid ring structure and tail conformation vary with lipid surface density. CPZ is found to promote ordering of the lipid tails for all surface densities because of a reduction in the effective surface area per lipid upon addition of the drug. Similar effects have been observed in previous studies of cholesterol in DPPC monolayers, in which lipid tails were seen to order around the solute. This feature, however, is quite distinct from what we observe for the most dense monolayer considered here (area per lipid of 50 A(2)), for which we find that CPZ promotes a local distortion of the lipid tails in its immediate vicinity and a concomitant ordering of lipid tails located further away from the solute. This view is further supported by the results obtained for an approximated nonlinear vibrational sum frequency generation susceptibility, which showed greater tail disorder close to CPZ.

Steered molecular dynamics simulations of ligand dissociation from Thyroid hormone receptors indicate that dissociation is favored via rearrangements in a mobile part of the LBD comprising H3, the loop between H1 and H2, and nearby beta-sheets, contrary to current models in which the H12 is mostly involved. Dissociation is facilitated in this path by the interaction of the hydrophilic part of the ligand with external water molecules, suggesting strategies to enhance ligand binding affinity.