Geometry of nonbonded interactions involving planar groups in proteins

Mutual influence of non-covalent interactions formed by imidazole: A systematic quantum-chemical study

Imidazole is a five-membered heterocycle that is part of a number of biologically important molecules such as the amino acid histidine and the hormone histamine. Imidazole has a unique ability to participate in a variety of non-covalent interactions involving the NH group, the pyridine-like nitrogen atom or the π-system. For many biologically active compounds containing the imidazole moiety, its participation in formation of hydrogen bond NH⋯O/N and following proton transfer is the key step of mechanism of their action. In this work a systematic study of the mutual influence of various paired combinations of non-covalent interactions (e.g., hydrogen bonds and π-interactions) involving the imidazole moiety was performed by means of quantum chemistry (PW6B95-GD3/def2-QZVPD) for a series of model systems constructed based on analysis of available x-ray data. It is shown that for considered complexes formation of additional non-covalent interactions can only enhance the proton-donating ability of imidazole. At the same time, its proton-accepting ability can be both enhanced and weakened, depending on what additional interactions are added to a given system. The mutual influence of non-covalent interactions involving imidazole can be classified as weak geometric and strong energetic cooperativity-a small change in the length of non-covalent interaction formed by imidazole can strongly influence its strength. The latter can be used to develop methods for controlling the rate and selectivity of chemical reactions involving the imidazole fragment in larger systems. It is shown that the strong mutual influence of non-covalent interactions involving imidazole is due to the unique ability of the imidazole ring to effectively redistribute electron density in non-covalently bound systems with its participation.

Understanding a point mutation signature D54K in the caspase activation recruitment domain of NOD1 capitulating concerted immunity via atomistic simulation

Point mutation D54K in the human N-terminal caspase recruitment domain (CARD) of nucleotide-binding oligomerization domain -1 (NOD1) abrogates an imperative downstream interaction with receptor-interacting protein kinase (RIPK2) that entails combating bacterial infections and inflammatory dysfunction. Here, we addressed the molecular details concerning conformational changes and interaction patterns (monomeric-dimeric states) of D54K by signature-based molecular dynamics simulation. Initially, the sequence analysis prioritized D54K as a pathogenic mutation, among other variants, based on a sequence signature. Since the mutation is highly conserved, we derived the distant ortholog to predict the sequence and structural similarity between native and mutant. This analysis showed the utility of 33 communal core residues associated with structural-functional preservation and variations, concurrently served to infer the cryptic hotspots Cys39, Glu53, Asp54, Glu56, Ile57, Leu74, and Lys78 determining the inter helical fold forming homodimers for putative receptor interaction. Subsequently, the atomistic simulations with free energy (MM/PB(GB)SA) calculations predicted structural alteration that takes place in the N-terminal mutant CARD where coils changed to helices (45 α3- L4-α4-L6- α683) in contrast to native (45T2-L4-α4-L6-T483). Likewise, the C-terminal helices 93T1-α7105 connected to the loops distorted compared to native 93α6-L7105 may result in conformational misfolding that promotes functional regulation and activation. These structural perturbations of D54K possibly destabilize the flexible adaptation of critical homotypic NOD1CARD-CARDRIPK2 interactions (α4Asp42-Arg488α5 and α6Phe86-Lys471α4) is consistent with earlier experimental reports. Altogether, our findings unveil the conformational plasticity of mutation-dependent immunomodulatory response and may aid in functional validation exploring clinical investigation on CARD-regulated immunotherapies to prevent systemic infection and inflammation.Communicated by Ramaswamy H. Sarma.

Molecular characterization of the PhiKo endolysin from Thermus thermophilus HB27 bacteriophage phiKo and its cryptic lytic peptide RAP-29

Introduction

In the era of increasing bacterial resistance to antibiotics, new bactericidal substances are sought, and lysins derived from extremophilic organisms have the undoubted advantage of being stable under harsh environmental conditions. The PhiKo endolysin is derived from the phiKo bacteriophage infecting Gram-negative extremophilic bacterium Thermus thermophilus HB27. This enzyme shows similarity to two previously investigated thermostable type-2 amidases, the Ts2631 and Ph2119 from Thermus scotoductus bacteriophages, that revealed high lytic activity not only against thermophiles but also against Gram-negative mesophilic bacteria. Therefore, antibacterial potential of the PhiKo endolysin was investigated in the study presented here.

Methods

Enzyme activity was assessed using turbidity reduction assays (TRAs) and antibacterial tests. Differential scanning calorimetry was applied to evaluate protein stability. The Collection of Anti-Microbial Peptides (CAMP) and Antimicrobial Peptide Calculator and Predictor (APD3) were used to predict regions with antimicrobial potential in the PhiKo primary sequence. The minimum inhibitory concentration (MIC) of the RAP-29 synthetic peptide was determined against Gram-positive and Gram-negative selected strains, and mechanism of action was investigated with use of membrane potential sensitive fluorescent dye 3,3'-Dipropylthiacarbocyanine iodide (DiSC3(5)).

Results and discussion

The PhiKo endolysin is highly thermostable with melting temperature of 91.70°C. However, despite its lytic effect against such extremophiles as: T. thermophilus, Thermus flavus, Thermus parvatiensis, Thermus scotoductus, and Deinococcus radiodurans, PhiKo showed moderate antibacterial activity against mesophiles. Consequently, its protein sequence was searched for regions with potential antibacterial activity. A highly positively charged region was identified and synthetized (PhiKo105-133). The novel RAP-29 peptide lysed mesophilic strains of staphylococci and Gram-negative bacteria, reducing the number of cells by 3.7-7.1 log units and reaching the minimum inhibitory concentration values in the range of 2-31 μM. This peptide is unstructured in an aqueous solution but forms an α-helix in the presence of detergents. Moreover, it binds lipoteichoic acid and lipopolysaccharide, and causes depolarization of bacterial membranes. The RAP-29 peptide is a promising candidate for combating bacterial pathogens. The existence of this cryptic peptide testifies to a much wider panel of antimicrobial peptides than thought previously.

An apical Phe-His pair defines the Orai1-coupling site and its occlusion within STIM1

