The importance of reactant positioning in enzyme catalysis: A hybrid quantum mechanics/molecular mechanics study of a haloalkane dehalogenase

Hybrid quantum mechanics/molecular mechanics calculations using Austin Model 1 system-specific parameters were performed to study the SN2 displacement reaction of chloride from 1,2-dichloroethane (DCE) by nucleophilic attack of the carboxylate of acetate in the gas phase and by Asp-124 in the active site of haloalkane dehalogenase from Xanthobacter autotrophicus GJ10. The activation barrier for nucleophilic attack of acetate on DCE depends greatly on the reactants having a geometry resembling that in the enzyme or an optimized gas-phase structure. It was found in the gas-phase calculations that the activation barrier is 9 kcal/mol lower when dihedral constraints are used to restrict the carboxylate nucleophile geometry to that in the enzyme relative to the geometries for the reactants without dihedral constraints. The calculated quantum mechanics/molecular mechanics activation barriers for the enzymatic reaction are 16.2 and 19.4 kcal/mol when the geometry of the reactants is in a near attack conformer from molecular dynamics and in a conformer similar to the crystal structure (DCE is gauche), respectively. This haloalkane dehalogenase lowers the activation barrier for dehalogenation of DCE by 2–4 kcal/mol relative to the single point energies of the enzyme's quantum mechanics atoms in the gas phase. SN2 displacements of this sort in water are infinitely slower than in the gas phase. The modest lowering of the activation barrier by the enzyme relative to the reaction in the gas phase is consistent with mutation experiments.

Large vortex-like structure of dipole field in computer models of liquid water and dipole-bridge between biomolecules

We propose a framework to describe the cooperative orientational motions of water molecules in liquid water and around solute molecules in water solutions. From molecular dynamics (MD) simulation a new quantity “site-dipole field” is defined as the averaged orientation of water molecules that pass through each spatial position. In the site-dipole field of bulk water we found large vortex-like structures of more than 10 Å in size. Such coherent patterns persist more than 300 ps although the orientational memory of individual molecules is quickly lost. A 1-ns MD simulation of systems consisting of two amino acids shows that the fluctuations of site-dipole field of solvent are pinned around the amino acids, resulting in a stable dipole-bridge between side-chains of amino acids. The dipole-bridge is significantly formed even for the side-chain separation of 14 Å, which corresponds to five layers of water. The way that dipole-bridge forms sensitively depends on the side-chain orientations and thereby explains the specificity in the solvent-mediated interactions between biomolecules.

Quantification of the hydrophobic interaction by simulations of the aggregation of small hydrophobic solutes in water

The hydrophobic interaction, the tendency for nonpolar molecules to aggregate in solution, is a major driving force in biology. In a direct approach to the physical basis of the hydrophobic effect, nanosecond molecular dynamics simulations were performed on increasing numbers of hydrocarbon solute molecules in water-filled boxes of different sizes. The intermittent formation of solute clusters gives a free energy that is proportional to the loss in exposed molecular surface area with a constant of proportionality of 45 ± 6 cal/mol⋅Å2. The molecular surface area is the envelope of the solute cluster that is impenetrable by solvent and is somewhat smaller than the more traditional solvent-accessible surface area, which is the area transcribed by the radius of a solvent molecule rolled over the surface of the cluster. When we apply a factor relating molecular surface area to solvent-accessible surface area, we obtain 24 cal/mol⋅Å2. Ours is the first direct calculation, to our knowledge, of the hydrophobic interaction from molecular dynamics simulations; the excellent qualitative and quantitative agreement with experiment proves that simple van der Waals interactions and atomic point-charge electrostatics account for the most important driving force in biology.

The arginine finger of RasGAP helps Gln-61 align the nucleophilic water in GAP-stimulated hydrolysis of GTP

