Modelagem molecular voltada para a construção de um nanobiossensor para detecção do herbicida Imazaquin

In this study, our goal was develop and describe a molecular model of the enzyme-inhibiting interaction which can be used for an optimized projection of a Microscope Force Atomic nanobiosensor to detect pesticides molecules, used in agriculture, to evaluate its accordance with limit levels stipulated in valid legislation for its use. The studied herbicide (imazaquin) is a typical member of imidazolinone family and is an inhibitor of the enzymatic activity of Acetohydroxiacid Synthase (AHAS) enzyme that is responsible for the first step of pathway for the synthesis of side-chains in amino acids. The analysis of this enzyme property in the presence of its cofactors was made to obtain structural information and charge distribution of the molecular surface to evaluate its capacity of became immobilized on the Microscopy Atomic Force tip. The computational simulation of the system, using Molecular Dynamics, was possible with the force-field parameters for the cofactor and the herbicides obtained by the online tool SwissParam and it was implemented in force-field CHARMM27, used by software GROMACS; then appropriated simulations were made to validate the new parameters. The molecular orientation of the AHAS was defined based on electrostatic map and the availability of the herbicide in the active site. Steered Molecular Dynamics (SMD) Simulations, followed by quantum mechanics calculations for more representative frames, according to the sequential QM/MM methodology, in a specific direction of extraction of the herbicide from the active site. Therefore, external harmonic forces were applied with similar force constants of AFM cantilever for to simulate herbicide detection experiments by the proposed nanobiosensor. Force value of 1391 pN and binding energy of -14048.52 kJ mol-1 were calculated.

The (Un)Folding of Multidomain Proteins Through the Lens of Single-molecule Force-spectroscopy and Computer Simulation

Proteins are specialized molecules that catalyze most of the reactions that can sustain life, and they become functional by folding into a specific 3D structure. Despite their importance, the question, "how do proteins fold?" - first pondered in in the 1930's - is still listed as one of the top unanswered scientific questions as of 2005, according to the journal Science. Answering this question would provide a foundation for understanding protein function and would enable improved drug targeting, efficient biofuel production, and stronger biomaterials. Much of what we currently know about protein folding comes from studies on small, single-domain proteins, which may be quite different from the folding of large, multidomain proteins that predominate the proteomes of all organisms.

In this thesis I will discuss my work to fill this gap in understanding by studying the unfolding and refolding of large, multidomain proteins using the powerful combination of single-molecule force-spectroscopy experiments and molecular dynamic simulations.

The three model proteins studied - Luciferase, Protein S, and Streptavidin - lend insight into the inter-domain dependence for unfolding and the subdomain stabilization of binding ligands, and ultimately provide new insight into atomistic details of the intermediate states along the folding pathway.

Studies of Protein-Protein and Protein-Water Interactions by Small Angle X-Ray Scattering, Terahertz Spectroscopy, ASMOS, And Computer Simulation

The protein folding problem has been one of the most challenging subjects in biological physics due to its complexity. Energy landscape theory based on statistical mechanics provides a thermodynamic interpretation of the protein folding process. We have been working to answer fundamental questions about protein-protein and protein-water interactions, which are very important for describing the energy landscape surface of proteins correctly. At first, we present a new method for computing protein-protein interaction potentials of solvated proteins directly from SAXS data. An ensemble of proteins was modeled by Metropolis Monte Carlo and Molecular Dynamics simulations, and the global X-ray scattering of the whole model ensemble was computed at each snapshot of the simulation. The interaction potential model was optimized and iterated by a Levenberg-Marquardt algorithm. Secondly, we report that terahertz spectroscopy directly probes hydration dynamics around proteins and determines the size of the dynamical hydration shell. We also present the sequence and pH-dependence of the hydration shell and the effect of the hydrophobicity. On the other hand, kinetic terahertz absorption (KITA) spectroscopy is introduced to study the refolding kinetics of ubiquitin and its mutants. KITA results are compared to small angle X-ray scattering, tryptophan fluorescence, and circular dichroism results. We propose that KITA monitors the rearrangement of hydrogen bonding during secondary structure formation. Finally, we present development of the automated single molecule operating system (ASMOS) for a high throughput single molecule detector, which levitates a single protein molecule in a 10 µm diameter droplet by the laser guidance. I also have performed supporting calculations and simulations with my own program codes.

