Molecular dynamics simulations:bring biomolecular structures alive on a computer

The molecular dynamics (MD) simulations play a very important role in science today. They have been used successfully in binding free-energy calculations and rational design of drugs and vaccines. MD simulations can help visualize and understand structures and dynamics at an atomistic level when combined with molecular graphics programs. The molecular and atomistic properties can be displayed on a computer in a time-dependent way, which opens a road toward a better understanding of the relationship of structure, dynamics, and function. In this chapter, the basics of MD are explained, together with a step-by-step description of setup and running an MD simulation.

Anisotropic surface chemistry of aspirin crystals

The wettability of the (001), (100), and (011) crystallographic facets of macroscopic aspirin crystals has been experimentally investigated using a sessile drop contact angle (θ) method. θ for a nonpolar liquid was very similar for all three facets, though significant θ differences were observed for three polar probe liquids. The observed hydrophobicity of the (001) and (100) facets is ascribed to a reduced hydrogen bonding potential at these surfaces, whilst the observed hydrophilicity of facet (011) may be attributed to presence of surface carboxylic functionalities as confirmed by X-ray photoelectron spectroscopy (XPS). The dispersive component of the surface free energy (γ) was similar for all three facets (35 ± 2 mJ/m). The total surface energy, γs varied between 46 and 60 mJ/m due to significant variations in the polar/acid-base components of γ for all facets. Surface polarity as determined by γ measurements and XPS data were in good agreement, linking the variations in wettability to the concentration of oxygen containing surface functional groups. In conclusion, the wettability and the surface energy of a crystalline organic solid, such as aspirin, was found to be anisotropic and facet dependant, and in this case, related to the presence of surface carboxylic functionalities. © 2007 Wiley-Liss, Inc. and the American Pharmacists Association.

An XPS study of pulsed plasma polymerised allyl alcohol film growth on polyurethane

The growth of highly functionalised poly allyl alcohol films by pulsed plasma polymerisation of CH2 double bond; length as m-dashCHCH2OH on biomedical grade polyurethane has been followed by X-ray photoelectron spectroscopy (XPS) and contact angle measurements. Film thickness is observed to increase approximately linearly with plasma modification time, suggesting a layer-by-layer growth mode of poly allyl alcohol. Water contact angle measurements reveal the change in the surface free energy of wetting decreases linearly with plasma modification up to the monolayer point after which a constant limiting value of −24 mJ m−2 was attained. Films prepared at 20 W plasma power with a duty cycle of 10 μs:500 μs exhibit a high degree of hydroxyl (single bondOH) retention with minimal fragmentation of the monomer observed. Increasing the plasma power up to 125 W is found to improve single bondOH retention at the expense of ether formation generating films close to the monomer stoichiometry. Duty cycle plays an important role in controlling both film composition and thickness, with longer off times increasing single bondOH retention, while longer on times enhance allyl alcohol film growth.

Typical performance of regular low-density parity-check codes over general symmetric channels

Typical performance of low-density parity-check (LDPC) codes over a general binary-input output-symmetric memoryless channel is investigated using methods of statistical mechanics. Relationship between the free energy in statistical-mechanics approach and the mutual information used in the information-theory literature is established within a general framework; Gallager and MacKay-Neal codes are studied as specific examples of LDPC codes. It is shown that basic properties of these codes known for particular channels, including their potential to saturate Shannon's bound, hold for general symmetric channels. The binary-input additive-white-Gaussian-noise channel and the binary-input Laplace channel are considered as specific channel models.

Variational mean-field algorithm for efficient inference in large systems of stochastic differential equations

This work introduces a Gaussian variational mean-field approximation for inference in dynamical systems which can be modeled by ordinary stochastic differential equations. This new approach allows one to express the variational free energy as a functional of the marginal moments of the approximating Gaussian process. A restriction of the moment equations to piecewise polynomial functions, over time, dramatically reduces the complexity of approximate inference for stochastic differential equation models and makes it comparable to that of discrete time hidden Markov models. The algorithm is demonstrated on state and parameter estimation for nonlinear problems with up to 1000 dimensional state vectors and compares the results empirically with various well-known inference methodologies.