The Ras family of GTPases is a collection of molecular switches that link receptors on the plasma membrane to signaling pathways that regulate cell proliferation and differentiation. The accessory GTPase-activating proteins (GAPs) negatively regulate the cell signaling by increasing the slow intrinsic GTP to GDP hydrolysis rate of Ras. Mutants of Ras are found in 25–30% of human tumors. The most dramatic property of these mutants is their insensitivity to the negative regulatory action of GAPs. All known oncogenic mutants of Ras map to a small subset of amino acids. Gln-61 is particularly important because virtually all mutations of this residue eliminate sensitivity to GAPs. Despite its obvious importance for carcinogenesis, the role of Gln-61 in the GAP-stimulated GTPase activity of Ras has remained a mystery. Our molecular dynamics simulations of the p21ras–p120GAP–GTP complex suggest that the local structure around the catalytic region can be different from that revealed by the x-ray crystal structure. We find that the carbonyl oxygen on the backbone of the arginine finger supplied in trans by p120GAP (Arg-789) interacts with a water molecule in the active site that is forming a bridge between the NH2 group of the Gln-61 and the γ-phosphate of GTP. Thus, Arg-789 may play a dual role in generating the nucleophile as well as stabilizing the transition state for P—O bond cleavage.

Intrinsic compressibility and volume compression in solvated proteins by molecular dynamics simulation at high pressure.

