Statistical and molecular dynamics studies of buried waters in globular proteins

Buried solvent molecules are common in the core of globular proteins and contribute to structural stability. Folding necessitates the burial of polar backbone atoms in the protein core, whose hydrogen-bonding capacities should be satisfied on average. Whereas the residues in alpha-helices and beta-sheets form systematic main-chain hydrogen bonds, the residues in turns, coils and loops often contain polar atoms that fail to form intramolecular hydrogen bonds. The statistical analysis of 842 high resolution protein structures shows that well-resolved, internal water molecules preferentially reside near residues without alpha-helical and beta-sheet secondary structures. These buried waters most often form primary hydrogen bonds to main-chain atoms not involved in intramolecular hydrogen bonds, providing strong evidence that hydrating main-chain atoms is a key structural role of buried water molecules. Additionally, the average B-factor of protein atoms hydrogen-bonded to waters is smaller than that of protein atoms forming intramolecular hydrogen bonds, and the average B-factor of water molecules involved in primary hydrogen bonds with main-chain atoms is smaller than the average B-factor of water molecules involved in secondary hydrogen bonds to protein atoms that form concurrent intramolecular hydrogen bonds. To study the structural coupling between internal waters and buried polar atoms in detail we simulated the dynamics of wild-type FKBP12, in which a buried water, Wat137, forms one side-chain and multiple main-chain hydrogen bonds. We mutated E60, whose side-chain hydrogen bonds with Wat137, to Q, N, S or A, to modulate the multiplicity and geometry of hydrogen bonds to the water. Mutating E60 to a residue that is unable to form a hydrogen bond with Wat137 results in reorientation of the water molecule and leads to a structural readjustment of residues that are both near and distant to the water. We predict that the E60A mutation will result in a significantly reduced affinity of FKBP12 for its ligand FK506. The propensity of internal waters to hydrogen bond to buried polar atoms suggests that ordered water molecules may constitute fundamental structural components of proteins, particularly in regions where alpha-helical or beta-sheet secondary structure is not present.     

Relationship between intramolecular hydrogen bonding and solvent accessibility of side-chain donors and acceptors in proteins

This study shows that intramolecular hydrogen bonding in proteins depends on the accessibility of donors and acceptors to water molecules. The frequency of occurrence of H-bonded side chains in proteins is inversely proportional to the solvent accessibility of their donors and acceptors. Estimates of the notional free energy of hydrogen bonding suggest that intramolecular hydrogen-bonding interactions of buried and half-buried donors and acceptors can contribute favorably to the stability of a protein, whereas those of solvent-exposed polar atoms become less favorable or unfavorable.     

Patterns in ionizable side chain interactions in protein structures

In a selected set of 44 high-resolution, non-homologous protein structures, the intramolecular hydrogen bonds or salt bridges formed by ionizable amino acid side chains were identified and analyzed. The analysis was based on the investigation of several properties of the involved residues such as their solvent exposure, their belonging to a certain secondary structural element, and their position relative to the N- and C-termini of their respective structural element. It was observed that two-thirds of the interactions made by basic or acidic side chains are hydrogen bonds to polar uncharged groups. In particular, the majority (78%) of the hydrogen bonds between ionizable side chains and main chain polar groups (sch:mch bonds) involved at least one buried atom, and in 42% of the cases both interacting atoms were buried. In α-helices, the sch:mch bonds observed in the proximity of the C- and N-termini show a clear preference for acidic and basic side chains, respectively. This appears to be due to the partial charges of peptide group atoms at the termini of α-helices, which establish energetically favorable electrostatic interactions with side chain carrying opposite charge, at distances even greater than 4.5 Å. The sch:mch interactions involving ionizable side chains that belong either to β-strands or to the central part of α-helices are based almost exclusively on basic residues. This results from the presence of main chain carbonyl oxygen atoms in the protein core which have unsatisfied hydrogen bonding capabilities.     

Comparative X-ray studies on the interaction of carotenoids with a model phosphatidylcholine membrane

