Computational and experimental studies of the catalytic mechanism of Thermobifida fusca cellulase Cel6A (E2)

Mutagenesis experiments suggest that Asp79 in cellulase Cel6A(E2) from Thermobifida fusca has a catalytic role, in spiteof the fact that this residue is more than 13 Å from thescissile bond in models of the enzyme–substrate complexbuilt upon the crystal structure of the protein. This suggeststhat there is a substantial conformational shift in the proteinupon substrate binding. Molecular mechanics simulations wereused to investigate possible alternate conformations of theprotein bound to a tetrasaccharide substrate, primarily involvingshifts of the loop containing Asp79, and to model the role ofwater in the active site complex for both the native conformationand alternative low-energy conformations. Several alternativeconformations of reasonable energy have been identified, includingone in which the overall energy of the enzyme–substratecomplex in solution is lower than that of the conformation inthe crystal structure. This conformation was found to be stablein molecular dynamics simulations with a cellotetraose substrateand water. In simulations of the substrate complexed with thenative protein conformation, the sugar ring in the –1binding site was observed to make a spontaneous transition fromthe ⁴C₁ conformation to a twist-boat conformer, consistent withgenerally accepted glycosidase mechanisms. Also, from thesesimulations Tyr73 and Arg78 were found to have important rolesin the active site. Based on the results of these various MDsimulations, a new catalytic mechanism is proposed. Using thismechanism, predictions about the effects of changes in Arg78were made which were confirmed by site-directed mutagenesis.     

Mechanism of the Conformational Change of the Protein Methyltransferase SMYD3: A Molecular Dynamics Simulation Study

SMYD3 is a SET-domain-containing methyltransferase that catalyzes the transfer of methyl groups onto lysine residues of substrate proteins. Methylation of MAP3K2 by SMYD3 has been implicated in Ras-driven tumorigenesis, which makes SMYD3 a potential target for cancer therapy. Of all SMYD family proteins, SMYD3 adopt a closed conformation in a crystal structure. Several studies have suggested that the conformational changes between the open and closed forms may regulate the catalytic activity of SMYD3. In this work, we carried out extensive molecular dynamics simulations on a series of complexes with a total of 21 μs sampling to investigate the conformational changes of SMYD3 and unveil the molecular mechanisms. Based on the C-terminal domain movements, the simulated models could be depicted in three different conformational states: the closed, intermediate and open states. Only in the case that both the methyl donor binding pocket and the target lysine-binding channel had bound species did the simulations show SMYD3 maintaining its conformation in the closed state, indicative of a synergetic effect of the cofactors and target lysine on regulating the conformational change of SMYD3. In addition, we performed analyses in terms of structure and energy to shed light on how the two regions might regulate the C-terminal domain movement. This mechanistic study provided insights into the relationship between the conformational change and the methyltransferase activity of SMYD3. The more complete understanding of the conformational dynamics developed here together with further work may lay a foundation for the rational drug design of SMYD3 inhibitors.     

Molecular modeling of interleukin-8 receptor beta and analysis of the receptor-ligand interaction

A structural model of interleukin-8 receptor type beta (IL-8R-beta) was constructed based on the structure of bacteriorhodopsin. High temperature molecular dynamics simulations were performed to search the possible conformations of loop regions in IL-8R-beta which recognize the ligand. The crystal structure of interleukin 8 (IL-8) was used as a geometric constraint of the extracellular loop regions of IL-8R-beta in the conformational search. 500 complex structures were extracted from the dynamics trajectory and five plausible models were selected based on the binding energy and known experimental data. To study further the interaction between IL-8R-beta and its ligands, the complex of IL-8R- beta and platelet factor 4 (PF4) C-terminal peptide was also modeled by molecular dynamics simulations. From these models, the N-terminus, extracellular domain 3 and extracellular domain 4 of IL-8R-beta were found to be important for ligand binding. Key residues of these regions involved in ligand binding were characterized. These models provide insight into the structural basis of biological activity of IL-8 and PF4 and may guide the design of potential therapeutic agents targeting IL-8 receptors. Furthermore, the approach developed from this study may have implications for the understanding of other chemokine receptor- ligand interactions that have been recently suggested to be involved in HIV infection.      

Conformational Solvation Studies of LIGNOLs with Molecular Dynamics and Conductor-Like Screening Model

