Homology‐based hydrogen bond information improves crystallographic structures in the PDB

The Protein Data Bank (PDB) is the global archive for structural information on macromolecules, and a popular resource for researchers, teachers, and students, amassing more than one million unique users each year. Crystallographic structure models in the PDB (more than 100,000 entries) are optimized against the crystal diffraction data and geometrical restraints. This process of crystallographic refinement typically ignored hydrogen bond (H‐bond) distances as a source of information. However, H‐bond restraints can improve structures at low resolution where diffraction data are limited. To improve low‐resolution structure refinement, we present methods for deriving H‐bond information either globally from well‐refined high‐resolution structures from the PDB‐REDO databank, or specifically from on‐the‐fly constructed sets of homologous high‐resolution structures. Refinement incorporating HOmology DErived Restraints (HODER), improves geometrical quality and the fit to the diffraction data for many low‐resolution structures. To make these improvements readily available to the general public, we applied our new algorithms to all crystallographic structures in the PDB: using massively parallel computing, we constructed a new instance of the PDB‐REDO databank ( https://pdb-redo.eu ). This resource is useful for researchers to gain insight on individual structures, on specific protein families (as we demonstrate with examples), and on general features of protein structure using data mining approaches on a uniformly treated dataset.     

The unusually strong hydrogen bond between the carbonyl of Q(A) and His M219 in the Rhodobacter sphaeroides reaction center is not essential for efficient electron transfer from Q(A)(-) to Q(B)

In native reaction centers (RCs) from photosynthetic purple bacteria the primary quinone (QA) and the secondary quinone (QB) are interconnected via a specific His-Fe-His bridge. In Rhodobacter sphaeroides RCs the C4=O carbonyl of QA forms a very strong hydrogen bond with the protonated Npi of His M219, and the Ntau of this residue is in turn coordinated to the non-heme iron atom. The second carbonyl of QA is engaged in a much weaker hydrogen bond with the backbone N-H of Ala M260. In previous work, a Trp side chain was introduced by site-directed mutagenesis at the M260 position in the RC of Rb. sphaeroides, resulting in a complex that is completely devoid of QA and therefore nonfunctional. A photochemically competent derivative of the AM260W mutant was isolated that contains a Cys side chain at the M260 position (denoted AM260(W-->C)). In the present work, the interactions between the carbonyl groups of QA and the protein in the AM260(W-->C) suppressor mutant have been characterized by light-induced FTIR difference spectroscopy of the photoreduction of QA. The QA-/QA difference spectrum demonstrates that the strong interaction between the C4=O carbonyl of QA and His M219 is lost in the mutant, and the coupled CO and CC modes of the QA- semiquinone are also strongly perturbed. In parallel, a band assigned to the perturbation of the C5-Ntau mode of His M219 upon QA- formation in the native RC is lacking in the spectrum of the mutant. Furthermore, a positive band between 2900 and 2400 cm-1 that is related to protons fluctuating within a network of highly polarizable hydrogen bonds in the native RC is reduced in amplitude in the mutant. On the other hand, the QB-/QB FTIR difference spectrum is essentially the same as for the native RC. The kinetics of electron transfer from QA- to QB were measured by the flash-induced absorption changes at 780 nm. Compared to native RCs the absorption transients are slowed by a factor of about 2 for both the slow phase (in the hundreds of microseconds range) and fast phase (microseconds to tens of microseconds range) in AM260(W-->C) RCs. We conclude that the unusually strong hydrogen bond between the carbonyl of QA and His M219 in the Rb. sphaeroides RC is not obligatory for efficient electron transfer from QA- to QB.     

Direct triazine herbicide detection using a self-assembled photosynthetic reaction center from purple bacterium

In this study, a direct detection system for triazine derivative herbicides was developed using the photosynthetic reaction center (RC) from the purple bacterium,Rhodobacter sphaeroides, and surface plasmon resonance (SPR) apparatus. The histidine-tagged RCs were immobilized on an SPR gold chip using nickel-nitrilotriacetic acid groups as a binder for one of the triazine herbicide, atrazine. The SPR responses were proportional to the sample concentrations of atrazine in the range 0.1–1 μg/mL. The sensitivity of the direct detection of atrazine using the RC-assembled sensor chip was higher than that using the antibody-immobilized chip. The other types of herbicides, DCMU or MCPP, were not detected with such high sensitivity. The results indicated the high binding selectivity of the RC complex.     

N-H...O, O-H...O, and C-H...O hydrogen bonds in protein-ligand complexes: strong and weak interactions in molecular recognition

The characteristics of N-H...O, O-H...O, and C-H...O hydrogen bonds are examined in a group of 28 high-resolution crystal structures of protein-ligand complexes from the Protein Data Bank and compared with interactions found in small-molecule crystal structures from the Cambridge Structural Database. It is found that both strong and weak hydrogen bonds are involved in ligand binding. Because of the prevalence of multifurcation, the restrictive geometrical criteria set up for hydrogen bonds in small-molecule crystal structures may need to be relaxed in macromolecular structures. For example, there are definite deviations from linearity for the hydrogen bonds in protein-ligand complexes. The formation of C-H...O hydrogen bonds is influenced by the activation of the C(alpha)-H atoms and by the flexibility of the side-chain atoms. In contrast to small-molecule structures, anticooperative geometries are common in the macromolecular structures studied here, and there is a gradual lengthening as the extent of furcation increases. C-H...O bonds formed by Gly, Phe, and Tyr residues are noteworthy. The numbers of hydrogen bond donors and acceptors agree with Lipinski's "rule of five" that predicts drug-like properties. Hydrogen bonds formed by water are also seen to be relevant in ligand binding. Ligand C-H...O(w) interactions are abundant when compared to N-H...O(w) and O-H...O(w). This suggests that ligands prefer to use their stronger hydrogen bond capabilities for use with the protein residues, leaving the weaker interactions to bind with water. In summary, the interplay between strong and weak interactions in ligand binding possibly leads to a satisfactory enthalpy-entropy balance. The implications of these results to crystallographic refinement and molecular dynamics software are discussed.     

HBAT: a complete package for analysing strong and weak hydrogen bonds in macromolecular crystal structures

The program HBAT is a tool to automate the analysis of potential hydrogen bonds and similar type of weak interactions like halogen bonds and non-canonical interactions in macromolecular structures, available in Brookhaven Protein Database (PDB) file format. HBAT is written using PERL and TK languages. The program generates an MSOFFICE Excel compatible output file for statistical analysis. HBAT identify potential interactions based on geometrical criteria. A series of analysis reports like frequency tables, geometry distribution tables, furcations list are generated. A user friendly GUI offers freedom to select several parameters and options. Graphviz based visualization of hydrogen bond networks in 2D helps to study the cooperativity and anticooperativity geometry in hydrogen bond. HBAT supports post docking interaction analysis between PDB files for any target/receptor (in PDB files) and docked ligands/poses (in SDF). This tool can be implemented in active site interaction analysis, structure based drug design and molecular dynamics simulations.     

Structural basis of the drastically increased initial electron transfer rate in the reaction center from a Rhodopseudomonas viridis mutant described at 2.00-A resolution

It has previously been shown that replacement of the residue His L168 with Phe (HL168F) in the Rhodopseudomonas viridis reaction center (RC) leads to an unprecedented drastic acceleration of the initial electron transfer rate. Here we describe the determination of the x-ray crystal structure at 2.00-A resolution of the HL168F RC. The electron density maps confirm that a hydrogen bond from the protein to the special pair is removed by this mutation. Compared with the wild-type RC, the acceptor of this hydrogen bond, the ring I acetyl group of the "special pair" bacteriochlorophyll, D(L), is rotated, and its acetyl oxygen is found 1.1 A closer to the bacteriochlorophyll-Mg(2+) of the other special pair bacteriochlorophyll, D(M). The rotation of this acetyl group and the increased interaction between the D(L) ring I acetyl oxygen and the D(M)-Mg(2+) provide the structural basis for the previously observed 80-mV decrease in the D(+)/D redox potential and the drastically increased rate of initial electron transfer to the accessory bacteriochlorophyll, B(A). The high quality of the electron density maps also allowed a reliable discussion of the mode of binding of the triazine herbicide terbutryn at the binding site of the secondary quinone, Q(B).     

Protein regulation of carotenoid binding; gatekeeper and locking amino acid residues in reaction centers of Rhodobacter sphaeroides

X-ray diffraction was used to determine high-resolution structures of the reaction center (RC) complex from the carotenoidless mutant, Rb. sphaeroides R-26.1, without or reconstituted with carotenoids. The results are compared with the structure of the RC from a semiaerobically grown Rb. sphaeroides strain 2.4.1. The investigation reveals the structure of the carotenoid in the different protein preparations, the nature of its binding site, and a plausible mechanism by which the carotenoid is incorporated unidirectionally in its characteristic geometric configuration. The structural data suggest that the accessibility of the carotenoid to the binding site is controlled by a specific "gatekeeper" residue which allows the carotenoid to approach the binding site from only one direction. Carotenoid binding to the protein is secured by hydrogen bonding to a separate "locking" amino acid. The study reveals the specific molecular interactions that control how the carotenoid protects the photosynthetic apparatus against photo-induced oxidative destruction.     

Metabolism of the herbicide atrazine by Rhodococcus strains.

Rhodococcus strains were screened for their ability to degrade the herbicide atrazine. Only rhodococci that degrade the herbicide EPTC (s-ethyl-dipropylthiocarbamate) metabolized atrazine. Rhodococcus strain TE1 metabolized atrazine under aerobic conditions to produce deethyl- and deisopropylatrazine, which were not degraded further and which accumulated in the incubation medium. The bacterium also metabolized the other s-triazine herbicides propazine, simazine, and cyanazine. The N dealkylation of triazine herbicides by Rhodococcus strain TE1 was associated with a 77-kb plasmid previously shown to be required for EPTC degradation.     

Discovering protein−ligand chalcogen bonding in the protein data bank using endocyclic sulfur-containing heterocycles as ligand search subsets

The chalcogen bond, the noncovalent, electrostatic attraction between covalently bonded atoms in group 16 and Lewis bases, is present in protein?ligand interactions based on X-ray structures deposited in the Protein Data Bank (PDB). Discovering protein?ligand chalcogen bonding in the PDB employed a strategy that focused on searching the database for protein complexes of five-membered, heterocyclic ligands containing endocyclic sulfur with endo electron-withdrawing groups (isothiazoles; thiazoles; 1,2,3-, 1,2.4-, 1,2,5-, 1,3,4-thiadiazoles) and thiophenes with exo electron-withdrawing groups, e.g., 2-chloro, 2-bromo, 2-amino, 2-alkylthio. Out of 930 ligands investigated, 33 or 3.5% have protein?ligand S---O interactions of which 31 are chalcogen bonds and two appear to be S---HO hydrogen bonds. The bond angles for some of the chalcogen bonds found in the PDB are less than 90°, and an electrostatic model is proposed to explain this phenomenon.     

Oxidative stress response in two representative bacteria exposed to atrazine

Bacteria are present extensively in the environment. Investigation of their antioxidant properties will be useful for further study on atrazine stress tolerance of bacteria and the defense mechanism of antioxidant enzymes against atrazine or other triazine herbicides. Superoxide dismutase (SOD), catalase (CAT), glutathione S-transferase (GST) and total antioxidant capacity (T-AOC) from one Gram-negative representative strain Escherichia coli K12 and one Gram-positive representative strain Bacillus subtilis B19, respectively, were tested for response to atrazine stress. The results indicated that SOD, CAT, GST and T-AOC were induced upon exposure to atrazine. The growth of two bacteria was better in the absence than in the presence of atrazine, indicating that atrazine can decrease bacterial growth. The changes of enzyme activities indicate the presence of oxidative stress. Oxidative stress induced by atrazine may be due to imbalance of redox potential in bacterial cells, which leads to bacterial metabolic disorder.     

Temporary Stabilization of Electron on Quinone Acceptor Side of Reaction Centers from the Bacterium Rhodobacter sphaeroides Wild Type and Mutant SA(L223) Depending on Duration of Light Activation

The dark reduction of photooxidized bacteriochlorophyll (P+) by photoreduced secondary quinone acceptor (QB-) in isolated reaction centers (RC) from the bacterium Rhodobacter sphaeroides wild type and mutant strain SA(L223) depending on the duration of light activation of RC was studied. The kinetics of the dark reduction of P+ decreased with increasing light duration, which is probably due to conformational changes occurring under prolonged light activation in RC from the wild type bacterium. In RC from bacteria of the mutant strain in which protonatable amino acid Ser L223 near QB is substituted by Ala, the dependence of reduction kinetics of P+ on duration of light was not observed. Such dependence, however, became observable after addition of cryoprotectors, namely glycerol and dimethylsulfoxide, to the RC samples from the mutant strain. It was concluded that substitution of Ser L223 with Ala disturbs the native mechanism of electrostatic stabilization of the electron in the RC quinone acceptor site. At the same time, an additional modification of RC hydrogen bonds by glycerol and dimethylsulfoxide probably includes various possibilities for more effective time delay of the electron on QB.     

Crystal structure of glycerophosphodiester phosphodiesterase (GDPD) from Thermoanaerobacter tengcongensis, a metal ion-dependent enzyme: insight into the catalytic mechanism

Glycerophosphodiester phosphodiesterase (GDPD; EC 3.1.4.46) catalyzes the hydrolysis of a glycerophosphodiester to an alcohol and glycerol 3-phosphate in glycerol metabolism. It has an important role in the synthesis of a variety of products that participate in many biochemical pathways. We report the crystal structure of the Thermoanaerobacter tengcongensis GDPD (ttGDPD) at 1.91 A resolution, with a calcium ion and glycerol as a substrate mimic coordinated at this calcium ion (PDB entry 2pz0). The ttGDPD dimer with an intermolecular disulfide bridge and two hydrogen bonds is considered as the potential functional unit. We used site-directed mutagenesis to characterize ttGDPD as a metal ion-dependent enzyme, identified a cluster of residues involved in substrate binding and the catalytic reaction, and we propose a possible general acid-base catalytic mechanism for ttGDPD. Superposing the active site with the homologous structure GDPD from Agrobacterium tumefaciens (PDB entry 1zcc), which binds a sulfate ion in the active site, the sulfate ion can represent the phosphate moiety of the substrate, simulating the binding mode of the true substrate of GDPD.     

Rapid and specific detection of herbicides using a self-assembled photosynthetic reaction center from purple bacterium on an SPR chip

In this study, a direct detection system for herbicides inhibiting photosynthetic electron transfer was developed using the photosynthetic reaction center (RC) from the purple bacterium, Rhodobacter sphaeroides, and surface plasmon resonance (SPR) apparatus. The heavy-subunit-histidine-tagged RCs (HHisRCs) were immobilized on an SPR sensor chip via nickel chelation chemistry as a binder for one of the triazine herbicides, atrazine. Immediately after injection of atrazine solution on the HHisRCs-immobilized chip, the SPR responses increased and reached plateaus within 1 min. The SPR signals were proportional to the sample concentrations of atrazine in the range 1-100 microg/ml. To evaluate the binding specificity to atrazine, chlorinated aromatic herbicides, DCMU and MCPP, were investigated using the HHisRCs-immobilized chip. An RC inhibitor, DCMU, could also be detected with a higher detection limit of 20 microg/ml than atrazine (1 microg/ml). MCPP showed no signals because its inhibition mechanism against plants is different from that of atrazine and DCMU. These results indicated that the sensor chip immobilized RCs could be used for the specific detection of photosynthetic inhibitors.     

Exploring the primary electron acceptor (QA)-site of the bacterial reaction center from Rhodobacter sphaeroides. Binding mode of vitamin K derivatives.

The functional replacement of the primary ubiquinone (QA) in the photosynthetic reaction center (RC) from Rhodobacter sphaeroides with synthetic vitamin K derivatives has provided a powerful tool to investigate the electron transfer mechanism. To investigate the binding mode of these quinones to the QA binding site we have determined the binding free energy and charge recombination rate from QA(-) to D+ (kAD) of 29 different 1,4-naphthoquinone derivatives with systematically altered structures. The most striking result was that none of the eight tested compounds carrying methyl groups in both positions 5 and 8 of the aromatic ring exhibited functional binding. To understand the binding properties of these quinones on a molecular level, the structures of the reaction center-naphthoquinone complexes were predicted with ligand docking calculations. All protein--ligand structures show hydrogen bonds between the carbonyl oxygens of the quinone and AlaM260 and HisM219 as found for the native ubiquinone-10 in the X-ray structure. The center-to-center distance between the naphthoquinones at QA and the native ubiquinone-10 at QB (the secondary electron acceptor) is essentially the same, compared to the native structure. A detailed analysis of the docking calculations reveals that 5,8-disubstitution prohibits binding due to steric clashes of the 5-methyl group with the backbone atoms of AlaM260 and AlaM249. The experimentally determined binding free energies were reproduced with an rmsd of approximately 4 kJ x mol(-1) in most cases providing a valuable tool for the design of new artificial electron acceptors and inhibitors.     

atz gene expressions during atrazine degradation in the soil drilosphere

One of the various ecosystemic services sustained by soil is pollutant degradation mediated by adapted soil bacteria. The pathways of atrazine biodegradation have been elucidated but in situ expression of the genes involved in atrazine degradation has yet to be demonstrated in soil. Expression of the atzA and atzD genes involved in atrazine dechlorination and s‐triazine ring cleavage, respectively, was investigated during in situ degradation of atrazine in the soil drilosphere and bulked samples from two agricultural soils that differed in their ability to mineralize atrazine. Interestingly, expression of the atzA gene, although present in both soils, was not detected. Atrazine mineralization was greatest in Epoisses soil, where a larger pool of atzD mRNA was consistently measured 7 days after atrazine treatment, compared with Vezin soil (146 vs. 49 mRNA per 10 ⁶ 16S rRNA, respectively). Expression of the atzD gene varied along the degradation time course and was profoundly modified in soil bioturbated by earthworms. The atzD mRNA pool was the highest in the soil drilosphere (casts and burrow‐linings) and it was significantly different in burrow‐linings compared with bulk soil (e.g. 363 vs. 146 mRNA per 10 ⁶ 16S rRNA, 7 days after atrazine treatment in Epoisses soil). Thus, consistent differences in atrazine mineralization were demonstrated between the soil drilosphere and bulk soil. However, the impact of bioturbation on atrazine mineralization depended on soil type. Mineralization was enhanced in casts, compared with bulk soil, from Epoisses soil but in burrow‐linings from Vezin soil. This study is the first to report the effects of soil bioturbation by earthworms on s‐triazine ring cleavage and its spatial variability in soil.     

Influence of carotenoid molecules on the structure of the bacteriochlorophyll binding site in peripheral light-harvesting proteins from Rhodobacter sphaeroides

In the LH2 proteins from Rhodobacter (Rb.) sphaeroides, the hydrogen bonds between the bacteriochlorophyll (Bchl) molecules and their proteic binding sites exhibit a strong variance with respect to carotenoid content and type. In the absence of the carotenoid molecule, such as in the LH2 from Rb. sphaeroides R26.1, the void in the protein structure induces a significant reorganization of the binding site of both Bchl molecules responsible for the 850 nm absorption, which is not observed when the 800 nm absorbing Bchl is selectively removed from these complexes. FT Raman spectra of LH2 complexes from Rb. sphaeroides show that the strength of the hydrogen bond between the 850 nm absorbing Bchl bound to the alpha polypeptide and the tyrosine alpha(45) depends precisely on the chemical nature of the bound carotenoid. These results suggest that the variable extremity of the carotenoid is embedded in these LH2 complexes, lying close to the interacting Bchl molecules. In the LH2 from Rhodopseudomonas acidophila, the equivalent part of the rhodopin glucoside, which bears the glucose group, lies close to the amino terminal of the antenna polypeptide. This contrast suggests that the structure of the carotenoid binding site in LH2 complexes strongly depends on the bacterial species and/or on the chemical nature of the bound carotenoid.     

On the molecular discrimination between adenine and guanine by proteins

The molecular recognition and discrimination of adenine and guanine ligand moieties in complexes with proteins have been studied using empirical observations on carefully selected crystal structures. The distribution of protein folds that bind these purines has been found to differ significantly from that across the whole PDB, but the most populated architectures and folds are also the most common in three genomes from the three different domains of life. The protein environments around the two nucleic acid bases were significantly different, in terms of the propensities of amino acid residues to be in the binding site, as well as their propensities to form hydrogen bonds to the bases. Plots of the distribution of protein atoms around the two purines clearly show different clustering of hydrogen bond donors and acceptors opposite complimentary acceptors and donors in the rings, with hydrophobic areas below and above the rings. However, the clustering pattern is fuzzy, reflecting the variety of ways that proteins have evolved to recognise the same molecular moiety. Furthermore, an analysis of the conservation of residues in the protein chains binding guanine shows that residues in contact with the base are in general better conserved than the rest of the chain.     

Refinement of a low-resolution crystal structure to better understand erythromycin interactions on large ribosomal subunit

Macrolides are a group of diverse class of naturally occurring and synthetic antibiotics made of macrocyclic-lactone ring carrying one or more sugar moieties linked to various atoms of the lactone ring. These macrolides selectively bind to a single high affinity site on the prokaryotic 50S ribosomal subunit, making them highly effective towards a wide range of bacterial pathogens. The understanding of binding between macrolides and ribosome serves a good basis in elucidating how they work at the molecular level and these findings would be important in rational drug design. Here, we report refinement of reconstructed PDB structure of erythromycin-ribosome system using molecular dynamics (MD) simulation. Interesting findings were observed in this refinement stage that could improve the understanding of the binding of erythromycin A (ERYA) onto the 50S subunit. The results showed ERYA was highly hydrated and water molecules were found to be important in bridging hydrogen bond at the binding pocket during the simulation time. ERYA binding to ribosome was also strengthened by hydrogen bond network and hydrophobic interactions between the antibiotic and the ribosome. Our MD simulation also demonstrated direct interaction of ERYA with Domains II, V and with C1773 (U1782EC), a residue in Domain IV that has yet been described of its role in ERYA binding. It is hoped that this refinement will serve as a starting model for a further enhancement of our understanding towards the binding of ERYA to ribosome.     

Allosteric Regulation of 3CL Protease of SARS-CoV-2 and SARS-CoV Observed in the Crystal Structure Ensemble

The 3C-like protease (3CL^pro) of SARS-CoV-2 is a potential therapeutic target for COVID-19. Importantly, it has an abundance of structural information solved as a complex with various drug candidate compounds. Collecting these crystal structures (83 Protein Data Bank (PDB) entries) together with those of the highly homologous 3CL^pro of SARS-CoV (101 PDB entries), we constructed the crystal structure ensemble of 3CL^pro to analyze the dynamic regulation of its catalytic function. The structural dynamics of the 3CL^pro dimer observed in the ensemble were characterized by the motions of four separate loops (the C-loop, E-loop, H-loop, and Linker) and the C-terminal domain III on the rigid core of the chymotrypsin fold. Among the four moving loops, the C-loop (also known as the oxyanion binding loop) causes the order (active)–disorder (collapsed) transition, which is regulated cooperatively by five hydrogen bonds made with the surrounding residues. The C-loop, E-loop, and Linker constitute the major ligand binding sites, which consist of a limited variety of binding residues including the substrate binding subsites. Ligand binding causes a ligand size dependent conformational change to the E-loop and Linker, which further stabilize the C-loop via the hydrogen bond between the C-loop and E-loop. The T285A mutation from SARS-CoV 3CL^pro to SARS-CoV-2 3CL^pro significantly closes the interface of the domain III dimer and allosterically stabilizes the active conformation of the C-loop via hydrogen bonds with Ser1 and Gly2; thus, SARS-CoV-2 3CL^pro seems to have increased activity relative to that of SARS-CoV 3CL^pro.