Pressure affects the structure and the dynamics of the D-galactose/D-glucose-binding protein from Escherichia coli by perturbing the C-terminal domain of the protein

The effect of the pressure on the structure and stability of the D-galactose/D-glucose binding protein from Escherichia coli in the absence (GGBP) and in the presence (GGBP/Glc) of glucose was studied by Fourier transform infrared (FT-IR) spectroscopy and molecular dynamic (MD) simulations. FT-IR spectroscopy experiments showed that the protein beta-structures are more resistant than alpha-helices structures to pressure value increases. In addition, the infrared data indicated that the binding of glucose stabilizes the protein structure against high pressure values, and the protein structure does not completely unfold up to pressure values close to 9000 bar. MD simulations allow a prediction of the most probable configuration of the protein, consistent with the increasing pressures on the two systems. The detailed analysis of the structures at molecular level confirms that, among secondary structures, alpha-helices are more sensitive than beta-structures to the destabilizing effect of high pressure and that glucose is able to preserve the structure of the protein in the complex. Moreover, the evidence of the different resistance of the two domains of this protein to high pressure is investigated and explained at a molecular level, indicating the importance of aromatic amino acid in protein stabilization.     

Folding propensity and biological activity of peptides: New insights from conformational properties of a novel peptide derived from Vitreoscilla haemoglobin

The synthetic peptide Vitr-p-13 (YPIVGQELLGAIK-NH(2)), derived from the bacterial dimeric Vitreoscilla haemoglobin (VHb) in the position 95-107, is characterized by a pre-eminent "statistical coil" conformation in water as demonstrated by CD experiments and long time-scale MD simulations. In particular, Vitr-p-13 does not spontaneously adopt an alpha-helix folding in water, but it is rather preferentially found in beta-hairpin-like conformations. Long time-scale MD simulations have also shown that Vitr-p-13 displays a "topological-trigger" which initiates alpha-helix folding within residues 7-10, exactly like seen in the temporins, a group of linear, membrane-active antimicrobial peptides of similar length. At variance with temporins, in Vitr-p-13 such a process is energetically very demanding (+10 kJ/mol) in water at 300 K, and the peptide was found to be unable to bind model membranes in vitro and was devoid of antimicrobial activity. The present results, compared with previous studies on similar systems, strengthen the hypothesis of the requirement of a partial folding when still in aqueous environment to allow a peptide to interact with cell-membranes and eventually exert membrane perturbation-related antibiotic effects on target microbial cells.     

Conformational equilibria of valine studied by dynamics simulation.

The conformational probability distribution of a valine residue in the valine dipeptide and of the valine side chain in an alpha-helix, as well as the change in helix stability for replacing alanine with valine, has been calculated by molecular dynamics simulations of explicitly hydrated systems: dipeptide, tetrapeptide and 10-, 14- and 18-residue oligoalanine helices. All computed free-energy differences are means from at least eight separate slow-growth simulations, four in each direction and are reported with their root-mean-square deviations. Different values for the change in free energy of folding (delta delta G degrees) have been calculated with the use of forcefields having an all-atom and a central-atom representation of methyl groups, etc. The value obtained with the all-atom forcefield agrees well with new experimental values (3 kJ/mol = 0.7 kcal/mol). Furthermore, the most stable valine side-chain rotamer in the helix is different for these two representations. The most stable rotamer for the all atom conformation is the same one that predominates for valines in alpha-helices in proteins of known conformation. The lower conformational freedom of the valine side chain in the helix contributes 1 kJ/mol to the difference in stability computed with the all-atom potential; unfavorable interactions of the side chain with helix, even in the most stable conformation, further increase delta delta G degrees.     

Water penetration and escape in proteins

