Repulsive interactions between polar and apolar atoms in globular proteins

A significant number of tetrahedral carbon-backbone nitrogen pairs in ferrocytochrome c have interatomic distances that are more than 0.4 Å shorter than the sum of the atomic van der Waals' radii. The non-bonded repulsions of these pairs destabilize the native structure by as much as 50 to 60 kcal/mol. A detailed examination of the close contacts suggests that these result from packing forces associated with the formation of secondary and tertiary structure in the protein.     

Refinement of the structure of carp muscle calcium-binding parvalbumin by model building and difference Fourier analysis.

The structure of carp muscle calcium-binding parvalbumin has been refined to an overall residual, ($\text{Σ¦F}{}_0\text{− F}{}_{\mathrm{c}}\text{¦}\text{Σ¦F}{}_0\text{¦}$), of 0.25 by a combination of model building and difference Fourier analyses. The atomic positions were allowed to vary up to 0.25 Å from their idealized positions; individual temperature factors ranged from 2 to 150; and 138 solvent peaks were included in this final structure factor calculation. The effects of varying these parameters were analyzed. The treatment of low-order reflections and the influence of solvent have been analyzed in terms of Babinet's principle. These procedures and results are applicable to other protein refinement problems. The errors in co-ordinates are estimated to range from 0.15 Å for the Ca²⁺ ions to 0.30 Å for internal side chain atoms.A van der Waals' radii study indicates that 42% of the crystal volume is occupied by solvent. There are 16 intermolecular contacts. Of the 108 main chain peptide hydrogen atoms, 15 are neither in contact with solvent nor involved in hydrogen bonds. These “lost” hydrogen bonds are considered to be important to the proposed model of function.     

Quantification of helix-helix binding affinities in micelles and lipid bilayers

A theoretical approach for estimating association free energies of alpha-helices in nonpolar media has been developed. The parameters of energy functions have been derived from DeltaDeltaG values of mutants in water-soluble proteins and partitioning of organic solutes between water and nonpolar solvents. The proposed approach was verified successfully against three sets of published data: (1) dissociation constants of alpha-helical oligomers formed by 27 hydrophobic peptides; (2) stabilities of 22 bacteriorhodopsin mutants, and (3) protein-ligand binding affinities in aqueous solution. It has been found that coalescence of helices is driven exclusively by van der Waals interactions and H-bonds, whereas the principal destabilizing contributions are represented by side-chain conformational entropy and transfer energy of atoms from a detergent or lipid to the protein interior. Electrostatic interactions of alpha-helices were relatively weak but important for reproducing the experimental data. Immobilization free energy, which originates from restricting rotational and translational rigid-body movements of molecules during their association, was found to be less than 1 kcal/mole. The energetics of amino acid substitutions in bacteriorhodopsin was complicated by specific binding of lipid and water molecules to cavities created in certain mutants.     

van der Waals forces influencing adhesion of cells

Adhesion molecules, often thought to be acting by a ‘lock and key’ mechanism, have been thought to control the adhesion of cells. While there is no doubt that a coating of adhesion molecules such as fibronectin on a surface affects cell adhesion, this paper aims to show that such surface contamination is only one factor in the equation. Starting from the baseline idea that van der Waals force is a ubiquitous attraction between all molecules, and thereby must contribute to cell adhesion, it is clear that effects from geometry, elasticity and surface molecules must all add on to the basic cell attractive force. These effects of geometry, elasticity and surface molecules are analysed. The adhesion force measured between macroscopic polymer spheres was found to be strongest when the surfaces were absolutely smooth and clean, with no projecting protruberances. Values of the measured surface energy were then about 35 mJ m⁻², as expected for van der Waals attractions between the non-polar molecules. Surface projections such as abrasion roughness or dust reduced the molecular adhesion substantially. Water cut the measured surface energy to 3.4 mJ m⁻². Surface active molecules lowered the adhesion still further to less than 0.3 mJ m⁻². These observations do not support the lock and key concept.     