The interaction of structurally different carotenoids with a membrane molecular model was examined by X-ray diffraction. The selected compounds were beta-carotene, lycopene, lutein, violaxanthin, zeaxanthin, and additionally carotane, a fully saturated derivative of beta-carotene. They present similarities and differences in their rigidity, the presence of terminal ionone rings and hydroxy and epoxy groups bound to the rings. The membrane models were multibilayers of dipalmitoylphosphatidylcholine (DPPC), chosen for this investigation because the 3 nm thickness of the hydrophobic core of its bilayer coincides with the thickness of the hydrophobic core of thylakoid membranes and the length of the carotenoid molecules. Results indicate that the six compounds induced different types and degrees of structural perturbations to DPPC bilayers in aqueous media. They were interpreted in terms of the molecular characteristics of DPPC and the carotenoids. Lycopene and violaxanthin induced the highest structural damage to the acyl chain and polar headgroup regions of DPPC bilayers, respectively.     

Influence of hydrogen bonds on edge-to-face interactions between pyridine molecules

Edge-to-face interactions between two pyridine molecules and the influence of simultaneous hydrogen bonding of one or both of the pyridines to water on those interactions were studied by analyzing data from ab initio calculations. The results show that the edge-to-face interactions of pyridine dimers that are hydrogen bonded to water are generally stronger than those of non-H-bonded pyridine dimers, especially when the donor pyridine forms a hydrogen bond. The binding energy of the most stable edge-to-face interacting H-bonded pyridine dimer is ?5.05 kcal/mol, while that for the most stable edge-to-face interacting non-H-bonded pyridine dimer is ?3.64 kcal/mol. The interaction energy data obtained in this study cannot be explained solely by the differences in electrostatic potential between pyridine and the pyridine–water dimer. However, the calculated cooperative effect can be predicted using electrostatic potential maps.     

The role of hydrogen bonding via interfacial water molecules in antigen-antibody complexation. The HyHEL-10-HEL interaction

To study the role of hydrogen bonding via interfacial water molecules in protein-protein interactions, we examined the interaction between hen egg white lysozyme (HEL) and its HyHEL-10 variable domain fragment (Fv) antibody. We constructed three antibody mutants (l-Y50F, l-S91A, and l-S93A) and investigated the interactions between the mutant Fvs and HEL. Isothermal titration calorimetry indicated that the mutations significantly decreased the negative enthalpy change (8-25 kJ mol(-1)), despite some offset by a favorable entropy change. X-ray crystallography demonstrated that the complexes had nearly identical structures, including the positions of the interfacial water molecules. Taken together, the isothermal titration calorimetric and x-ray crystallographic results indicate that hydrogen bonding via interfacial water enthalpically contributes to the Fv-HEL interaction despite the partial offset because of entropy loss, suggesting that hydrogen bonding stiffens the antigen-antibody complex.     

Intramolecular compactization of circular DNA in interaction with dansyl hydrazide trivaline leading to a formation of a new type of structure--triple rings

Compactization of supercoiled circular plasmid pBR322 caused by interaction with synthetic oligopeptide dansyl hydrazide trivaline capable of beta-structure formation was studied by electron microscopy. The results show that at rising input peptide concentration circular DNA molecules undergo intramolecular structural transition with the formation of compact ring structures. The compact ring structures are formed by the fiber having the thickness of 60 A. The analysis of morphology of intermediate structures and the contour length measurements enable us to conclude that 60 A-fiber contains three lying side-by-side and interwound double-stranded DNA segments. Thus, the compact ring structures are addressed to as triple rings. The triple ring have one special point, where the triple region ends are locked by a duplex DNA segment. The mechanisms responsible for the triple ring formation may be of importance for DNA and chromatin compactization processes in vivo.     

Structural features which distinguish estrogen agonists and antagonists