The kinetics of water penetration and escape in cytochrome c (cyt c) is studied by molecular dynamics (MD) simulations at various temperatures. Water molecules that penetrate the protein interior during the course of an MD simulation are identified by monitoring the number of water molecules in the first coordination shell (within 3.5 A) of each water molecule in the system. Water molecules in the interior of cyt c have 0-3 water molecules in their first hydration shell and this coordination number persists for extended periods of time. At T = 300 K we identify over 200 events in which water molecules penetrate the protein and reside inside for at least 5 picoseconds (ps) within a 1.5 nanoseconds (ns) time period. Twenty-seven (27) water molecules reside for at least 300 ps, 17 water molecules reside in the protein interior for times longer than 500 ps, and two interior water molecules do not escape; at T = 360 K one water molecule does not escape; at 430 K all water molecules exchange. Some of the internal water molecules show mean square displacements (MSD) of 1 A2 characteristic of structural waters. Others show MSD as large as 12 A2, suggesting that some of these water molecules occupy transient cavities and diffuse extensively within the protein. Motions of protein-bound water molecules are rotationally hindred, but show large librations. Analysis of the kinetics of water escape in terms of a survival time correlation function shows a power law behavior in time that can be interpreted in terms of a broad distribution of energy barriers, relative to kappa BT, for water exchange. At T = 300 K estimates of the roughness of the activation energy distribution is 4-10 kJ/mol (2-4 kappa BT). Activation enthalpies for water escape are 6-23 kJ/mol. The difference in activation entropies between fast exchanging (0.01 ns) and slow exchanging (0.1-1 ns) water molecules is -27 J/K/mol. Dunitz (Science 1997;264:670.) has estimated the maximum entropy loss of a water molecule due to binding to be 28 J/K/mol. Therefore, our results suggest that the entropy of interior water molecules is similar to entropy of bulk water.     

Interpretation of H-NMR Experiments on the Orientation of the Transmembrane Helix WALP23 by Computer Simulations

Orientation, dynamics, and packing of transmembrane helical peptides are important determinants of membrane protein structure, dynamics, and function. Because it is difficult to investigate these aspects by studying real membrane proteins, model transmembrane helical peptides are widely used. NMR experiments provide information on both orientation and dynamics of peptides, but they require that motional models be interpreted. Different motional models yield different interpretations of quadrupolar splittings (QS) in terms of helix orientation and dynamics. Here, we use coarse-grained (CG) molecular dynamics (MD) simulations to investigate the behavior of a well-known model transmembrane peptide, WALP23, under different hydrophobic matching/mismatching conditions. We compare experimental ²H-NMR QS (directly measured in experiments), as well as helix tilt angle and azimuthal rotation (not directly measured), with CG MD simulation results. For QS, the agreement is significantly better than previously obtained with atomistic simulations, indicating that equilibrium sampling is more important than atomistic details for reproducing experimental QS. Calculations of helix orientation confirm that the interpretation of QS depends on the motional model used. Our simulations suggest that WALP23 can form dimers, which are more stable in an antiparallel arrangement. The origin of the preference for the antiparallel orientation lies not only in electrostatic interactions but also in better surface complementarity. In most cases, a mixture of monomers and antiparallel dimers provides better agreement with NMR data compared to the monomer and the parallel dimer. CG MD simulations allow predictions of helix orientation and dynamics and interpretation of QS data without requiring any assumption about the motional model.     

Hydrophobic Interactions Are a Key to MDM2 Inhibition by Polyphenols as Revealed by Molecular Dynamics Simulations and MM/PBSA Free Energy Calculations

p53, a tumor suppressor protein, has been proven to regulate the cell cycle, apoptosis, and DNA repair to prevent malignant transformation. MDM2 regulates activity of p53 and inhibits its binding to DNA. In the present study, we elucidated the MDM2 inhibition potential of polyphenols (Apigenin, Fisetin, Galangin and Luteolin) by MD simulation and MM/PBSA free energy calculations. All polyphenols bind to hydrophobic groove of MDM2 and the binding was found to be stable throughout MD simulation. Luteolin showed the highest negative binding free energy value of -173.80 kJ/mol followed by Fisetin with value of -172.25 kJ/mol. It was found by free energy calculations, that hydrophobic interactions (vdW energy) have major contribution in binding free energy.     

The human prion protein alpha2 helix: a thermodynamic study of its conformational preferences

We have synthesized both free and terminally-blocked peptide corresponding to the second helical region of the globular domain of normal human prion protein, which has recently gained the attention of structural biologists because of a possible role in the nucleation process and fibrillization of prion protein. The profile of the circular dichroism spectrum of the free peptide was that typical of alpha-helix, but was converted to that of beta-structure in about 16 h. Instead, below 2.1 x 10(-5) M, the spectrum of the blocked peptide exhibited a single band centered at 200 nm, unequivocally associated to random conformations, which did not evolve even after 24 h. Conformational preferences of this last peptide have been investigated as a function of temperature, using trifluoroethanol or low-concentration sodium dodecyl sulfate as alpha- or beta-structure inducers, respectively. Extrapolation of free energy data to zero concentration of structuring agent highlighted that the peptide prefers alpha-helical to beta-type organization, in spite of results from prediction algorithms. However, the free energy difference between the two forms, as obtained by a thermodynamic cycle, is subtle (roughly 5-8 kJ mol(-1) at any temperature from 280 K to 350 K), suggesting conformational ambivalence. This result supports the view that, in the prion protein, the structural behavior of the peptide is governed by the cellular microenvironment.     