The thermodynamics and statistical mechanics of protein folding

The physics of self-organization and complexity is manifested on a variety of biological scales, from large ecosystems to the molecular level. Protein molecules exhibit characteristics of complex systems in terms of their structure, dynamics, and function. Proteins have the extraordinary ability to fold to a specific functional three-dimensional shape, starting from a random coil, in a biologically relevant time. How they accomplish this is one of the secrets of life. In this work, theoretical research into understanding this remarkable behavior is discussed. Thermodynamic and statistical mechanical tools are used in order to investigate the protein folding dynamics and stability. Theoretical analyses of the results from computer simulation of the dynamics of a four-helix bundle show that the excluded volume entropic effects are very important in protein dynamics and crucial for protein stability. The dramatic effects of changing the size of sidechains imply that a strategic placement of amino acid residues with a particular size may be an important consideration in protein engineering. Another investigation deals with modeling protein structural transitions as a phase transition. Using finite size scaling theory, the nature of unfolding transition of a four-helix bundle protein was investigated and critical exponents for the transition were calculated for various hydrophobic strengths in the core. It is found that the order of the transition changes from first to higher order as the strength of the hydrophobic interaction in the core region is significantly increased. Finally, a detailed kinetic and thermodynamic analysis was carried out in a model two-helix bundle. The connection between the structural free-energy landscape and folding kinetics was quantified. I show how simple protein engineering, by changing the hydropathy of a small number of amino acids, can enhance protein folding by significantly changing the free energy landscape so that kinetic traps are removed. The results have general applicability in protein engineering as well as understanding the underlying physical mechanisms of protein folding. ^

Isolation and characterization of the U1 and U2 snRNAs in the silk gland of the 5th instar larval stage silk moth Bombyx mori

In isolation and characterization studies, expression level U1 and U2 snRNA isoforms were obtained from the 5th instar larval stage silk gland (SG). The DNA content of the SG cells is approximately 200,000-fold higher compared to the usual (2N) somatic cells of B. mori due to endoreduplication. In this study, the existence of U1 and U2 snRNA isoforms in the SG of the organism is investigated. Bombyx mori U1 and U2-specific RT-PCR libraries from the silk gland were generated. Five U1 and eight U2 isoforms were isolated and characterized. Nucleotide differences, structural alterations, as well as protein and RNA interaction sites were analyzed in these variants. For the U1 snRNA variants, they were compared to the previously reported BmN isoforms. In all these U-snRNA variants, polymorphic sites do not predominate at the core of known functional sequences, which were interspecifically conserved. Variant sites and inter-species differences are located in moderately conserved regions. Free energy (ΔG) values for the entire U1 and U2 snRNA secondary structures and for the individual stem/loops domains of the isoforms were generated and compared to determine their structural stability. This will be the first time that U1 and U2 variants are shown specific for a development stage (larval) other than embryonic or adult. ^ Using phylogenetic analysis, evolutionary trees were generated for the U1 and U2 snRNAs using animal, plant, protista and fungal species. The resulting trees were boostrapped for robustness and rooted with the self-splicing RNA group II intron sequence from the cyanobacterium Calothrix. Using phylogenetic analyses, possible structural and functional evolutionary interdependence between the U1 and U2 snRNAs was investigated. ^

Structural flexibility and oxygen diffusion pathways in monomeric fluorescent proteins