Molecular dynamics (MD) simulations were performed on sterically hindered α-conidendrin-based chiral 1,4-diols (LIGNOLs) from the naturally occurring lignan hydroxymatairesinol (HMR) using the GROMACS software. The aim of this study was to explore the conformational behaviour of the LIGNOLs in aqueous solution adopting the TIP4P model. The topologies of the LIGNOLs were constructed manually and they were modeled with the OPLS-AA force field implemented in GROMACS. The four most relevant torsional angles in the LIGNOLs were properly analyzed during the simulations. The determining property for the conformation preferred in aqueous solution was found to be the lowest energy in gas phase. The solvation effects on the LIGNOLs were also studied by quantum chemical calculations applying the COnductor-like Screening MOdel (COSMO). The hydration studies of the MD simulations showed that several of these LIGNOLs, produced from a renewable source, have a great potential of acting as chiral catalysts.     

Analysis of Binding Determinants for Different Classes of Competitive and Noncompetitive Inhibitors of Glycine Transporters

Glycine transporters are interesting therapeutic targets as they play significant roles in glycinergic and glutamatergic systems. The search for new selective inhibitors of particular types of glycine transporters (GlyT-1 and GlyT-2) with beneficial kinetics is hampered by limited knowledge about the spatial structure of these proteins. In this study, a pool of homology models of GlyT-1 and GlyT-2 in different conformational states was constructed using the crystal structures of related transporters from the SLC6 family and the recently revealed structure of GlyT-1 in the inward-open state, in order to investigate their binding sites. The binding mode of the known GlyT-1 and GlyT-2 inhibitors was determined using molecular docking studies, molecular dynamics simulations, and MM-GBSA free energy calculations. The results of this study indicate that two amino acids, Gly373 and Leu476 in GlyT-1 and the corresponding Ser479 and Thr582 in GlyT-2, are mainly responsible for the selective binding of ligands within the S1 site. Apart from these, one pocket of the S2 site, which lies between TM3 and TM10, may also be important. Moreover, selective binding of noncompetitive GlyT-1 inhibitors in the intracellular release pathway is affected by hydrophobic interactions with Ile399, Met382, and Leu158. These results can be useful in the rational design of new glycine transporter inhibitors with desired selectivity and properties in the future.     

Synthesis, biological evaluation, and molecular modeling studies of rebeccamycin analogues modified in the carbohydrate moiety

A new series of indolocarbazole glycosides containing disaccharides were synthesized and their in vitro antiproliferative activity was evaluated against three human cancer cell lines (A2780, H460, and GLC4). Cytotoxicity appeared to be remarkably affected by the regio- and stereochemical features of the disaccharide moiety. In vivo antitumor activity of the compounds studied, two of which having IC(50)<100 nm, was determined using ovarian cancer cell line A2780 xenografted on nude mice. One compound showed an efficacy similar to that of the reference compound edotecarin, though with a lower long-lasting activity. The topoisomerase I inhibitory properties of some compounds were also examined. Molecular dynamics simulations of the ternary topoisomerase I-DNA-ligand complexes were performed to analyze the structural features of topoisomerase I poisoning with this class of indolocarbazoles. A plausible explanation of their biological behavior was provided. These theoretical results were compared with the recently published crystal structure of an indolocarbazole monosaccharide bound to the covalent human topoisomerase I-DNA complex.     

Computer Simulation and the Design of New Biological Molecules

Molecular dynamics simulations have been carried out for the enzyme bovine trypsin, for the small inhibitor benzamidine, and for the complex formed by these two molecules, all in the presence of enough water to approximate dilute solution conditions. The simulations have been analyzed to characterize the structure and dynamics of the enzyme and its solvent surroundings. Using the recently developed thermodynamic cycle-perturbation method, the simulations have also been analyzed to determine how changes in the structure of the enzyme and the inhibitor alter the free energy of complex formation. The results suggest that such simulations may be useful in the design of new enzymes and enzyme inhibitors.     

Study on the Application of the Combination of TMD Simulation and Umbrella Sampling in PMF Calculation for Molecular Conformational Transitions

Free energy calculations of the potential of mean force (PMF) based on the combination of targeted molecular dynamics (TMD) simulations and umbrella samplings as a function of physical coordinates have been applied to explore the detailed pathways and the corresponding free energy profiles for the conformational transition processes of the butane molecule and the 35-residue villin headpiece subdomain (HP35). The accurate PMF profiles for describing the dihedral rotation of butane under both coordinates of dihedral rotation and root mean square deviation (RMSD) variation were obtained based on the different umbrella samplings from the same TMD simulations. The initial structures for the umbrella samplings can be conveniently selected from the TMD trajectories. For the application of this computational method in the unfolding process of the HP35 protein, the PMF calculation along with the coordinate of the radius of gyration (R_g) presents the gradual increase of free energies by about 1 kcal/mol with the energy fluctuations. The feature of conformational transition for the unfolding process of the HP35 protein shows that the spherical structure extends and the middle α-helix unfolds firstly, followed by the unfolding of other α-helices. The computational method for the PMF calculations based on the combination of TMD simulations and umbrella samplings provided a valuable strategy in investigating detailed conformational transition pathways for other allosteric processes.     