Ca2+ signal-generation through inter-membrane junctional coupling between endoplasmic reticulum (ER) STIM proteins and plasma membrane (PM) Orai channels, remains a vital but undefined mechanism. We identify two unusual overlapping Phe-His aromatic pairs within the STIM1 apical helix, one of which (F394-H398) mediates important control over Orai1-STIM1 coupling. In resting STIM1, this locus is deeply clamped within the folded STIM1-CC1 helices, likely near to the ER surface. The clamped environment in holo-STIM1 is critical-positive charge replacing Phe-394 constitutively unclamps STIM1, mimicking store-depletion, negative charge irreversibly locks the clamped-state. In store-activated, unclamped STIM1, Phe-394 mediates binding to the Orai1 channel, but His-398 is indispensable for transducing STIM1-binding into Orai1 channel-gating, and is spatially aligned with Phe-394 in the exposed Sα2 helical apex. Thus, the Phe-His locus traverses between ER and PM surfaces and is decisive in the two critical STIM1 functions-unclamping to activate STIM1, and conformational-coupling to gate the Orai1 channel.

Creating an Amyloid 'Kaleidoscope' Using Short Iodinated Peptides

Amyloid fibrils formed by peptides with different sequences exhibit diversified morphologies, material properties and activities, making them valuable for developing functional bionanomaterials. However, the molecular understanding underlying the structural diversity of peptide fibrillar assembly at atomic level is still lacking. In this study, by using cryogenic electron microscopy, we first revealed the structural basis underlying the highly reversible assembly of 1 GFGGNDNFG9 (referred to as hnRAC1) peptide fibril. Furthermore, by installing iodine at different sites of hnRAC1, we generated a collection of peptide fibrils with distinct thermostability. By determining the atomic structures of the iodinated fibrils, we discovered that iodination at different sites of the peptide facilitates the formation of diverse halogen bonds and triggers the assembly of entirely different structures of iodinated fibrils. Finally, based on this structural knowledge, we designed an iodinated peptide that assembles into new atomic structures of fibrils, exhibiting superior thermostability, that aligned with our design. Our work provides an in-depth understanding of the atomic-level processes underlying the formation of diverse peptide fibril structures, and paves the way for creating an amyloid "kaleidoscope" by employing various modifications and peptide sequences to fine-tune the atomic structure and properties of fibrillar nanostructures.

Categorization of hotspots into three types - weak, moderate and strong to distinguish protein-protein versus protein-peptide interactions

Protein-protein and protein-peptide interactions (PPI and PPepI) belong to a similar category of interactions, yet seemingly subtle differences exist among them. To characterize differences between protein-protein (PP) and protein-peptide (PPep) interactions, we have focussed on two important classes of residues-hotspot and anchor residues. Using implicit solvation-based free energy calculations, a very large-scale alanine scanning has been performed on benchmark datasets, consisting of over 5700 interface residues. The differences in the two categories are more pronounced, if the data were divided into three distinct types, namely - weak hotspots (having binding free energy loss upon Ala mutation, ΔΔG, ∼2-10 kcal/mol), moderate hotspots (ΔΔG, ∼10-20 kcal/mol) and strong hotspots (ΔΔG ≥ ∼20 kcal/mol). The analysis suggests that for PPI, weak hotspots are predominantly populated by polar and hydrophobic residues. The distribution shifts towards charged and polar residues for moderate hotspot and charged residues (principally Arg) are overwhelmingly present in the strong hotspot. On the other hand, in the PPepI dataset, the distribution shifts from predominantly hydrophobic and polar (in the weak type) to almost similar preference for polar, hydrophobic and charged residues (in moderate type) and finally the charged residue (Arg) and Trp are mostly occupied in the strong type. The preferred anchor residues in both categories are Arg, Tyr and Leu, possessing bulky side chain and which also strike a delicate balance between side chain flexibility and rigidity. The present knowledge should aid in effective design of biologics, by augmentation or disruption of PPIs with peptides or peptidomimetics.Communicated by Ramaswamy H. Sarma.

Noncovalent π Interactions in Mutated Aquomet-Myoglobin Proteins: A QM/MM and Local Vibrational Mode Study

Protein dynamics and function is strongly connected to the energy flow taking place. Myoglobin (Mb) and its mutations are ideal systems to study the process of vibrational energy transfer (VET) at the molecular level. Anti-Stokes ultraviolet resonance Raman studies using a tryptophan (Trp) probe, introduced at different Mb positions by amino acid replacement, have suggested that the amount of VET depends on the position of the Trp probe relative to the heme group. Inspired by this experimental work, we explored the strength of noncovalent π interactions, as well as covalent interactions for both the axial and distal ligands bound to iron in aquomet-Mb with the local vibrational mode analysis (LMA), originally developed by Konkoli and Cremer. Two sets of noncovalent interactions were investigated: (1) the interaction between the water ligand and Trp rings and (2) the interaction between the Trp and the porphyrin rings of the heme group. We assessed the strength of these noncovalent interactions via a special local mode force constant. Various Trp-modified water-bound ferric Mb proteins in the ground state were studied (6 in total) using gas-phase and QM/MM calculations followed by LMA. Our results disclose that VET is indeed dependent on the position of the Trp probe relative to the heme group but also on the tautomeric nature of distal histidine. They provide new guidelines on how to assess noncovalent π interactions in proteins utilizing LMA and how to use these data to explore VET, and more generally protein dynamics and function.

Structures and Supramolecular Properties of Inclusion Complexes of Anthracene-Triptycene Nanocages with Fullerene Guests and Their Dynamic Motion as Molecular Gyroscopes

Three derivatives of macrocyclic cage compounds consisting of diarylanthracene and triptycene units were synthesized. These nanocages formed host-guest complexes with C₆₀ and other fullerene guests as confirmed by ¹ H NMR and fluorescence spectroscopy. The association constant of the mesityl and 2,4,6-tributoxyphenyl derivatives with C₆₀ was determined to be 2.2 × 10⁴  L mol^-1 , which was larger than that of the pentafluorophenyl derivative. Direct experimental evidence of the complexation was obtained by X-ray diffraction analysis: the guest C₆₀ molecule was included in the cavity via multipoint CH⋅⋅⋅π interactions. Dynamic disorders of the included C₆₀ molecule in variable-temperature X-ray analysis indicated uniaxial motion, such as gyroscopic motion. The unique dynamic behavior of the spherical C₆₀ rotor anchored by the cage stator via CH⋅⋅⋅π interactions in the crystal, as well as substituent effects on the association properties, are discussed with the aid of DFT calculations.