The molecular structures of compounds that compete for the estrogen binding site suggest that when estradiol is bound to the receptor, there is a close fit only at the A-ring end of the steroid. The most potent antagonists have phenolic rings capable of mimicking the estradiol A-ring in promoting high affinity binding to the receptor. They fail to stabilize the conformational change or molecular interaction needed to achieve hormonal response because they either lack an essential functional group (hydrogen bond donor) or they present a steric block or topological feature incompatible with transformation or interaction subsequent to the initiation of binding by the A-ring or A-ring analogues. X-Ray analyses of 12 triphenylethylene derivatives demonstrate the stability of a specific conformation that is unaffected by crystal packing forces. This conformation has a pinwheel-like orientation of phenyl rings, the direction of which is correlated with a 10 degrees twist about the central double bond and appears to depend upon the orientation of the non-phenyl substituent relative to the double bond. Empirical energy calculations are insensitive to this intramolecular structural dependence and incorrectly predict that the pinwheel of the opposite direction is of lower energy.     

Structure of Sphingomyelin Bilayers: A Simulation Study

We have carried out a molecular dynamics simulation of a hydrated 18:0 sphingomyelin lipid bilayer. The bilayer contained 1600 sphingomyelin (SM) molecules, and 50,592 water molecules. After construction and initial equilibration, the simulation was run for 3.8 ns at a constant temperature of 50°C and a constant pressure of 1 atm. We present properties of the bilayer calculated from the simulation, and compare with experimental data and with properties of dipalmitoyl phosphatidylcholine (DPPC) bilayers. The SM bilayers are significantly more ordered and compact than DPPC bilayers at the same temperature. SM bilayers also exhibit significant intramolecular hydrogen bonding between phosphate ester oxygen and hydroxyl hydrogen atoms. This results in a decreased hydration in the polar region of the SM bilayer compared with DPPC. Since our simulation system is very large we have calculated the power spectrum of bilayer undulation and peristaltic modes, and we compare these data with similar calculations for DPPC bilayers. We find that the SM bilayer has significantly larger bending modulus and area compressibility compared to DPPC.     

<Emphasis Type="Italic">Z</Emphasis>-Isomers of (4α→6″, 2α→<Emphasis Type="Italic">O</Emphasis>→1″)-phenylflavan substituted with R′=R=OH. Conformational properties,electronic structure and aqueous solvent effects

Procyanidins are highly hydroxylated polymers known as antioxidant compounds, thereby exhibiting beneficial effects. These compounds are protective agents against oxidative stress and the damage induced by free radicals in membranes and nucleic acids. This paper describes a study of the conformational space of (4α→6″, 2α→O→1″)-phenylflavan substituted with R′=R=OH as part of a larger study of similar structures with different substitutions. The relationships between aqueous solution–vacuum variations of some properties were studied, as well as the stabilization and reactivity of (4α→6″, 2α→O→1″)-phenylflavan substituted with R′=R=H, R′=H, R=OH, R′=R=OH, and (+)-catechin. The variations in geometric parameters and electronic properties due to conformational changes, as well as the effects of substituents and polar solvents, were evaluated and analyzed. Bader’s theory of atoms in molecules was applied to characterize intramolecular interactions, along with a natural bond orbital analysis for each conformer described. The molecular electrostatic potential was rationalized by charge delocalization mechanisms and interatomic intramolecular interactions, relating them to the structural changes and topological properties of the electron charge density. Molecular polarizability and permanent electric dipole moment values were estimated. The results show the importance of a knowledge of the conformational space, and values for each conformer. Based on our previous results, we showed the existence of electron charge delocalization mechanisms acting cooperatively as “delocalization routes”, showing interactions between different rings not even sharing the same plane. These “delocalization routes” were more effective for (4α→6″, 2α→O→1″)-phenylflavan substituted with R′=R=OH than for (+)-catechin, and are proposed as adding insight into the structure–antioxidant activity relationship of flavans.     

Crystal and molecular structures of trichloro-cobalt(II) complexes of epiquinine, epiquinidine, and epidihydrocinchonine