Comparative thermodynamic analysis of DNA--protein interactions using surface plasmon resonance and fluorescence correlation spectroscopy

We report a kinetic and thermodynamic analysis of interactions between ssDNA and replication protein A (RPA) using surface plasmon resonance (SPR) and fluorescence correlation spectroscopy (FCS) at variable temperature. The two methods yield different values for the Gibbs free energy but nearly the same value for the reaction enthalpy of ssDNA-RPA complex formation. The Gibbs free energy was determined by SPR and FCS to be -62.6 and -54.7 kJ/mol, respectively. The values for the reaction enthalpy are -64.4 and -66.5 kJ/mol. It is concluded that the difference in Gibbs free energy measured by the two methods is due to different reaction entropies. The entropic contribution to the free energy at 25 degrees C is -1.8 kJ/mol for SPR and -11.8 kJ/mol for FCS. In SPR, the reaction is restricted to two dimensions because of immobilization of the DNA molecules to the sensor surface. In contrast, FCS is able to follow complex formation without spatial restrictions. In consequence, the reaction entropy determined from SPR experiments is lower than for FCS experiments.     

Salt bridges in the hyperthermophilic protein Ssh10b are resilient to temperature increases

A double mutant cycle (DMC) approach was employed to estimate the effect of temperature on the contribution of two highly conserved salt bridges to protein stability in the hyperthermophilic protein Ssh10b. The coupling free energy were 2.4 +/- 0.4 kJ/mol at 298 K and 2.2 +/- 0.4 kJ/mol at 353 K for Glu-54/Arg-57, and 6.0 +/- 0.2 kJ/mol at 298 K and 5.9 +/- 0.6 kJ/mol at 353 K for Glu-36/Lys-68. The stability free energy of Ssh10b decrease greatly with increasing temperature, while the direct contribution of these two salt bridges to protein stability remain almost constant, providing evidence supporting the theoretical prediction that salt bridges are extremely resilient to temperature increases and thus are specially suited to improving protein stability at high temperatures. The reason for the difference in coupling free energy between salt bridges Glu-54/Arg-57 and Glu-36/Lys-68 is discussed. Comparing our results with published DMC data for the contribution of salt bridges to stability in other proteins, we found that the energy contribution of a salt bridge formed by two charged residues far apart in the primary sequence is higher than that of those formed between two very close ones. Implications of this finding are useful for engineering proteins with enhanced thermostability.     

Interaction of lipid-free apolipoprotein A-I with cholesterol revealed by molecular modeling

We report the modeling of the interaction of differently self-associated lipid-free apoA-I with cholesterol monomer and tail-to-tail (TT) or face-to-face (FF) cholesterol dimer. Cholesterol dimerization is exploited to reconcile the existing experimental data on cholesterol binding to apoA-I with extremely low critical micelle concentration of cholesterol. Two crystal structures of 1–43 N-truncated apolipoprotein Δ(1-43)A-I tetramer (PDB ID: 1AV1, structure B), 185–243 C-truncated apolipoprotein Δ(185-243)A-I dimer (PDB ID: 3R2P, structure M) were analyzed. Cholesterol monomers bind to multiple binding sites in apoA-I monomer, dimer and tetramer with low, moderate and high energy (?10 to ?28 kJ/mol with Schrödinger package), still insufficient to overcome the thermodynamic restriction by cholesterol micellization (?52.8 kJ/mol). The binding sites partially coincide with the putative cholesterol-binding motifs. However, apoA-I monomer and dimer existing in structure B, that contain nonoverlapping and non-interacting pairs of binding sites with high affinity for TT and FF cholesterol dimers, can bind in common 14 cholesterol molecules that correspond to existing values. ApoA-I monomer and dimer in structure M can bind in common 6 cholesterol molecules. The values of respective total energy of cholesterol binding up to 64.5 and 67.0 kJ/mol for both B and M structures exceed the free energy of cholesterol micellization. We hypothesize that cholesterol dimers may simultaneously interact with extracellular monomer and dimer of lipid-free apoA-I, that accumulate at acid pH in atheroma. The thermodynamically allowed apolipoprotein-cholesterol interaction outside the macrophage may represent a new mechanism of cholesterol transport by apoA-I from atheroma, in addition to ABCA1-mediated cholesterol efflux.     