Engineering Silk Protein to Modulate Polymorphic Transitions for Green Lithography Resists

Silk protein is being increasingly introduced as a prospective material for biomedical devices. However, a limited locus to intervene in nature-oriented silk protein makes it challenging to implement on-demand functions to silk. Here, we report how polymorphic transitions are related with molecular structures of artificially synthesized silk protein and design principles to construct a green-lithographic and high-performative protein resist. The repetition number and ratio of two major building blocks in synthesized silk protein are essential to determine the size and content of β-sheet crystallites, and radicals resulting from tyrosine cleavages by the 193 nm laser irradiation induce the β-sheet to α-helix transition. Synthesized silk is designed to exclusively comprise homogeneous building blocks and exhibit high crystallization and tyrosine-richness, thus constituting an excellent basis for developing a high-performance deep-UV photoresist. Additionally, our findings can be conjugated to design an electron-beam resist governed by the different irradiation-protein interaction mechanisms. All synthesis and lithography processes are fully water-based, promising green lithography. Using the engineered silk, a nanopatterned planar color filter showing the reduced angle dependence can be obtained. Our study provides insights into the industrial scale production of silk protein with on-demand functions.

Molecular duplexes featuring NH···N, CH···O and CH···π interactions in solid-state self-assembly of triazine-based compounds

Synthetic supramolecular structures constructed through the cooperative action of numerous non-covalent forces are highly desirable as models to unravel and understand the complexity of systems created in nature via self-assembly. Taking advantage of the low cost of 2,4,6-trichloro-1,3,5-triazine (cyanuric chloride) and the sequential nucleophilic substitution reactions with almost all types of nucleophiles, a series of six structurally related novel s-triazine derivatives 1-6 were synthesized and structurally characterized based on their physical, spectral and crystallographic data. The solid-state structures of all the six compounds showed intriguing and unique molecular duplexes featuring NH···N, CH···O and CH···π interactions. Careful analysis of different geometric parameters of the involved H-bonds indicates that they are linear, significant and are therefore responsible for guiding the three-dimensional structure of these compounds in the solid state. The prevalence of sextuple hydrogen bond array-driven molecular duplexes and the possibility of structural modifications on the s-triazine ring render these novel triazine derivatives 1-6 attractive as a platform to create heteroduplex constructs and their subsequent utility in the field of supramolecular chemistry and crystal engineering.

Mining anion-aromatic interactions in the Protein Data Bank

Mutual positioning and non-covalent interactions in anion–aromatic motifs are crucial for functional performance of biological systems. In this context, regular, comprehensive Protein Data Bank (PDB) screening that involves various scientific points of view and individual critical analysis is of utmost importance. Analysis of anions in spheres with radii of 5 Å around all 5- and 6-membered aromatic rings allowed us to distinguish 555 259 unique anion–aromatic motifs, including 92 660 structures out of the 171 588 structural files in the PDB. The use of a scarcely exploited (x, h) coordinate system led to (i) identification of three separate areas of motif accumulation: A – over the ring, B – over the ring-substituent bonds, and C – roughly in the plane of the aromatic ring, and (ii) unprecedented simultaneous comparative description of various anion–aromatic motifs located in these areas. Of the various residues considered, i.e. aminoacids, nucleotides, and ligands, the latter two exhibited a considerable tendency to locate in region Avia archetypal anion–π contacts. The applied model not only enabled statistical quantitative analysis of space around the ring, but also enabled discussion of local intermolecular arrangements, as well as detailed sequence and secondary structure analysis, e.g. anion–π interactions in the GNRA tetraloop in RNA and protein helical structures. As a purely practical issue of this work, the new code source for the PDB research was produced, tested and made freely available at https://github.com/chemiczny/PDB_supramolecular_search.

The comprehensive analysis of non-redundant PDB macromolecular structures investigating anion distributions around all aromatic molecules in available biosystems is presented.

Interlayer Interactions as Design Tool for Large-Pore COFs

Covalent organic frameworks (COFs) with a pore size beyond 5 nm are still rarely seen in this emerging field. Besides obvious complications such as the elaborated synthesis of large linkers with sufficient solubility, more subtle challenges regarding large-pore COF synthesis, including pore occlusion and collapse, prevail. Here we present two isoreticular series of large-pore imine COFs with pore sizes up to 5.8 nm and correlate the interlayer interactions with the structure and thermal behavior of the COFs. By adjusting interlayer interactions through the incorporation of methoxy groups acting as pore-directing “anchors”, different stacking modes can be accessed, resulting in modified stacking polytypes and, hence, effective pore sizes. A strong correlation between stacking energy toward highly ordered, nearly eclipsed structures, higher structural integrity during thermal stress, and a novel, thermally induced phase transition of stacking modes in COFs was found, which sheds light on viable design strategies for increased structural control and stability in large-pore COFs.

Controlling the self-assembly of human calcitonin: a theoretical approach using molecular dynamics simulations

Human calcitonin (hCT) is a 32-residue amino acid poly-peptide hormone which is secreted by the C-cells (also known as parafollicular cells) of thyroid glands. It acts to inhibit osteoclast cell hormones by reducing the cell function and regulating calcium and phosphate in blood. hCT has a high tendency to assemble into protofilaments with β-sheet conformations. Amyloid fibril formation of hCT reduces its bio-activity and limits its application as a therapeutic drug. Salmon calcitonin (sCT), which also carries the same disulfide bridge at the N and C-terminus, but differs at the 16 residue position from hCT, has less propensity to aggregate than hCT. Human calcitonin has much higher bio-activity than sCT if its aggregation propensity is reduced. Substituting the key residues which are responsible for the aggregation of hCT, is one of the ways to reduce its aggregation and fibril formation. hCT analogues with less aggregation tendency can be exploited as therapeutic drugs. In this work, we study the amyloidogenic behavior of hCT and its peptide based derivatives i.e., sCT, phCT, N17H hCT, Y12L hCT and DM hCT, through classical molecular dynamics (MD) simulations. Our study reveals that sCT is the least aggregation prone derivative, and the double mutation at position 12 and 17 can reduce the aggregation propensity of this peptide. Also, we have applied these mutant variants of hCT as peptide inhibitors in the self-aggregation of hCT. This study could help in understanding and preparing peptide-based inhibitors for hCT fibrillation and their applications as therapeutic drugs.