Cinchona alkaloids are very well known antimalarials but the mechanism of their biological action still remains to be elucidated. The structural studies of active erythro and inactive threo alkaloid complexes are an important step to this aim. In this paper results of crystal structure analysis of three cobalt complexes of threo alkaloids are presented: (epiquininium)trichlorocobalt(II) (EpiQnCoCl3), (epiquinidinium)trichlorocobalt(II) (EpiQdCoCl3) and (epidihydrocinchoninium)trichlorocobalt(II) (EpiCnCoCl3). The complexes are zwitterions in which trichlorocobalt substituents are coordinated to quinoline nitrogen atoms and quinuclidine nitrogen atoms are protonated. EpiQnCoCl3 adopts uncommon conformation with quinoline moiety oriented in the opposite direction in comparison to the analogous uncomplexed alkaloid. The packing in the crystal structures is determined mainly by the hydrogen bonds, in which the chlorine atoms of substituents and solvent molecules contribute. Atoms participating in hydrogen bonds in EpiQnCoCl3 and EpiQdCoCl3 form large rings, while in EpiCnCoCl3 only chains are present. Solvent molecules are very important for the packing mode. In contrast to most erythro alkaloids, the hydroxyl oxygen atom in the title complexes forms weak or not well defined hydrogen bonds. The contribution of very weak intramolecular interactions N1--H1...O12 cannot be excluded. Such "trace" interactions can be considered a relic of the unprotonated status of an epi alkaloid.     

Interactions of d(10) metal ions with hippuric acid and cytosine. X-ray structure of the first cadmium (II)-amino acid derivative-nucleobase ternary compound.

The interactions of Zn(II), Cd(II) and Hg(II) with hippuric acid (hipH) were studied and several novel compounds were synthesized and studied by NMR. Some new metal-hippuric-cytosine ternary compounds were formed and the structure of the [Cd(hip)(2)(cyt)(H(2)O)](2) ternary complex resolved. Each cadmium (II) atom has a distorted trigonal bipyramid coordination which is linked to a water molecule, a cytosine via N(3), a carboxylic oxygen atom of a hippurate moiety and two bridging dicoordinated hippurates bound through the carboxylic oxygen atoms. To these five main bonds, two longer ancillary interactions can be observed: the second oxygen of the monocoordinated hippurate group and the carboxylic oxygen of the cytosine ligand. The compound is stabilized by an intramolecular stacking between the benzene and cytosine rings and by the hydrogen bonds between the coordinated water molecules and the ligands. This is, to our knowledge, the first structure of a cadmium-amino acid derivative-natural nucleobase compound described so far.     

Novicidin's membrane permeabilizing activity is driven by membrane partitioning but not by helicity: A biophysical study of the impact of lipid charge and cholesterol

We have investigated the interactions between the antimicrobial peptide Novicidin (Nc) and vesicles containing the phospholipid DOPC, with various amounts of DOPG and cholesterol using circular dichroism spectroscopy, calcein release, equilibrium dialysis and isothermal titration calorimetry. Nc adopts a random coil structure in the absence of lipids and in the presence of vesicles containing 100% DOPC. Lipids with 25–40% DOPG induce the highest level of helicity in Nc; higher DOPG levels lead to lower helicity levels and an altered tertiary arrangement of the peptide. However, the ability of Nc to permeabilize vesicles correlates not with helicity but rather with its overall membrane affinity, which is enthalpically favorable but opposed by entropy. Permeabilization declines with increasing mole percentage PG. Changes in helicity correlate with changes in enthalpy, reflecting the enthalpy of helix formation, but not with affinity. There is also a large favorable enthalpic interaction between Nc and lipids in the absence of negative charge and structural changes. Cholesterol slightly reduces membrane permeabilization but has little effect on Nc affinity and secondary structure, and probably protects the membrane by inducing the liquid ordered state. We conclude that helicity is not a prerequisite for activity, and charge–charge interactions are not the only major driving force for AMP interactions with membranes. Our data are compatible with a model in which a superficial binding mode with a large membrane surface binding area per peptide is more efficient than a more intimate embedding within the membrane environment.     

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.     

Sequence dependence of beta-hairpin structure: comparison of a salt bridge and an aromatic interaction