Relative binding free energy calculations of inhibitors to two mutants (Glu46----Ala/Gln) of ribonuclease T1 using molecular dynamics/free energy perturbation approaches.

We present free energy perturbation calculations on the complexes of Glu46----Ala46 (E46A) and Glu46----Gln46 (E46Q) mutants of ribonuclease T1 (RNaseT1) with inhibitors 2'-guanosine monophosphate (GMP) and 2'-adenosine monophosphate (AMP) by a thermodynamic perturbation method implemented with molecular dynamics (MD). Using the available crystal structure of the RNaseT1-GMP complex, the structures of E46A-GMP and E46Q-GMP were model built and equilibrated with MD simulations. The structures of E46A-AMP and E46Q-AMP were obtained as a final structure of the GMP----AMP perturbation calculation respectively. The calculated difference in the free energy of binding (delta delta Gbind) was 0.31 kcal/mol for the E46A system and -1.04 kcal/mol for the E46Q system. The resultant free energies are much smaller than the experimental and calculated value of approximately 3 kcal/mol for the native RNaseT1, which suggests that both mutants have greater relative adenine affinities than native RNaseT1. Especially E46Q is calculated to have a larger affinity for adenine than guanine, as we suggested previously from the calculation on the native RNaseT1. Thus, the molecular dynamics/free energy perturbation method may be helpful in protein engineering, directed toward increasing or changing the substrate specificity of enzymes.     

Anisotropic elastic network modeling of entire microtubules

Microtubules are supramolecular structures that make up the cytoskeleton and strongly affect the mechanical properties of the cell. Within the cytoskeleton filaments, the microtubule (MT) exhibits by far the highest bending stiffness. Bending stiffness depends on the mechanical properties and intermolecular interactions of the tubulin dimers (the MT building blocks). Computational molecular modeling has the potential for obtaining quantitative insights into this area. However, to our knowledge, standard molecular modeling techniques, such as molecular dynamics (MD) and normal mode analysis (NMA), are not yet able to simulate large molecular structures like the MTs; in fact, their possibilities are normally limited to much smaller protein complexes. In this work, we developed a multiscale approach by merging the modeling contribution from MD and NMA. In particular, MD simulations were used to refine the molecular conformation and arrangement of the tubulin dimers inside the MT lattice. Subsequently, NMA was used to investigate the vibrational properties of MTs modeled as an elastic network. The coarse-grain model here developed can describe systems of hundreds of interacting tubulin monomers (corresponding to up to 1,000,000 atoms). In particular, we were able to simulate coarse-grain models of entire MTs, with lengths up to 350 nm. A quantitative mechanical investigation was performed; from the bending and stretching modes, we estimated MT macroscopic properties such as bending stiffness, Young modulus, and persistence length, thus allowing a direct comparison with experimental data.     

Structure,stability and water permeation of ([D-Leu-L-Lys-(D-Gln-L-Ala)3]) cyclic peptide nanotube: a molecular dynamics study

The structural stability of 8 × ([D-Leu-L-Lys-(D-Gln-L-Ala)₃]) cyclic peptide nanotube (CPN) in water and different phospholipid bilayers were explored by 100 ns independent molecular dynamics (MD) simulations. The role of non-bonded interaction energy between the side and main chains of cyclic peptide rings in different membrane environments assessed, wherein the repulsive electrostatic interaction energy between neighbouring cyclic peptide rings was found adequate to break hydrogen bond energy thereby to crumple CPN. Further, the water permeation across the CPN channel was studied in four types of phospholipid bilayers- DMPG (1,2-Dimyristoyl-sn-glycero-3-phosphorylglycerol), DMPS (1,2-Dimyristoyl-sn-glycero-3-phosphoserine), POPC (1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) and POPE (1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine) from MD simulations. DMPS membrane shows higher non-bonded interaction energies (?1913.06 kJ/mol of electrostatic interaction energy and ?994.13 kJ/mol of van der Waals interaction energy) with CPN due to the presence of polar molecules in lipid structure. Thusly, the non-bonded interaction energies were essential towards the stability of CPN than hydrogen bonds between the nearby cyclic peptides. The result also reveals the role of side chains, hydrogen bonds and non-bonded interaction energies in an aqueous environment. The diffusion coefficient of water obtained from means square deviation calculation shows similar coefficients irrespective of the lipid surroundings. However, the permeation coefficients demonstrate water flow in the channel relies upon the environment.     

Antifreeze proteins at the ice/water interface: three calculated discriminating properties for orientation of type I proteins

