Enzyme Microheterogeneous Hydration and Stabilization in Supercritical Carbon Dioxide

Supercritical carbon dioxide is a promising green-chemistry solvent for many enzyme-catalyzed chemical reactions, yet the striking stability of some enzymes in such unconventional environments is not well understood. Here, we investigate the stabilization of the Candida antarctica Lipase B (CALB) in supercritical carbon dioxide-water biphasic systems using molecular dynamics simulations. The preservation of the enzyme structure and optimal activity depend on the presence of small amounts of water in the supercritical dispersing medium. When the protein is at least partially hydrated, water molecules bind to specific sites on the enzyme surface and prevent carbon dioxide from penetrating its catalytic core. Strikingly, water and supercritical carbon dioxide cover the protein surface quite heterogeneously. In the first solvation layer, the hydrophilic residues at the surface of the protein are able to pin down patches of water, whereas carbon dioxide solvates preferentially hydrophobic surface residues. In the outer solvation shells, water molecules tend to cluster predominantly on top of the larger water patches of the first solvation layer instead of spreading evenly around the remainder of the protein surface. For CALB, this exposes the substrate-binding region of the enzyme to carbon dioxide, possibly facilitating diffusion of nonpolar substrates into the catalytic funnel. Therefore, by means of microheterogeneous solvation, enhanced accessibility of hydrophobic substrates to the active site can be achieved, while preserving the functional structure of the enzyme. Our results provide a molecular picture on the nature of the stability of proteins in nonaqueous media.

The boundedness of penalty parameters in an augmented Lagrangian method with constrained subproblems

Augmented Lagrangian methods are effective tools for solving large-scale nonlinear programming problems. At each outer iteration, a minimization subproblem with simple constraints, whose objective function depends on updated Lagrange multipliers and penalty parameters, is approximately solved. When the penalty parameter becomes very large, solving the subproblem becomes difficult; therefore, the effectiveness of this approach is associated with the boundedness of the penalty parameters. In this paper, it is proved that under more natural assumptions than the ones employed until now, penalty parameters are bounded. For proving the new boundedness result, the original algorithm has been slightly modified. Numerical consequences of the modifications are discussed and computational experiments are presented.

Platinum-Based Nanoalloys PtnTM55-n (TM=Co, Rh, Au): A Density Functional Theory Investigation

Structural and electronic properties of the PtnTM55-n (TM = Co, Rh, Au) nanoalloys are investigated using density functional theory within the generalized gradient approximation and employing the all-electron projected augmented wave method. For TM = Co and Rh, the excess energy, which measures the relative energy stability of the nanoalloys, is negative for all Pt compositions. We found that the excess energy has similar values for a wide range of Pt compositions, i.e., n = 20-42 and n = 28-42 for Co and Rh, respectively, with the core shell icosahedron-like configuration (n = 42) being slightly more stable for both Co and Rh systems because of the larger release of the strain energy due to the smaller atomic size of the Co and Rh atoms. For TM = Au, the excess energy is positive for all compositions, except for n = 13, which is energetically favorable due to the formation of the core-shell structure (Pt in the core and Au atoms at the surface). Thus, our calculations confirm that the formation of core-shell structures plays an important role to increase the stability of nanoalloys. The center of gravity of the occupied d-states changes almost linearly as a function of the Pt composition, and hence, based on the d-band model, the magnitude of the adsorption energy of an adsorbate can be tuned by changing the Pt composition. The magnetic moments of PtnCo55-n decrease almost linearly as a function of the Pt composition; however, the same does not hold for PtRh and PtAu. We found an enhancement of the magnetic moments of PtRh by a few times by increasing Pt composition, which we explain by the compression effects induced by the large size of the Pt atoms compared with the Rh atoms.

Looking over Toxin-K+ Channel Interactions. Clues from the Structural and Functional Characterization of alpha-KTx Toxin Tc32, a Kv1.3 Channel Blocker

alpha-KTx toxin Tc32, from the Amazonian scorpion Tityus cambridgei, lacks the dyad motif; including Lys27, characteristic of the family and generally associated with channel blockage. The toxin has been cloned and expressed for the first time. Electrophysiological experiments, by showing that the recombinant form blocks Kv1.3 channels of olfactory bulb periglomerular cells like the natural Tc32 toxin, when tested on the Kv1.3 channel of human T lymphocytes, confirmed it is in an active fold. The nuclear magnetic resonance-derived structure revealed it exhibits an alpha/beta scaffold typical of the members of the alpha-KTx family. TdK2 and TdK3, all belonging to the same alpha-KTx 18 subfamily, share significant sequence identity with Tc32 but diverse selectivity and affinity for Kv1.3 and Kv1.1 channels. To gain insight into the structural features that may justify those differences, we used the recombinant Tc32 nuclear magnetic resonance-derived structure to model the other two toxins, for which no experimental structure is available. Their interaction with Kv1.3 and Kv1.1 has been investigated by means of docking simulations. The results suggest that differences in the electrostatic features of the toxins and channels, in their contact surfaces, and in their total dipole moment orientations govern the affinity and selectivity of toxins. In addition, we found that, regardless of whether the dyad motif is present, it is always a Lys side chain that physically blocks the channels, irrespective of its position in the toxin sequence.