Quantitative determination of the conformation of cyclic 3',5'-adenosine monophosphate in solution using lanthanide ions as nuclear magnetic resonance probes

The conformation of cyclic 3′,5′-adenosine monophosphate in deuterium oxide has been determined at pH 2.0 and pH 5.5, using lanthanide ions as paramagnetic nuclear magnetic resonance probes.The lanthanide ion-induced shifts in the nuclear magnetic resonance energy for a given nucleus are dependent on the geometric position of that nucleus relative to the bound lanthanide ion. As expected, these shifts are pseudocontact in origin and are consistent with axial symmetry. Analysis of the concentration dependence of the shift shows that the lanthanide ion is bound to the phosphate entity giving a 1:1 complex. Further, base stacking and other intermolecular interactions are negligible.To confirm the conformation, which is found from a computer search with the above shift data, we have measured the changes in relaxation times, T₁ and T₂, induced by binding of Gd³⁺. The geometric dependence of these relaxation effects is different from that of shifts, being dependent only on distance. The agreement of these data with the computer “shift” conformation is satisfactory.Some ³¹P nuclear magnetic resonance experiments were done to confirm the metal co-ordination position although, here, there are contact contributions to both shift and relaxation.The computer program finds the conformations that have the correct geometry to account for the shift data, by searching all possible conformations. Non-bond rotations were used as a method of changing the pucker of the phosphate and ribose rings, the position of the base being defined by a single bond rotation. The nuclear magnetic resonance data and minimum van der Waals' distances were used as “active filters” in the computer search.At both values of the pH we have found closely related families of solutions, with the pucker of the phosphate and ribose rings roughly similar to those in an approximate X-ray study of cyclic AMP. The orientation of the base varies with pH.     

Structural response to mutation at a protein-protein interface

We have crystallised three mutants of the barnase-barstar complex in which interactions across the interface have been deleted by simultaneous mutation of both residues involved in the interaction. Each mutant deletes a different type of interaction at the interface: the first complex bnHis102-->Ala-bsTyr29-->Phe (bn, barnase; bs, barstar), deletes a van der Waals packing interaction; the second complex, bnLys27-->Ala-bsThr42-->Ala, deletes a hydrogen bond; the third, bnLys27-->Ala-bsAsp35-->Ala, deletes a long-range charge-charge interaction. The contribution of each of these side-chains to the stability of the complex is known; the coupling energy between the deleted side-chains is also known. Despite each of the double mutants being significantly destabilised compared with the wild-type, the effects of mutation are local. Only small movements in the main-chain surrounding the sites of mutation and some larger movements of neighbouring side-chains are observed in the mutant complexes. The exact response to mutation is context-dependent and for the same mutant can vary depending upon the environment within the crystal. In some double mutant complexes, interfacial pockets, which are accessible to bulk solvent are formed, whereas interfacial cavities which are isolated from bulk solvent, are formed in others. In all double mutants, water molecules fill the created pockets and cavities. These water molecules mimic the deleted side-chains by occupying positions close to the non-carbon atoms of truncated side-chains and re-making many hydrogen bonds made by the truncated side-chains in the wild-type. It remains extremely difficult, however, to correlate energetic and structural responses to mutation because of unknown changes in entropy and entropy-enthalpy compensation.     

Molecular Dynamics Simulation of Model Lipid Membranes: Structural Effects of Impurities

Abstract

The gel to fluid phase transition or ordered to disordered phase transition observed in biological membranes are simulated by using constant energy Molecular Dynamics. The surface part of the membrane is modelled as a two-dimensional matrix formed by the head groups of the phospholipid molecules. Head molecules which are modelled as three spheres fused with three force centers, interact with each other via van der Waals and Coulomb type interactions. The -so called- impurity or foreign molecule embedded in the surface represents the protein type molecule which is present in biological membranes and control its activity. It is modelled as a pentagon having one force centers in each corner. It also interacts with the surface molecules again via van der Waals and Coulomb type interactions. The surface density is kept constant in the simulations of the systems with or without impurity. Structural and orientational changes due to impurity were observed and proved by monitoring two-dimensional order parameter. It has been shown that melting of the surface or breakage of the ordering of the surface molecules becomes easier and ordered to disordered phase transition temperature was lowered by 100 K if the impurity is present.     