Liquid Mixtures Involving Hydrogenated and Fluorinated Alcohols: Thermodynamics, Spectroscopy, and Simulation

This article reports a combined thermodynamic, spectroscopic, and computational study on the interactions and structure of binary mixtures of hydrogenated and fluorinated substances that simultaneously interact through strong hydrogen bonding. Four binary mixtures of hydrogenated and fluorinated alcohols have been studied, namely, (ethanol + 2,2,2-trifluoroethanol (TFE)), (ethanol + 2,2,3,3,4,4,4-heptafluoro-1-butanol), (1-butanol (BuOH) + TFE), and (BuOH + 2,2,3,3,4,4,4-heptafluoro-1-butanol). Excess molar volumes and vibrational spectra of all four binary mixtures have been measured as a function of composition at 298 K, and molecular dynamics simulations have been performed. The systems display a complex behavior when compared with mixtures of hydrogenated alcohols and mixtures of alkanes and perfluoroalkanes. The combined analysis of the results from different approaches indicates that this results from a balance between preferential hydrogen bonding between the hydrogenated and fluorinated alcohols and the unfavorable dispersion forces between the hydrogenated and fluorinated chains. As the chain length increases, the contribution of dispersion increases and overcomes the contribution of H-bonds. In terms of the liquid structure, the simulations suggest the possibility of segregation between the hydrogenated and fluorinated segments, a hypothesis corroborated by the spectroscopic results. Furthermore, a quantitative analysis of the infrared spectra reveals that the presence of fluorinated groups induces conformational changes in the hydrogenated chains from the usually preferred all-trans to more globular arrangements involving gauche conformations. Conformational rearrangements at the CCOH dihedral angle upon mixing are also disclosed by the spectra.

Vapor pressure and liquid density of fluorinated alcohols:Experimental, simulation and GC-SAFT-VR predictions

The vapor pressure of four liquid 1H,1H-perfluoroalcohols (CF3(CF2)n(CH2)OH, n ¼ 1, 2, 3, 4), often called odd-fluorotelomer alcohols, was measured as a function of temperature between 278 K and 328 K. Liquid densities were also measured for a temperature range between 278 K and 353 K. Molar enthalpies of vaporization were calculated from the experimental data. The results are compared with data from the literature for other perfluoroalcohols as well as with the equivalent hydrogenated alcohols. The results were modeled and interpreted using molecular dynamics simulations and the GC-SAFT-VR equation of state.

Exploiting molecular simulations as an atomic-level microscope to elucidate clinically relevant cation chloride cotransporters

The cation chloride cotransporters (CCCs) represent a vital family of ion transporters, with several members implicated in significant neurological disorders. Specifically, conditions such as cerebrospinal fluid accumulation, epilepsy, Down’s syndrome, Asperger’s syndrome, and certain cancers have been attributed to various CCCs. This thesis delves into these pharmacological targets using advanced computational methodologies. I primarily employed GPU-accelerated all-atom molecular dynamics simulations, deep learning-based collective variables, enhanced sampling methods, and custom Python scripts for comprehensive simulation analyses. Our research predominantly centered on KCC1 and NKCC1 transporters. For KCC1, I examined its equilibrium dynamics in the presence/absence of an inhibitor and assessed the functional implications of different ion loading states. In contrast, our work on NKCC1 revealed its unique alternating access mechanism, termed the rocking-bundle mechanism. I identified a previously unobserved occluded state and demonstrated the transporter's potential for water permeability under specific conditions. Furthermore, I confirmed the actual water flow through its permeable states. In essence, this thesis leverages cutting-edge computational techniques to deepen our understanding of the CCCs, a family of ion transporters with profound clinical significance.

Mechanical Properties And Fracture Dynamics Of Silicene Membranes.

As graphene has become one of the most important materials, there is renewed interest in other similar structures. One example is silicene, the silicon analogue of graphene. It shares some of the remarkable graphene properties, such as the Dirac cone, but presents some distinct ones, such as a pronounced structural buckling. We have investigated, through density functional based tight-binding (DFTB), as well as reactive molecular dynamics (using ReaxFF), the mechanical properties of suspended single-layer silicene. We calculated the elastic constants, analyzed the fracture patterns and edge reconstructions. We also addressed the stress distributions, unbuckling mechanisms and the fracture dependence on the temperature. We analysed the differences due to distinct edge morphologies, namely zigzag and armchair.