Molecular dynamics simulation of polymers adsorbed onto an alumina surface

This work presents a relatively simple simulation procedure to demonstrate the effects of polymers on an alumina surface. The procedure employs molecular dynamics (MD) techniques to execute real-time simulations on the interactions of polyolefin, polyacrylate, polyoxide, polyol, and polyphenyl linkages with an idealized alumina surface. According to the technique, the adsorption energy is dependent on the geometrical structure of the monomers and decreases for polymer chains with alkyl side-groups in the backbone, but increases for those with functional groups. The results from this simulation procedure indicate that polymer chains with more -CH₂- or functional groups in the framework can markedly increase the adsorption energy. In addition, polyphenyl linkages reveal a wide range of the low-energy region in the rotations of torsional angles. The result is a favorable deformation of the polymer chains with phenylenes in the backbone, thereby leading to a large adsorption energy.     

The Odd Faces of Oligomers: The Case of TRAF2-C,A Trimeric C-Terminal Domain of TNF Receptor-Associated Factor

TNF Receptor Associated Factor 2 (TRAF2) is a trimeric protein that belongs to the TNF receptor associated factor family (TRAFs). The TRAF2 oligomeric state is crucial for receptor binding and for its interaction with other proteins involved in the TNFR signaling. The monomer-trimer equilibrium of a C- terminal domain truncated form of TRAF2 (TRAF2-C), plays also a relevant role in binding the membrane, causing inward vesiculation. In this study, we have investigated the conformational dynamics of TRAF2-C through circular dichroism, fluorescence, and dynamic light scattering, performing temperature-dependent measurements. The data indicate that the protein retains its oligomeric state and most of its secondary structure, while displaying a significative increase in the heterogeneity of the tyrosines signal, increasing the temperature from ≈15 to ≈35 °C. The peculiar crowding of tyrosine residues (12 out of 18) at the three subunit interfaces and the strong dependence on the trimer concentration indicate that such conformational changes mainly involve the contact areas between each pair of monomers, affecting the oligomeric state. Molecular dynamic simulations in this temperature range suggest that the interfaces heterogeneity is an intrinsic property of the trimer that arises from the continuous, asymmetric approaching and distancing of its subunits. Such dynamics affect the results of molecular docking on the external protein surface using receptor peptides, indicating that the TRAF2-receptor interaction in the solution might not involve three subunits at the same time, as suggested by the static analysis obtainable from the crystal structure. These findings shed new light on the role that the TRAF2 oligomeric state might have in regulating the protein binding activity in vivo.     

Towards an improved understanding of the β-TCP crystal structure by means of “checkerboard” atomistic simulations

The large variability of biological responses to β-TCP implants reported in the literature could possibly be related to subtle differences in the β-TCP crystal structure. The structure contains one partially occupied site Ca(4). In order to better understand the ordering of this site, 12 pairs of unit cells with different Ca(4) site occupancies were combined in different checkerboard patterns with an average occupancy of 3. Atomistic simulations were conducted to identify the lowest energy configurations. The previously published low energy configuration is not the most stable one when considering a larger supercell. Plotting the 662 simulation outputs by lattice parameters a or c versus relative lattice energy revealed clusters of high density which are composed of configurations with predominant motifs of Ca(4) occupancy. The tools introduced in this study can be applied in future simulation studies to better explain the Ca(4) site occupancy in the β-TCP crystal structure.     

Energetics for displacing a single chain from the surface of microcrystalline cellulose into the active site of Acidothermus cellulolyticus Cel5A

A series of molecular mechanics calculations were used to analyzethe energetics for moving a single polysaccharide chain fromthe surface of microcrystalline cellulose into the binding cleftof the Cel5A cellulase from Acidothermus cellulolyticus. A build-upprocedure was used to model the placement of a 12-residue oligosaccharidechain along the surface of the enzyme, using as a guide thefour residues of the tetrasaccharide substrate co-crystallizedwith the protein in the crystallographic structure determination.The position of this 12-residue oligosaccharide was used toorient the enzyme properly above two different surfaces of cellulose1ß, the (1,0,0) and the (1,1,0) faces of the crystal.Constrained molecular dynamics simulations were then used topull a target chain directly below the enzyme up out of thecrystal surface and into the binding groove. The energeticsfor this process were favorable for both faces, with the stepface being more favorable than the planar face, implying thatthis surface could be hydrolyzed more readily. Received January 21, 2003; revised October 6, 2003; accepted October 9, 2003     