Atomic resolution analysis of a 2:1 complex of CpG and acridine orange.

Cytidylyl-3', 5'-guanosine and acridine orange crystallize in a highly-ordered triclinic lattice which diffracts X-rays to 0.85 angstrom resolution. The crystal structure has been solved and refined to a residual factor of 9.5%. The two dinucleoside phosphate molecules form an antiparallel double helix with the acridine orange intercalated between them. The two base pairs of the double helical fragment have a twist angle of 10 degrees and it is found to have a C3' endo-(3', 5')-C2' endo mixed sugar puckering along the nucleotide backbone as has been observed for other simple intercalator complexes. Twenty-five water molecules have been located in the lattice together with a sodium ion. The intercalator double helical fragments form sheets which are held together by van der Waals interactions in one direction and hydrogen bonding interactions in the other. The crystal lattice contains aqueous channels in which sixteen water molecules are hydrogen bonded to the nucleotide, none to the intercalator, five water molecules are coordinated about the sodium ion and four water molecules bind solely to other water molecules. The bases in the base pairs have a dihedral angle of 7 to 8 degrees between them.     

Three-dimensional structures of complexes of Lathyrus ochrus isolectin I with glucose and mannose: fine specificity of the monosaccharide-binding site

The structure of the methyl-alpha-D-mannopyranoside-LOL I complex has been solved by the molecular replacement method using the refined saccharide-free LOL I coordinates as starting model. The methyl-alpha-D-mannopyranoside-LOL I complex was refined by simulated annealing using the program X-PLOR. The final R-factor value is 0.182 [Fo greater than 1 sigma(Fo)]. The isostructural methyl-alpha-D-glucopyranoside-LOL I complex was refined by X-Ray coupled energy minimization using the methyl-alpha-D-mannopyranoside-LOL I structure as a starting model to an R factor of 0.179 (all data). In both crystal forms, each dimer binds two molecules of sugar in pockets found near the calcium ions. The two saccharide moieties, which are in the C1 chair conformation, establish the same hydrogen bond pattern with the lectin. However, the van der Waals contacts are different between the O2, C2, C6, and O6 atoms of the two molecules and the backbone atoms of residues 208-211. Mannose, due to its axial C2 conformation, encloses the backbone atoms of the protein in a clamplike way. Van der Waals energy calculations suggest that this better complementarity of the mannoside molecule with the lectin could explain its higher affinity for isolectin I.     

IS50-mediated inverse transposition. Discrimination between the two ends of an IS element

The crystal and molecular structure of l-pyroglutamyl-β-(2-thienyl)-l-alanyl-l-prolinamide, < Glu-Thi-Pro-NH₂(Thi²-TRH), C₁₇H₂₂N₄O₄S, has been determined from X-ray diffraction data. Thi²-TRH is a highly active analogue of thyroliberin, a thyrotropin-releasing hormone (TRH), in which the imidazole ring of the central histidine moiety in the natural hormone has been replaced by a 2-thienyl group. Thi²-TRH crystallizes from water in the monoclinic space group P2₁, $\text{a = 9.340(1)}\text{A}\text{?}\text{, b = 21.961(3)}\text{A}\text{?}\text{, c = 9.449(1)}\text{A}\text{?}$ and β = 109.58(1) °, with two molecules per asymmetric unit. These independent molecules, A and B, have the same general backbone conformation with the φ₂, ψ₂ and ψ₃ torsional angles close to ?90 °, +120 ° and +150 °, respectively, but they show different magnitudes of rotational disorder in the thiophene ring as well as a certain disorder in the pyrrolidine ring. A and B are cross-linked by four interchain hydrogen bonds, forming a two-stranded antiparallel β-pleated sheet structure. The molecules in these dimer fragments are further hydrogen-bonded to successive translated molecules along the a and c axes, forming a pronounced two-dimensional predominantly hydrophobic layer structure. These layers, in which the atoms are almost equally arranged on both sides, are separated by ordinary van der Waals' distances. A close correlation between the molecular conformation in the solid state and the preferential conformation in solution is found. It is concluded that the crystalline structure of Thi²-TRH possesses structural features which may be of relevance in the hormone-receptor interaction process.     