An "up" oriented methionine-aromatic structural motif in SUMO is critical for its stability and activity

Protein structural bioinformatic analyses suggest preferential associations between methionine and aromatic amino acid residues in proteins. Ab initio energy calculations highlight a conformation-dependent stabilizing interaction between the interacting sulfur-aromatic molecular pair. However, the relevance of buried methionine-aromatic motifs to protein folding and function is relatively unexplored. The Small Ubiquitin-Like Modifier (SUMO) is a β-grasp fold protein and a common posttranslational modifier that affects diverse cellular processes, including transcriptional regulation, chromatin remodeling, metabolic regulation, mitosis, and meiosis. SUMO is a member of the Ubiquitin-Like (UBL) protein family. Herein, we report that a highly conserved and buried methionine-phenylalanine motif is a unique signature of SUMO proteins but absent in other homologous UBL proteins. We also detect that a specific "up" conformation between the methionine-phenylalanine pair of interacting residues in SUMO is critical to its β-grasp fold. The noncovalent interactions of SUMO with its ligands are dependent on the methionine-phenylalanine pair. MD simulations, NMR, and biophysical and biochemical studies suggest that perturbation of the methionine-aromatic motif disrupts native contacts, modulates noncovalent interactions, and attenuates SUMOylation activity. Our results highlight the importance of conserved orientations of Met-aromatic structural motifs inside a protein core for its structure and function.

3D Interaction Homology: Hydropathic Analyses of the "π-Cation" and "π-π" Interaction Motifs in Phenylalanine, Tyrosine, and Tryptophan Residues

Three-dimensional (3D) maps of the hydropathic environments of protein amino acid residues are information-rich descriptors of preferred conformations, interaction types and energetics, and solvent accessibility. The interactions made by each residue are the primary factor for rotamer selection and the secondary, tertiary, and even quaternary protein structure. Our evolving basis set of environmental data for each residue type can be used to understand the protein structure. This work focuses on the aromatic residues phenylalanine, tyrosine, and tryptophan and their structural roles. We calculated and analyzed side chain-to-environment 3D maps for over 70,000 residues of these three types that reveal, with respect to hydrophobic and polar interactions, the environment around each. After binning with backbone ϕ/ψ and side chain χ₁, we clustered each bin by 3D similarities between map-map pairs. For each of the three residue types, four bins were examined in detail: one in the β-pleat, two in the right-hand α-helix, and one in the left-hand α-helix regions of the Ramachandran plot. For high degrees of side chain burial, encapsulation of the side chain by hydrophobic interactions is ubiquitous. The more solvent-exposed side chains are more likely to be involved in polar interactions with their environments. Evidence for π-π interactions was observed in about half of the residues surveyed [phenylalanine (PHE): 53.3%, tyrosine (TYR): 34.1%, and tryptophan (TRP): 55.7%], but on an energy basis, this contributed to only ∼4% of the total. Evidence for π-cation interactions was observed in 14.1% of PHE, 8.3% of TYR, and 26.8% of TRP residues, but on an energy basis, this contributed to only ∼1%. This recognition of even these subtle interactions in the 3D hydropathic environment maps is key support for our interaction homology paradigm of protein structure elucidation and possibly prediction.

Functionally specialized human CD4+ T-cell subsets express physicochemically distinct TCRs

The organizational integrity of the adaptive immune system is determined by functionally discrete subsets of CD4⁺ T cells, but it has remained unclear to what extent lineage choice is influenced by clonotypically expressed T-cell receptors (TCRs). To address this issue, we used a high-throughput approach to profile the αβ TCR repertoires of human naive and effector/memory CD4⁺ T-cell subsets, irrespective of antigen specificity. Highly conserved physicochemical and recombinatorial features were encoded on a subset-specific basis in the effector/memory compartment. Clonal tracking further identified forbidden and permitted transition pathways, mapping effector/memory subsets related by interconversion or ontogeny. Public sequences were largely confined to particular effector/memory subsets, including regulatory T cells (Tregs), which also displayed hardwired repertoire features in the naive compartment. Accordingly, these cumulative repertoire portraits establish a link between clonotype fate decisions in the complex world of CD4⁺ T cells and the intrinsic properties of somatically rearranged TCRs.

Crystallographic binding studies of rat peroxisomal multifunctional enzyme type 1 with 3-ketodecanoyl-CoA: capturing active and inactive states of its hydratase and dehydrogenase catalytic sites

The peroxisomal multifunctional enzyme type 1 (MFE1) catalyzes two successive reactions in the β-oxidation cycle: the 2E-enoyl-CoA hydratase (ECH) and NAD+-dependent 3S-hydroxyacyl-CoA dehydrogenase (HAD) reactions. MFE1 is a monomeric enzyme that has five domains. The N-terminal part (domains A and B) adopts the crotonase fold and the C-terminal part (domains C, D and E) adopts the HAD fold. A new crystal form of MFE1 has captured a conformation in which both active sites are noncompetent. This structure, at 1.7 Å resolution, shows the importance of the interactions between Phe272 in domain B (the linker helix; helix H10 of the crotonase fold) and the beginning of loop 2 (of the crotonase fold) in stabilizing the competent ECH active-site geometry. In addition, protein crystallographic binding studies using optimized crystal-treatment protocols have captured a structure with both the 3-ketodecanoyl-CoA product and NAD+bound in the HAD active site, showing the interactions between 3-ketodecanoyl-CoA and residues of the C, D and E domains. Structural comparisons show the importance of domain movements, in particular of the C domain with respect to the D/E domains and of the A domain with respect to the HAD part. These comparisons suggest that the N-terminal part of the linker helix, which interacts tightly with domains A and E, functions as a hinge region for movement of the A domain with respect to the HAD part.

DSS1 allosterically regulates the conformation of the tower domain of BRCA2 that has dsDNA binding specificity for homologous recombination

