Molecular modelling studies of binding of DACD derivatives into G-Quadruplex DNA: comparison of force field and quantum polarized ligand docking methods

The DNA G-qadruplexes are one of the targets being actively explored for anti-cancer therapy by inhibiting them through small molecules. This computational study was conducted to predict the binding strengths and orientations of a set of novel dimethyl-amino-ethyl-acridine (DACA) analogues that are designed and synthesized in our laboratory, but did not diffract in Synchrotron light.Thecrystal structure of DNA G-Quadruplex(TGGGGT)4(PDB: 1O0K) was used as target for their binding properties in our studies.We used both the force field (FF) and QM/MM derived atomic charge schemes simultaneously for comparing the predictions of drug binding modes and their energetics. This study evaluates the comparative performance of fixed point charge based Glide XP docking and the quantum polarized ligand docking schemes. These results will provide insights on the effects of including or ignoring the drug-receptor interfacial polarization events in molecular docking simulations, which in turn, will aid the rational selection of computational methods at different levels of theory in future drug design programs. Plenty of molecular modelling tools and methods currently exist for modelling drug-receptor or protein-protein, or DNA-protein interactionssat different levels of complexities.Yet, the capasity of such tools to describevarious physico-chemical propertiesmore accuratelyis the next step ahead in currentresearch.Especially, the usage of most accurate methods in quantum mechanics(QM) is severely restricted by theirtedious nature. Though the usage of massively parallel super computing environments resulted in a tremendous improvement in molecular mechanics (MM) calculations like molecular dynamics,they are still capable of dealing with only a couple of tens to hundreds of atoms for QM methods. One such efficient strategy that utilizes thepowers of both MM and QM are the QM/MM hybrid methods. Lately, attempts have been directed towards the goal of deploying several different QM methods for betterment of force field based simulations, but with practical restrictions in place. One of such methods utilizes the inclusion of charge polarization events at the drug-receptor interface, that is not explicitly present in the MM FF.

GSHMC: an efficient method for molecular simulation

The hybrid Monte Carlo (HMC) method is a popular and rigorous method for sampling from a canonical ensemble. The HMC method is based on classical molecular dynamics simulations combined with a Metropolis acceptance criterion and a momentum resampling step. While the HMC method completely resamples the momentum after each Monte Carlo step, the generalized hybrid Monte Carlo (GHMC) method can be implemented with a partial momentum refreshment step. This property seems desirable for keeping some of the dynamic information throughout the sampling process similar to stochastic Langevin and Brownian dynamics simulations. It is, however, ultimate to the success of the GHMC method that the rejection rate in the molecular dynamics part is kept at a minimum. Otherwise an undesirable Zitterbewegung in the Monte Carlo samples is observed. In this paper, we describe a method to achieve very low rejection rates by using a modified energy, which is preserved to high-order along molecular dynamics trajectories. The modified energy is based on backward error results for symplectic time-stepping methods. The proposed generalized shadow hybrid Monte Carlo (GSHMC) method is applicable to NVT as well as NPT ensemble simulations.

Molecular insights into aggregates made of amphiphilic Fmoc-tetrapeptides

Pre-assembled aggregates made of Fmoc-conjugated RGDS and GRDS peptides, where Fmoc refers to fluorenylmethoxycarbonyl, have been investigated using atomistic molecular dynamics simulations. The structural characteristics of twelve different models involving two sheets packed with the Fmoc-aligned or with the charged side groups oriented face-to-face, each one containing seven explicit peptide molecules arranged in parallel or antiparallel, have been evaluated for each Fmoc-tetrapeptide. Structural criteria have been used to select the preferred assembly for each Fmoc-tetrapeptide. The two peptides have been found to prefer b-sheet assemblies with a parallel configuration under simulated low concentration conditions. Furthermore, the assembly is dominated by the interactions among Fmoc units. The overall results provide a complete atomistic view of the interactions between Fmoc-peptide molecules comprised within the same sheet or in different sheets that was not achieved experimentally.

Study of the mechanism of action of anoplin, a helical antimicrobial decapeptide with ion channel-like activity, and the role of the amidated C-terminus