The Equilibrium Structure of Lithium Salt Solutions in Ether-Functionalized Ammonium Ionic Liquids

Molecular dynamics simulations have been performed for ionic liquids based on a ternary mixture of lithium and ammonium cations and a common anion, bis(trifluoromethylsulfonyl)imide, [Tf2N](-). We address structural changes resulting from adding Li+ in ionic liquids with increasing length of an ether-functionalized chain in the ammonium cation. The calculation of static structure factors reveals the lithium effect on charge ordering and intermediate range order in comparison with the neat ionic liquids. The charge ordering is modified in the lithium solution because the coordination of [Tf2N](-) toward Li+ is much stronger than ammonium cations. Intermediate range order is observed in neat ionic liquids based on ammonium cations with a long chain, but in the lithium solutions, there is also a nonhomogenous distribution of Li+ cations. The presence of Li+ enhances interactions between the ammonium cations due to correlations between the oxygen atom of the ether chain and the nitrogen atom of another ammonium cation.

Photoisomerization for a model protonated Schiff base in solution: Sloped/peaked conical intersection perspective

The topographical character of conical intersections (CIs)-either sloped or peaked-has played a fundamental and important role in the discussion of the efficiency of CIs as photochemical "funnels." Here this perspective is employed in connection with a recent study of a model protonated Schiff base (PSB) cis to trans photoisomerization in solution [Malhado et al., J. Phys. Chem. A 115, 3720 (2011)]. In that study, the calculated reduced photochemical quantum yield for the successful production of trans product versus cis reactant in acetonitrile solvent compared to water was interpreted in terms of a dynamical solvent effect related to the dominance, for the acetonitrile case, of S-1 to S-0 nonadiabatic transitions prior to the reaching the seam of CIs. The solvent influence on the quantum yield is here re-examined in the sloped/peaked CI topographical perspective via conversion of the model's two PSB internal coordinates and a nonequilibrium solvent coordinate into an effective branching space description, which is then used to re-analyze the generalized Langevin equation/surface hopping results. The present study supports the original interpretation and enriches it in terms of topographical detail. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4754505]

Toll-like receptor 1 N248S single-nucleotide polymorphism is associated with leprosy risk and regulates immune activation during mycobacterial infection

Conflicting findings about the association between leprosy and TLR1 variants N248S and I602S have been reported. Here, we performed case-control and family based studies, followed by replication in 2 case-control populations from Brazil, involving 3162 individuals. Results indicated an association between TLR1 248S and leprosy in the case-control study (SS genotype odds ratio [OR], 1.81; P = .004) and the family based study (z = 2.02; P = .05). This association was consistently replicated in other populations (combined OR, 1.51; P < .001), corroborating the finding that 248S is a susceptibility factor for leprosy. Additionally, we demonstrated that peripheral blood mononuclear cells (PBMCs) carrying 248S produce a lower tumor necrosis factor/interleukin-10 ratio when stimulated with Mycobacterium leprae but not with lipopolysaccharide or PAM3cysK4. The same effect was observed after infection of PBMCs with the Moreau strain of bacillus Calmette-Guerin but not after infection with other strains. Finally, molecular dynamics simulations indicated that the Toll-like receptor 1 structure containing 248S amino acid is different from the structure containing 248N. Our results suggest that TLR1 248S is associated with an increased risk for leprosy, consistent with its hypoimmune regulatory function.

Theoretical Study of the Molecules in Solution and Phospholipid Bilayers

Using Molecular Dynamics simulations and Quantum Mechanical calculations, we study the behavior of molecules with biophysical and pharmacological interest in solution and in phospholipid bilayer.

Movement of NH3 through the human urea transporter B: a new gas channel