DSS1 is an evolutionary conserved, small intrinsically disordered protein that regulates various cellular functions. Although several studies have elucidated the role of DSS1 in stabilizing BRCA2 and its importance in homologous recombination repair (HRR), yet the structural mechanism behind the stability and HRR remains elusive. In this study, using molecular dynamics simulation we show that DSS1 stabilizes linearly arranged DNA/DSS1 binding domains of BRCA2 with many native contacts. These contacts are absent in the complexes with two missense DSS1 mutants associated with germline breast cancer and somatic mouth carcinoma. Most importantly, our protein energy-based network models show DSS1 allosterically regulates the conformation of the distant tower domain of BRCA2 that has dsDNA binding specificity for HRR. We further postulate that the unique conformation of the tower domain with kinked-helices might be responsible for DNA strand invasion and initiation of HRR. Induced conformation of the tower domain by the kinked-helices is absent in the unbound BRCA2, as well as in the two mutant DSS1-BRCA2 complexes. This suggests that DSS1 allosterically regulates the tower domain conformations of BRCA2 that affects dsDNA binding, essential for HRR. Our results add a new dimension to the function of DSS1 and its role in regulating HRR.

The susceptibility of disulfide bonds towards radiation damage may be explained by S⋯O interactions

Radiation-induced damage to protein crystals during X-ray diffraction data collection is a major impediment to obtaining accurate structural information on macromolecules. Some of the specific impairments that are inflicted upon highly brilliant X-ray irradiation are metal-ion reduction, disulfide-bond cleavage and a loss of the integrity of the carboxyl groups of acidic residues. With respect to disulfide-bond reduction, previous results have indicated that not all disulfide bridges are equally susceptible to damage. A careful analysis of the chemical environment of disulfide bonds in the structures of elastase, lysozyme, acetylcholinesterase and other proteins suggests that S-S bonds which engage in a close contact with a carbonyl O atom along the extension of the S-S bond vector are more susceptible to reduction than the others. Such an arrangement predisposes electron transfer to occur from the O atom to the disulfide bond, leading to its reduction. The interaction between a nucleophile and an electrophile, akin to hydrogen bonding, stabilizes protein structures, but it also provides a pathway of electron transfer to the S-S bond, leading to its reduction during exposure of the protein crystal to an intense X-ray beam. An otherwise stabilizing interaction can thus be the cause of destabilization under the condition of radiation exposure.

Relevance of Electrostatic Charges in Compactness, Aggregation, and Phase Separation of Intrinsically Disordered Proteins

The abundance of intrinsic disorder in the protein realm and its role in a variety of physiological and pathological cellular events have strengthened the interest of the scientific community in understanding the structural and dynamical properties of intrinsically disordered proteins (IDPs) and regions (IDRs). Attempts at rationalizing the general principles underlying both conformational properties and transitions of IDPs/IDRs must consider the abundance of charged residues (Asp, Glu, Lys, and Arg) that typifies these proteins, rendering them assimilable to polyampholytes or polyelectrolytes. Their conformation strongly depends on both the charge density and distribution along the sequence (i.e., charge decoration) as highlighted by recent experimental and theoretical studies that have introduced novel descriptors. Published experimental data are revisited herein in the frame of this formalism, in a new and possibly unitary perspective. The physicochemical properties most directly affected by charge density and distribution are compaction and solubility, which can be described in a relatively simplified way by tools of polymer physics. Dissecting factors controlling such properties could contribute to better understanding complex biological phenomena, such as fibrillation and phase separation. Furthermore, this knowledge is expected to have enormous practical implications for the design, synthesis, and exploitation of bio-derived materials and the control of natural biological processes.

Cooperativity in alkane-1,2- and 1,3-polyols: NMR, QTAIM, and IQA study of O─H… OH and C─H… OH bonding interactions

Proton nuclear magnetic resonance chemical shifts and atom-atom interaction energies for alkanepolyols with 1,2-diol and 1,3-diol repeat units, and for their 1:1 pyridine complexes, are computed by density functional theory calculations. In the 1,3-polyols, based on a tG'Gg' repeat unit, the only important intramolecular hydrogen bonding interactions are O─H^… OH. By quantum theory of atoms in molecules analysis of the electron density, unstable bond and ring critical points are found for such interactions in 1,2-polyols with tG'g repeat units, from butane-1,2,3,4-tetrol onwards and in their pyridine complexes from propane-1,2,3-triol onwards. Several features (OH proton shifts and charges, and interaction energies computed by the interacting quantum atoms approach) are used to monitor the dependence of cooperativity on chain length: This is much less regular in 1,2-polyols than in 1,3-polyols and by most criteria has a higher damping factor. Well defined C─H^… OH interactions are found in butane-1,2,3,4-tetrol and higher members of the 1,2-polyol series, as well as in their pyridine complexes: There is no evidence for cooperativity with O─H^… OH bonding. For the 1,2-polyols, there is a tenuous empirical relationship between the existence of a bond critical point for O─H^… OH hydrogen bonding and the interaction energies of competing exchange channels, but the primary/secondary ratio is always less than unity.

Electrostatics does not dictate the slip-stacked arrangement of aromatic π-π interactions

Benzene dimer has long been an archetype for π-stacking. According to the Hunter–Sanders model, quadrupolar electrostatics favors an edge-to-face CH⋯π geometry but competes with London dispersion that favors cofacial π-stacking, with a compromise “slip-stacked” structure emerging as the minimum-energy geometry. This model is based on classical electrostatics, however, and neglects charge penetration. A fully quantum-mechanical analysis, presented here, demonstrates that electrostatics actually exerts very little influence on the conformational landscape of (C₆H₆)₂. Electrostatics also cannot explain the slip-stacked arrangement of C₆H₆⋯C₆F₆, where the sign of the quadrupolar interaction is reversed. Instead, the slip-stacked geometry emerges in both systems due to competition between dispersion and Pauli repulsion, with electrostatics as an ambivalent spectator. This revised interpretation helps to rationalize the persistence of offset π-stacking in larger polycyclic aromatic hydrocarbons and across the highly varied electrostatic environments that characterize π–π interactions in proteins.

According to the Hunter–Sanders model, geometries in π–π systems arise from competition between quadrupolar electrostatics (favoring an edge-to-face geometry) and London dispersion (favoring stacking), but this model misrepresents the molecular physics.

Revisiting high-resolution crystal structure of Phormidium rubidum phycocyanin