Antifreeze proteins (AFPs) protect many plants and organisms from freezing in low temperatures. Of the different AFPs, the most studied AFP Type I from winter flounder is used in the current computational studies to gain molecular insight into its adsorption at the ice/water interface. Employing molecular dynamics simulations, we calculate the free energy difference between the hydrophilic and hydrophobic faces of the protein interacting with ice. Furthermore, we identify three properties of Type I "antifreeze" proteins that discriminate among these two orientations of the protein at the ice/water interface. The three properties are: the "surface area" of the protein; a measure of the interaction of the protein with neighboring water molecules as determined by the number of hydrogen bond count, for example; and the side-chain orientation angles of the threonine residues. All three discriminants are consistent with our free energy results, which clearly show that the hydrophilic protein face orientations toward the ice/water interface, as hypothesized from experimental and ice/vacuum simulations, are incorrect and support the hypothesis that the hydrophobic face is oriented toward the ice/water interface. The adsorption free energy is calculated to be 2-3 kJ/mol.     

Free energy surface for rotamers of cis-enol malonaldehyde in aqueous solution studied by molecular dynamics calculations

Two-dimensional free energy surfaces for four rotamers of cis-enol malonaldehyde in water have been investigated by umbrella sampling molecular dynamics (MD) calculations. Biasing potential used in the umbrella sampling calculation was adopted to be the minus of conformational free energy preliminary obtained by the thermodynamic integration MD calculations for the rigid malonaldehyde whose stretching and bending were all fixed. The calculated free energy surface shows that, in water, a rotamer that has an intramolecular hydrogen bond is most stable among the rotamers. This is the same as that in vacuum, while order of relative stability of the other three rotamers is different in water and in vacuum. Inclusion of intramolecular vibrations changed the free energy surface little, i.e. at most 2.6 kJ/mol, which is much smaller than the solvation free energy. Free energy barriers from the most stable intramolecular hydrogen bonded rotamer to the others are lowered by hydration but they are still very high, >50 kJ/mol, such that the malonaldehyde molecule spends most of its time in water taking this conformation. Thus, reaction coordinate for intramolecular proton transfer reaction in water may be constructed assuming this rotamer.     

Molecular dynamics and atomic charge calculations in the study of heparin conformation in aqueous solution

HF/6-31G** and molecular dynamics (MD) simulations were used to evaluate the performance of different atomic charge basis sets (i.e., Mulliken, Lowdin, and Electrostatic Potential Derived Charges--ESP) in heparin simulations. HF/3-21 G calculations were also used to study the NMR conformation of the IdoA residue. The results thus obtained indicated that ESP and Lowdin charges gave the better results in heparin simulations, followed by Mulliken charges, and that the minimum-energy conformation of IdoA can be different from that observed by NMR spectroscopy by less than 1 Angstrom. However, it was found that this small conformational modification is capable of inducing a change of almost 200 kJ/mol in the interactions of heparin with the surrounding environment, which is a meaningful amount of energy in the context of ligand-receptor interactions. This information can be potentially of great relevance in the design of heparin-derived antithrombotic compounds.     

Conformation and stability of the Streptococcus pyogenes pSM19035-encoded site-specific beta recombinase, and identification of a folding intermediate

Solution properties of beta recombinase were studied by circular dichroism and fluorescence spectroscopy, size exclusion chromatography, analytical ultracentrifugation, denaturant-induced unfolding and thermal unfolding experiments. In high ionic strength buffer (1 M NaCl) beta recombinase forms mainly dimers, and strongly tends to aggregate at ionic strength lower than 0.3 M NaCl. Urea and guanidinium chloride denaturants unfold beta recombinase in a two-step process. The unfolding curves have bends at approximately 5 M and 2.2 M in urea and guanidinium chloride-containing buffers. Assuming a three-state unfolding model (N2-->2I-->2U), the total free energy change from 1 mol of native dimers to 2 mol of unfolded monomers amounts to deltaG(tot) = 17.9 kcal/mol, with deltaG(N2-->2I) = 4.2 kcal/mol for the first transition and deltaG(I-->U) = 6.9 kcal/mol for the second transition. Using sedimentation-equilibrium analytical ultracentrifugation, the presence of beta recombinase monomers was indicated at 5 M urea, and the urea dependence of the circular dichroism at 222 nm strongly suggests that folded monomers represent the unfolding intermediate.     

Insight into ribonuclease A domain swapping by molecular dynamics unfolding simulations