Aquaporins and Rh proteins can function as gas (CO2 and NH3) channels. The present study explores the urea, H2O, CO2, and NH3 permeability of the human urea transporter B (UT-B) (SLC14A1), expressed in Xenopus oocytes. We monitored urea uptake using [14C]urea and measured osmotic water permeability (Pf) using video microscopy. To obtain a semiquantitative measure of gas permeability, we used microelectrodes to record the maximum transient change in surface pH (∆pHS) caused by exposing oocytes to 5% CO2/33 mM HCO3- (pHS increase) or 0.5 mM NH3/NH4+ (pHS decrease). UT-B expression increased oocyte permeability to urea by >20-fold, and Pf by 8-fold vs. H2O-injected control oocytes. UT-B expression had no effect on the CO2-induced ∆pHS but doubled the NH3-induced ∆pHS. Phloretin reduced UT-B-dependent urea uptake (Jurea * ) by 45%, Pf * by 50%, and (- ∆pHS * )NH3 by 70%. p-Chloromercuribenzene sulfonate reduced Jurea * by 25%, Pf * by 30%, and (∆pHS * )NH3 by 100%. Molecular dynamics (MD) simulations of membrane-embedded models of UT-B identified the monomeric UT-B pores as the main conduction pathway for both H2O and NH3 and characterized the energetics associated with permeation of these species through the channel. Mutating each of two conserved threonines lining the monomeric urea pores reduced H2O and NH3 permeability. Our data confirm that UT-B has significant H2O permeability and for the first time demonstrate significant NH3 permeability. Thus the UTs become the third family of gas channels. Inhibitor and mutagenesis studies and results of MD simulations suggest that NH3 and H2O pass through the three monomeric urea channels in UT-B.

Atomistic Study of Water Confined in Silica.

In this work, we have used a combined of atomistic simulation methods to explore the effects of confinement of water molecules between silica surfaces. Firstly, the mechanical properties of water severe confined (~3A) between two silica alpha-quartz was determined based on first principles calculations within the density functional theory (DFT). Simulated annealing methods were employed due to the complex potential energry surface, and the difficulties to avoid local minima. Our results suggest that much of the stiffness of the material (46%) remains, even after the insertion of a water monolayer in the silica. Secondly, in order to access typical time scales for confined systems, classical molecular dynamics was used to determine the dynamical properties of water confined in silica cylindrical pores, with diameters varying from 10 to 40A. in this case we have varied the passivation of the silica surface, from 13% to 100% of SiOH, and the other terminations being SiOH2 and SiOH3, the distribution of the different terminations was obtained with a Monte Carlo simulation. The simulations indicates a lowering of the diffusion coefficientes as the diameter decreases, due to the structuration of hydrogen bonds of water molecules; we have also obtained the density profiles of the confined water and the interfacial tension.

Structure-based drug design studies on a series of aldolase inhibitors

Human African trypanosomiasis, also known as sleeping sickness, is a major cause of death in Africa, and for which there are no safe and effective treatments available. The enzyme aldolase from Trypanosoma brucei is an attractive, validated target for drug development. A series of alkyl‑glycolamido and alkyl-monoglycolate derivatives was studied employing a combination of drug design approaches. Three-dimensional quantitative structure-activity relationships (3D QSAR) models were generated using the comparative molecular field analysis (CoMFA). Significant results were obtained for the best QSAR model (r2 = 0.95, non-cross-validated correlation coefficient, and q2 = 0.80, cross-validated correlation coefficient), indicating its predictive ability for untested compounds. The model was then used to predict values of the dependent variables (pKi) of an external test set,the predicted values were in good agreement with the experimental results. The integration of 3D QSAR, molecular docking and molecular dynamics simulations provided further insight into the structural basis for selective inhibition of the target enzyme.

Synthesis, structure and conformation of oligo- and polysaccharides

Carbohydrates are a complex group of biomolecules with a high structural diversity. Their almost omnipresent occurrence has generated a broad field of research in both biology and chemistry. This thesis focuses on three different aspects of carbohydrate chemistry, synthesis, structure elucidation and the conformational analysis of carbohydrates. The first paper describes the synthesis of a penta- and a tetrasaccharide related to the highly branched capsular polysaccharide from Streptococcus pneumoniae type 37. In the second paper, the structure of the O-antigenic repeating unit from the lipopolysaccharide of E. coli 396/C1 was determined along with indications of the structure of the biological repeating unit. In addition, its structural and immunological relationship with E. coli O126 is discussed. In the third paper, partially protected galactopyranosides were examined to clarify the origin of an intriguing 4JHO,H coupling, where a W-mediated coupling pathway was found to operate. In the fourth paper, the conformation of methyl a-cellobioside is studied with a combination of molecular dynamics simulations and NMR spectroscopy. In addition to the expected syn-conformation, detection and quantification of anti-ø and anti-ψ conformers was also possible.

Computational chemistry studies of UV induced processes in human skin

