Local dynamics of proteins and DNA evaluated from crystallographic B factors

Differential Responses in the Core, Active Site and Peripheral Regions of Cytochrome c Peroxidase to Extreme Pressure and Temperature

In consideration of life in extreme environments, the effects of hydrostatic pressure on proteins at the atomic level have drawn substantial interest. Large deviations of temperature and pressure from ambient conditions can shift the free energy landscape of proteins to reveal otherwise lowly populated structural states and even promote unfolding. We report the crystal structure of the heme-containing peroxidase, cytochrome c peroxidase (CcP) at 1.5 and 3.0 kbar and make comparisons to structures determined at 1.0 bar and cryo-temperatures (100 K). Pressure produces anisotropic changes in CcP, but compressibility plateaus after 1.5 kbar. CcP responds to pressure with volume declines at the periphery of the protein where B-factors are relatively high but maintains nearly intransient core structure, hydrogen bonding interactions and active site channels. Changes in active-site solvation and heme ligation reveal pressure sensitivity to protein-ligand interactions and a potential docking site for the substrate peroxide. Compression at the surface affects neither alternate side-chain conformers nor B-factors. Thus, packing in the core, which resembles a crystalline solid, limits motion and protects the active site, whereas looser packing at the surface preserves side-chain dynamics. These data demonstrate that conformational dynamics and packing densities are not fully correlated in proteins and that encapsulation of cofactors by the polypeptide can provide a precisely structured environment resistant to change across a wide range of physical conditions.

RNAfcg: RNA Flexibility Prediction Based on Topological Centrality and Global Features

The dynamics of RNAs are related intimately to their functions. Molecular flexibility, as a starting point for understanding their dynamics, has been utilized to predict many characteristics associated with their functions. Since the experimental measurement methods are time-consuming and labor-intensive, it is urgently needed to develop reliable theoretical methods to predict RNA flexibility. In this work, we develop an effective machine learning method, RNAfcg, to predict RNA flexibility, where the Random Forest (RF) is trained by features including the topological centralities, flexibility-rigidity index, and global characteristics first introduced by us, as well as some traditional sequence and structural features. The analyses show that the three types of features introduced first have significant contributions to RNA flexibility prediction, among which the topological type contributes the most, which indicates the importance of structural topology in determining RNA flexibility. The performance comparison indicates that RNAfcg outperforms the state-of-the-art machine learning methods and the commonly used Gaussian Network Model (GNM) models, achieving a much higher Pearson correlation coefficient (PCC) of 0.6619 on the test data set. This work is helpful for understanding RNA dynamics and can be used to predict RNA function information. The source code is available at https://github.com/ChunhuaLab/RNAfcg/.

Discovery of Nonretinoid Inhibitors of CRBP1: Structural and Dynamic Insights for Ligand-Binding Mechanisms

The dysregulation of retinoid metabolism has been linked to prevalent ocular diseases including age-related macular degeneration and Stargardt disease. Modulating retinoid metabolism through pharmacological approaches holds promise for the treatment of these eye diseases. Cellular retinol-binding protein 1 (CRBP1) is the primary transporter of all-trans-retinol (atROL) in the eye, and its inhibition has recently been shown to protect mouse retinas from light-induced retinal damage. In this report, we employed high-throughput screening to identify new chemical scaffolds for competitive, nonretinoid inhibitors of CRBP1. To understand the mechanisms of interaction between CRBP1 and these inhibitors, we solved high-resolution X-ray crystal structures of the protein in complex with six selected compounds. By combining protein crystallography with hydrogen/deuterium exchange mass spectrometry, we quantified the conformational changes in CRBP1 caused by different inhibitors and correlated their magnitude with apparent binding affinities. Furthermore, using molecular dynamic simulations, we provided evidence for the functional significance of the "closed" conformation of CRBP1 in retaining ligands within the binding pocket. Collectively, our study outlines the molecular foundations for understanding the mechanism of high-affinity interactions between small molecules and CRBPs, offering a framework for the rational design of improved inhibitors for this class of lipid-binding proteins.

Knowledge-based prediction of DNA hydration using hydrated dinucleotides as building blocks