Bovine pancreatic ribonuclease (RNase A) deserves a special place among the numerous proteins that form oligomers by three-dimensional domain swapping. In fact, under destabilizing conditions and at high protein concentrations, it can swap two different domains, the N-terminal alpha-helix or the C-terminal beta-strand, leading to dimers with different quaternary structures. With the change in the unfolding conditions, the relative abundance of the two dimers varies, and the prevalence of one dimer over the other is inverted. To investigate the dynamic behavior of the termini, four independent 10 ns high-temperature molecular dynamics simulations of RNase A were carried out at two different pH values in an attempt to reproduce the experimental conditions of neutral and very low pH that favor the formation of the N- and C-terminal domain-swapped dimers, respectively. In agreement with experimental data, under mild unfolding conditions, a partial or complete opening of the N-terminal arm is observed, whereas the dislocation of the C-terminus away from the core of the structure occurs only during the low-pH simulations. Furthermore, the picture emerging from this study indicates that the same protein can have different pathways for domain swapping. Indeed, in RNase A the C-terminal swapping requires a substantial unfolding of the monomers, whereas the N-terminal swapping can occur through only partial unfolding.     

Substituting leucine for alanine-86 in the tether region of the iron-sulfur protein of the cytochrome bc1 complex affects the mobility of the [2Fe2S] domain

Mutating three conserved alanine residues in the tether region of the iron-sulfur protein of the yeast cytochrome bc(1) complex resulted in 22-56% decreases in enzymatic activity [Obungu et al. (2000) Biochim. Biophys. Acta 1457, 36-44]. The activity of the cytochrome bc(1) complex isolated from A86L was decreased 60% compared to the wild-type without loss of heme or protein and without changes in the 2Fe2S cluster or proton-pumping ability. The activity of the bc(1) complex from mutant A92R was identical to the wild-type, while loss of both heme and activity was observed in the bc(1) complex isolated from mutant A90I. Computer simulations indicated that neither mutation A86L nor mutation A92R affects the alpha-helical backbone in the tether region; however, the side chain of the leucine substituted for Ala-86 interacts with the side chain of Leu-89. The Arrhenius plot for mutant A86L was apparently biphasic with a transition observed at 17-19 degrees C and an activation energy of 279.9 kJ/mol below 17 degrees C and 125.1 kJ/mol above 17 degrees C. The initial rate of cytochrome c(1) reduction was lowered 33% in mutant A86L; however, the initial rate of cytochrome b reduction was unaffected, suggesting that movement of the tether region of the iron-sulfur protein is necessary for maximum rates of enzymatic activity. Substituting a leucine for Ala-86 impedes the unwinding of the alpha-helix and hence movement of the tether.     

Behavior of Aromatic Molecules in Silicalite by the Direct Integration of the Configurational Integral

Abstract

Adsorption data of aromatic molecules adsorbed in silicalite show highly unusual characteristics which were attributed to structural effects caused by the comparable size of molecules and pores. In this study, the interaction of aromatic compounds with silicalite are examined on the molecular level. The interactions are calculated by atom-atom approximation using Lennard-Jones potentials. The constants are calculated, without fitting, from Kirkwood-Muller formulas. Benzene and p-xylene are represented as a rigid structure of 12 and 18 atom centers. The model is anisotropic.

The diffusional behavior of molecules is examined by minimizing the potential energy in the channels which requires less computational time than Molecular Dynamics. The activation energy for the diffusion of benzene, 27.6 kJ/mol, is in excellent agreement with data, 28.8 kJ/mol. The results indicate that both molecules can enter the smaller zig-zag channels. The energetically most favorable location in the main channels is the mid-point between intersections. All rotations are restricted in the channels but the molecules can rotate in any direction (with some movement of the center) at intersections.

The Henry's law constant and internal energy of adsorption at zero coverage are calculated by direct integration of the configurational integral. Direct integration is more efficient than Monte Carlo and Molecular Dynamics simulations since the molecules are highly restricted in the pores. The predicted internal energy of adsorption, ? 54.86 and ? 75.30 kJ/mol for benzene and p-xylene is in good agreement with data of ? 50.92 and ? 62.15 kJ/mol respectively. There is appreciable difference between the predicted and experimental Henry's law constants. The agreement can be improved by fitting the Lennard-Jones constants which has not been attempted.

Although the calculations are performed at infinite dilution and entropy effects are not included, the results bring insight to the behavior of molecules in highly restricted environments such as in tight pores. Similar simplified calculations can be used to close the gap between highly idealized molecular simulations and complicated systems common in real applications.