Anoplin, an antimicrobial, helical decapeptide from wasp venom, looses its biological activities by mere deamidation of its C-terminus. Secondary structure determination, by circular dichroism spectroscopy in amphipathic environments, and lytic activity in zwitterionic and anionic vesicles showed quite similar results for the amidated and the carboxylated forms of the peptide. The deamidation of the C-terminus introduced a negative charge at an all-positive charged peptide, causing a loss of amphipathicity, as indicated by molecular dynamics simulations in TFE/water mixtures and this subtle modification in a peptide`s primary structure disturbed the interaction with bilayers and biological membranes. Although being poorly lytic, the amidated form, but not the carboxylated, presented ion channel-like activity on anionic bilayers with a well-defined conductance step; at approximately the same concentration it showed antimicrobial activity. The pores remain open at trans-negative potentials, preferentially conducting cations, and this situation is equivalent to the interaction of the peptide with bacterial membranes that also maintain a high negative potential inside. Copyright (C) 2007 European Peptide Society and John Wiley & Sons, Ltd.

Theoretical investigation of Hf and Zr defects in c-Ge

The introduction of high-permittivity gate dielectric materials into complementary metal oxide semiconductor technology has reopened the interest in Ge as a channel material mainly due to its high hole mobility. Since HfO(2) and ZrO(2) are two of the most promising dielectric candidates, it is important to investigate if Hf and Zr may diffuse into the Ge channel. Therefore, using ab initio density functional theory calculations, we have studied substitutional and interstitial Hf and Zr impurities in c-Ge, looking for neutral defects. We find that (i) substitutional Zr and Hf defects are energetically more favorable than interstitial defects; (ii) under oxygen-rich conditions, neither Zr nor Hf migration towards the channel is likely to occur; (iii) either under Hf- or Zr-rich conditions it is very likely, particularly for Zr, that defects will be incorporated in the channel.

A Multitechnique Study of Structure and Dynamics of Polyfluorene Cast Films and the Influence on Their Photoluminescence

This article describes the microstructure and dynamics in the solid state of polyfluorene-based polymers, poly(9,)-dioctylfluorenyl-2,7-diyl) (PFO), a semicrystalline polymer, and poly [(9,9-dioctyl- 2,7-divinylene-fluorenylene)-alt-co-{2-methoxy-5-(2-ethyl-hexyloxy)- 1,4-phenylene vinylene}, a copolymer with mesomorphic phase properties. These Structures were determined by wide-angle X-ray scattering (WAXS) measurements, Assuming a packing model for the copolymer structure, where the planes of the phenyl rings are stacked and separated by an average distance of similar to 4.5 angstrom and laterally spaced by about similar to 16 angstrom, we followed the evolution of these distances as a function of temperature using WAXS and associated the changes observed to the polymer relaxation processes identified by dynamical mechanical thermal analysis. Specific molecular motions were studied by solid-state nuclear magnetic resonance. The onset of the side-chain motion at about 213 K (beta-relaxation) produced a small increase in the lateral spacing and in the stacking distance of the phenyl rings in them aggregated Structures, Besides, at about 383 K (alpha-relaxation) there occurs a significant increase in the amplitude of the torsion motion in the backbone, producing a greater increase in the stacking distance of the phenyl rings. Similar results were observed in the semicrystalline phase of PFO, but in this case the presence of the crystalline structure affects considerably the overall dynamics, which tends to be more hindered. Put together, Our data explain many features of the temperature dependence of the photoluminescence of these two polymers.

The Mobility of Chondroitin Sulfate in Articular and Artificial Cartilage Characterized by (13)C Magic-Angle Spinning NMR Spectroscopy

We have studied the molecular dynamics of one of the major macromolecules in articular cartilage, chondroitin sulfate. Applying (13)C high-resolution magic-angle spinning NMR techniques, the NMR signals of all rigid macromolecules in cartilage can be suppressed, allowing the exclusive detection of the highly mobile chondroitin sulfate. The technique is also used to detect the chondroitin sulfate in artificial tissue-engineered cartilage. The tissue-engineered material that is based on matrix producing chondrocytes cultured in a collagen gel should provide properties as close as possible to those of the natural cartilage. Nuclear relaxation times of the chondroitin sulfate were determined for both tissues. Although T(1) relaxation times are rather similar, the T(2) relaxation in tissue-engineered cartilage is significantly shorter. This suggests that the motions of chondroitin sulfate in data:rat and artificial cartilage different. The nuclear relaxation times of chondroitin sulfate in natural and tissue-engineered cartilage were modeled using a broad distribution function for the motional correlation times. Although the description of the microscopic molecular dynamics of the chondroitin sulfate in natural and artificial cartilage required the identical broad distribution functions for the correlation times of motion, significant differences in the correlation times of motion that are extracted from the model indicate that the artificial tissue does not fully meet the standards of the natural ideal. This could also be confirmed by macroscopic biomechanical elasticity measurements. Nevertheless, these results suggest that NMR is a useful tool for the investigation of the quality of artificially engineered tissue. (C) 2010 Wiley Periodicals, Inc. Biopolymers 93: 520-532, 2010.

Dynamics of meso and thermo citrate synthases with implicit solvation

The dynamics of hydration of meso and thermo citrate synthases has been investigated using the EEF1 methodology implemented with the CHARNM program. The native enzymes are composed of two identical subunits, each divided into a small and large domain. The dynamics behavior of both enzymes at 30 degrees C and 60 degrees C has been compared. The results of simulations show that during the hydration process, each subunit follows a different pathway of hydration, in spite of the identical sequence. The hydrated structures were compared with the crystalline structure, and the root mean square deviation (RMSD) of each residue along the trajectory was calculated. The regions with larger and smaller mobility were identified. In particular, helices belonging to the small domain are more mobile than those of the large domain. In contrast, the residues that constitute the active site show a much lower displacement compared with the crystalline structure. Hydration free energy calculations point out that Thermoplasma acidophilum citrate synthase (TCS) is more stable than chicken citrate synthase (CCS), at high temperatures. Such result has been ascribed to the higher number of superficial charges in the thermophilic homologue, which stabilizes the enzyme, while the mesophilic homologue denatures. These results are in accord with the experimental found that TCS keeps activity at temperatures farther apart from the catalysis regular temperature than the CCS. (c) 2005 Wiley Periodicals, Inc.

Molecular characterization and phylogenetic analysis of JussuMP-I: A RGD-P-III class hemorrhagic metalloprotease from Bothrops jararacussu snake venom

Snake venom metalloproteases (SVMPs) embody zinc-dependent multidomain enzymes responsible for a relevant pathophysiology in envenomation. including local and systemic hemorrhage. The molecular features responsible for hemorrhagic potency of SVMPs have been associated with their multidomains structures which can target these proteins them to several receptors of different tissues and cellular types. BjussuMP-I. a SVMP isolated from the Bothrops jararacussu venom, has been characterized as a P-III hemorrhagic metalloprotease. The complete cDNA sequence of BjussuMP-I with 1641bp encodes open reading frames of 547 amino acid residues, which conserve the common domains of P-III high molecular weight hemorrhagic metalloproteases: (i) pre-pro-peptide, (ii) metalloprotease, (iii) disintegrin-like and (iv) rich cysteine domain. BjussuMP-I induced lyses in fibrin clots and inhibited collagen- and ADP-induced platelet aggregation. We are reporting, for the first time, the primary structure of an RGD-P-III class snake venom metalloprotease. A phylogenetic analysis of the BjussuMP-1 metalloprotease/catalytic domain was performed to get new insights into the molecular evolution of the metalloproteases. A theoretical molecular model of this domain was built through folding recognition (threading) techniques and refined by molecular dynamics simulation. Then, the final BjussuMP-I catalytic domain model was compared to other SVMPs and Reprolysin family proteins in order to identify eventual structural differences, which could help to understand the biochemical activities of these enzymes. The presence of large hydrophobic areas and some conserved surface charge-positive residues were identified as important features of the SVMPs and other metalloproteases. (C) 2006 Elsevier B.V. All rights reserved.

Cupin: A candidate molecular structure for the Nep1-like protein family

Background: NEP1-like proteins (NLPs) are a novel family of microbial elicitors of plant necrosis. Some NLPs induce a hypersensitive-like response in dicot plants though the basis for this response remains unclear. In addition, the spatial structure and the role of these highly conserved proteins are not known.Results: We predict a 3d-structure for the beta-rich section of the NLPs based on alignments, prediction tools and molecular dynamics. We calculated a consensus sequence from 42 NLPs proteins, predicted its secondary structure and obtained a high quality alignment of this structure and conserved residues with the two Cupin superfamily motifs. The conserved sequence GHRHDWE and several common residues, especially some conserved histidines, in NLPs match closely the two cupin motifs. Besides other common residues shared by dicot Auxin-Binding Proteins (ABPs) and NLPs, an additional conserved histidine found in all dicot ABPs was also found in all NLPs at the same position.Conclusion: We propose that the necrosis inducing protein class belongs to the Cupin superfamily. Based on the 3d-structure, we are proposing some possible functions for the NLPs.

Molecular cloning and biochemical characterization of a myotoxin inhibitor from Bothrops alternatus snake plasma

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

How C-terminal carboxyamidation alters the biological activity of peptides from the venom of the eumenine solitary wasp

Inflammatory peptides display different types of post-transcriptional modifications, such as C-terminal amidation, that alter their biological activity. Here we describe the structural and molecular dynamics features of the mast cell degranulating peptide, eumenine mastoparan-AF (EMP-AF-NH2), found in the venom of the solitary wasp, and of its carboxyl-free C-terminal form (EMP-AF-COO-) characterized by a reduced activity. Circular dichroism indicates that both peptides switch from a random coil conformation in water to a helical structure in TFE and SDS micelles. NMR data, in 30% TFE, reveal that the two peptides fold into an alpha-helix spanning most of their length, while they differ in terms of molecular rigidity. To understand the origins of the conformational flexibility observed in the case of EMP-AF-COO-, a 5 ns MD simulation was carried out for each peptide, in an explicit water/TFE environment. The results show that the two peptides differ in an H-bond between Leu14 NH2 and the backbone carbonyl of Ile11. The loss of that H-bond in EMP-AF-COO- leads to a significant modification of its structural dynamics. In fact, as evidenced by essential dynamics analysis, while EMP-AF-NH2 exists mainly as a rigid structure, EMP-AF-COO- presents two helical stretches that fluctuate in some sort of independent fashion. We conclude that the diverse biological activity of the two peptides is not simply due to the reduction of the net positive charge, as generally suggested, but also to a structural perturbation of the amphipathic alpha-helix that affects their ability to perturb the cell membrane.

SMase II, a new sphingomyelinase D from Loxosceles laeta venom gland: Molecular cloning, expression, function and structural analysis

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Coordinate and time-dependent diffusion dynamics in protein folding

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Study of the mechanism of action of anoplin, a helical antimicrobial decapeptide with ion channel-like activity, and the role of the amidated C-terminus

Anoplin, an antimicrobial, helical decapeptide from wasp venom, looses its biological activities by mere deamidation of its C-terminus. Secondary structure determination, by circular dichroism spectroscopy in amphipathic environments, and lytic activity in zwitterionic and anionic vesicles showed quite similar results for the amidated and the carboxylated forms of the peptide. The deamidation of the C-terminus introduced a negative charge at an all-positive charged peptide, causing a loss of amphipathicity, as indicated by molecular dynamics simulations in TFE/water mixtures and this subtle modification in a peptide's primary structure disturbed the interaction with bilayers and biological membranes. Although being poorly lytic, the amidated form, but not the carboxylated, presented ion channel-like activity on anionic bilayers with a well-defined conductance step; at approximately the same concentration it showed antimicrobial activity. The pores remain open at trans-negative potentials, preferentially conducting cations, and this situation is equivalent to the interaction of the peptide with bacterial membranes that also maintain a high negative potential inside. Copyright (C) 2007 European Peptide Society and John Wiley & Sons, Ltd.