Monte Carlo Study of the Buckingham Exponential-six Fluid

Abstract

In this paper we report the results of extensive Monte Carlo simulations of a pure fluid of Buckingham modified exponential-six molecules. Data are presented for the configurational energy and pressure covering a wide range of temperatures and densities. These data are interpreted using the generalized van der Waals partition function with a novel separation into free volume and mean potential terms. We find, surprisingly, that the Buckingham fluid is described by a simple van der Waals-like equation of state provided that the b parameter is temperature dependent and chosen in a theoretically correct manner.     

Filling small, empty protein cavities: structural and energetic consequences

Most proteins contain small cavities that can be filled by replacing cavity-lining residues by larger ones. Since shortening mutations in hydrophobic cores tend to destabilize proteins, it is expected that cavity-filling mutations may conversely increase protein stability. We have filled three small cavities in apoflavodoxin and determined by NMR and equilibrium unfolding analysis their impact in protein structure and stability. The smallest cavity (14 A3) has been filled, at two different positions, with a variety of residues and, in all cases, the mutant proteins are locally unfolded, their structure and energetics resembling those of an equilibrium intermediate of the thermal unfolding of the wild-type protein. In contrast, two slightly larger cavities of 20 A3 and 21 A3 have been filled with Val to Ile or Val to Leu mutations and the mutants preserve both the native fold and the equilibrium unfolding mechanism. From the known relationship, observed in shortening mutations, between stability changes and the differential hydrophobicity of the exchanged residues and the volume of the cavities, the filling of these apoflavodoxin cavities is expected to stabilize the protein by approximately 1.5 kcal mol(-1). However, both urea and thermal denaturation analysis reveal much more modest stabilizations, ranging from 0.0 kcal mol(-1) to 0.6 kcal mol(-1), which reflects that the accommodation of single extra methyl groups in small cavities requires some rearrangement, necessarily destabilizing, that lowers the expected theoretical stabilization. As the size of these cavities is representative of that of the typical small, empty cavities found in most proteins, it seems unlikely that filling this type of cavities will give rise to large stabilizations.     

Side-chain contributions to membrane protein structure and stability

The molecular forces that stabilize membrane protein structure are poorly understood. To investigate these forces we introduced alanine substitutions at 24 positions in the B helix of bacteriorhodopsin and examined their effects on structure and stability. Although most of the results can be rationalized in terms of the folded structure, there are a number of surprises. (1) We find a remarkably high frequency of stabilizing mutations (17%), indicating that membrane proteins are not highly optimized for stability. (2) Helix B is kinked, with the kink centered around Pro50. The P50A mutation has no effect on stability, however, and a crystal structure reveals that the helix remains bent, indicating that tertiary contacts dominate in the distortion of this helix. (3) We find that the protein is stabilized by about 1kcal/mol for every 38A(2) of surface area buried, which is quite similar to soluble proteins in spite of their dramatically different environments. (4) We find little energetic difference, on average, in the burial of apolar surface or polar surface area, implying that van der Waals packing is the dominant force that drives membrane protein folding.     

Measurements of conformational changes during adhesion of lipid and protein (polylysine and S-layer) surfaces

The adhesion forces between various surfaces were measured using the "surface forces apparatus" technique. This technique allows for the thickness of surface layers and the adhesion force between them to be directly measured in controlled vapor or liquid environments. Three types of biological surfaces were prepared by depositing various lipid-protein monolayers (with thicknesses ranging from 1 to 4 nm) on the inert, molecularly smooth mica surface: (i) hydrophobic lipid monolayers; (ii) amphiphilic polyelectrolyte surfaces of adsorbed polylysine; and (iii) deposited bacterial S-layer proteins. The adhesion, swelling, and wetting properties of these surfaces was measured as a function of relative humidity and time. Initial adhesion is due mainly to the van der Waals forces arising from nonpolar (hydrophobic) contacts. Following adhesive contact, significant molecular rearrangements can occur which alter their hydrophobic-hydrophilic balance and increase their adhesion with time. Increased adhesion is generally enhanced by (i) increased relative humidity (or degree of hydration); (ii) increased contact time; and (iii) increased rates of separation. The results are likely to be applicable to the adhesion of many other biosurfaces, and show that the hydrophobicity of a lipid or protein surface is not an intrinsic property of that surface but depends on its environment (e.g., on whether it is in aqueous solution or exposed to the atmosphere), and on the relative humidity of the atmosphere. It also depends on whether the surface is in adhesive contact with another surface and-when considering dynamic (nonequilibrium) conditions-on the time and previous history of its interaction with that surface. (c) 1993 John Wiley & Sons, Inc.     

Sampling the Initial Conditions for Quasiclassical Trajectory Studies of Vibrational Predissociation Dynamics

A method is proposed for selecting initial conditions in order to study the classical dynamics of vibrational predissociation of van der Waals (vdW) clusters. The method starts from the quantum initial state of the system, which is used to sample and weight the initial positions of the different modes. The initial values of the associated momenta are calculated in such a way that they correspond with the total energy and angular momentum of the system initial state, at which the classical trajectory simulation is to be carried out. An application to the case of Cl2–-He2 is presented and discussed.     

