Développement d’une méthode de simulation multi-échelle pour l’étude des grandes transformations dans les protéines

Les films de simulations qui accompagnent le document ont été réalisés avec Pymol.

Molecular computations for reactions and phase transitions: applications to protein stabilization, hydrates and catalysis

In this work we have made significant contributions in three different areas of interest: therapeutic protein stabilization, thermodynamics of natural gas clathrate-hydrates, and zeolite catalysis. In all three fields, using our various computational techniques, we have been able to elucidate phenomena that are difficult or impossible to explain experimentally. More specifically, in mixed solvent systems for proteins we developed a statistical-mechanical method to model the thermodynamic effects of additives in molecular-level detail. It was the first method demonstrated to have truly predictive (no adjustable parameters) capability for real protein systems. We also describe a novel mechanism that slows protein association reactions, called the “gap effect.” We developed a comprehensive picture of methioine oxidation by hydrogen peroxide that allows for accurate prediction of protein oxidation and provides a rationale for developing strategies to control oxidation. The method of solvent accessible area (SAA) was shown not to correlate well with oxidation rates. A new property, averaged two-shell water coordination number (2SWCN) was identified and shown to correlate well with oxidation rates. Reference parameters for the van der Waals Platteeuw model of clathrate-hydrates were found for structure I and structure II. These reference parameters are independent of the potential form (unlike the commonly used parameters) and have been validated by calculating phase behavior and structural transitions for mixed hydrate systems. These calculations are validated with experimental data for both structures and for systems that undergo transitions from one structure to another. This is the first method of calculating hydrate thermodynamics to demonstrate predictive capability for phase equilibria, structural changes, and occupancy in pure and mixed hydrate systems. We have computed a new mechanism for the methanol coupling reaction to form ethanol and water in the zeolite chabazite. The mechanism at 400°C proceeds via stable intermediates of water, methane, and protonated formaldehyde.

Mechanism of action of cyclic antimicrobial peptides

This PhD thesis is the result of the combination of experimental and computational techniques with the aim of understanding the mechanism of action of de novo cyclic decapeptides with high antimicrobial activity. By experimental techniques the influence of the replacement of the phenylalanine for tryptophan residue in their antimicrobial activity was tested and the stability in human serum was also analyzed, in order to evaluate their potential therapeutic application as antitumor agents. On the other hand, the interaction amongst the peptide BPC194 c(KKLKKFKKLQ), the best candidate from the whole library of cyclic peptides, and a model anionic membrane was simulated. The results showed a structure-function relationship derived from the stable conformation of the peptides involved in the membrane permeabilization. As a result, a rational design was performed being BPC490 the peptide with best antimicrobial activity compared with the best active peptide from the original library.

Computer simulations of structures, energetics and dynamics of myoglobin center dot center dot center dot ligand complexes

Myoglobin has been studied in considerable detail using different experimental and computational techniques over the past decades. Recent developments in time-resolved spectroscopy have provided experimental data amenable to detailed atomistic simulations. The main theme of the present review are results on the structures, energetics and dynamics of ligands ( CO, NO) interacting with myoglobin from computer simulations. Modern computational methods including free energy simulations, mixed quantum mechanics/molecular mechanics simulations, and reactive molecular dynamics simulations provide insight into the dynamics of ligand dynamics in confined spaces complementary to experiment. Application of these methods to calculate and understand experimental observations for myoglobin interacting with CO and NO are presented and discussed.

Twists and turns: a tale of two shikimate-pathway enzymes

We are studying two enzymes from the shikimate pathway, dehydroquinate synthase (DHQS) and 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS). Both enzymes have been the subject of numerous studies to elucidate their reaction mechanisms. Crystal structures of DHQS and EPSPS in the presence and absence of substrates, cofactors and/or inhibitors are now available. These structures reveal movements of domains, rearrangements of loops and changes in side-chain positions necessary for the formation of a catalytically competent active site. The potential for using complementary small-angle X-ray scattering (SAXS) studies to confirm the presence of these structural differences in solution has also been explored. Comparative analysis of crystal structures, in the presence and absence of ligands, has revealed structural features critical for substrate-binding and catalysis. We have also analysed these structures by generating GRID energy maps to detect favourable binding sites. The combination of X-ray crystallography, SAXS and computational techniques provides an enhanced analysis of structural features important for the function of these complex enzymes.

Standardized test to evaluate numerical weather prediction algorithms

In order to assist in comparing the computational techniques used in different models, the authors propose a standardized set of one-dimensional numerical experiments that could be completed for each model. The results of these experiments, with a simplified form of the computational representation for advection, diffusion, pressure gradient term, Coriolis term, and filter used in the models, should be reported in the peer-reviewed literature. Specific recommendations are described in this paper.

The interplay of metal-atom ordering, Fermi leveltuning and thermoelectric properties in cobalt shandites Co3M2S2 (M = Sn, In)