This thesis presents and uses the techniques of computational chemistry to explore two different processes induced in human skin by ultraviolet light. The first is the transformation of urocanic acid into a immunosuppressing agent, and the other is the enzymatic action of the 8-oxoguanine glycosylase enzyme. The photochemistry of urocanic acid is investigated by time-dependent density functional theory. Vertical absorption spectra of the molecule in different forms and environments is assigned and candidate states for the photochemistry at different wavelengths are identified. Molecular dynamics simulations of urocanic acid in gas phase and aqueous solution reveals considerable flexibility under experimental conditions, particularly for for the cis isomer where competition between intra- and inter-molecular interactions increases flexibility. A model to explain the observed gas phase photochemistry of urocanic acid is developed and it is shown that a reinterpretation in terms of a mixture between isomers significantly enhances the agreement between theory and experiment , and resolves several peculiarities in the spectrum. A model for the photochemistry in the aqueous phase of urocanic acid is then developed, in which two excited states governs the efficiency of photoisomerization. The point of entrance into a conical intersection seam is shown to explain the wavelength dependence of photoisomerization quantum yield. Finally some mechanistic aspects of the DNA repair enzyme 8-oxoguanine glycosylase is investigated with density functional theory. It is found that the critical amino acid of the active site can provide catalytic power in several different manners, and that a recent proposal involving a SN1 type of mechanism seems the most efficient one.

Local Structure of Hydrogen-Bonded Liquids

Ordinary yet unique, water is the substance on which life is based. Water seems, at first sight, to be a very simple molecule, consisting of two hydrogen atoms attached to one oxygen. Its small size belies the complexity of its action and its numerous anomalies, central to a broad class of important phenomena, ranging from global current circulation, terrestrial water and CO2 cycles to corrosion and wetting. The explanation of this complex behavior comes from water's unique ability to form extensive three-dimensional networks of hydrogen-bonds, whose nature and structures, in spite of a great deal of efforts involving a plethora of experimental and theoretical techniques, still lacks a complete scientific understanding. This thesis is devoted to the study of the local structure of hydrogen-bonded liquids, with a particular emphasis on water, taking advantage of a combination of core-level spectroscopies and density functional theory spectra calculations. X-ray absorption, in particular, is found to be sensitive to the local hydrogen-bond environment, thus offering a very promising tool for spectroscopic identification of specific structural configurations in water, alcohols and aqueous solutions. More specifically, the characteristic spectroscopic signature of the broken hydrogen-bond at the hydrogen side is used to analyze the structure of bulk water, leading to the finding that most molecules are arranged in two hydrogen-bond configurations, in contrast to the picture provided by molecular dynamics simulations. At the liquid-vapor interface, an interplay of surface sensitive measurements and theoretical calculations enables us to distinguish a new interfacial species in equilibrium with the gas. In a similar approach the cluster form of the excess proton in highly concentrated acid solutions and the different coordination of methanol at the vacuum interface and in the bulk can also be clearly identified. Finally the ability of core-level spectroscopies, aided by sophisticated density functional theory calculations, to directly probe the valence electronic structure of a system is used to observe the nature of the interaction between water molecules and solvated ions in solution. Water around transition metal ions is found to interact with the solute via orbital mixing with the metal d-orbitals. The hydrogen-bond between water molecules is explained in terms of electrostatic interactions enhanced by charge rehybridization in which charge transfer between connecting molecules is shown to be fundamental.

Probing unoccupied electronic states in aqueous solutions

Water is one of the most common compounds on earth and is essential for all biological activities. Water has, however, been a mystery for many years due to the large number of unusual chemical and physical properties, e.g. decreased volume during melting and maximum density at 4 °C. The origin of the anomalies behavior is the nature of the hydrogen bond. This thesis will presented an x-ray absorption spectroscopy (XAS) study to reveal the hydrogen bond structure in liquid water. The x-ray absorption process is faster than a femtosecond and thereby reflects the molecular orbital structure in a frozen geometry locally around the probed water molecules. The results indicate that the electronic structure of liquid water is significantly different from that of the solid and gaseous forms. The molecular arrangement in the first coordination shell of liquid water is actually very similar as the two-hydrogen-bonded configurations at the surface of ice. This discovery suggests that most molecules in liquid water have two-hydrogen-bonded configurations with one donor and one acceptor hydrogen bond compared to the four-hydrogen-bonded tetrahedral structure in ice. This result is controversial since the general picture is that the structure of liquid water is very similar to the structure of ice. The results are, however, consistent with x-ray and neutron diffraction data but reveals serious discrepancies with structures based on current molecular dynamics simulations. The two-hydrogen-bond configuration in liquid water is rigid and heating from 25 °C to 90 °C introduce a minor change in the hydrogen-bonded configurations. Furthermore, XAS studies of water in aqueous solutions show that ion hydration does not affect the hydrogen bond configuration of the bulk. Only water molecules in the close vicinity to the ions show changes in the hydrogen bond formation. XAS data obtained with fluorescence yield are sensitive enough to resolved electronic structure of water molecules in the first hydration sphere and to distinguish between different protonated species. Hence, XAS is a useful tool to provide insight into the local electronic structure of a hydrogen-bonded liquid and it is applied for the first time on water revealing unique information of high importance.