B-factor accuracy in protein crystal structures

ATLAS: protein flexibility description from atomistic molecular dynamics simulations

Dynamical behaviour is one of the most crucial protein characteristics. Despite the advances in the field of protein structure resolution and prediction, analysis and prediction of protein dynamic properties remains a major challenge, mostly due to the low accessibility of data and its diversity and heterogeneity. To address this issue, we present ATLAS, a database of standardised all-atom molecular dynamics simulations, accompanied by their analysis in the form of interactive diagrams and trajectory visualisation. ATLAS offers a large-scale view and valuable insights on protein dynamics for a large and representative set of proteins, by combining data obtained through molecular dynamics simulations with information extracted from experimental structures. Users can easily analyse dynamic properties of functional protein regions, such as domain limits (hinge positions) and residues involved in interaction with other biological molecules. Additionally, the database enables exploration of proteins with uncommon dynamic properties conditioned by their environment such as chameleon subsequences and Dual Personality Fragments. The ATLAS database is freely available at https://www.dsimb.inserm.fr/ATLAS.

Using graphlet degree vectors to predict atomic displacement parameters in protein structures

In structural biology, atomic displacement parameters, commonly used in the form of B values, describe uncertainties in atomic positions. Their distribution over the structure can provide hints on local structural reliability and mobility. A spatial macromolecular model can be represented by a graph whose nodes are atoms and whose edges correspond to all interatomic contacts within a certain distance. Small connected subgraphs, called graphlets, provide information about the wiring of a particular atom. The multiple linear regression approach based on this information aims to predict a distribution of values of isotropic atomic displacement parameters (B values) within a protein structure, given the atomic coordinates and molecular packing. By modeling the dynamic component of atomic uncertainties, this method allows the B values obtained from experimental crystallographic or cryo-electron microscopy studies to be reproduced relatively well.

How ambiguity codes specify molecular descriptors and information flow in Code Biology

The article presents IUPAC ambiguity codes for incomplete nucleic acid specification, and their use in Code Biology. It is shown how to use this nomenclature in order to extract accurate information on different properties of the biological systems. We investigated the use of ambiguity codes, as mathematical and logical operators and truth table elements, for the encoding of amino acids by means of the Standard Genetic Code. It is explained how to use ambiguity codes and truth functions in order to obtain accurate information on different properties of the biological systems. Nucleotide ambiguity codes could be applied to: 1. encoding descriptive information of nucleotides, amino acids and proteins (e.g., of polarity, relative solvent accessibility, atom depth, etc.), and 2. system modelling ranging from standard bioinformatics tools to classic evolutionary models (i.e. from Miyazawa-Jernigan statistical potential to Kimura three-substitution-type model, respectively). It is shown that the algorithms based on IUPAC ambiguity codes, Boolean functions and truth table, Probabilistic Square of Opposition/Semiotic Square and Klein 4-groups-could be used for the bioinformatics analyses and Relational data modelling in natural science. Underlying mathematical, logical and semiotic concepts of interest are presented and addressed.

Linking B-factor and temperature-induced conformational transition

The crystallographic B-factor, also called temperature factor or Debye-Waller factor, has long been used as a surrogate for local protein flexibility. However, the use of the absolute B-factor as a probe for protein motion requires reproducible validation against conformational changes against chemical and physical variables. Here we report the investigation of the thermal dependence of the crystallographic B-factor and its correlation with conformational changes of the protein. We obtained the crystal protein structure coordinates and B-factors at high resolution (1.5 Å) over a broad temperature range (100 K to 325 K). The exponential thermal dependence of B-factor as a function of temperature was equal for both the diffraction intensity data (Wilson B-factor) and for all modeled atoms of the system (protein and non-protein atoms), with a thermal diffusion constant of about 0.0045 K^-1, similar for all atoms. The extrapolated B-factor at zero Kelvin (or zero-point fluctuation) varies among the atoms, although with no apparent correlation with temperature-dependent protein conformational changes. These data suggest that the thermal vibration of the atom does not necessarily correlate with the conformational dynamics of the protein.

Molecular docking and molecular dynamics simulation study on the toxicity mechanism of bongkrekic acid