A combination of structural, physical and computational techniques including powder X-ray and neutron diffraction, SQUID magnetometry, electrical and thermal transport measurements, DFT calculations and 119Sn Mössbauer and X-ray photoelec-tron spectroscopies has been applied to Co3Sn2-xInxS2 (0 ≤ x ≤ 2) in an effort to understand the relationship between metal-atom ordering and physical properties as the Fermi level is systematically varied. Whilst solid solution behavior is found throughout the composition region, powder neutron diffraction reveals that indium preferentially occupies an inter-layer site over an alternative kagome-like intra-layer site. DFT calculations indicate that this ordering, which leads to a lowering of energy, is related to the dif-fering bonding properties of tin and indium. Spectroscopic data suggest that throughout the composition range 0 ≤ x ≤ 2, all ele-ments adopt oxidation states that are significantly reduced from expectations based on formal charges. Chemical substitution ena-bles the electrical transport properties to be controlled through tuning of the Fermi level within a region of the density of states, which comprises narrow bands of predominantly Co d-character. This leads to a compositionally-induced double metal-to-semiconductor-to-metal transition. The marked increase in the Seebeck coefficient as the semiconducting region is approached leads to a substantial improvement in the thermoelectric figure of merit, ZT, which exhibits a maximum of ZT = 0.32 at 673 K. At 425 K, the figure of merit for phases in the region 0.8 ≤ x ≤ 0.85 is amongst the highest reported for sulphide phases, suggesting these materials may have applications in low-grade waste heat recovery.

Electronic Structure and Spectro-Structural Correlations of Fe(III)Zn(II) Biomimetics for Purple Acid Phosphatases: Relevance to DNA Cleavage and Cytotoxic Activity

Purple acid phosphatases (PAPs) are a group of metallohydrolases that contain a dinuclear Fe(II)M(II) center (M(II) = Fe, Mn, Zn) in the active site and are able to catalyze the hydrolysis of a variety of phosphoric acid esters. The dinuclear complex [(H(2)O)Fe(III)(mu-OH)Zn(II)(L-H)](CIO(4))(2) (2) with the ligand 2-[N-bis(2-pyridylmethyl)aminomethyl]-4-methyl-6-[N-(2-pyridylmethyl)(2-hydroxybenzyl) aminomethyl]phenol (H(2)L-H) has recently been prepared and is found to closely mimic the coordination environment of the Fe(III)Zn(II) active site found in red kidney bean PAP (Neves et al. J. Am. Chem. Soc. 2007, 129, 7486). The biomimetic shows significant catalytic activity in hydrolytic reactions. By using a variety of structural, spectroscopic, and computational techniques the electronic structure of the Fe(III) center of this biomimetic complex was determined. In the solid state the electronic ground state reflects the rhombically distorted Fe(III)N(2)O(4) octahedron with a dominant tetragonal compression align ad along the mu-OH-Fe-O(phenolate) direction. To probe the role of the Fe-O(phenolate) bond, the phenolate moiety was modified to contain electron-donating or -withdrawing groups (-CH(3), -H, -Br, -NO(2)) in the 5-position. Tie effects of the substituents on the electronic properties of the biomimetic complexes were studied with a range of experimental and computational techniques. This study establishes benchmarks against accurate crystallographic struck ral information using spectroscopic techniques that are not restricted to single crystals. Kinetic studies on the hydrolysis reaction revealed that the phosphodiesterase activity increases in the order -NO(2)<- Br <- H <- CH(3) when 2,4-bis(dinitrophenyl)phosphate (2,4-bdnpp) was used as substrate, and a linear free energy relationship is found when log(k(cat)/k(0)) is plotted against the Hammett parameter a. However, nuclease activity measurements in the cleavage of double stranded DNA showed that the complexes containing the electron-withdrawing -NO(2) and electron-donating CH3 groups are the most active while the cytotoxic activity of the biomimetics on leukemia and lung tumoral cells is highest for complexes with electron-donating groups.

Integration of Ligand- and Target-Based Virtual Screening for the Discovery of Cruzain Inhibitors

A myriad of methods are available for virtual screening of small organic compound databases. In this study we have successfully applied a quantitative model of consensus measurements, using a combination of 3D similarity searches (ROCS and EON), Hologram Quantitative Structure Activity Relationships (HQSAR) and docking (FRED, FlexX, Glide and AutoDock Vina), to retrieve cruzain inhibitors from collected databases. All methods were assessed individually and then combined in a Ligand-Based Virtual Screening (LBVS) and Target-Based Virtual Screening (TBVS) consensus scoring, using Receiving Operating Characteristic (ROC) curves to evaluate their performance. Three consensus strategies were used: scaled-rank-by-number, rank-by-rank and rank-by-vote, with the most thriving the scaled-rank-by-number strategy, considering that the stiff ROC curve appeared to be satisfactory in every way to indicate a higher enrichment power at early retrieval of active compounds from the database. The ligand-based method provided access to a robust and predictive HQSAR model that was developed to show superior discrimination between active and inactive compounds, which was also better than ROCS and EON procedures. Overall, the integration of fast computational techniques based on ligand and target structures resulted in a more efficient retrieval of cruzain inhibitors with desired pharmacological profiles that may be useful to advance the discovery of new trypanocidal agents.

Modelos de dispositivos de microondas e ópticos através de redes neurais artificiais de alimentação direta

