Estudos da interação de peptídeos antimicrobianos com modelo de membrana por simulações de dinâmica molecular

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

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)

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)

Application of a new reverse Monte Carlo algorithm to polyatomic molecular systems. I. Liquid water

Using a new reverse Monte Carlo algorithm, we present simulations that reproduce very well several structural and thermodynamic properties of liquid water. Both Monte Carlo, molecular dynamics simulations and experimental radial distribution functions used as input are accurately reproduced using a small number of molecules and no external constraints. Ad hoc energy and hydrogen bond analysis show the physical consistency and limitations of the generated RMC configurations. (C) 2001 American Institute of Physics.

Generalized simulated annealing: Application to silicon clusters

We have compared the recently introduced generalized simulated annealing (GSA) with conventional simulated annealing (CSA). GSA was tested as a tool to obtain the ground-state geometry of molecules. We have used selected silicon clusters (Sin, n=4-7,10) as test cases. Total energies were calculated through tight-binding molecular dynamics. We have found that the replacement of Boltzmann statistics (CSA) by Tsallis's statistics (GSA) has the potential to speed up optimizations with no loss of accuracy. Next, we applied the GSA method to study the ground-state geometry of a 20-atom silicon cluster. We found an original geometry, apparently lower in energy than those previously described in the literature.

Dynamics of a collapsing and exploding Bose-Einstein condensed vortex state

The dynamics of small repulsive Bose-Einstein condensed vortex states of 85Rb atoms in a cylindrical traps with low angular momentum was studied. The time-dependent mean-field Gross-Pitaevskii equation was used for the study. The condensates collapsed and atoms ejected via explosion and a remnant condensate with a smaller number of atoms emerges that survived for a long time.

Electromagnetic form factor of the pion in the space- and time-like regions within the front-form dynamics

The pion electromagnetic form factor is calculated in the space- and time-like regions from -10 (GeV/c)2 up to 10 (GeV/c)2, within a front-form model. The dressed photon vertex where a photon decays in a quark-antiquark pair is depicted generalizing the vector meson dominance ansatz, by means of the vector meson vertex functions. An important feature of our model is the description of the on-mass-shell vertex functions in the valence sector, for the pion and the vector mesons, through the front-form wave functions obtained within a realistic quark model. The theoretical results show an excellent agreement with the data in the space-like region, while in the time-like region the description is quite encouraging. © 2003 Elsevier B.V. All rights reserved.

Multiple minima hypersurfaces procedures for biomimetic ligands screening

In this work the interaction of the pesticide carbaryl with two groups of biomimetic ligands, peptides and MIPs was screened by multiple minima hypersurfaces (MMH) procedures, through the AM1 semiempirical method. Data related to the properties of the molecular association of the complex biomimetic ligand-pesticide were obtained and compared with another molecular modeling algorithm named Leapfrog, as included in the Sybyl software package, and experimental results from the literature, remarking good correlation between them. All important MMH program parameters (cells number, box size, conformers) were studied and optimized with the aim of getting the minimum computation time without losing the correlation with experimental data. The data demonstrated that MMH approach can be used as a fast biomimetic ligand screening tool for MIPs. In the case of peptides the computation time was not comparable with the molecular dynamics methods conventionally used for this approach. © 2011 Springer Science+Business Media B.V.

Kinetic parameters for thermal decomposition of supramolecular polymers derived from diclofenac-meglumine supramolecular adducts

Meglumine is an aminocarbohydrate able to form supramolecular adducts with organic acids. The recognition is based on hydrogen bonds and the structures resulting from the complexation have high solubility in water. This property has been exploited by the pharmaceutical industry in the improvement of existing drugs, and the successful example of this approach involves the poorly soluble non-steroidal anti-inflammatory drugs (NSAIDs). Investigation of the thermal behavior of adduct obtained from meglumine and the NSAID diclofenac revealed that a polymer-like material is formed from the self-assembly of diclofenac-meglumine adducts in the melt. This polymer showed a high molecular weight around 2.0×105kDa. The kinetic parameters for the thermal decomposition step of the polymer were determined by the Capela-Ribeiro non-linear isoconversional method. From data for the TG curves in nitrogen atmosphere and heating rates of 5, 10, 15 and 20°Cmin-1, the Eα and Bα terms could be determined, and consequently the pre-exponential factor, Aα, as well as the kinetic model, g(α). © 2012 Elsevier B.V.

Cooperative RNA Polymerase Molecules Behavior on a Stochastic Sequence-Dependent Model for Transcription Elongation

The transcription process is crucial to life and the enzyme RNA polymerase (RNAP) is the major component of the transcription machinery. The development of single-molecule techniques, such as magnetic and optical tweezers, atomic-force microscopy and single-molecule fluorescence, increased our understanding of the transcription process and complements traditional biochemical studies. Based on these studies, theoretical models have been proposed to explain and predict the kinetics of the RNAP during the polymerization, highlighting the results achieved by models based on the thermodynamic stability of the transcription elongation complex. However, experiments showed that if more than one RNAP initiates from the same promoter, the transcription behavior slightly changes and new phenomenona are observed. We proposed and implemented a theoretical model that considers collisions between RNAPs and predicts their cooperative behavior during multi-round transcription generalizing the Bai et al. stochastic sequence-dependent model. In our approach, collisions between elongating enzymes modify their transcription rate values. We performed the simulations in Mathematica® and compared the results of the single and the multiple-molecule transcription with experimental results and other theoretical models. Our multi-round approach can recover several expected behaviors, showing that the transcription process for the studied sequences can be accelerated up to 48% when collisions are allowed: the dwell times on pause sites are reduced as well as the distance that the RNAPs backtracked from backtracking sites. © 2013 Costa et al.

Modelling polymers with side chains: MEH-PPV and P3HT

Modelling polymers with side chains is always a challenge once the degrees of freedom are very high. In this study, we present a successful methodology to model poly[2-methoxy-5-(2′-ethyl-hexyloxy)-p-phenylenevinylene] (MEH-PPV) and poly[3-hexylthiophene] (P3HT) in solutions, taking into account the influence of side chains on the polymer conformation. Molecular dynamics and semi-empirical quantum mechanical methods were used for structure optimisation and evaluation of optical properties. The methodology allows to describe structural and optical characteristics of the polymers in a satisfactory way, as well as to evaluate some usual simplifications adopted for modelling these systems. Effective conjugation lengths of 8-14.6 and 21 monomers were obtained for MEH-PPV and P3HT, respectively, in accordance with experimental findings. In addition, anti/syn conformations of these polymers could be predicted based on intrinsic interactions of the lateral branches. © 2013 Copyright Taylor and Francis Group, LLC.