Nanobody Paratope Ensembles in Solution Characterized by MD Simulations and NMR

Variable domains of camelid antibodies (so-called nanobodies or V_HH) are the smallest antibody fragments that retain complete functionality and therapeutic potential. Understanding of the nanobody-binding interface has become a pre-requisite for rational antibody design and engineering. The nanobody-binding interface consists of up to three hypervariable loops, known as the CDR loops. Here, we structurally and dynamically characterize the conformational diversity of an anti-GFP-binding nanobody by using molecular dynamics simulations in combination with experimentally derived data from nuclear magnetic resonance (NMR) spectroscopy. The NMR data contain both structural and dynamic information resolved at various timescales, which allows an assessment of the quality of protein MD simulations. Thus, in this study, we compared the ensembles for the anti-GFP-binding nanobody obtained from MD simulations with results from NMR. We find excellent agreement of the NOE-derived distance maps obtained from NMR and MD simulations and observe similar conformational spaces for the simulations with and without NOE time-averaged restraints. We also compare the measured and calculated order parameters and find generally good agreement for the motions observed in the ps–ns timescale, in particular for the CDR3 loop. Understanding of the CDR3 loop dynamics is especially critical for nanobodies, as this loop is typically critical for antigen recognition.     

The Influences of Sulphation,Salt Type,and Salt Concentration on the Structural Heterogeneity of Glycosaminoglycans

The increasing recognition of the biochemical importance of glycosaminoglycans (GAGs) has in recent times made them the center of attention of recent research investigations. It became evident that subtle conformational factors play an important role in determining the relationship between the chemical composition of GAGs and their activity. Therefore, a thorough understanding of their structural flexibility is needed, which is addressed in this work by means of all-atom molecular dynamics (MD) simulations. Four major GAGs with different substitution patterns, namely hyaluronic acid as unsulphated GAG, heparan-6-sulphate, chondroitin-4-sulphate, and chondroitin-6-sulphate, were investigated to elucidate the influence of sulphation on the dynamical features of GAGs. Moreover, the effects of increasing NaCl and KCl concentrations were studied as well. Different structural parameters were determined from the MD simulations, in combination with a presentation of the free energy landscape of the GAG conformations, which allowed us to unravel the conformational fingerprints unique to each GAG. The largest effects on the GAG structures were found for sulphation at position 6, as well as binding of the metal ions in the absence of chloride ions to the carboxylate and sulphate groups, which both increase the GAG conformational flexibility.     

Dispersion of a bundle of carbon nanotubes by mechanical torsional energy

Dispersion of a bundle of carbon nanotubes by applying torsional energy is realized by molecular dynamics simulations. The torsional energy applied on the two ends of the bundle leads to a local buckling of the tube structures. The impulse between carbon nano-tubes owing to the local buckling of the tubes forms the driving force through a process of releasing the energy for a successful dispersion. The critical dispersion energy between two carbon nanotubes in a bundle in different solutions and the critical torsional energy for a successful dispersion are calculated by molecular dynamics simulations, and thus the dispersion efficiency is obtained. Effects of different parameters, such as the length and the diameter of the carbon nano-tubes and the temperature of the solution, on the dispersion efficiency are discussed. The dispersion of a bundle of five carbon nano-tubes in an aqueous solution is realized to further prove the feasibility of the proposed dispersion method by torsional energy. This paper provides an effective mechanical approach for dispersion of carbon nano-tubes from their bundle structure.     

Protein dynamics determined by backbone conformation and atom packing

To study the factors determining the collective motions in thermal, conformational fluctuations of a globular protein, molecular dynamics simulations were performed with a backbone model and an atomic-level model. In the backbone model, only the C alpha atoms were explicitly treated with two types of pairwise interactions assigned between the C alpha atoms; atom-packing interactions to take into account the effect of tight atom packing in the protein interior and chain-restoring interactions to maintain the backbone around the native conformation. A quasi-harmonic method was used to decompose the overall fluctuations into independent, collective modes. The modes assigned to large conformational fluctuations showed a good correlation between the backbone and atomic-level models. From this study, it was suggested that the collective modes were motions in which a protein fluctuates, keeping the tertiary structure around the native one and avoiding backbone overlap and, hence, rough aspects of the collective modes can be derived without details of the atomic interactions. The backbone model is useful in obtaining the overall backbone motions of a protein without heavy simulations, even though the simulation starts from a poorly determined conformation of experiments and in sampling main chain conformations, from which the side chain conformations may be predicted.