This dissertation contributes for the development of methodologies through feed forward artificial neural networks for microwave and optical devices modeling. A bibliographical revision on the applications of neuro-computational techniques in the areas of microwave/optical engineering was carried through. Characteristics of networks MLP, RBF and SFNN, as well as the strategies of supervised learning had been presented. Adjustment expressions of the networks free parameters above cited had been deduced from the gradient method. Conventional method EM-ANN was applied in the modeling of microwave passive devices and optical amplifiers. For this, they had been proposals modular configurations based in networks SFNN and RBF/MLP objectifying a bigger capacity of models generalization. As for the training of the used networks, the Rprop algorithm was applied. All the algorithms used in the attainment of the models of this dissertation had been implemented in Matlab

Aplicação da modelagem estrutural computacional direta na análise da geometria de falhas: exemplo da falha de Baixa Grande (Bacia Potiguar - NE do Brasil)

The Baixa grande fault is located on the edge of the S-SW Potiguar Rift. It limits the south part of Umbuzeiro Graben and the Apodi Graben. Although a number of studies have associated the complex deformation styles in the hanging wall of the Baixa Grande Fault with geometry and displacement variations, none have applied the modern computational techniques such as geometrical and kinematic validations to address this problem. This work proposes a geometric analysis of the Baixa Fault using seismic interpretation. The interpretation was made on 3D seismic data of the Baixa Grande fault using the software OpendTect (dGB Earth Sciences). It was also used direct structural modeling, such as Analog Direct Modeling know as Folding Vectors and, 2D and 3D Direct Computational Modeling. The Folding Vectors Modeling presented great similarity with the conventional structural seismic interpretations of the Baixa Grande Fault, thus, the conventional interpretation was validated geometrically. The 2D direct computational modeling was made on some sections of the 3D data of the Baixa Grande Fault on software Move (Midland Valley Ltd) using the horizon modeling tool. The modeling confirms the influence of fault geometry on the hanging wall. The Baixa Grande Fault ramp-flat-ramp geometry generates synform on the concave segments of the fault and antiform in the convex segments. On the fault region that does not have segments angle change, the beds are dislocated without deformation, and on the listric faults occur rollover. On the direct 3D computational modeling, structural attributes were obtained as horizons on the hanging wall of the main fault, after the simulation of several levels of deformation along the fault. The occurrence of structures that indicates shortening in this modeling, also indicates that the antiforms on the Baixa Grande Fault were influenced by fault geometry

Influence of nearest-neighbor Coulomb interactions on the phase diagram of the ferromagnetic Kondo model

The influence of a nearest-neighbor Coulomb repulsion of strength V on the properties of the ferromagnetic Kondo model is analyzed using computational techniques. The Hamiltonian studied here is defined on a chain using localized S = 1/2 spins, and one orbital per site. Special emphasis is given to the influence of the Coulomb repulsion on the regions of phase separation recently discovered in this family of models, as well as on the double-exchange-induced ferromagnetic ground state. When phase separation dominates at V= 0, the Coulomb interaction breaks the large domains of the two competing phases into small islands of one phase embedded into the other. This is in agreement with several experimental results, as discussed in the text. Vestiges of the original phase separation regime are found in the spin structure factor as incommensurate peaks, even at large values of V. In the ferromagnetic regime close to density n = 0.5, the Coulomb interaction induces tendencies to charge ordering without altering the fully polarized character of the state. This regime of charge-ordered ferromagnetism may be related with experimental observations of a similar phase by Chen and Cheong [Phys. Rev. Lett. 76, 4042 (1996)]. Our results reinforce the recently introduced notion [see, e.g., S. Yunoki et al., Phys. Rev. Lett. 80, 845 (1998)] that in realistic models for manganites analyzed with unbiased many-body techniques, the ground state properties arise from a competition between ferromagnetism and phase-separation - charge-ordering tendencies. ©1999 The American Physical Society.

Chronological Monte Carlo-based assessment of distribution system reliability

Regulatory authorities in many countries, in order to maintain an acceptable balance between appropriate customer service qualities and costs, are introducing a performance-based regulation. These regulations impose penalties, and in some cases rewards, which introduce a component of financial risk to an electric power utility due to the uncertainty associated with preserving a specific level of system reliability. In Brazil, for instance, one of the reliability indices receiving special attention by the utilities is the Maximum Continuous Interruption Duration per customer (MCID). This paper describes a chronological Monte Carlo simulation approach to evaluate probability distributions of reliability indices, including the MCID, and the corresponding penalties. In order to get the desired efficiency, modern computational techniques are used for modeling (UML -Unified Modeling Language) as well as for programming (Object- Oriented Programming). Case studies on a simple distribution network and on real Brazilian distribution systems are presented and discussed. © Copyright KTH 2006.

Abordagens de descoberta de conhecimento em bases de dados aplicadas ao Cadastro Territorial Multifinalitário

Pós-graduação em Ciências Cartográficas - FCT

Desenvolvimento de um sistema dinâmico para predição de cargas elétricas por redes neurais através do paradigma de programação orientada a objeto sob a linguagem JAVA

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