Ice-binding mechanism of winter flounder antifreeze proteins

We have studied the winter flounder antifreeze protein (AFP) and two of its mutants using molecular dynamics simulation techniques. The simulations were performed under four conditions: in the gas phase, solvated by water, adsorbed on the ice (2021) crystal plane in the gas phase and in aqueous solution. This study provided details of the ice-binding pattern of the winter flounder AFP. Simulation results indicated that the Asp, Asn, and Thr residues in the AFP are important in ice binding and that Asn and Thr as a group bind cooperatively to the ice surface. These ice-binding residues can be collected into four distinct ice-binding regions: Asp-1/Thr-2/Asp-5, Thr-13/Asn-16, Thr-24/Asn-27, and Thr-35/Arg-37. These four regions are 11 residues apart and the repeat distance between them matches the ice lattice constant along the (1102) direction. This match is crucial to ensure that all four groups can interact with the ice surface simultaneously, thereby, enhancing ice binding. These Asx (x = p or n)/Thr regions each form 5-6 hydrogen bonds with the ice surface: Asn forms about three hydrogen bonds with ice molecules located in the step region while Thr forms one to two hydrogen bonds with the ice molecules in the ridge of the (2021) crystal plane. Both the distance between Thr and Asn and the ordering of the two residues are crucial for effective ice binding. The proper sequence is necessary to generate a binding surface that is compatible with the ice surface topology, thus providing a perfect "host/guest" interaction that simultaneously satisfies both hydrogen bonding and van der Waals interactions. The results also show the relation among binding energy, the number of hydrogen bonds, and the activity. The activity is correlated to the binding energy, and in the case of the mutants we have studied the number of hydrogen bonds. The greater the number of the hydrogen bonds the greater the antifreeze activity. The roles van der Waals interactions and the hydrophobic effect play in ice binding are also highlighted. For the latter it is demonstrated that the surface of ice has a clathratelike structure which favors the partitioning of hydrophobic groups to the surface of ice. It is suggested that mutations that involve the deletion of hydrophobic residues (e.g., the Leu residues) will provide insight into the role the hydrophobic effect plays in partitioning these peptides to the surface of ice.     

Accurate representation of B-DNA double helical structure with implicit solvent and counterions

High-resolution nuclear magnetic resonance (NMR) and crystallographic data have been taken to refine the force field used in the torsion angle space nucleic acids molecular mechanics program DUPLEX. The population balance deduced from NMR studies of two carcinogen-modified DNA conformers in equilibrium was used to fine tune a sigmoidal, distance-dependent dielectric function so that reasonable relative energies could be obtained. In addition, the base-pair and backbone geometry from high-resolution crystal structures of the Dickerson-Drew dodecamer was used to re-evaluate the deoxyribose pseudorotation profile and the Lennard-Jones nonbonded energy terms. With a modified dielectric function that assumes a very steep distance-dependent form, a deoxyribose pseudorotation profile with reduced energy barriers between C2'- and C3'-endo minima, and a shift of the Lennard-Jones potential energy minimum to a distance approximately 0.4 A greater than the sum of the van der Waals' radii, the sequence-dependent conformational features of the Dickerson-Drew dodecamer in both the solid state and the aqueous liquid crystalline phase are well reproduced. The robust performance of the revised force field, in conjunction with its efficiency through implicit treatment of solvent and counterions, provides a valuable tool for elucidating conformations and structure-function relationships of DNA, including those of molecules modified by carcinogens and other ligands.     

