Recent developments in research on water oxidation by photosystem II

Photosynthetic water oxidation chemistry at the unique manganese-calcium complex of photosystem II (PSII) is of fundamental importance and serves as a paragon in the development of efficient synthetic catalysts. A recent crystal structure of PSII shows the atoms of the water-oxidizing complex; its Mn4CaO5 core resembles inorganic manganese-calcium oxides. Merging of crystallographic and spectroscopic information reverses radiation-induced modifications at the Mn-complex in silico and facilitates discussion of the O-O bond chemistry. Coordinated proton movements are promoted by a water network connecting the Mn4CaO5 core with the oxidant, a tyrosine radical and one possibly mobile chloride ion. A basic reaction-cycle model predicts an alternating proton and electron removal from the catalytic site, which facilitates energetically efficient water oxidation.     

Influence of stereochemistry on proton transfer in protonated tripeptide models

Vectorial proton transfer among carbonyl oxygen atoms was studied in two models of tripeptide via quantum chemical calculations using the hybrid B3LYP functional and the 6-31++G** basis set. Two principal proton transfer pathways were found: a first path involving isomerization of the proton around the double bond of the carbonyl group, and a second based on the large conformational flexibility of the tripeptide model where all carbonyl oxygen atoms cooperate. The latter pathway has a rate-determining step energy barrier that is only around half of that for the first pathway. As conformational flexibility plays a crucial role in second pathway, the effect of attaching methyl groups to the alpha carbon atoms was studied. The results obtained are presented for all four possible stereochemical configurations.     

Mechanical property of a TIM-barrel protein

The mechanical response of a TIM-barrel protein to an applied pressure has been studied. We generated structures under an applied pressure by assuming the volume change to be a linear function of normal mode variables. By Delaunay tessellation, the space occupied by protein atoms is divided uniquely into tetrahedra, whose four vertices correspond to atomic positions. Based on the atoms that define them, the resulting Delaunay tetrahedra are classified as belonging to various secondary structures in the protein. The compressibility of various regions identified with respect to secondary structural elements in this protein is obtained from volume changes of respective regions in two structures with and without an applied pressure. We found that the β barrel region located at the core of the protein is quite soft. The interior of the β barrel, occupied by side chains of β strands, is the softest. The helix, strand, and loop segments themselves are extremely rigid, while the regions existing between these secondary structural elements are soft. These results suggest that the regions between secondary structural elements play an important role in protein dynamics. Another aspect of tetrahedra, referred to as bond distance, is introduced to account for rigidities of the tetrahedra. Bond distance is a measure of separation of the atoms of a tetrahedron in terms of number of bonds along the polypeptide chain or side chains. Tetrahedra with longer bond distances are found to be softer on average. From this behavior, we derive a simple empirical equation, which well describes the compressibilities of various regions. © 1997 Wiley-Liss Inc.     

Analysis of conformational parameters in nucleic acid fragments. II. Co-crystal complexes of nucleic acid bases.