Fluorescent proteins are valuable tools as biochemical markers for studying cellular processes. Red fluorescent proteins (RFPs) are highly desirable for in vivo applications because they absorb and emit light in the red region of the spectrum where cellular autofluorescence is low. The naturally occurring fluorescent proteins with emission peaks in this region of the spectrum occur in dimeric or tetrameric forms. The development of mutant monomeric variants of RFPs has resulted in several novel FPs known as mFruits. Though oxygen is required for maturation of the chromophore, it is known that photobleaching of FPs is oxygen sensitive, and oxygen-free conditions result in improved photostabilities. Therefore, understanding oxygen diffusion pathways in FPs is important for both photostabilites and maturation of the chromophores. We used molecular dynamics calculations to investigate the protein barrel fluctuations in mCherry, which is one of the most useful monomeric mFruit variants, and its GFP homolog citrine. We employed implicit ligand sampling and locally enhanced sampling to determine oxygen pathways from the bulk solvent into the mCherry chromophore in the interior of the protein. The pathway contains several oxygen hosting pockets, which were identified by the amino acid residues that form the pocket. We calculated the free-energy of an oxygen molecule at points along the path. We also investigated an RFP variant known to be significantly less photostable than mCherry and find much easier oxygen access in this variant. We showed that oxygen pathways can be blocked or altered, and barrel fluctuations can be reduced by strategic amino acid substitutions. The results provide a better understanding of the mechanism of molecular oxygen access into the fully folded mCherry protein barrel and provide insight into the photobleaching process in these proteins.

Estudo teórico de intermediários tetraédricos acidez / basicidade e estereosseletividade enzimáticos

The present work aimed first, the theoretical study of tetrahedral intermediate stability formed from carbonyl addition reactions using the second (MP2) and third (MP3) order Møller–Plesset perturbation theory. Linear correlations between electronic energy difference of reactions with Wiberg Indexes and C-O bond lengths were obtained, and was observed that the stability of adducts formed depends directly of electronic density involved between these atoms. The knowing of electronic parameters of these structures has an important hole due to the large use on reactions that in his course forms this tetrahedral intermediate. Employing the ONIOM (B3LYP:AMBER) methodology, was evaluated the stereoselectivity of a enzymatic reaction between CAL B enzyme and a long chain ester. In this study, were obtained the electronic energies of ground state and intermediate state of transesterification rate-determing step from two possible proquirals faces Re and Si. The objective was study the enantioselectivity of CAL B and rationalizes it using quantum theory of atoms in molecules (QTAIM). A theoretical study employing inorganic compounds was performed using ab initio CBS-QB3 method aiming to find a link between thermodynamic and equilibrium involving acids and bases. The results observed showed an excellent relationship between difference in Gibbs free energy, ΔG of acid dissociation reaction and ΔG of hydrolysis reaction of the corresponding conjugate base. It was also observed, a relationship between ΔG of hydrolysis reaction of conjugate acids and their corresponding atomic radius showing that stability plays an important role in hydrolysis reactions. The importance of solvation in acid/base behavior when compared to theoretical and experimental ΔG´s also was evaluated.

Síntese e caracterização de materiais mesoporosos para a captura de CO2: influência do óxido de níquel