A priori crystal structure prediction of native celluloses

The packing of beta-1,4-glucopyranose chains has been modeled to further elaborate the molecular structures of native cellulose microfibrils. A chain pairing procedure was implemented that evaluates the optimal interchain distance and energy for all possible settings of the two chains. Starting with a rigid model of an isolated chain, its interaction with a second chain was studied at various helix-axis translations and mutual rotational orientations while keeping the chains at van der Waals separation. For each setting, the sum of the van der Waals and hydrogen-bonding energy was calculated. No energy minimization was performed during the initial screening, but the energy and interchain distances were mapped to a three-dimensional grid, with evaluation of parallel settings of the cellulose chains. The emergence of several energy minima suggests that parallel chains of cellulose can be paired in a variety of stable orientations. A further analysis considered all possible parallel arrangements occurring between a cellulose chain pair and a further cellulose chain. Among all the low-energy three-chain models, only a few of them yield closely packed three-dimensional arrangements. From these, unit-cell dimensions as well as lattice symmetry were derived; interestingly two of them correspond closely to the observed allomorphs of crystalline native cellulose. The most favorable structural models were then optimized using a minicrystal procedure in conjunction with the MM3 force field. The two best crystal lattice predictions were for a triclinic (P(1)) and a monoclinic (P2(1)) arrangement with unit cell dimensions a = 0.63, b = 0.69, c = 1.036 nm, alpha = 113.0, beta = 121.1, gamma = 76.0 degrees, and a = 0.87, b = 0.75, c = 1.036 nm, gamma = 94.1 degrees, respectively. They correspond closely to the respective lattice symmetry and unit-cell dimensions that have been reported for cellulose Ialpha and cellulose Ibeta allomorphs. The suitability of the modeling protocol is endorsed by the agreement between the predicted and experimental unit-cell dimensions. The results provide pertinent information toward the construction of macromolecular models of microfibrils.     

Role of 3'-end sequences in infectivity of poliovirus transcripts made in vitro

It has been shown by van der Werf et al. (S. van der Werf, J. Bradley, E. Wimmer, F. W. Studier, and J. Dunn, Proc. Natl. Acad. Sci. USA 83:2330-2334, 1986) that in vitro synthesis of poliovirus RNA by T7 RNA polymerase gives rise to infectious RNA molecules; however, these molecules are only 5% as infectious as RNA isolated from virions. A plasmid, T7D-polio, was constructed that allows the in vitro synthesis of full-length RNA molecules with two additional guanine residues at the 5' end. However, T7D-polio differed from the construct of van der Werf et al. in that RNA transcribed from T7D-polio has an authentic 3' end, ending with only a polyadenine nucleotide sequence. Transfection of these RNA molecules into mammalian cells produced wild-type poliovirus with an efficiency similar to that of virion RNA. The use of this vector in the characterization of viral mutants in vivo and in vitro is discussed.     

Temperature dependence of kinetic interactions between progesterone and the uterine cytoplasmic receptor

The temperature dependence of the rates of dissociation and association for progesterone-receptor interactions was measured over the temperature range of 0–20°C. The dissociation process is biphasic indicating that either two forms of receptor are present or that the binding of progesterone to the receptor is a concatenated reaction.The enthalpy of activation for the dissociation of progesterone from the receptor is about 26–28 kcal/mol and the entropic energy of activation is about ?5 kcal/mol. The enthalpy of activation for the association of these molecules is about 3 kcal/mol and the entropic energy of activation is about 6 kcal/mol. These data are consistent with a model of progesterone binding to the receptor that includes hydrogen bonds between each of the two ketone groups and hydrogen donors on the receptor protein and involves van der Waals' interactions, due to the close proximity of the receptor binding site to a large fraction of the progesterone surface.     

Influence of cation-pi interactions in different folding types of membrane proteins