ProteinVolume: calculating molecular van der Waals and void volumes in proteins

Background

Voids and cavities in the native protein structure determine the pressure unfolding of proteins. In addition, the volume changes due to the interaction of newly exposed atoms with solvent upon protein unfolding also contribute to the pressure unfolding of proteins. Quantitative understanding of these effects is important for predicting and designing proteins with predefined response to changes in hydrostatic pressure using computational approaches. The molecular surface volume is a useful metric that describes contribution of geometrical volume, which includes van der Waals volume and volume of the voids, to the total volume of a protein in solution, thus isolating the effects of hydration for separate calculations.

Results

We developed ProteinVolume, a highly robust and easy-to-use tool to compute geometric volumes of proteins. ProteinVolume generates the molecular surface of a protein and uses an innovative flood-fill algorithm to calculate the individual components of the molecular surface volume, van der Waals and intramolecular void volumes. ProteinVolume is user friendly and is available as a web-server or a platform-independent command-line version.

Conclusions

ProteinVolume is a highly accurate and fast application to interrogate geometric volumes of proteins. ProteinVolume is a free web server available on http://gmlab.bio.rpi.edu. Free-standing platform-independent Java-based ProteinVolume executable is also freely available at this web site.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0531-2) contains supplementary material, which is available to authorized users.     

Methyl 6-deoxy-6-iodo-α-d-glucopyranoside: crystal structure and high-resolution NMR spectroscopy

Single-crystal X-ray diffraction and high-resolution ¹H and ¹³C NMR spectral data for methyl 6-deoxy-6-iodo-α-d-glucopyranoside are reported. The ⁴C₁ conformation was found to be the preferred form for this compound, both in the crystal lattice and in solution. The rotational preferences of all the groups bound to the pyranose ring are presented. The stabilization of the crystal structure by a network of O-H···O intra- and intermolecular interactions as well as the short contacts of the iodine atoms is discussed.     

On the existence of MHn species with M=Al, Ga and n=4, 5, 6. Computational study of structures, stabilities and bonding

Based on second-order perturbation theory (MP2) predictions with large 6-311++G(3df, 3pd) basis set we have reviewed the possible structures and stabilities of a series of neutral MH_n(M=Al, Ga; n=4, 5, 6) species. For AlH₄ and AlH₅, our results confirm the previous theoretical findings, which indicate the dihydrogen C_s complexes (²A′) AlH₂(H₂) and (¹A′) AlH₃(H₂), respectively, as the lowest energy isomers. We found, similarly, C_s (²A′) GaH₂(H₂) and (¹A′) GaH₃(H₂) van der Waals complexes as the most stable species of the gallium analogues GaH₄ and GaH₅. The calculated H₂ dissociation energies (D_e) for AlH₂(H₂) and AlH₃(H₂) are of the order 1.8–2.5 kcalmol¹, whereas this range of values for GaH₂(H₂) and GaH₃(H₂) is 1.4–1.8 kcalmol¹ . Symmetry-adapted perturbation theory (SAPT) was used to analyze the interaction energies of these dihydrogen complexes (for n=5) to determine why the Ga species show a smaller binding energy than the Al species. The SAPT partitioning of the interaction energy showed significant differences between AlH₃(H₂) and GaH₃(H₂), resulting from the much stronger “hydride” character of the aluminum species. The experimental observation of AlH₂(H₂) and AlH₃(H₂), and likely GaH₃(H₂), via low-temperature matrix isolation has been reported recently by Pullumbi et al. and Andrews et al., supporting the theoretical predictions. For n=6, we found the degenerate C₂(²A) and C_s(²A′) MH₂(H₂)₂ “double H₂” type van der Waals complexes as the lowest energy species for both M=Al and Ga.Electronic Supplementary Material Supplementary material is available for this article at     

PRI-Modeler: extracting RNA structural elements from PDB files of protein-RNA complexes

A complete understanding of protein and RNA structures and their interactions is important for determining the binding sites in protein-RNA complexes. Computational approaches exist for identifying secondary structural elements in proteins from atomic coordinates. However, similar methods have not been developed for RNA, due in part to the very limited structural data so far available. We have developed a set of algorithms for extracting and visualizing secondary and tertiary structures of RNA and for analyzing protein-RNA complexes. These algorithms have been implemented in a web-based program called PRI-Modeler (protein-RNA interaction modeler). Given one or more protein data bank files of protein-RNA complexes, PRI-Modeler analyzes the conformation of the RNA, calculates the hydrogen bond (H bond) and van der Waals interactions between amino acids and nucleotides, extracts secondary and tertiary RNA structure elements, and identifies the patterns of interactions between the proteins and RNAs. This paper presents PRI-Modeler and its application to the hydrogen bond and van der Waals interactions in the most representative set of protein-RNA complexes. The analysis reveals several interesting interaction patterns at various levels. The information provided by PRI-Modeler should prove useful for determining the binding sites in protein-RNA complexes. PRI-Modeler is accessible at http://wilab.inha.ac.kr/primodeler/, and supplementary materials are available in the analysis results section at http://wilab.inha.ac.kr/primodeler/.