The crystal structure of phycocyanin (pr-PC) isolated from Phormidium rubidum A09DM (P. rubidum) is described at a resolution of 1.17 Å. Electron density maps derived from crystallographic data showed many clear differences in amino acid sequences when compared with the previously obtained gene-derived sequences. The differences were found in 57 positions (30 in α-subunit and 27 in β-subunit of pr-PC), in which all residues except one (β145Arg) are not interacting with the three phycocyanobilin chromophores. Highly purified pr-PC was then sequenced by mass spectrometry (MS) using LC-MS/MS. The MS data were analyzed using two independent proteomic search engines. As a result of this analysis, complete agreement between the polypeptide sequences and the electron density maps was obtained. We attribute the difference to multiple genes in the bacterium encoding the phycocyanin apoproteins and that the gene sequencing sequenced the wrong ones. We are not implying that protein sequencing by mass spectrometry is more accurate than that of gene sequencing. The final 1.17 Å structure of pr-PC allows the chromophore interactions with the protein to be described with high accuracy.

cis alkenes stabilized by intramolecular sulphurπ interactions

A series of alkenes with bistable isomers were obtained containing a thiophene/azoheteroaryl backbone. Visible light and heat-induced reversible cis ⇌ trans isomerizations were evidenced by UV-Vis and 1H NMR spectra. The stabilization of cis alkenes was attributed to intramolecular sulphurπ (Sπ) interactions, which were further supported by theoretical calculations.

Delineation of a new structural motif involving NHN γ-turn

Macromolecules are characterized by distinctive arrangement of hydrogen bonds. Different patterns of hydrogen bonds give rise to distinct and stable structural motifs. An analysis of 4114 non-redundant protein chains reveals the existence of a three-residue, (i - 1) to (i + 1), structural motif, having two hydrogen-bonded five-membered pseudo rings (the first, an NH···OC involving the first residue, and the second being NH∙∙∙N involving the last two residues), separated by a peptide bond. There could be an additional hydrogen bond between the side-chain at (i-1) and the main-chain NH of (i + 1). The average backbone torsion angles of -76(±21)° and - 12(±17)° at i creates a tight turn in the polypeptide chain, akin to a γ-turn. Indeed, a search of three-residue fragments with restriction on the terminal C^α ···C^α distance and the existence of the two pseudo rings on either side revealed the presence 14 846 cases of a variant, termed NHN γ-turn, distinct from the NHO γ-turn (2032 cases) that has traditionally been characterized by the presence of NHO hydrogen bond linking the terminal main-chain atoms. As in the latter, the newly identified γ-turns are also of two types-classical and inverse, occurring in the ratio of 1:6. The propensities of residues to occur in these turns and their secondary structural features have been enumerated. An understanding of these turns would be useful for structure prediction and loop modeling, and may serve as models to represent some of the unfolded state or disordered region in proteins.

An ab initio molecular dynamics study of benzene in water at supercritical conditions: Structure, dynamics, and polarity of hydration shell water and the solute

Anisotropic structure and dynamics of the hydration shell of a benzene solute in supercritical water are investigated by means of ab initio molecular dynamics simulations. The polarity and structural distortion of the benzene solute in supercritical water are also investigated in this study. Calculations are done at 673 K for three different densities of the solvent. The simulations are carried out using the Becke-Lee-Yang-Parr (BLYP) and also the Becke-Lee-Yang-Parr functional including dispersion corrections of Grimme (BYLP-D). The structural anisotropy is found to exist even at supercritical conditions as elucidated by the radial distribution functions of different conical regions and also by angular and spatial distribution functions. The benzene-water πH-bond and also the water-water hydrogen bonds are found to exist even at the supercritical temperature of 673 K. However, the numbers of these hydrogen bonds are reduced substantially with a decrease in water density. The water molecules in the axial region of benzene are found to be preferably oriented with one OH vector pointing toward the benzene ring, whereas the water molecules located in the equatorial region are found to orient their dipoles mostly parallel to the ring plane. The orientational distributions, however, are found to be rather broad at the supercritical temperature due to thermal fluctuations. Although the water molecules have faster dynamics at these supercritical conditions, a slight difference is observed in the dynamics of the solvation shell and bulk molecules. The conformational flexibility of the ring is found to be enhanced which causes an increase in polarity of the benzene solute in water under supercritical conditions.

Multiple RNA-RNA tertiary interactions are dispensable for formation of a functional U2/U6 RNA catalytic core in the spliceosome

The active 3D conformation of the spliceosome's catalytic U2/U6 RNA core is stabilised by a network of secondary and tertiary RNA interactions, but also depends on spliceosomal proteins for its formation. To determine the contribution towards splicing of specific RNA secondary and tertiary interactions in the U2/U6 RNA core, we introduced mutations in critical U6 nucleotides and tested their effect on splicing using a yeast in vitro U6 depletion/complementation system. Elimination of selected RNA tertiary interactions involving the U6 catalytic triad, or deletions of the bases of U6-U80 or U6-A59, had moderate to no effect on splicing, showing that the affected secondary and tertiary interactions are not required for splicing catalysis. However, removal of the base of U6-G60 of the catalytic triad completely blocked splicing, without affecting assembly of the activated spliceosome or its subsequent conversion into a B*-like complex. Our data suggest that the catalytic configuration of the RNA core that allows catalytic metal M1 binding can be maintained by Protein–RNA contacts. However, RNA stacking interactions in the U2/U6 RNA core are required for productive coordination of metal M2. The functional conformation of the U2/U6 RNA core is thus highly buffered, with overlapping contributions from RNA–RNA and Protein–RNA interactions.

Modulating Integrin αIIbβ3 Activity through Mutagenesis of Allosterically Regulated Intersubunit Contacts

Integrin αIIbβ3, a transmembrane heterodimer, mediates platelet aggregation when it switches from an inactive to an active ligand-binding conformation following platelet stimulation. Central to regulating αIIbβ3 activity is the interaction between the αIIb and β3 extracellular stalks, which form a tight heterodimer in the inactive state and dissociate in the active state. Here, we demonstrate that alanine replacements of sensitive positions in the heterodimer stalk interface destabilize the inactive conformation sufficiently to cause constitutive αIIbβ3 activation. To determine the structural basis for this effect, we performed a structural bioinformatics analysis and found that perturbing intersubunit contacts with favorable interaction geometry through substitutions to alanine quantitatively accounted for the degree of constitutive αIIbβ3 activation. This mutational study directly assesses the relationship between favorable interaction geometry at mutation-sensitive positions and the functional activity of those mutants, giving rise to a simple model that highlights the importance of interaction geometry in contributing to the stability between protein-protein interactions.