Cation-pi interactions play an important role to the stability of protein structures. In this work, we analyze the influence of cation-pi interactions in three-dimensional structures of membrane proteins. We found that transmembrane strand (TMS) proteins have more number of cation-pi interactions than transmembrane helical (TMH) proteins. In TMH proteins, both the positively charged residues Lys and Arg equally experience favorable cation-pi interactions whereas in TMS proteins, Arg is more likely than Lys to be in such interactions. There is no relationship between number of cation-pi interactions and number of residues in TMH proteins whereas a good correlation was observed in TMS proteins. The average cation-pi interaction energy for TMH proteins is -16 kcal/mol and that for TMS proteins is -27 kcal/mol. The pair-wise cation-pi interaction energy between aromatic and positively charged residues showed that Lys-Trp energy is stronger in TMS proteins than TMH proteins; Arg-Phe, Arg-Tyr and Lys-Phe have higher energy in TMH proteins than TMS proteins. The decomposition of energies into electrostatic and van der Waals revealed that the contribution from electrostatic energy is twice as that from van der Waals energy in both TMH and TMS proteins. The results obtained in the present study would be helpful to understand the contribution of cation-pi interactions to the stability of membrane proteins.     

Long-distance interactions between Ehrlich ascites tumour cells.

In previous work, electron micrographs were made of adjacent surfaces of aldehyde-fixed, Ehrlich ascites tumour cells cultured on coverslips, after reacting some of their negatively charged surface sites with colloidal iron hydroxide (CIH) particles. It was observed that microvilli from one cell were aligned with intermicrovillus regions on another, where the density of the adsorbed CIH particles was significantly lower than in adjacent regions. Alignment, which was considered to represent interactions between the two peripheral cellular regions, took place when these regions were apparently separated by more than 200 nm, in an environment of physiologic ionic strength ( 0·145 m NaCl).In this communication we attempt to find feasible mechanisms for the alignment phenomenon in physical terms, in cases where the observed separation of 200 nm is correct, and in cases where the distances are overestimated due to preparative artifacts.It is concluded, that at distances of separation in excess of 200 nm, one feasible mechanism for alignment is that net negatively charged macromolecules diffusing out of cells in the region of their microvilli, electrostatically repel CIH-binding anionic sites in the lipid-rich “fluid” matrix of the periphery of the opposed cell, causing gaps in their distribution. The role of electrostatic and electrodynamic (van der Waals') forces in causing alignment is also discussed in terms of distance of separation.This communication is concerned with the interpretation in terms of various interactions, of electron micrographs showing evidence of alignment between microvilli from one cell with specific areas of another.     

Crystalline A-dna: the X-ray analysis of the fragment d(G-G-T-A-T-A-C-C)

An A-DNA type double helical conformation was observed in the single crystal X-ray structure of the octamer d(G-G-T-A-T-A-C-C), 1, and its 5-bromouracil-containing analogue, 2. The structure of the isomorphous crystals (space group P61) was solved by a search technique based on packing criteria and R-factor calculations, with use of only low order data. At the present stage of refinement the R factors are 31% for 1 and 28% for 2 at a resolution of 2.25 A (0.225 nm). The molecules interact through their minor grooves by hydrogen bonding and base to sugar van der Waals contacts. The stable A conformation observed in the crystal may have some structural relevance to promoter regions where the T-A-T-A sequence is frequently found.     

The phage 434 Cro/OR1 complex at 2.5 A resolution

The crystal structure of phage 434 Cro protein in complex with a 20 base-pair DNA fragment has been determined to 2.5 A resolution. The DNA fragment contains the sequence of the OR1 operator site. The structure shows a bent conformation for the DNA, straighter at the center and more bent at the ends. The central base-pairs adopt conformations with significant deviations from coplanarity. The two molecules interact extensively along their common interface, both through hydrogen bonds and van der Waals interactions. The significance of these interactions for operator binding and recognition is discussed.     

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.