Studies have been made of conformational parameters in co-crystal complexes and compounds of nucleic acid bases in which there is the possibility of formation of hetero-base-pairs. Using published data extracted from the Cambridge structural database, a total of 37 base-pairs were found, of which 25 were hetero-pairs and 12 homo-pairs. These base-pairs were subject to analysis to reveal hydrogen bond parameters, propeller twist, buckle and C1'-C1' separation (or a similar parameter if C1' atoms were not present). Hetero-pairs were found to show larger twists than homo-pairs, the magnitude of twist being unrelated to hydrogen bond parameters or buckle value. The propeller twisting is less pronounced in these nucleic acid bases than in nucleosides, but still has a significant magnitude. Propeller twisting in hetero-pairs is found to be larger than in homo-pairs. Hetero-pairs appear to be formed preferentially in competitive situations.     

Relaxivities of paramagnetic liposomes: on the importance of the chain type and the length of the amphiphilic complex

Nuclear magnetic relaxation dispersion (NMRD) profiles of unilamellar DPPC liposomes incorporating Gd-DTPA-bisamides with alkyl chains of 12 to 18 C atoms in their external and internal layers were recorded in order to study the influence that the chain length and structure of Gd-bisamides incorporated in the liposomal membrane have on their proton relaxivity. The NMRD profiles recorded at 310 K show that the relaxivity reaches a minimum value when the carbon chain lengths of the phospholipid and of the Gd complex match and is at a maximum in the presence of a carbon-carbon double bond. For these DPPC paramagnetic liposomes, the longer the aliphatic chains of the complex, the larger will be its immobilization in the membrane. In addition, the presence of an unsaturated carbon-carbon bond in the alkyl chain of the Gd complex induces an increase of its mobility and of its water exchange rate with, as a result, a much greater efficiency as an MRI contrast agent.     

The role of hydrogen bonding in the enzymatic reaction catalyzed by HIV-1 protease

The hydrogen-bond network in various stages of the enzymatic reaction catalyzed by HIV-1 protease was studied through quantum-classical molecular dynamics simulations. The approximate valence bond method was applied to the active site atoms participating directly in the rearrangement of chemical bonds. The rest of the protein with explicit solvent was treated with a classical molecular mechanics model. Two possible mechanisms were studied, general-acid/general-base (GA/GB) with Asp 25 protonated at the inner oxygen, and a direct nucleophilic attack by Asp 25. Strong hydrogen bonds leading to spontaneous proton transfers were observed in both reaction paths. A single-well hydrogen bond was formed between the peptide nitrogen and outer oxygen of Asp 125. The proton was diffusely distributed with an average central position and transferred back and forth on a picosecond scale. In both mechanisms, this interaction helped change the peptide-bond hybridization, increased the partial charge on peptidyl carbon, and in the GA/GB mechanism, helped deprotonate the water molecule. The inner oxygens of the aspartic dyad formed a low-barrier, but asymmetric hydrogen bond; the proton was not positioned midway and made a slightly elongated covalent bond, transferring from one to the other aspartate. In the GA/GB mechanism both aspartates may help deprotonate the water molecule. We observed the breakage of the peptide bond and found that the protonation of the peptidyl amine group was essential for the peptide-bond cleavage. In studies of the direct nucleophilic mechanism, the peptide carbon of the substrate and oxygen of Asp 25 approached as close as 2.3 A.     

Short-strong hydrogen bonds and a low barrier transition state for the proton transfer reaction in RNase A catalysis: a quantum chemical study

There is growing evidence that some enzymes catalyze reactions through the formation of short-strong hydrogen bonds as first suggested by Gerlt and Gassman. Support comes from several experimental and quantum chemical studies that include correlation energies on model systems. In the present study, the process of proton transfer between hydroxyl and imidazole groups, a model of the crucial step in the hydrolysis of RNA by the enzymes of the RNase A family, is investigated at the quantum mechanical level of density functional theory and perturbation theory at the MP2 level. The model focuses on the nature of the formation of a complex between the important residues of the protein and the hydroxyl group of the substrate. We have also investigated different configurations of the ground state that are important in the proton transfer reaction. The nature of bonding between the catalytic unit of the enzyme and the substrate in the model is investigated by Bader's atoms in molecule theory. The contributions of solvation and vibrational energies corresponding to the reactant, the transition state and the product configurations are also evaluated. Furthermore, the effect of protein environment is investigated by considering the catalytic unit surrounded by complete proteins--RNase A and Angiogenin. The results, in general, indicate the formation of a short-strong hydrogen bond and the formation of a low barrier transition state for the proton transfer model of the enzyme.     

Refined crystal structure of ascorbate oxidase at 1.9 A resolution.

The crystal structure of the fully oxidized form of ascorbate oxidase (EC 1.10.3.3) from Zucchini has been refined at 1.90 A (1 A = 0.1 nm) resolution, using an energy-restrained least-squares refinement procedure. The refined model, which includes 8764 protein atoms, 9 copper atoms and 970 solvent molecules, has a crystallographic R-factor of 20.3% for 85,252 reflections between 8 and 1.90 A resolution. The root-mean-square deviation in bond lengths and bond angles from ideal values is 0.011 A and 2.99 degrees, respectively. The subunits of 552 residues (70,000 Mr) are arranged as tetramers with D2 symmetry. One of the dyads is realized by the crystallographic axis parallel to the c-axis giving one dimer in the asymmetric unit. The dimer related about this crystallographic axis is suggested as the dimer present in solution. Asn92 is the attachment site for one of the two N-linked sugar moieties, which has defined electron density for the N-linked N-acetyl-glucosamine ring. Each subunit is built up by three domains arranged sequentially on the polypeptide chain and tightly associated in space. The folding of all three domains is of a similar beta-barrel type and related to plastocyanin and azurin. An analysis of intra- and intertetramer hydrogen bond and van der Waals interactions is presented. Each subunit has four copper atoms bound as mononuclear and trinuclear species. The mononuclear copper has two histidine, a cysteine and a methionine ligand and represents the type-1 copper. It is located in domain 3. The bond lengths of the type-1 copper centre are comparable to the values for oxidized plastocyanin. The trinuclear cluster has eight histidine ligands symmetrically supplied from domain 1 and 3. It may be subdivided into a pair of copper atoms with histidine ligands whose ligating N-atoms (5 NE2 atoms and one ND1 atom) are arranged trigonal prismatic. The pair is the putative type-3 copper. The remaining copper has two histidine ligands and is the putative spectroscopic type-2 copper. Two oxygen atoms are bound to the trinuclear species as OH- or O2- and bridging the putative type-3 copper pair and as OH- or H2O bound to the putative type-2 copper trans to the copper pair. The bond lengths within the trinuclear copper site are similar to comparable binuclear model compounds. The putative binding site for the reducing substrate is close to the type-1 copper.(ABSTRACT TRUNCATED AT 400 WORDS)     

Structural basis of DNA flexibility

It has been recently shown by us, on the basis of crystal structure database that the flexibility of B-DNA double helices depends significantly on their base sequence. Our model building studies further indicated that the existence of bifurcated cross-strand hydrogen bonds between successive base pairs is possibly the main factor behind the sequence directed DNA flexibility. These cross-strand hydrogen bonds are, of course, weaker than the usual Watson-Crick hydrogen bonds and their bond geometry is characterized by relatively larger bond lengths and smaller bond angles. We have tried to improve our model structures by incorporating non-planarity of the amino groups in DNA bases due to the presence of lone pair electrons at the nitrogen atoms. Energy minimization studies have been carried out by using different quantum chemical methods, whereby it is found that in all cases of N-H....O type cross-strand hydrogen bonds, the bond geometry improves significantly. In the cases of N-H....N type hydrogen bonds, however, no such consistent improvements can be noticed. Perhaps the true picture would emerge only if all the other interactions present in the DNA macromolecule could be appropriately taken into account.     

Mechanisms for ligand binding to GluR0 ion channels: crystal structures of the glutamate and serine complexes and a closed apo state

High-resolution structures of the ligand binding core of GluR0, a glutamate receptor ion channel from Synechocystis PCC 6803, have been solved by X-ray diffraction. The GluR0 structures reveal homology with bacterial periplasmic binding proteins and the rat GluR2 AMPA subtype neurotransmitter receptor. The ligand binding site is formed by a cleft between two globular alpha/beta domains. L-Glutamate binds in an extended conformation, similar to that observed for glutamine binding protein (GlnBP). However, the L-glutamate gamma-carboxyl group interacts exclusively with Asn51 in domain 1, different from the interactions of ligand with domain 2 residues observed for GluR2 and GlnBP. To address how neutral amino acids activate GluR0 gating we solved the structure of the binding site complex with L-serine. This revealed solvent molecules acting as surrogate ligand atoms, such that the serine OH group makes solvent-mediated hydrogen bonds with Asn51. The structure of a ligand-free, closed-cleft conformation revealed an extensive hydrogen bond network mediated by solvent molecules. Equilibrium centrifugation analysis revealed dimerization of the GluR0 ligand binding core with a dissociation constant of 0.8 microM. In the crystal, a symmetrical dimer involving residues in domain 1 occurs along a crystallographic 2-fold axis and suggests that tetrameric glutamate receptor ion channels are assembled from dimers of dimers. We propose that ligand-induced conformational changes cause the ion channel to open as a result of an increase in domain 2 separation relative to the dimer interface.     

Analysis of pH-dependent elements in proteins: geometry and properties of pairs of hydrogen-bonded carboxylic acid side-chains

A rather frequent but so far little discussed observation is that pairs of carboxylic acid side-chains in proteins can share a proton in a hydrogen bond. In the present article, quantum chemical calculations of simple model systems for carboxyl-carboxylate interactions are compared with structural observations from proteins. A detailed structural analysis of the proteins deposited in the PDB revealed that, in a subset of proteins sharing less than 90% sequence identity, 19% (314) contain at least one pair of carboxylic acids with their side-chain oxygen atoms within hydrogen-bonding distance. As the distance between those interacting oxygen atoms is frequently very short ( approximately 2.55 A), many of these carboxylic acids are suggested to share a proton in a strong hydrogen bond. When situated in an appropriate structural environment (low dielectric constant), some might even form a low barrier hydrogen bond. The quantum chemical studies show that the most frequent geometric features of carboxyl-carboxylate pairs found in proteins, and no or symmetric ligation, are also the most stable arrangements at low dielectric constants, and they also suggest at medium and low pH a higher stability than for isosteric amide-carboxylate pairs. The presence of these pairs in 119 different enzymes found in the BRENDA database is set in relation to their properties and functions. This analysis shows that pH optima of enzymes with carboxyl-carboxylate pairs are shifted to lower than average values, whereas temperature optima seem to be increased. The described structural principles can be used as guidelines for rational protein design (e.g., in order to improve pH or temperature stability).     

High-precision measurement of hydrogen bond lengths in proteins by nuclear magnetic resonance methods.

We have compared hydrogen bond lengths on enzymes derived with high precision (< or = +/- 0.05 A) from both the proton chemical shifts (delta) and the fractionation factors (phi) of the proton involved with those obtained from protein X-ray crystallography. Hydrogen bond distances derived from proton chemical shifts were obtained from a correlation of 59 O--H....O hydrogen bond lengths, measured by small molecule high-resolution X-ray crystallography, with chemical shifts determined by solid-state nuclear magnetic resonance (NMR) in the same crystals (McDermott A, Ridenour CF, Encyclopedia of NMR, Sussex, U.K.: Wiley, 1996:3820-3825). Hydrogen bond distances were independently obtained from fractionation factors that yield distances between the two proton wells in quartic double minimum potential functions (Kreevoy MM, Liang TM, J Am Chem Soc, 1980;102:3315-3322). The high-precision hydrogen bond distances derived from their corresponding NMR-measured proton chemical shifts and fractionation factors agree well with each other and with those reported in protein X-ray structures within the larger errors (+/-0.2-0.8 A) in distances obtained by protein X-ray crystallography. The increased precision in measurements of hydrogen bond lengths by NMR has provided insight into the contributions of short, strong hydrogen bonds to catalysis for several enzymatic reactions.     

Abundance and arrangement of single,double and bifurcated hydrogen bonds in protein α-helices: Statistical analysis of PDB-select

Statistics are collected and analyzed for the possibility of hydrogen bonding in the secondary structures of globular proteins, based on geometric criteria. Double and bifurcated bonds are considered as pairs of admissible H-bonds with two proton donors or two proton acceptors, respectively. Most of such bonds belong to peptide groups in α-helices, with O _i …N _{i + 3} nearly as frequent as O _i …N _{i + 4}; in contrast, most of the 3/10-helical segments are too short to have any. Alternating double and bifurcated bonds in α-helices form an apparently cooperative network structure. A typical α-helical segment perhaps carries two stretches of the H-bond network broken in the middle. The constituent H-bonds are nonlinear: the hydrogen atom is off the straight line connecting the proton donor and proton acceptor atoms. This deflection is larger for H _{i + 3} vs. bond line O _i −N _{i + 3} than for H _{i + 4} vs. O _i −N _{i + 4}, and though the two kinds of bond have about the same length (exceeding those typical of low-molecular compounds), O _i …N _{i + 4} must be stronger than O _i …N _{i + 3}. Double/bifurcated bonds are also not coplanar, i.e., hydrogen atoms are beyond the N…O…N (or O…N…O) plane. The text was submitted by the authors in English.     

Viscosity dependence of protein dynamics

The influence of solvent viscosity on protein dynamics was investigated with molecular dynamics simulations of factor Xa in two solvents differing only in viscosity, by a factor of 10. We obtained this viscosity change by changing the masses of the solvent atoms by a factor of 100. Equilibrium properties of the protein, that is, the average structure, its fluctuations, and the secondary structure, show no significant dependence on the solvent viscosity. The dynamic properties of the protein, that is, the atom-positional correlation times and torsional angle transitions, however, depend on the solvent viscosity. The protein appears to be much more mobile in the solvent of lower viscosity. It feels the influence of the solvent not only on the surface but even in its core. With increasing solvent viscosity, the positional relaxation times of atoms in the protein core increase as much as those of atoms on the protein surface, and the relative increase in the core is even larger than on the surface.     

Conformational difference in cation complexes of tetranactin in solution

The conformational changes and binding behavior of tetranactin on complexation with sodium, potassium, rubidium, cesium, and ammonium ions were investigated by the measurements of proton magnetic resonance, ir, and Raman spectra. It has been clearly shown that alkali cations coordinate to the oxygen atoms of both the carbonyl group and the tetra-hydrofuran ring, but the ammonium ion coordinates only to the oxygen atom of the tetrahydrofuran. Among the alkali cations the potassium ion most strongly coordinates to the tetrahydrofuran oxygen atoms. The complexation with larger cations induces an expansion of the cavity of the macrocyclic ring of tetranactin and smaller cations contract the cavity. The evidence is revealed by the coupling constants of the methylene protons and the frequency separation between the carbonyl stretching vibrations of the ir- and Raman-active modes. The conformations of the cation complexes in the solid are maintained in solution but that of the cation free form is not.     

Quantum-dynamical picture of a multistep enzymatic process: reaction catalyzed by phospholipase A(2)

A quantum-classical molecular dynamics model (QCMD), applying explicit integration of the time-dependent Schr?dinger equation (QD) and Newtonian equations of motion (MD), is presented. The model is capable of describing quantum dynamical processes in complex biomolecular systems. It has been applied in simulations of a multistep catalytic process carried out by phospholipase A(2) in its active site. The process includes quantum-dynamical proton transfer from a water molecule to histidine localized in the active site, followed by a nucleophilic attack of the resulting OH(-) group on a carbonyl carbon atom of a phospholipid substrate, leading to cleavage of an adjacent ester bond. The process has been simulated using a parallel version of the QCMD code. The potential energy function for the active site is computed using an approximate valence bond (AVB) method. The dynamics of the key proton is described either by QD or classical MD. The coupling between the quantum proton and the classical atoms is accomplished via Hellmann-Feynman forces, as well as the time dependence of the potential energy function in the Schr?dinger equation (QCMD/AVB model). Analysis of the simulation results with an Advanced Visualization System revealed a correlated rather than a stepwise picture of the enzymatic process. It is shown that an sp(2)--> sp(3) configurational change at the substrate carbonyl carbon is mostly responsible for triggering the activation process.     

Periodic conformations of deoxyribonucleic acids

The conformation of a deoxyribonucleotide unit in a deoxyribonucleic acid molecule can be defined by six angles, each of which specifies the relative orientations between two groups of atoms adjacent to a covalent bond. With the assumption that these atoms are hard spheres with fixed van der Waal radii, conformations are sought to minimize their overlapping. The additional requirement of the polymer having a periodic structure further reduces the allowable conformations.     

Quantum chemical determination of Young's modulus of lignin. Calculations on a beta-O-4' model compound

The calculation of Young's modulus of lignin has been examined by subjecting a dimeric model compound to strain, coupled with the determination of energy and stress. The computational results, derived from quantum chemical calculations, are in agreement with available experimental results. Changes in geometry indicate that modifications in dihedral angles occur in response to linear strain. At larger levels of strain, bond rupture is evidenced by abrupt changes in energy, structure, and charge. Based on the current calculations, the bond scission may be occurring through a homolytic reaction between aliphatic carbon atoms. These results may have implications in the reactivity of lignin especially when subjected to processing methods that place large mechanical forces on the structure.     

Molecular dynamics simulation of proton transport near the surface of a phospholipid membrane

The structural and dynamical properties of a hydrated proton near the surface of DMPC membrane were studied using a molecular dynamics simulation. The proton transport between water molecules was modeled using the second generation multistate empirical valence bond model. The proton diffusion was found to be inhibited at the membrane surface. The potential of mean force for the proton adsorption to the membrane surface and its release back into the bulk water was also determined, yielding a small barrier in each direction. An efficient algorithm for Ewald summation calculations for the multistate empirical valence bond model is also introduced.     

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.