A comparison of the contributions and position dependence of cross-strand electrostatic and aromatic side-chain interactions to beta-sheet stability has been performed by using nuclear magnetic resonance in a well-folded beta-hairpin peptide of the general sequence XRTVXVdPGOXITQX. Phe-Phe and Glu-Lys pairs were varied at the internal and terminal non-hydrogen-bonded position, and the resulting stability was measured by the effects on alpha-hydrogen and aromatic hydrogen chemical shifts. It was determined that the introduction of a Phe-Phe pair resulted in a more folded peptide, regardless of position, and a more tightly folded core. Substitution of the Glu-Lys pair at the internal position results in a less folded peptide and increased fraying at the terminal residues. Upfield shifting of the aromatic hydrogens provided evidence for an edge-face aromatic interaction, regardless of position of the Phe-Phe pair. In peptides with two Phe-Phe pairs, substitution with Glu-Lys at either position resulted in a weakening of the aromatic interaction and a subsequent decrease in peptide stability. Thermal denaturation of the peptides containing Phe-Phe indicates that the aromatic interaction is enthalpically favored, whereas the folding of hairpins with cross-strand Glu-Lys pairs was less enthalpically favorable but entropically more favorable.     

Models of haloadaptation in bacterial membranes

Abstract Cell membranes consist of a complex assortment of amphipathic lipids. These lipids exist in one of three phases in aqueous systems at the growth temperature of the organism: namely, lamellar gel, lamellar liquid-crystalline or hexagonal-II. The phase behaviour is modified by interaction of the lipids with other membrane components and electrolytes. A stable membrane structure is achieved when the polar and non-polar interactions are balanced such that a durable bilayer arrangement is formed into which the various membrane proteins are integrated. The effect of surface charge on phase domain behaviour of the membrane lipids and the modulation by electrolytes is crucial to understanding how halophiles adapt to high-salt environments.     

Prudent modeling of core polar residues in computational protein design

Hydrogen bond interactions were surveyed in a set of protein structures. Compared to surface positions, polar side-chains at core positions form a greater number of intra-molecular hydrogen bonds. Furthermore, the majority of polar side-chains at core positions form at least one hydrogen bond to main-chain atoms that are not involved in hydrogen bonds to other main-chain atoms. Based on this structural survey, hydrogen bond rules were generated for each polar amino acid for use in protein core design. In the context of protein core design, these prudent polar rules were used to eliminate from consideration polar amino acid rotamers that do not form a minimum number of hydrogen bonds. As an initial test, the core of Escherichia coli thioredoxin was selected as a design target. For this target, the prudent polar strategy resulted in a minor increase in computational complexity compared to a strategy that did not allow polar residues. Dead-end elimination was used to identify global minimum energy conformations for the prudent polar and no polar strategies. The prudent polar strategy identified a protein sequence that was thermodynamically stabilized by 2.5 kcal/mol relative to wild-type thioredoxin and 2.2 kcal/mol relative to a thioredoxin variant whose core was designed without polar residues.     

Molecular dynamics simulation of solvated azurin: correlation between surface solvent accessibility and water residence times

A system containing the globular protein azurin and 3,658 water molecules has been simulated to investigate the influence on water dynamics exerted by a protein surface. Evaluation of water mean residence time for elements having different secondary structure did not show any correlation. Identically, comparison of solvent residence time for atoms having different charge and polarity did not show any clear trend. The main factor influencing water residence time in proximity to a specific site was found to be its solvent accessibility. In detail for atoms belonging to lateral chains and having solvent-accessible surface lower than approximately 16 A(2)a relation is found for which charged and polar atoms are surrounded by water molecules characterized by residence times longer than the non polar ones. The involvement of the low accessible protein atom in an intraprotein hydrogen bond further modulates the length of the water residence time. On the other hand for surfaces having high solvent accessibility, all atoms, independently of their character, are surrounded by water molecules which rapidly exchange with the bulk solvent. Proteins 2000;39:56-67.