π-Hydrogen bonding and aromaticity: a systematic interplay study

Quantum DFT calculations, corrected for long-range interactions, have been carried out on complex models formed between HF as a proton donor and 2-methylene-2H-indene derivatives as proton acceptors. Using various exocyclic X substitutions, mutual effects of the aromaticity and the strength of the resulting π-hydrogen bond (after its evaluation by AIM methodology) have been investigated. The results show that the aromaticity of 6-membered rings and the hydrogen bond strength increase upon increasing the electron-donating character of the X-substituents. Based on some aromaticity indices (HOMA, FLU, SA and NICS(1)zz), it has been shown that the formation of a π-hydrogen bond causes an increase of aromaticity of the 6-membered ring. Also, the strength of the resulting π-hydrogen bond (with an energy of about 4.0 to 7.0 kcal mol-1) depends on the aromaticity of the 6-membered ring and increases with an increase in the aromaticity. In addition, a linear relationship was found between the most negative value of the molecular electrostatic potential (Vmin) and the HOMA, which confirms that the Vmin in the region of the studied ring could be used as a new index to estimate the amount of aromaticity. The electronic properties of the complexes have also been evaluated by means of the molecular electrostatic potential (MEP), the atoms in molecules (AIM) and the natural bond orbital (NBO) analyses.

The Kinetics, Thermodynamics and Mechanisms of Short Aromatic Peptide Self-Assembly

The self-assembly of short aromatic peptides and peptide derivatives into a variety of different nano- and microstructures (fibrillar gels, crystals, spheres, plates) is a promising route toward the creation of bio-compatible materials with often unexpected and useful properties. Furthermore, such simple self-assembling systems have been proposed as model systems for the self-assembly of longer peptides, a process that can be linked to biological function and malfunction. Much effort has been made in the last 15 years to explore the space of peptide sequences, chemical modifications and solvent conditions in order to maximise the diversity of assembly morphologies and properties. However, quantitative studies of the corresponding mechanisms of, and driving forces for, peptide self-assembly have remained relatively scarce until recently. In this chapter we review the current state of understanding of the thermodynamic driving forces and self-assembly mechanisms of short aromatic peptides into supramolecular structures. We will focus on experimental studies of the assembly process and our perspective will be centered around diphenylalanine (FF), a key motif of the amyloid β sequence and a paradigmatic self-assembly building block. Our main focus is the basic physical chemistry and key structural aspects of such systems, and we will also compare the mechanism of dipeptide aggregation with that of longer peptide sequences into amyloid fibrils, with discussion on how these mechanisms may be revealed through detailed analysis of growth kinetics, thermodynamics and other fundamental properties of the aggregation process.

Structural motif, topi and its role in protein function and fibrillation

A protein chain is arranged into regions in which the backbone is organized into regular patterns (of conformation and hydrogen bonding) to form the most common secondary structures, α-helix and β-sheet, which are interspersed by turns and more irregular loop regions. A structural motif, topi, is discussed in which a pair of 2-residue segments, each containing hydrogen-bonded five-membered fused-ring motifs, distant in sequence are linked to each other by a hydrogen bond. Though a small motif, it appears to be important in the context of local folding patterns of proteins and occurs near protein active sites. The motif shows quite significant residue preference, and a Cys (or Ser) occupying the second position may further stabilize the motif by forming an additional hydrogen bond across it. Remarkably, topi is found within disease causing misfolded proteins, such as the fibrilled form of Aβ42, and also across the interface between two protein chains. This motif may be an important component of fibrillation and useful for modeling loop regions.

Aromatic interactions in β-hairpin scaffold stability: A historical perspective

Non-covalent interactions between naturally occurring aromatic residues have been widely exploited as scaffold stabilizing agents in de novo designed peptides and in Nature - inspired structures. Our understanding of the factors driving aromatic interactions and their observed interaction geometries have advanced remarkably with improvements in conventional structural studies, availability of novel molecular methods and in silico studies, which have together provided atomistic information on aromatic interactions and interaction strengths. This review attempts to recapitulate the early advances in our understanding of aromatic interactions as stabilizing agents of peptide β-hairpins.

Tipping the Balance between S-π and O-π Interactions

A comprehensive experimental survey consisting of 36 molecular balances was conducted to compare 18 pairs of S-π versus O-π interactions over a wide range of structural, geometric, and solvent parameters. A strong linear correlation was observed between the folding energies of the sulfur and oxygen balances across the entire library of balance pairs. The more stable interaction systematically switched from the O-π to S-π interaction. Computational studies of bimolecular PhSCH₃-arene and PhOCH₃-arene complexes were able to replicate the experimental trends in the molecular balances. The change in preference for the O-π to S-π interaction was due to the interplay of stabilizing (dispersion and solvophobic) and destabilizing (exchange-repulsion) terms arising from the differences in size and polarizability of the oxygen and sulfur atoms.

The Pathways of the iRFP713 Unfolding Induced by Different Denaturants

Near-infrared fluorescent proteins (NIR FPs) based on the complexes of bacterial phytochromes with their natural biliverdin chromophore are widely used as genetically encoded optical probes for visualization of cellular processes and deep-tissue imaging of cells and organs in living animals. In this work, we show that the steady-state and kinetic dependencies of the various spectral characteristics of iRFP713, developed from the bacterial phytochrome RpBphP2 and recorded at protein unfolding induced by guanidine hydrochloride (GdnHCl), guanidine thiocyanate (GTC), and urea, differ substantially. A study of the unfolding of three single-tryptophan mutant forms of iRFP713 expectedly revealed that protein unfolding begins with the dissociation of the native dimer, while the monomers remain compact. A further increase in the denaturant concentration leads to the formation of an intermediate state of iRFP713 having hydrophobic areas exposed on the protein surface (I). The total surface charge of iRFP713 (pI 5.86) changes from negative to positive with an increase in the concentration of GdnHCl and GTC because the negative charge of glutamic and aspartic acids is neutralized by forming salt bridges between the carboxyl groups and GdnH⁺ ions and because the guanidinium cations bind to amide groups of glutamines and asparagines. The coincidence of both the concentration of the denaturants at which the intermediate state of iRFP713 accumulates and the concentration of GdnH⁺ ions at which the neutralization of the surface charge of the protein in this state is ensured results in strong protein aggregation. This is evidently realized by iRFP713 unfolding by GTC. At the unfolding of the protein by GdnHCl, an intermediate state is populated at higher denaturant concentrations and a strong aggregation is not observed. As expected, protein aggregates are not formed in the presence of the urea. The aggregation of the protein upon neutralization of the charge on the macromolecule surface is the main indicator of the intermediate state of protein. The unfolded state of iRFP713, whose formation is accompanied by a significant decrease in the parameter A, was found to have a different residual structure in the denaturants used.