Propriedades dinâmicas de fluidos por simulação computacional: métodos híbridos atomístico-contínuo

Computational methods for the calculation of dynamical properties of fluids might consider the system as a continuum or as an assembly of molecules. Molecular dynamics (MD) simulation includes molecular resolution, whereas computational fluid dynamics (CFD) considers the fluid as a continuum. This work provides a review of hybrid methods MD/CFD recently proposed in the literature. Theoretical foundations, basic approaches of computational methods, and dynamical properties typically calculated by MD and CFD are first presented in order to appreciate the similarities and differences between these two methods. Then, methods for coupling MD and CFD, and applications of hybrid simulations MD/CFD, are presented.

Study of the optical and magnetic properties of pyrimidine in water combining PCM and QM/MM methodologies

The solvent effects on the low-lying absorption spectrum and on the (15)N chemical shielding of pyrimidine in water are calculated using the combined and sequential Monte Carlo simulation and quantum mechanical calculations. Special attention is devoted to the solute polarization. This is included by an iterative procedure previously developed where the solute is electrostatically equilibrated with the solvent. In addition, we verify the simple yet unexplored alternative of combining the polarizable continuum model (PCM) and the hybrid QM/MM method. We use PCM to obtain the average solute polarization and include this in the MM part of the sequential QM/MM methodology, PCM-MM/QM. These procedures are compared and further used in the discrete and the explicit solvent models. The use of the PCM polarization implemented in the MM part seems to generate a very good description of the average solute polarization leading to very good results for the n-pi* excitation energy and the (15)N nuclear chemical shield of pyrimidine in aqueous environment. The best results obtained here using the solute pyrimidine surrounded by 28 explicit water molecules embedded in the electrostatic field of the remaining 472 molecules give the statistically converged values for the low lying n-pi* absorption transition in water of 36 900 +/- 100 (PCM polarization) and 36 950 +/- 100 cm(-1) (iterative polarization), in excellent agreement among one another and with the experimental value observed with a band maximum at 36 900 cm(-1). For the nuclear shielding (15)N the corresponding gas-water chemical shift obtained using the solute pyrimidine surrounded by 9 explicit water molecules embedded in the electrostatic field of the remaining 491 molecules give the statistically converged values of 24.4 +/- 0.8 and 28.5 +/- 0.8 ppm, compared with the inferred experimental value of 19 +/- 2 ppm. Considering the simplicity of the PCM over the iterative polarization this is an important aspect and the computational savings point to the possibility of dealing with larger solute molecules. This PCM-MM/QM approach reconciles the simplicity of the PCM model with the reliability of the combined QM/MM approaches.

Hydrogen bond interactions between acetone and supercritical water

Hydrogen bond interactions between acetone and supercritical water are investigated using a combined and sequential Monte Carlo/quantum mechanics (S-MC/QM) approach. Simulation results show a dominant presence of con. gurations with one hydrogen bond for different supercritical states, indicating that this specific interaction plays an important role on the solvation properties of acetone in supercritical water. Using QM MP2/aug-cc-pVDZ the calculated average interaction energy reveals that the hydrogen-bonded acetone-water complex is energetically more stable under supercritical conditions than ambient conditions and its stability is little affected by variations of temperature and/or pressure. All average results reported here are statistically converged.

The isotropic nuclear magnetic shielding constants of acetone in supercritical water: A sequential Monte Carlo/quantum mechanics study including solute polarization