BKA belongs to gram-negative brevibacterium. It can cause poisoning in humans or animals and can be fatal in severe cases. There are few investigations on toxic mechanisms of BKA because of foodborne factors. MD simulations were used to study the stability and intermolecular interactions of BKA and ANT complexes to reveal the mechanism of BKA in this paper. BKA blocked ANT protein translocation mainly through Van der Waals force, hydrophobic and hydrogen bonding interactions by the MD simulations. The conformational flexibility of the complex system during different simulation times indicated that BKA affected the conformational changes of ANT through strong interactions of hydrogen bonds with active domain residues Gln-93, Tyr-196, Arg-287 and Arg-245. The results of binding free energy, principal component analysis, hydrophobic interactions and root-mean-square fluctuation showed that the prominent binding force of Tyr-196 with C₂₆ of BKA was significant to the toxicity. The active site interactions analysis indicated that the essential positively charged polar amino acids which play a crucial role within the active site of the ANT protein undergo conformational changes with BKA as the branch point.

Thermodynamic characterization of a macrocyclic Zika virus NS2B/NS3 protease inhibitor and its acyclic analogs

Cyclization of small molecules is a widely applied strategy in drug design for ligand optimization to improve affinity, as it eliminates the putative need for structural preorganization of the ligand before binding, or to improve pharmacokinetic properties. In this work, we provide a deeper insight into the binding thermodynamics of a macrocyclic Zika virus NS2B/NS3 protease inhibitor and its linear analogs. Characterization of the thermodynamic binding profiles by isothermal titration calorimetry experiments revealed an unfavorable entropy of the macrocycle compared to the open linear reference ligands. Molecular dynamic simulations and X-ray crystal structure analysis indicated only minor benefits from macrocyclization to fixate a favorable conformation, while linear ligands retained some flexibility even in the protein-bound complex structure, possibly explaining the initially surprising effect of a higher entropic penalty for the macrocyclic ligand.

Soft disorder modulates the assembly path of protein complexes

The relationship between interactions, flexibility and disorder in proteins has been explored from many angles over the years: folding upon binding, flexibility of the core relative to the periphery, entropy changes, etc. In this work, we provide statistical evidence for the involvement of highly mobile and disordered regions in complex assembly. We ordered the entire set of X-ray crystallographic structures in the Protein Data Bank into hierarchies of progressive interactions involving identical or very similar protein chains, yielding 40205 hierarchies of protein complexes with increasing numbers of partners. We then examine them as proxies for the assembly pathways. Using this database, we show that upon oligomerisation, the new interfaces tend to be observed at residues that were characterised as softly disordered (flexible, amorphous or missing residues) in the complexes preceding them in the hierarchy. We also rule out the possibility that this correlation is just a surface effect by restricting the analysis to residues on the surface of the complexes. Interestingly, we find that the location of soft disordered residues in the sequence changes as the number of partners increases. Our results show that there is a general mechanism for protein assembly that involves soft disorder and modulates the way protein complexes are assembled. This work highlights the difficulty of predicting the structure of large protein complexes from sequence and emphasises the importance of linking predictors of soft disorder to the next generation of predictors of complex structure. Finally, we investigate the relationship between the Alphafold2’s confidence metric pLDDT for structure prediction in unbound versus bound structures, and soft disorder. We show a strong correlation between Alphafold2 low confidence residues and the union of all regions of soft disorder observed in the hierarchy. This paves the way for using the pLDDT metric as a proxy for predicting interfaces and assembly paths.

Toward Real Real-Space Refinement of Atomic Models

High-quality atomic models providing structural information are the results of their refinement versus diffraction data (reciprocal-space refinement), or versus experimental or experimentally based maps (real-space refinement). A proper real-space refinement can be achieved by comparing such a map with a map calculated from the atomic model. Similar to density distributions, the maps of a limited and even inhomogeneous resolution can also be calculated as sums of terms, known as atomic images, which are three-dimensional peaky functions surrounded by Fourier ripples. These atomic images and, consequently, the maps for the respective models, can be expressed analytically as functions of coordinates, atomic displacement parameters, and the local resolution. This work discusses the practical feasibility of such calculation for the real-space refinement of macromolecular atomic models.