Water plays an important role in stabilizing the structure of DNA and mediating its interactions. Here, the hydration of DNA was analyzed in terms of dinucleotide fragments from an ensemble of 2727 nonredundant DNA chains containing 41 853 dinucleotides and 316 265 associated first-shell water molecules. The dinucleotides were classified into categories based on their 16 sequences and the previously determined structural classes known as nucleotide conformers (NtCs). The construction of hydrated dinucleotide building blocks allowed dinucleotide hydration to be calculated as the probability of water density distributions. Peaks in the water densities, known as hydration sites (HSs), uncovered the interplay between base and sugar-phosphate hydration in the context of sequence and structure. To demonstrate the predictive power of hydrated DNA building blocks, they were then used to predict hydration in an independent set of crystal and NMR structures. In ten tested crystal structures, the positions of predicted HSs and experimental waters were in good agreement (more than 40% were within 0.5 Å) and correctly reproduced the known features of DNA hydration, for example the `spine of hydration' in B-DNA. Therefore, it is proposed that hydrated building blocks can be used to predict DNA hydration in structures solved by NMR and cryo-EM, thus providing a guide to the interpretation of experimental data and computer models. The data for the hydrated building blocks and the predictions are available for browsing and visualization at the website https://watlas.datmos.org/watna/.

Conformational variability of loops in the SARS-CoV-2 spike protein

The SARS‐CoV‐2 spike (S) protein facilitates viral infection, and has been the focus of many structure determination efforts. Its flexible loop regions are known to be involved in protein binding and may adopt multiple conformations. This article identifies the S protein loops and studies their conformational variability based on the available Protein Data Bank structures. While most loops had essentially one stable conformation, 17 of 44 loop regions were observed to be structurally variable with multiple substantively distinct conformations based on a cluster analysis. Loop modeling methods were then applied to the S protein loop targets, and the prediction accuracies discussed in relation to the characteristics of the conformational clusters identified. Loops with multiple conformations were found to be challenging to model based on a single structural template.

Decline of protein structure rigidity with interatomic distance

BACKGROUND

Protein structural rigidity was analyzed in a non-redundant ensemble of high-resolution protein crystal structures by means of the Hirshfeld test, according to which the components (uX and uY) of the B-factors of two atoms (X and Y) along the interatomic direction is related to their degree of rigidity: the atoms may move as a rigid body if uX = uY and they cannot if uX ≠ uY.

RESULTS

It was observed that the rigidity degree diminishes if the number of covalent bonds intercalated between the two atoms (d_seq) increases, while it is rather independent on the Euclidean distance between the two atoms (d): for a given value of d_seq, the difference between uX and uY does not depend on d. No additional rigidity decline is observed when d_seq ≥  ~ 30 and this upper limit is very modest, close to 0.015 Å.

CONCLUSIONS

This suggests that protein flexibility is not fully described by B-factors that capture only partially the wide range of distortions that proteins can afford.

The β-hairpin from the Thermus thermophilus HB27 laccase works as a pH-dependent switch to regulate laccase activity

The multi-copper oxidase from the hyper-thermophilic bacteria Thermus thermophilus (Tth-MCO), has been previously characterized and described as an example of a laccase with low catalytic properties, especially when it is compared with the activity of fungal laccases, but it is active at high temperatures. Structurally, Tth-MCO has a unique feature: a β-hairpin near the T1Cu site, which is not present in any other laccases deposited at the PDB. This β-hairpin has an expected crystallographic behavior in solvent-exposed areas of a crystallized protein: lack of electron density, high B-values and several crystalline contacts with neighboring crystallographic copies; however, its dynamical behavior in solution and its biological implications have not been described. Here, we describe four new Tth-MCO crystallographic structures, and the β-hairpin behavior has been analyzed by molecular dynamics simulations, considering the effect of pH and temperature. The β-hairpin new crystallographic conformations described here, together with their dynamics, were used to understand the pH-restrained laccase activity of Tth-MCO against substrates as syringaldazine. Remarkably, there are insertions in laccases from Thermus and Meiothermus genus, sharing the same position and a methionine-rich composition of the Tth-MCO β-hairpin. This unique high methionine content of the Tth-MCO β-hairpin is responsible to coordinate, Ag⁺¹ and Hg⁺¹ in oxidative conditions, but Cu⁺¹ and Cu⁺² are not coordinated in crystallographic experiments, regardless of the redox conditions; however, Ag⁺¹ addition does not affect Tth-MCO laccase activity against syringaldazine. Here, we propose that the pH-dependent β-hairpin dynamical behavior could explain, at least in part, the inefficient laccase activity displayed by Tth-MCO in acidic pH values.

Conformational Strain Indicated by Ramachandran Angles for the Protein Backbone Is Only Weakly Related to the Flexibility

Studies on energy associated with free dipeptides have shown that conformers with unfavorable (ϕ,ψ) torsion angles have higher energy compared to conformers with favorable (ϕ,ψ) angles. It is expected that higher energy confers higher dynamics and flexibility to that part of the protein. Here, we explore a potential relationship between conformational strain in a residue due to unfavorable (ϕ,ψ) angles and its flexibility and dynamics in the context of protein structures. We compared flexibility of strained and relaxed residues, which are recognized based on outlier/allowed and favorable (ϕ,ψ) angles respectively, using normal-mode analysis (NMA). We also performed in-depth analysis on flexibility and dynamics at catalytic residues in protein kinases, which exhibit different strain status in different kinase structures using NMA and molecular dynamics simulations. We underline that strain of a residue, as defined by backbone torsion angles, is almost unrelated to the flexibility and dynamics associated with it. Even the overall trend observed among all high-resolution structures in which relaxed residues tend to have slightly higher flexibility than strained residues is counterintuitive. Consequently, we propose that identifying strained residues based on (ϕ,ψ) values is not an effective way to recognize energetic strain in protein structures.

Molecular basis for the interaction of cellular retinol binding protein 2 (CRBP2) with nonretinoid ligands

Present in the small intestine, cellular retinol binding protein 2 (CRBP2) plays an important role in the uptake, transport, and metabolism of dietary retinoids. However, the recent discovery of the interactions of CRBP2 with 2-arachidonoylglycerol and other monoacylglycerols (MAGs) suggests the broader involvement of this protein in lipid metabolism and signaling. To better understand the physiological role of CRBP2, we determined its protein-lipid interactome using a fluorescence-based retinol replacement assay adapted for a high-throughput screening format. By examining chemical libraries of bioactive lipids, we provided evidence for the selective interaction of CRBP2 with a subset of nonretinoid ligands with the highest affinity for sn-1 and sn-2 MAGs that contain polyunsaturated C18-C20 acyl chains. We also elucidated the structure-affinity relationship for nonretinoid ligands of this protein. We further dissect the molecular basis for this ligand's specificity by analyzing high-resolution crystal structures of CRBP2 in complex with selected derivatives of MAGs. Finally, we identify T51 and V62 as key amino acids that enable the broadening of ligand selectivity to MAGs in CRBP2 as compared with retinoid-specific CRBP1. Thus, our study provides the molecular framework for understanding the lipid selectivity and diverse functions of CRBPs in controlling lipid homeostasis.

Structural Insights into the Interaction of Clinically Relevant Phosphorothioate mRNA Cap Analogs with Translation Initiation Factor 4E Reveal Stabilization via Electrostatic Thio-Effect

mRNA-based therapies and vaccines constitute a disruptive technology with the potential to revolutionize modern medicine. Chemically modified 5′ cap structures have provided access to mRNAs with superior translational properties that could benefit the currently flourishing mRNA field. Prime examples of compounds that enhance mRNA properties are antireverse cap analog diastereomers that contain an O-to-S substitution within the β-phosphate (β-S-ARCA D1 and D2), where D1 is used in clinically investigated mRNA vaccines. The compounds were previously found to have high affinity for eukaryotic translation initiation factor 4E (eIF4E) and augment translation in vitro and in vivo. However, the molecular basis for the beneficial “thio-effect” remains unclear. Here, we employed multiple biophysical techniques and captured 11 cap analog-eIF4E crystallographic structures to investigate the consequences of the β-O-to-S or -Se substitution on the interaction with eIF4E. We determined the S_P/R_P configurations of β-S-ARCA and related compounds and obtained structural insights into the binding. Unexpectedly, in both stereoisomers, the β-S/Se atom occupies the same binding cavity between Lys162 and Arg157, indicating that the key driving force for complex stabilization is the interaction of negatively charged S/Se with positively charged amino acids. This was observed for all structural variants of the cap and required significantly different conformations of the triphosphate for each diastereomer. This finding explains why both β-S-ARCA diastereomers have higher affinity for eIF4E than unmodified caps. Binding affinities determined for di-, tri-, and oligonucleotide cap analogs suggested that the “thio-effect” was preserved in longer RNAs. Our observations broaden the understanding of thiophosphate biochemistry and enable the rational design of translationally active mRNAs and eIF4E-targeting drugs.

BANΔIT: B'-Factor Analysis for Drug Design and Structural Biology

The analysis of B‐factor profiles from X‐ray protein structures can be utilized for structure‐based drug design since protein mobility changes have been associated with the quality of protein‐ligand interactions. With the BANΔIT (B’‐factor analysis and ΔB’ interpretation toolkit), we have developed a JavaScript‐based browser application that provides a graphical user interface for the normalization and analysis of B’‐factor profiles. To emphasize the usability for rational drug design applications, we have analyzed a selection of crystallographic protein‐ligand complexes and have given exemplary conclusions for further drug optimization including the development of a B’‐factor‐supported pharmacophore model for SARS CoV‐2 main protease inhibitors. BANΔIT is available online at https://bandit.uni‐mainz.de. The source code can be downloaded from https://github.com/FBarthels/BANDIT.

Essential site scanning analysis: A new approach for detecting sites that modulate the dispersion of protein global motions

Despite the wealth of methods developed for exploring the molecular basis of allostery in biomolecular systems, there is still a need for structure-based predictive tools that can efficiently detect susceptible sites for triggering allosteric responses. Toward this goal, we introduce here an elastic network model (ENM)-based method, Essential Site Scanning Analysis (ESSA). Essential sites are here defined as residues that would significantly alter the protein's global dynamics if bound to a ligand. To mimic the crowding induced upon substrate binding, the heavy atoms of each residue are incorporated as additional network nodes into the α-carbon-based ENM, and the resulting shifts in soft mode frequencies are used as a metric for evaluating the essentiality of each residue. Results on a dataset of monomeric proteins indicate the enrichment of allosteric and orthosteric binding sites, as well as global hinge regions among essential residues, highlighting the significant role of these sites in controlling the overall structural dynamics. Further integration of ESSA with information on predicted pockets and their local hydrophobicity density enables successful predictions of allosteric pockets for both ligand-bound and -unbound structures. ESSA can be efficiently applied to large multimeric systems. Three case studies, namely (i) G-protein binding to a GPCR, (ii) heterotrimeric assembly of the Ser/Thr protein phosphatase PP2A, and (iii) allo-targeting of AMPA receptor, demonstrate the utility of ESSA for identifying essential sites and narrowing down target allosteric sites identified by druggability simulations.

Ligand-Induced Conformational Dynamics of A Tyramine Receptor from Sitophilus oryzae

Tyramine receptor (TyrR) is a biogenic amine G protein-coupled receptor (GPCR) associated with many important physiological functions in insect locomotion, reproduction, and pheromone response. Binding of specific ligands to the TyrR triggers conformational changes, relays the signal to G proteins, and initiates an appropriate cellular response. Here, we monitor the binding effect of agonist compounds, tyramine and amitraz, to a Sitophilus oryzae tyramine receptor (SoTyrR) homology model and their elicited conformational changes. All-atom molecular dynamics (MD) simulations of SoTyrR-ligand complexes have shown varying dynamic behavior, especially at the intracellular loop 3 (IL3) region. Moreover, in contrast to SoTyrR-tyramine, SoTyrR-amitraz and non-liganded SoTyrR shows greater flexibility at IL3 residues and were found to be coupled to the most dominant motion in the receptor. Our results suggest that the conformational changes induced by amitraz are different from the natural ligand tyramine, albeit being both agonists of SoTyrR. This is the first attempt to understand the biophysical implication of amitraz and tyramine binding to the intracellular domains of TyrR. Our data may provide insights into the early effects of ligand binding to the activation process of SoTyrR.

Discrete analyses of protein dynamics

Protein structures are highly dynamic macromolecules. This dynamics is often analysed through experimental and/or computational methods only for an isolated or a limited number of proteins. Here, we explore large-scale protein dynamics simulation to observe dynamics of local protein conformations using different perspectives. We analysed molecular dynamics to investigate protein flexibility locally, using classical approaches such as RMSf, solvent accessibility, but also innovative approaches such as local entropy. First, we focussed on classical secondary structures and analysed specifically how β-strand, β-turns, and bends evolve during molecular simulations. We underlined interesting specific bias between β-turns and bends, which are considered as the same category, while their dynamics show differences. Second, we used a structural alphabet that is able to approximate every part of the protein structures conformations, namely protein blocks (PBs) to analyse (i) how each initial local protein conformations evolve during dynamics and (ii) if some exchange can exist among these PBs. Interestingly, the results are largely complex than simple regular/rigid and coil/flexible exchange. AbbreviationsN_eqnumber of equivalentPBProtein BlocksPDBProtein DataBankRMSfroot mean square fluctuationsCommunicated by Ramaswamy H. Sarma.

Structural changes in DNA-binding proteins on complexation

Characterization and prediction of the DNA-biding regions in proteins are essential for our understanding of how proteins recognize/bind DNA. We analyze the unbound (U) and the bound (B) forms of proteins from the protein–DNA docking benchmark that contains 66 binary protein–DNA complexes along with their unbound counterparts. Proteins binding DNA undergo greater structural changes on complexation (in particular, those in the enzyme category) than those involved in protein–protein interactions (PPI). While interface atoms involved in PPI exhibit an increase in their solvent-accessible surface area (ASA) in the bound form in the majority of the cases compared to the unbound interface, protein–DNA interactions indicate increase and decrease in equal measure. In 25% structures, the U form has missing residues which are located in the interface in the B form. The missing atoms contribute more toward the buried surface area compared to other interface atoms. Lys, Gly and Arg are prominent in disordered segments that get ordered in the interface on complexation. In going from U to B, there may be an increase in coil and helical content at the expense of turns and strands. Consideration of flexibility cannot distinguish the interface residues from the surface residues in the U form.

DNA-protein interaction: identification, prediction and data analysis

Life in living organisms is dependent on specific and purposeful interaction between other molecules. Such purposeful interactions make the various processes inside the cells and the bodies of living organisms possible. DNA-protein interactions, among all the types of interactions between different molecules, are of considerable importance. Currently, with the development of numerous experimental techniques, diverse methods are convenient for recognition and investigating such interactions. While the traditional experimental techniques to identify DNA-protein complexes are time-consuming and are unsuitable for genome-scale studies, the current high throughput approaches are more efficient in determining such interaction at a large-scale, but they are clearly too costly to be practice for daily applications. Hence, according to the availability of much information related to different biological sequences and clearing different dimensions of conditions in which such interactions are formed, with the developments related to the computer, mathematics, and statistics motivate scientists to develop bioinformatics tools for prediction the interaction site(s). Until now, there has been much progress in this field. In this review, the factors and conditions governing the interaction and the laboratory techniques for examining such interactions are addressed. In addition, developed bioinformatics tools are introduced and compared for this reason and, in the end, several suggestions are offered for the promotion of such tools in prediction with much more precision.

Are crystallographic B-factors suitable for calculating protein conformational entropy?

It is not possible to obtain reliable entropy estimates from crystallographic B-factors even with re-refined or room-temperature crystal structures.

Cysteine 180 Is a Redox Sensor Modulating the Activity of Human Pyridoxal 5'-Phosphate Histidine Decarboxylase

Histidine decarboxylase is a pyridoxal 5'-phosphate enzyme catalyzing the conversion of histidine to histamine, a bioactive molecule exerting its role in many modulatory processes. The human enzyme is involved in many physiological functions, such as neurotransmission, gastrointestinal track function, cell growth, and differentiation. Here, we studied the functional properties of the human enzyme and, in particular, the effects exerted at the protein level by two cysteine residues: Cys-180 and Cys-418. Surprisingly, the enzyme exists in an equilibrium between a reduced and an oxidized form whose extent depends on the redox state of Cys-180. Moreover, we determined that (i) the two enzymatic redox species exhibit modest structural changes in the coenzyme microenvironment and (ii) the oxidized form is slightly more active and stable than the reduced one. These data are consistent with the model proposed by bioinformatics analyses and molecular dynamics simulations in which the Cys-180 redox state could be responsible for a structural transition affecting the C-terminal domain reorientation leading to active site alterations. Furthermore, the biochemical properties of the purified C180S and C418S variants reveal that C180S behaves like the reduced form of the wild-type enzyme, while C418S is sensitive to reductants like the wild-type enzyme, thus allowing the identification of Cys-180 as the redox sensitive switch. On the other hand, Cys-418 appears to be a residue involved in aggregation propensity. A possible role for Cys-180 as a regulatory switch in response to different cellular redox conditions could be suggested.

The crystal structure of XdpB, the bacterial old yellow enzyme, in an FMN-free form

Old Yellow Enzymes (OYEs) are NAD(P)H dehydrogenases of not fully resolved physiological roles that are widespread among bacteria, plants, and fungi and have a great potential for biotechnological applications. We determined the apo form crystal structure of a member of the OYE class, glycerol trinitrate reductase XdpB, from Agrobacterium bohemicum R89-1 at 2.1 Å resolution. In agreement with the structures of the related bacterial OYEs, the structure revealed the TIM barrel fold with an N-terminal β-hairpin lid, but surprisingly, the structure did not contain its cofactor FMN. Its putative binding site was occupied by a pentapeptide TTSDN from the C-terminus of a symmetry related molecule. Biochemical experiments confirmed a specific concentration-dependent oligomerization and a low FMN content. The blocking of the FMN binding site can exist in vivo and regulates enzyme activity. Our bioinformatic analysis indicated that a similar self-inhibition could be expected in more OYEs which we designated as subgroup OYE C1. This subgroup is widespread among G-bacteria and can be recognized by the conserved sequence GxxDYP in proximity of the C termini. In proteobacteria, the C1 subgroup OYEs are typically coded in one operon with short-chain dehydrogenase. This operon is controlled by the tetR-like transcriptional regulator. OYEs coded in these operons are unlikely to be involved in the oxidative stress response as the other known members of the OYE family because no upregulation of XdpB was observed after exposing A. bohemicum R89-1 to oxidative stress.

RABDAM: quantifying specific radiation damage in individual protein crystal structures

Radiation damage remains one of the major limitations to accurate structure determination in protein crystallography (PX). Despite the use of cryo-cooling techniques, it is highly probable that a number of the structures deposited in the Protein Data Bank (PDB) have suffered substantial radiation damage as a result of the high flux densities of third generation synchrotron X-ray sources. Whereas the effects of global damage upon diffraction pattern reflection intensities are readily detectable, traditionally the (earlier onset) site-specific structural changes induced by radiation damage have proven difficult to identify within individual PX structures. More recently, however, development of the BDamage metric has helped to address this problem. BDamage is a quantitative, per-atom metric identifies potential sites of specific damage by comparing the atomic B-factor values of atoms that occupy a similar local packing density environment in the structure. Building upon this past work, this article presents a program, RABDAM, to calculate the BDamage metric for all selected atoms within any standard-format PDB or mmCIF file. RABDAM provides several useful outputs to assess the extent of damage suffered by an input PX structure. This free and open-source software will allow assessment and improvement of the quality of PX structures both previously and newly deposited in the PDB.

QStatin, a Selective Inhibitor of Quorum Sensing in Vibrio Species

Pathogenic Vibrio species cause diseases in diverse marine animals reared in aquaculture. Since their pathogenesis, persistence, and survival in marine environments are regulated by quorum sensing (QS), QS interference has attracted attention as a means to control these bacteria in aquatic settings. A few QS inhibitors of Vibrio species have been reported, but detailed molecular mechanisms are lacking. Here, we identified a novel, potent, and selective Vibrio QS inhibitor, named QStatin [1-(5-bromothiophene-2-sulfonyl)-1H-pyrazole], which affects Vibrio harveyi LuxR homologues, the well-conserved master transcriptional regulators for QS in Vibrio species. Crystallographic and biochemical analyses showed that QStatin binds tightly to a putative ligand-binding pocket in SmcR, the LuxR homologue in V. vulnificus, and changes the flexibility of the protein, thereby altering its transcription regulatory activity. Transcriptome analysis revealed that QStatin results in SmcR dysfunction, affecting the expression of SmcR regulon required for virulence, motility/chemotaxis, and biofilm dynamics. Notably, QStatin attenuated representative QS-regulated phenotypes in various Vibrio species, including virulence against the brine shrimp (Artemia franciscana). Together, these results provide molecular insights into the mechanism of action of an effective, sustainable QS inhibitor that is less susceptible to resistance than other antimicrobial agents and useful in controlling the virulence of Vibrio species in aquacultures.IMPORTANCE Yields of aquaculture, such as penaeid shrimp hatcheries, are greatly affected by vibriosis, a disease caused by pathogenic Vibrio infections. Since bacterial cell-to-cell communication, known as quorum sensing (QS), regulates pathogenesis of Vibrio species in marine environments, QS inhibitors have attracted attention as alternatives to conventional antibiotics in aquatic settings. Here, we used target-based high-throughput screening to identify QStatin, a potent and selective inhibitor of V. harveyi LuxR homologues, which are well-conserved master QS regulators in Vibrio species. Structural and biochemical analyses revealed that QStatin binds tightly to a putative ligand-binding pocket on SmcR, the LuxR homologue in V. vulnificus, and affects expression of QS-regulated genes. Remarkably, QStatin attenuated diverse QS-regulated phenotypes in various Vibrio species, including pathogenesis against brine shrimp, with no impact on bacterial viability. Taken together, the results suggest that QStatin may be a sustainable antivibriosis agent useful in aquacultures.

A DNA Structural Alphabet Distinguishes Structural Features of DNA Bound to Regulatory Proteins and in the Nucleosome Core Particle

We analyzed the structural behavior of DNA complexed with regulatory proteins and the nucleosome core particle (NCP). The three-dimensional structures of almost 25 thousand dinucleotide steps from more than 500 sequentially non-redundant crystal structures were classified by using DNA structural alphabet CANA (Conformational Alphabet of Nucleic Acids) and associations between ten CANA letters and sixteen dinucleotide sequences were investigated. The associations showed features discriminating between specific and non-specific binding of DNA to proteins. Important is the specific role of two DNA structural forms, A-DNA, and BII-DNA, represented by the CANA letters AAA and BB2: AAA structures are avoided in non-specific NCP complexes, where the wrapping of the DNA duplex is explained by the periodic occurrence of BB2 every 10.3 steps. In both regulatory and NCP complexes, the extent of bending of the DNA local helical axis does not influence proportional representation of the CANA alphabet letters, namely the relative incidences of AAA and BB2 remain constant in bent and straight duplexes.

DNA conformational transitions inferred from re-evaluation of m|Fo| - D|Fc| electron-density maps

Conformational flexibility of DNA plays important roles in biological processes such as transcriptional regulation and DNA packaging etc. To understand the mechanisms of these processes, it is important to analyse when, where and how DNA shows conformational variations. Recent analyses have indicated that conventional refinement methods do not always provide accurate models of crystallographic heterogeneities and that some information on polymorphism has been overlooked in previous crystallographic studies. In the present study, the m|Fo| - D|Fc| electron-density maps of double-helical DNA crystal structures were calculated at a resolution equal to or better than 1.5 Å and potential conformational transitions were found in 27% of DNA phosphates. Detailed analyses of the m|Fo| - D|Fc| peaks indicated that some of these unassigned densities correspond to ZI ↔ ZII or A/B → BI conformational transitions. A relationship was also found between ZI/ZII transitions and metal coordination in Z-DNA from the detected peaks. The present study highlights that frequent transitions of phosphate backbones occur even in crystals and that some of these transitions are affected by the local molecular environment.

Crystal structure of a DNA aptamer bound to PvLDH elucidates novel single-stranded DNA structural elements for folding and recognition

Structural elements are key elements for understanding single-stranded nucleic acid folding. Although various RNA structural elements have been documented, structural elements of single-stranded DNA (ssDNA) have rarely been reported. Herein, we determined a crystal structure of PvLDH in complex with a DNA aptamer called pL1. This aptamer folds into a hairpin-bulge contact by adopting three novel structural elements, viz, DNA T-loop-like motif, base-phosphate zipper, and DNA G·G metal ion zipper. Moreover, the pL1:PvLDH complex shows unique properties compared with other protein:nucleic acid complexes. Generally, extensive intermolecular hydrogen bonds occur between unpaired nucleotides and proteins for specific recognitions. Although most protein-interacting nucleotides of pL1 are unpaired nucleotides, pL1 recognizes PvLDH by predominant shape complementarity with many bridging water molecules owing to the combination of three novel structural elements making protein-binding unpaired nucleotides stable. Moreover, the additional set of Plasmodium LDH residues which were shown to form extensive hydrogen bonds with unpaired nucleotides of 2008s does not participate in the recognition of pL1. Superimposition of the pL1:PvLDH complex with hLDH reveals steric clashes between pL1 and hLDH in contrast with no steric clashes between 2008s and hLDH. Therefore, specific protein recognition mode of pL1 is totally different from that of 2008s.

Molecular Dynamics of Photosystem II Embedded in the Thylakoid Membrane

Photosystem II (PSII) is one of the key protein complexes in photosynthesis. We introduce a coarse grained model of PSII and present the analysis of 60 μs molecular dynamics simulations of PSII in both monomeric and dimeric form, embedded in a thylakoid membrane model that reflects its native lipid composition. We describe in detail the setup of the protein complex and the many natural cofactors and characterize their mobility. Overall we find that the protein subunits and cofactors are more flexible toward the periphery of the complex as well as near the PLQ exchange cavity and at the dimer interface. Of all cofactors, β-carotenes show the highest mobility. Some of the β-carotenes diffuse in and out of the protein complex via the thylakoid membrane. In contrast with the PSII dimer, the monomeric form adopts a tilted conformation in the membrane, with strong interactions between the soluble PsbO subunit and the glycolipid headgroups. Interestingly, the tilted conformation causes buckling of the membrane. Together, our results provide an unprecedented view of PSII dynamics on a microsecond time scale. Our data may be used as basis for the interpretation of experimental data as well as for theoretical models describing exciton energy transfer.

Hydration of proteins and nucleic acids: Advances in experiment and theory. A review

BACKGROUND

Most biological processes involve water, and the interactions of biomolecules with water affect their structure, function and dynamics.

SCOPE OF REVIEW

This review summarizes the current knowledge of protein and nucleic acid interactions with water, with a special focus on the biomolecular hydration layer. Recent developments in both experimental and computational methods that can be applied to the study of hydration structure and dynamics are reviewed, including software tools for the prediction and characterization of hydration layer properties.

MAJOR CONCLUSIONS

In the last decade, important advances have been made in our understanding of the factors that determine how biomolecules and their aqueous environment influence each other. Both experimental and computational methods contributed to the gradually emerging consensus picture of biomolecular hydration.

GENERAL SIGNIFICANCE

An improved knowledge of the structural and thermodynamic properties of the hydration layer will enable a detailed understanding of the various biological processes in which it is involved, with implications for a wide range of applications, including protein-structure prediction and structure-based drug design.

Ligand Binding Induces Conformational Changes in Human Cellular Retinol-binding Protein 1 (CRBP1) Revealed by Atomic Resolution Crystal Structures

Important in regulating the uptake, storage, and metabolism of retinoids, cellular retinol-binding protein 1 (CRBP1) is essential for trafficking vitamin A through the cytoplasm. However, the molecular details of ligand uptake and targeted release by CRBP1 remain unclear. Here we report the first structure of CRBP1 in a ligand-free form as well as ultra-high resolution structures of this protein bound to either all-trans-retinol or retinylamine, the latter a therapeutic retinoid that prevents light-induced retinal degeneration. Superpositioning of human apo- and holo-CRBP1 revealed major differences within segments surrounding the entrance to the retinoid-binding site. These included α-helix II and hairpin turns between β-strands βC-βD and βE-βF as well as several side chains, such as Phe-57, Tyr-60, and Ile-77, that change their orientations to accommodate the ligand. Additionally, we mapped hydrogen bond networks inside the retinoid-binding cavity and demonstrated their significance for the ligand affinity. Analyses of the crystallographic B-factors indicated several regions with higher backbone mobility in the apoprotein that became more rigid upon retinoid binding. This conformational flexibility of human apo-CRBP1 facilitates interaction with the ligands, whereas the more rigid holoprotein structure protects the labile retinoid moiety during vitamin A transport. These findings suggest a mechanism of induced fit upon ligand binding by mammalian cellular retinol-binding proteins.

Decoding the Structural Bases of D76N ß2-Microglobulin High Amyloidogenicity through Crystallography and Asn-Scan Mutagenesis

D76N is the first natural variant of human β-2 microglobulin (β2m) so far identified. Contrary to the wt protein, this mutant readily forms amyloid fibres in physiological conditions, leading to a systemic and severe amyloidosis. Although the Asp76Asn mutant has been extensively characterized, the molecular bases of its instability and aggregation propensity remain elusive. In this work all Asp residues of human β2m were individually substituted to Asn; D-to-N mutants (D34N, D38N, D53N, D59N, D96N and D98N) were characterised in terms of thermodynamic stability and aggregation propensity. Moreover, crystal structures of the D38N, D53N, D59N and D98N variants were solved at high-resolution (1.24–1.70 Å). Despite showing some significant variations in their thermal stabilities, none showed the dramatic drop in melting temperature (relative to the wt protein) as observed for the pathogenic mutant. Consistently, none of the variants here described displayed any increase in aggregation propensity under the experimental conditions tested. The crystal structures confirmed that D-to-N mutations are generally well tolerated, and lead only to minor reorganization of the side chains in close proximity of the mutated residue. D38N is the only exception, where backbone readjustments and a redistribution of the surface electrostatic charges are observed. Overall, our results suggest that neither removing negative charges at sites 34, 38, 53, 59, 96 and 98, nor the difference in β2m pI, are the cause of the aggressive phenotype observed in D76N. We propose that the dramatic effects of the D76N natural mutation must be linked to effects related to the crucial location of this residue within the β2m fold.

Crystal structure of group II intron domain 1 reveals a template for RNA assembly

Although the importance of large noncoding RNAs is increasingly appreciated, our understanding of their structures and architectural dynamics remains limited. In particular, we know little about RNA folding intermediates and how they facilitate the productive assembly of RNA tertiary structures. Here, we report the crystal structure of an obligate intermediate that is required during the earliest stages of group II intron folding. Composed of domain 1 from the Oceanobacillus iheyensis group II intron (266 nucleotides), this intermediate retains native-like features but adopts a compact conformation in which the active site cleft is closed. Transition between this closed and the open (native) conformation is achieved through discrete rotations of hinge motifs in two regions of the molecule. The open state is then stabilized by sequential docking of downstream intron domains, suggesting a 'first come, first folded' strategy that may represent a generalizable pathway for assembly of large RNA and ribonucleoprotein structures.

Structure of the ordered hydration of amino acids in proteins: analysis of crystal structures

Crystallography provides unique information about the arrangement of water molecules near protein surfaces. Using a nonredundant set of 2818 protein crystal structures with a resolution of better than 1.8 Å, the extent and structure of the hydration shell of all 20 standard amino-acid residues were analyzed as function of the residue conformation, secondary structure and solvent accessibility. The results show how hydration depends on the amino-acid conformation and the environment in which it occurs. After conformational clustering of individual residues, the density distribution of water molecules was compiled and the preferred hydration sites were determined as maxima in the pseudo-electron-density representation of water distributions. Many hydration sites interact with both main-chain and side-chain amino-acid atoms, and several occurrences of hydration sites with less canonical contacts, such as carbon-donor hydrogen bonds, OH-π interactions and off-plane interactions with aromatic heteroatoms, are also reported. Information about the location and relative importance of the empirically determined preferred hydration sites in proteins has applications in improving the current methods of hydration-site prediction in molecular replacement, ab initio protein structure prediction and the set-up of molecular-dynamics simulations.

Two Loops Undergoing Concerted Dynamics Regulate the Activity of the ASH1L Histone Methyltransferase

ASH1L (absent, small, or homeotic-like 1) is a histone methyltransferase (HMTase) involved in gene activation that is overexpressed in multiple forms of cancer. Previous studies of ASH1L's catalytic SET domain identified an autoinhibitory loop that blocks access of histone substrate to the enzyme active site. Here, we used both nuclear magnetic resonance and X-ray crystallography to identify conformational dynamics in the ASH1L autoinhibitory loop. Using site-directed mutagenesis, we found that point mutations in the autoinhibitory loop that perturb the structure of the SET domain result in decreased enzyme activity, indicating that the autoinhibitory loop is not a simple gate to the active site but is rather a key feature critical to ASH1L function. We also identified a second loop in the SET-I subdomain of ASH1L that experiences conformational dynamics, and we trapped two different conformations of this loop using crystallographic studies. Mutation of the SET-I loop led to a large decrease in ASH1L enzymatic activity in addition to a significant conformational change in the SET-I loop, demonstrating the importance of the structure and dynamics of the SET-I loop to ASH1L function. Furthermore, we found that three C-terminal chromatin-interacting domains greatly enhance ASH1L enzymatic activity and that ASH1L requires native nucleosome substrate for robust activity. Our study illuminates the role of concerted conformational dynamics in ASH1L function and identifies structural features important for ASH1L enzymatic activity.

Protein flexibility in the light of structural alphabets

Protein structures are valuable tools to understand protein function. Nonetheless, proteins are often considered as rigid macromolecules while their structures exhibit specific flexibility, which is essential to complete their functions. Analyses of protein structures and dynamics are often performed with a simplified three-state description, i.e., the classical secondary structures. More precise and complete description of protein backbone conformation can be obtained using libraries of small protein fragments that are able to approximate every part of protein structures. These libraries, called structural alphabets (SAs), have been widely used in structure analysis field, from definition of ligand binding sites to superimposition of protein structures. SAs are also well suited to analyze the dynamics of protein structures. Here, we review innovative approaches that investigate protein flexibility based on SAs description. Coupled to various sources of experimental data (e.g., B-factor) and computational methodology (e.g., Molecular Dynamic simulation), SAs turn out to be powerful tools to analyze protein dynamics, e.g., to examine allosteric mechanisms in large set of structures in complexes, to identify order/disorder transition. SAs were also shown to be quite efficient to predict protein flexibility from amino-acid sequence. Finally, in this review, we exemplify the interest of SAs for studying flexibility with different cases of proteins implicated in pathologies and diseases.

Redesigning protein cavities as a strategy for increasing affinity in protein-protein interaction: interferon- γ receptor 1 as a model

Combining computational and experimental tools, we present a new strategy for designing high affinity variants of a binding protein. The affinity is increased by mutating residues not at the interface, but at positions lining internal cavities of one of the interacting molecules. Filling the cavities lowers flexibility of the binding protein, possibly reducing entropic penalty of binding. The approach was tested using the interferon-γ receptor 1 (IFNγR1) complex with IFNγ as a model. Mutations were selected from 52 amino acid positions lining the IFNγR1 internal cavities by using a protocol based on FoldX prediction of free energy changes. The final four mutations filling the IFNγR1 cavities and potentially improving the affinity to IFNγ were expressed, purified, and refolded, and their affinity towards IFNγ was measured by SPR. While individual cavity mutations yielded receptor constructs exhibiting only slight increase of affinity compared to WT, combinations of these mutations with previously characterized variant N96W led to a significant sevenfold increase. The affinity increase in the high affinity receptor variant N96W+V35L is linked to the restriction of its molecular fluctuations in the unbound state. The results demonstrate that mutating cavity residues is a viable strategy for designing protein variants with increased affinity.