The nuclear isotropic shielding constants sigma((17)O) and sigma((13)C) of the carbonyl bond of acetone in water at supercritical (P=340.2 atm and T=673 K) and normal water conditions have been studied theoretically using Monte Carlo simulation and quantum mechanics calculations based on the B3LYP/6-311++G(2d,2p) method. Statistically uncorrelated configurations have been obtained from Monte Carlo simulations with unpolarized and in-solution polarized solute. The results show that solvent effects on the shielding constants have a significant contribution of the electrostatic interactions and that quantitative estimates for solvent shifts of shielding constants can be obtained modeling the water molecules by point charges (electrostatic embedding). In supercritical water, there is a decrease in the magnitude of sigma((13)C) but a sizable increase in the magnitude of sigma((17)O) when compared with the results obtained in normal water. It is found that the influence of the solute polarization is mild in the supercritical regime but it is particularly important for sigma((17)O) in normal water and its shielding effect reflects the increase in the average number of hydrogen bonds between acetone and water. Changing the solvent environment from normal to supercritical water condition, the B3LYP/6-311++G(2d,2p) calculations on the statistically uncorrelated configurations sampled from the Monte Carlo simulation give a (13)C chemical shift of 11.7 +/- 0.6 ppm for polarized acetone in good agreement with the experimentally inferred result of 9-11 ppm. (C) 2008 American Institute of Physics.

Adiabatic intramolecular movements for water systems

An effective treatment of the intramolecular degrees of freedom is presented for water, where these modes are decoupled from the intermolecular ones, ""adiabatically"" allowing these coordinates to be positioned at their local minimum of the potential energy surface. We perform ab initio Monte Carlo simulations with the configurational energies obtained via density functional theory. We study a water dimer as a prototype system, and even in this simple case the intramolecular relaxations are very important to properly describe properties such as the dipole moment. We show that rigid simulations do not correctly sample the phase space, resulting in an average dipole moment smaller than the one obtained with the adiabatic model, which is closer to the experimental result. (c) 2008 American Institute of Physics.

Contrasting LH-HH subband splitting of strained quantum wells grown along [001] and [113] directions

Contrasting responses for the temperature tuning of the electronic structure in semiconductor quantum wells are discussed for heterolayered structures grown along (001) and (113) directions. The temperature affects the strain modulation of the deformation potentials and the effective optical gap is tuned along with the intersub-band splitting in the valence band. A multiband theoretical model accounts for the characterization of the electronic structure, highlighting the main qualitative and quantitative differences between the two systems under study. The microscopic source of strain fields and the detailed mapping of their distribution are provided by a simulation using classical molecular-dynamics technics.

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

Phospholipases A(2) (PLA(2)s) are important components of Bothrops snake venoms, that can induce several effects on envenomations such as myotoxicity, inhibition or induction of platelet aggregation and edema. It is known that venomous and non-venomous snakes present PLA(2) inhibitory proteins (PLIs) in their blood plasma. An inhibitory protein that neutralizes the enzymatic and toxic activities of several PLA2s from Bothrops venoms was isolated from Bothrops alternatus snake plasma by affinity chromatography using the immobilized myotoxin BthTX-I on CNBr-activated Sepharose. Biochemical characterization of this inhibitory protein, denominated alpha BaltMIP, showed it to be a glycoprotein with Mr of similar to 24,000 for the monomeric subunit. CD spectra of the PLA(2)/inhibitor complexes are considerably different from those corresponding to the individual proteins and data deconvolution suggests that the complexes had a relative gain of helical structure elements in comparison to the individual protomers, which may indicate a more compact structure upon complexation. Theoretical and experimental structural studies performed in order to obtain insights into the structural features of aBaltMIP indicated that this molecule may potentially trimerize in solution, thus strengthening the hypothesis previously raised by other authors about snake PLIs oligomerization. (C) 2010 Elsevier Masson SAS. All rights reserved.

Determination of chemical shielding tensor of an indole carbon and application of tryptophan orientation of a membrane peptide

Using tryptophan C-13-enriched at the C-4 (C epsilon(3)) of the indole, the orientation of the C epsilon(3) chemical shift tensor relative to the C epsilon(3)-H dipolar axis was determined from the C-13 chemical shift/C-13-H-1 dipolar 2D NMR powder pattern. The principal values obtained were 208, 137 and 15 ppm with sigma(33) perpendicular to the indole plane, and sigma(11) (least shielded direction) 5 degrees off the C epsilon(3)-H bond toward C xi(3). The side off the C epsilon(3)-H bond was determined by comparing the reduced chemical shift anisotropies obtained by solid-state NMR and from molecular dynamics calculations of [4-C-13] tryptophans in gramicidin A aligned in phospholipid membranes. (C) 1999 Elsevier Science B.V. All rights reserved.