Intramolecular O-H⋯O and C-H⋯O hydrogen bond cooperativity in D-glucopyranose and D-galactopyranose-A DFT/GIAO, QTAIM/IQA, and NCI approach

Density functional theory calculations are used to compute proton nuclear magnetic resonance (NMR) chemical shifts, interatomic distances, atom-atom interaction energies, and atomic charges for partial structures and conformers of α-D-glucopyranose, β-D-glucopyranose, and α-D-galactopyranose built up by introducing OH groups into 2-methyltetrahydropyran stepwisely. For the counterclockwise conformers, the most marked effects on the NMR shift and the charge on the OH1 proton are produced by OH2, those of OH3 and OH4 being somewhat smaller. This argues for a diminishing cooperative effect. The effect of OH6 depends on the configuration of the hydroxymethyl group and the position, axial or equatorial, of OH4, which controls hydrogen bonding in the 1,3-diol motif. Variations in the interaction energies reveal that a "new" hydrogen bond is sometimes formed at the expense of a preexisting one, probably due to geometrical constraints. Whereas previous work showed that complexing a conformer with pyridine affects only the nearest neighbour, successive OH groups increase the interaction energy of the N⋯H1 hydrogen bond and reduce its length. Analogous results are obtained for the clockwise conformers. The interaction energies for C-H⋯OH hydrogen bonding between axial CH protons and OH groups in certain conformers are much smaller than for O-H⋯OH bonds but they are largely covalent, whereas those of the latter are predominantly coulombic. These interactions are modified by complexation with pyridine in the same way as O-H⋯OH interactions: the computed NMR shifts of the CH protons increase, the atom-atom distances are shorter, and interaction energies are enhanced.

The Changing Landscape of Naive T Cell Receptor Repertoire With Human Aging

Human aging is associated with a profound loss of thymus productivity, yet naïve T lymphocytes still maintain their numbers by division in the periphery for many years. The extent of such proliferation may depend on the cytokine environment, including IL-7 and T-cell receptor (TCR) “tonic” signaling mediated by self pMHCs recognition. Additionally, intrinsic properties of distinct subpopulations of naïve T cells could influence the overall dynamics of aging-related changes within the naïve T cell compartment. Here, we investigated the differences in the architecture of TCR beta repertoires for naïve CD4, naïve CD8, naïve CD4⁺CD25⁻CD31⁺ (enriched with recent thymic emigrants, RTE), and mature naïve CD4⁺CD25⁻CD31⁻ peripheral blood subsets between young and middle-age/old healthy individuals. In addition to observing the accumulation of clonal expansions (as was shown previously), we reveal several notable changes in the characteristics of T cell repertoire. We observed significant decrease of CDR3 length, NDN insert, and number of non-template added N nucleotides within TCR beta CDR3 with aging, together with a prominent change of physicochemical properties of the central part of CDR3 loop. These changes were similar across CD4, CD8, RTE-enriched, and mature CD4 subsets of naïve T cells, with minimal or no difference observed between the latter two subsets for individuals of the same age group. We also observed an increase in “publicity” (fraction of shared clonotypes) of CD4, but not CD8 naïve T cell repertoires. We propose several explanations for these phenomena built upon previous studies of naïve T-cell homeostasis, and call for further studies of the mechanisms causing the observed changes and of consequences of these changes in respect of the possible holes formed in the landscape of naïve T cell TCR repertoire.

Dynamics of water in conical solvation shells around a benzene solute under different thermodynamic conditions

Water molecules in different parts of the anisotropic hydration shell of an aromatic molecule experience different interactions. In the present study, we investigate the anisotropic dynamics of water molecules in different non-overlapping conical shells around a benzene solute at sub- and supercritical conditions by means of molecular dynamics simulations using both non-polarizable and polarizable models. In addition to the dynamical properties, the effects of polarizability on the hydration structure of benzene at varying thermodynamic conditions are also investigated in the current study. The presence of πH-bonding in the solvation shell is found to be an important factor that influences the anisotropic dynamics of the benzene hydration shell. The water molecules located axial to the benzene plane are found to be maximally influenced by the πH-bonding. The extent of πH-bonding is found to be somewhat reduced on inclusion of polarizability. The πH-bonded water molecules are found to reorient through large-amplitude angular jumps where the jump-angle amplitude increases at higher temperatures and lower densities. For both non-polarizable and polarizable models, it is found that the water molecules in the axial conical shells possess faster orientational and hydrogen bond dynamics compared to those in the equatorial plane. Water molecules in the axial conical shells are also found to diffuse at a faster rate than bulk molecules due to the relatively weaker benzene-water πH-bonding interactions in the axial region of the hydration shell. The residence dynamics of water molecules in different conical solvation shells around the solute is also investigated in the current study.

Charge-Transfer Knowledge Graph among Amino Acids Derived from High-Throughput Electronic Structure Calculations for Protein Database

The charge-transfer coupling is an important component in tight-binding methods. Because of the highly complex chemical structure of biomolecules, the anisotropic feature of charge-transfer couplings in realistic proteins cannot be ignored. In this work, we have performed the first large-scale quantitative assessment of charge-transfer preference by calculating the charge-transfer couplings in all 20 × 20 possible amino acid side-chain combinations, which are extracted from available high-quality structures of thousands of protein complexes. The charge-transfer database quantitatively shows distinct features of charge-transfer couplings among millions of amino acid side-chain combinations. The overall distribution of charge-transfer couplings reveals that only one average or representative structure cannot be regarded as the typical charge-transfer preference in realistic proteins. This work provides us an alternative route to comprehensively understand the charge-transfer couplings for the overall distribution of realistic proteins in the foreseen big data scenario.

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.

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.