Several materials are currently under study for the CO2 capture process, like the metal oxides and mixed metal oxides, zeolites, carbonaceous materials, metal-organic frameworks (MOF's) organosilica and modified silica surfaces. In this work, evaluated the adsorption capacity of CO2 in mesoporous materials of different structures, such as MCM-48 and SBA- 15 without impregnating and impregnated with nickel in the proportions 5 %, 10 % and 20 % (m/m), known as 5Ni-MCM-48, 10Ni-MCM-48, 20Ni-MCM-48 and 5Ni-SBA-15, 10NiSBA-15, 20Ni-SBA-15. The materials were characterized by means of X-ray diffraction (XRD), thermal analysis (TG and DTG), Fourier transform infrared spectroscopy (FT-IR), N2 adsorption and desorption (BET) and scanning electron microscopy (SEM) with EDS. The adsorption process was performed varying the pressure of 100 - 4000 kPa and keeping the temperature constant and equal to 298 K. At a pressure of 100 kPa, higher concentrations of adsorption occurred for the materials 5Ni-MCM-48 (0.795 mmol g-1 ) and SBA-15 (0.914 mmol g-1 ) is not impregnated, and at a pressure of 4000 kPa for MCM-48 materials (14.89 mmol g-1) and SBA-15 (9.97 mmol g-1) not impregnated. The results showed that the adsorption capacity varies positively with the specific area, however, has a direct dependency on the type and geometry of the porous structure of channels. The data were fitted using the Langmuir and Freundlich models and were evaluated thermodynamic parameters Gibbs free energy and entropy of the adsorption system

Contribuição ao estudo de treliças não lineares considerando efeitos térmicos

This work presents the numerical analysis of nonlinear trusses summited to thermomechanical actions with Finite Element Method (FEM). The proposed formulation is so-called positional FEM and it is based on the minimum potential energy theorem written according to nodal positions, instead of displacements. The study herein presented considers the effects of geometric and material nonlinearities. Related to dynamic problems, a comparison between different time integration algorithms is performed. The formulation is extended to impact problems between trusses and rigid wall, where the nodal positions are constrained considering nullpenetration condition. In addition, it is presented a thermodynamically consistent formulation, based on the first and second law of thermodynamics and the Helmholtz free-energy for analyzing dynamic problems of truss structures with thermoelastic and thermoplastic behavior. The numerical results of the proposed formulation are compared with examples found in the literature.

Electrostatic Contributions to the Thermodynamics of Ribonuclease P Protein Folding

Electrostatic interactions are of fundamental importance in determining the structure and stability of macromolecules. For example, charge-charge interactions modulate the folding and binding of proteins and influence protein solubility. Electrostatic interactions are highly variable and can be both favorable and unfavorable. The ability to quantify these interactions is challenging but vital to understanding the detailed balance and major roles that they have in different proteins and biological processes. Measuring pKa values of ionizable groups provides a sensitive method for experimentally probing the electrostatic properties of a protein.

pKa values report the free energy of site-specific proton binding and provide a direct means of studying protein folding and pH-dependent stability. Using a combination of NMR, circular dichroism, and fluorescence spectroscopy along with singular value decomposition, we investigated the contributions of electrostatic interactions to the thermodynamic stability and folding of the protein subunit of Bacillus subtilis ribonuclease P, P protein. Taken together, the results suggest that unfavorable electrostatics alone do not account for the fact that P protein is intrinsically unfolded in the absence of ligand because the pKa differences observed between the folded and unfolded state are small. Presumably, multiple factors encoded in the P protein sequence account for its IUP property, which may play an important role in its function.

Phase coexistence of cluster crystals: beyond the Gibbs phase rule

We report a study of the phase behavior of multiple-occupancy crystals through simulation. We argue that in order to reproduce the equilibrium behavior of such crystals it is essential to treat the number of lattice sites as a constraining thermodynamic variable. The resulting free-energy calculations thus differ considerably from schemes used for single-occupancy lattices. Using our approach, we obtain the phase diagram and the bulk modulus for a generalized exponential model that forms cluster crystals at high densities. We compare the simulation results with existing theoretical predictions. We also identify two types of density fluctuations that can lead to two sound modes and evaluate the corresponding elastic constants.

Numerical detection of the Gardner transition in a mean-field glass former.

Recent theoretical advances predict the existence, deep into the glass phase, of a novel phase transition, the so-called Gardner transition. This transition is associated with the emergence of a complex free energy landscape composed of many marginally stable sub-basins within a glass metabasin. In this study, we explore several methods to detect numerically the Gardner transition in a simple structural glass former, the infinite-range Mari-Kurchan model. The transition point is robustly located from three independent approaches: (i) the divergence of the characteristic relaxation time, (ii) the divergence of the caging susceptibility, and (iii) the abnormal tail in the probability distribution function of cage order parameters. We show that the numerical results are fully consistent with the theoretical expectation. The methods we propose may also be generalized to more realistic numerical models as well as to experimental systems.