Constant pressure and temperature molecular dynamics techniques have been employed to investigate the changes in structure and volumes of two globular proteins, superoxide dismutase and lysozyme, under pressure. Compression (the relative changes in the proteins' volumes), computed with the Voronoi technique, is closely related with the so-called protein intrinsic compressibility, estimated by sound velocity measurements. In particular, compression computed with Voronoi volumes predicts, in agreement with experimental estimates, a negative bound water contribution to the apparent protein compression. While the use of van der Waals and molecular volumes underestimates the intrinsic compressibilities of proteins, Voronoi volumes produce results closer to experimental estimates. Remarkably, for two globular proteins of very different secondary structures, we compute identical (within statistical error) protein intrinsic compressions, as predicted by recent experimental studies. Changes in the protein interatomic distances under compression are also investigated. It is found that, on average, short distances compress less than longer ones. This nonuniform contraction underlines the peculiar nature of the structural changes due to pressure in contrast with temperature effects, which instead produce spatially uniform changes in proteins. The structural effects observed in the simulations at high pressure can explain protein compressibility measurements carried out by fluorimetric and hole burning techniques. Finally, the calculation of the proteins static structure factor shows significant shifts in the peaks at short wavenumber as pressure changes. These effects might provide an alternative way to obtain information concerning compressibilities of selected protein regions.

The human AQP4 gene: definition of the locus encoding two water channel polypeptides in brain.

The aquaporin family of membrane water transport proteins are expressed in diverse tissues, and in brain the predominant water channel protein is AQP4. Here we report the isolation and characterization of the human AQP4 cDNAs and genomic DNA. Two cDNAs were isolated corresponding to the two initiating methionines (M1 in a 323-aa polypeptide and M23 in a 301-aa polypeptide) previously identified in rat [Jung, J.S., Bhat, R.V., Preston, G.M., Guggino, W.B. & Agre, P. (1994) Proc. Natl. Acad. Sci. USA 91, 13052-13056]. Similar to other aquaporins, the AQP4 gene is composed of four exons encoding 127, 55, 27, and 92 amino acids separated by introns of 0.8, 0.3, and 5.2 kb. Unlike other aquaporins, an alternative coding initiation sequence (designated exon 0) was located 2.7 kb upstream of exon 1. When spliced together, M1 and the subsequent 10 amino acids are encoded by exon 0; the next 11 amino acids and M23 are encoded by exon 1. Transcription initiation sites have been mapped in the proximal promoters of exons 0 and 1. RNase protection revealed distinct transcripts corresponding to M1 and M23 mRNAs, and AQP4 immunoblots of cerebellum demonstrated reactive polypeptides of 31 and 34 kDa. Using a P1 and a lambda EMBL subclone, the chromosomal site of the human AQP4 gene was mapped to chromosome 18 at the junction of q11.2 and q12.1 by fluorescence in situ hybridization. These studies may now permit molecular characterization of AQP4 during human development and in clinical disorders.

Glass transition in DNA from molecular dynamics simulations.

Molecular dynamics simulations of the oligonucleotide duplex d(CGCGCG)2 in aqueous solution are used to investigate the glass transition phenomenon. The simulations were performed at temperatures in the 20 K to 340 K range. The mean square atomic fluctuations showed that the behavior of the oligonucleotide duplex was harmonic at low temperatures. A glass transition temperature at 223 K to 234 K was inferred for the oligonucleotide duplex, which is in agreement with experimental observations. The largest number of hydrogen bounds between the polar atoms of the oligonucleotide duplex and the water molecules was obtained at the glass transition temperature. With increasing temperature we observed a decrease in the average lifetime of the hydrogen bonds to water molecules.

Crystal and molecular structure of paclitaxel (taxol).

Paclitaxel (formerly called taxol), an important anticancer drug, inhibits cell replication by binding to and stabilizing microtubule polymers. As drug-receptor interactions are governed by the three-dimensional stereochemistries of both participants, we have determined the crystal structure of paclitaxel to identify its conformational preferences that may be related to biological activity. The monoclinic crystals contain two independent paclitaxel molecules in the asymmetric unit plus several water and dioxane solvent molecules. Taxane ring conformation is very similar in both paclitaxel molecules and is similar to the taxane ring conformation found in the crystal structure of the paclitaxel analogue docetaxel (formerly called taxotere). The two paclitaxel molecules have carbon-13 side-chain conformations that differ from each other and from that of the corresponding side chain in the docetaxel crystal structure. The carbon-13 side-chain conformation of one paclitaxel molecule is similar to what was proposed from NMR studies done in polar solvents, while that of the other paclitaxel molecule is different and hitherto unobserved. The paclitaxel molecules interact with each other and with solvent atoms through an extensive network of hydrogen bonds. Analysis of the hydrogen-bonding network together with structure-activity studies may suggest which atoms of paclitaxel are important for binding to microtubule receptors.

Immunolocalization of the mercurial-insensitive water channel and glycerol intrinsic protein in epithelial cell plasma membranes.

Two water channel homologs were cloned recently from rat kidney, mercurial-insensitive water channel (MIWC) and glycerol intrinsic protein (GLIP). Polyclonal antibodies were raised against synthetic C-terminal peptides and purified by affinity chromatography. MIWC and GLIP antibodies recognized proteins in rat kidney with an apparent molecular mass of 30 and 27 kDa, respectively, and did not cross-react. By immunofluorescence, MIWC and GLIP were expressed together on the basolateral plasma membrane of collecting duct principal cells in kidney. By immunohistochemistry, MIWC and GLIP were expressed on tracheal epithelial cells with greater expression of GLIP on the basal plasma membrane and MIWC on the lateral membrane; only MIWC was expressed in bronchial epithelia. In eye, GLIP was expressed in conjunctival epithelium, whereas MIWC was found in iris, ciliary body, and neural cell layers in retina. MIWC and GLIP colocalized on the basolateral membrane of villus epithelial cells in colon and brain ependymal cells. Expression of MIWC and GLIP was not detected in small intestine, liver, spleen, endothelia, and cells that express water channels CHIP28 or WCH-CD. These studies suggest water/solute transporting roles for MIWC and GLIP in the urinary concentrating mechanism, cerebrospinal fluid absorption, ocular fluid balance, fecal dehydration, and airway humidification. The unexpected membrane colocalization of MIWC and GLIP in several tissues suggests an interaction at the molecular and/or functional levels.

Interfacial activation-based molecular bioimprinting of lipolytic enzymes.

Interfacial activation-based molecular (bio)-imprinting (IAMI) has been developed to rationally improve the performance of lipolytic enzymes in nonaqueous environments. The strategy combinedly exploits (i) the known dramatic enhancement of the protein conformational rigidity in a water-restricted milieu and (ii) the reported conformational changes associated with the activation of these enzymes at lipid-water interfaces, which basically involves an increased substrate accessibility to the active site and/or an induction of a more competent catalytic machinery. Six model enzymes have been assayed in several model reactions in nonaqueous media. The results, rationalized in light of the present biochemical and structural knowledge, show that the IAMI approach represents a straightforward, versatile method to generate manageable, activated (kinetically trapped) forms of lipolytic enzymes, providing under optimal conditions nonaqueous rate enhancements of up to two orders of magnitude. It is also shown that imprintability of lipolytic enzymes depends not only on the nature of the enzyme but also on the "quality" of the interface used as the template.

Heterogeneity in molecular recognition by restriction endonucleases: osmotic and hydrostatic pressure effects on BamHI, Pvu II, and EcoRV specificity.

The cleavage specificity of the Pvu II and BamHI restriction endonucleases is found to be dramatically reduced at elevated osmotic pressure. Relaxation in specificity of these otherwise highly accurate and specific enzymes, previously termed "star activity," is uniquely correlated with osmotic pressure between 0 and 100 atmospheres. No other colligative solvent property exhibits a uniform correlation with star activity for all of the compounds tested. Application of hydrostatic pressure counteracts the effects of osmotic pressure and restores the natural selectivity of the enzymes for their canonical recognition sequences. These results indicate that water solvation plays an important role in the site-specific recognition of DNA by many restriction enzymes. Osmotic pressure did not induce an analogous effect on the specificity of the EcoRV endonuclease, implying that selective hydration effects do not participate in DNA recognition in this system. Hydrostatic pressure was found to have little effect on the star activity induced by changes in ionic strength, pH, or divalent cation, suggesting that distinct mechanisms may exist for these observed alterations in specificity. Recent evidence has indicated that BamHI and EcoRI share similar structural motifs, while Pvu II and EcoRV belong to a different structural family. Evidently, the use of hydration water to assist in site-specific recognition is a motif neither limited to nor defined by structural families.

Novel silica membrane material for molecular sieve applications

Development of new silica membranes properties, e.g., molecular sieving properties, has been increasingly gaining importance in the last few years. A novel unsupported silica membrane, referred to as hydrophobic metal-doped silica, was developed by cobalt-doping within the organic templated silica matrix. The novel material was prepared by the acid-catalyzed hydrolysis and condensation process of tetraethylorthosilicate (TEOS) and methyltriethoxysilane (MTES), which is the precursor for methyl ligand covalently bounded to the silica matrix. The synthesis and surface properties of the novel unsupported silica membrane as well as the unsupported blank silica and modified silica membranes were revealed by surface and microstructural techniques, such as water contact angle measurement, FTIR, X-ray, Solid-state 29Si MAS NMR, TGA and N2 and CO2 adsorption measurements. The results showed that the thermal stability of the organic templated silica matrix was enhanced by cobalt-doping process. A hydrophobic microporous silica membrane material with high thermal stability up to ∼560 °C in oxidizing atmosphere and a narrow pore size distribution centered at 1.1 nm was obtained. Therefore, a novel precursor material for molecular sieve silica membranes applications has been achieved and developed.

Water adsorption in hydrophilic zeolites: experiment and simulation

We have measured experimental adsorption isotherms of water in zeolite LTA4A, and studied the regeneration process by performing subsequent adsorption cycles after degassing at different temperatures. We observed incomplete desorption at low temperatures, and cation rearrangement at successive adsorption cycles. We also developed a new molecular simulation force field able to reproduce experimental adsorption isotherms in the range of temperatures between 273 K and 374 K. Small deviations observed at high pressures are attributed to the change in the water dipole moment at high loadings. The force field correctly describes the preferential adsorption sites of water at different pressures. We tested the influence of the zeolite structure, framework flexibility, and cation mobility when considering adsorption and diffusion of water. Finally, we performed checks on force field transferability between different hydrophilic zeolite types, concluding that classical, non-polarizable water force fields are not transferable.

Synthesis of denim waste-based adsorbents and their application in water defluoridation

This study evaluates the application of denim fiber scraps as a precursor for the synthesis of adsorbents for water treatment via pyrolysis and their application in water defluoridation. The best pyrolysis conditions for the synthesis of this novel adsorbent have been identified and a metal doping route with different salts of Al3 +, La3 + and Fe3 + was proposed to improve its fluoride adsorption behavior. Different spectroscopic and microscopic techniques (i.e., FTIR, XPS, XRF, SEM) were used to characterize the precursor and adsorbents, and to analyze the surface interactions involved in the fluoride removal mechanism. Experimental results showed that these adsorbents were effective for fluoride adsorption showing uptakes up to 4.25 mg/g. The Si-O–metal–F interactions appear to be highly relevant for the fluoride removal. This study highlights the potential of denim textile waste as a raw material for the production of added-value products, thus minimizing their associated disposal cost. It also shows the performance of denim textile waste as a precursor of adsorbents for addressing relevant environmental concerns such as fluoride pollution.