Electrostatic calculations and model-building suggest that DNA bound to CAP is sharply bent

Two observations suggest that DNA, upon binding to E. coli catabolite gene activator protein (CAP), is sharply bent by a total angle of at least 100-150 degrees: (1) The electrostatic potential field of CAP shows regions of positive potential that form a ramp on 3 sides of the protein. (2) The DNA binding site size as determined by DNA ethylation interference with binding, (Majors: "Control of the E. coli Lac Operon at the Molecular Level." Ph.D. Thesis, Harvard University, Cambridge, 1977) and by relative affinities of DNA fragments of various lengths (Liu-Johnson et al.: Cell 47:995-1005, 1986) requires severe bending of the DNA to maintain its favorable electrostatic contact with the protein.     

Optimization of the electrostatic interactions in proteins of different functional and folding type

The 3-dimensional optimization of the electrostatic interactions between the charged amino acid residues was studied by Monte Carlo simulations on an extended representative set of 141 protein structures with known atomic coordinates. The proteins were classified by different functional and structural criteria, and the optimization of the electrostatic interactions was analyzed. The optimization parameters were obtained by comparison of the contribution of charge-charge interactions to the free energy of the native protein structures and for a large number of randomly distributed charge constellations obtained by the Monte Carlo technique. On the basis of the results obtained, one can conclude that the charge-charge interactions are better optimized in the enzymes than in the proteins without enzymatic functions. Proteins that belong to the mixed αβ folding type are electrostatically better optimized than pure α-helical or β-strand structures. Proteins that are stabilized by disulfide bonds show a lower degree of electrostatic optimization. The electrostatic interactions in a native protein are effectively optimized by rejection of the conformers that lead to repulsive charge-charge interactions. Particularly, the rejection of the repulsive contacts seems to be a major goal in the protein folding process. The dependence of the optimization parameters on the choice of the potential function was tested. The majority of the potential functions gave practically identical results.     

RosettaDDGPrediction for high-throughput mutational scans: From stability to binding

Reliable prediction of free energy changes upon amino acid substitutions (ΔΔGs) is crucial to investigate their impact on protein stability and protein–protein interaction. Advances in experimental mutational scans allow high-throughput studies thanks to multiplex techniques. On the other hand, genomics initiatives provide a large amount of data on disease-related variants that can benefit from analyses with structure-based methods. Therefore, the computational field should keep the same pace and provide new tools for fast and accurate high-throughput ΔΔG calculations. In this context, the Rosetta modeling suite implements effective approaches to predict folding/unfolding ΔΔGs in a protein monomer upon amino acid substitutions and calculate the changes in binding free energy in protein complexes. However, their application can be challenging to users without extensive experience with Rosetta. Furthermore, Rosetta protocols for ΔΔG prediction are designed considering one variant at a time, making the setup of high-throughput screenings cumbersome. For these reasons, we devised RosettaDDGPrediction, a customizable Python wrapper designed to run free energy calculations on a set of amino acid substitutions using Rosetta protocols with little intervention from the user. Moreover, RosettaDDGPrediction assists with checking completed runs and aggregates raw data for multiple variants, as well as generates publication-ready graphics. We showed the potential of the tool in four case studies, including variants of uncertain significance in childhood cancer, proteins with known experimental unfolding ΔΔGs values, interactions between target proteins and disordered motifs, and phosphomimetics. RosettaDDGPrediction is available, free of charge and under GNU General Public License v3.0, at https://github.com/ELELAB/RosettaDDGPrediction .     

Treatment of NOE constraints involving equivalent or nonstereoassigned protons in calculations of biomacromolecular structures

Summary Two modifications to the commonly used protocols for calculating NMR structures are developed, relating to the treatment of NOE constraints involving groups of equivalent protons or nonstereoassigned diastereotopic protons. Firstly, a modified method is investigated for correcting for multiplicity, which is applicable whenever all NOE intensities are calibrated as a single set and categorised in broad intensity ranges. Secondly, a new set of values for pseudoatom corrections is proposed for use with calculations employing centre-averaging. The effect of these protocols on structure calculations is demonstrated using two proteins, one of which is well defined by the NOE data, the other less so. It is shown that failure to correct for multiplicity when using r^-6 averaging results in overly precise structures, higher NOE energies and deviations from geometric ideality, while failure to correct for multiplicity when using r^-6 summation can cause an avoidable degradation of precision if the NOE data are sparse. Conversely, when multiplicities are treated correctly, r^-6 averaging, r^-6 summation and centre averaging all give closely comparable results when the structure is well defined by the data. When the NOE data contain less information, r^-6 averaging or r^-6 summation offer a significant advantage over centre averaging, both in terms of precision and in terms of the proportion of calculations that converge on a consisten result.Abbreviations HMG high mobility group - NOE nuclear Overhauser enhancement - NOESY nuclear Overhauser enhancement spectroscopy - rmsd root-mean-square deviation - YASAP yet another simulated-annealing protocol     

Theoretical Study of the TiO2 and MgO Surface Acidity and the Adsorption of Acids and Bases

We present ab-initio periodic Hartree–Fock calculations (crystal program) of small molecules on TiO₂ and MgO. The adsorption of the molecules may be molecular or dissociative. This depends on their acid and basic properties in the gas phase. For the molecular adsorption, the molecules are adsorbed as bases on Ti(+IV) sites, the adsorption energies correlate with the proton affinities. The dissociations on the surface correlate with the gas phase cleavages: thus, the dissociation of MeOH leads to a preferential basic cleavage (the fragment HO– is adsorbed on a Ti⁺⁴ ion and the fragment Me⁺ is adsorbed on a O²– ion of the oxide). The opposite result is obtained with MeSH. Another important factor is the adsorbate–adsorbate interaction: favorable cases are a sequence of H-bonds for the hydroxyl groups resulting from the water dissociation and the mode of adsorption for the ammonium ions. Lateral interactions also force the adsorbed CO₂ molecules to bend over the surface so that their mutual orientation resembles the geometry of the CO₂ dimer. With respect to water adsorption, MgO appears to be a basic oxide. As experimentally observed, NH₃ adsorbs preferentially on TiO₂ and CO₂ on MgO. However, this difference of reactivity should not be expressed in terms of acid vs. basic behaviour but in terms of hard and soft acidity. The MgO surface is a 'soft' acidic surface that reacts preferentially with the soft base, CO₂.     

Chemical shifts and three-dimensional protein structures

Summary During the past three years it has become possible to compute ab initio the ¹³C, ¹⁵N and ¹⁹F NMR chemical shifts of many sites in native proteins. Chemical shifts are beginning to become a useful supplement to more established methods of solution structure determination, and may find utility in solid-state analysis as well. From ¹³C NMR, information on , and torsions can be obtained, permitting both assignment verification, and structure refinement and prediction. For ¹⁵N, both torsional and hydrogen-bonding effects are important, while for ¹⁹F, chemical shifts are primarily indicators of the local charge field. Chemical shift calculations are still slow, but shielding hypersurfaces — the shift as a function of the dihedral angles that define the molecular conformation — are becoming accessible. Over the next few years, theoretical and computer hardware improvements will enable more routine use of chemical shifts in structural studies, including the study of metal-ligand interactions, the analysis of drug and substrate binding and catalysis, the study of folding/unfolding pathways, as well as the characterization of conformational substates. Rather than simply being a necessary prerequisite for multidimensional NMR, chemical shifts and chemical shift non-equivalence due to folding are now beginning to be useful for structural characterization.     

Pseudo-disesquiterpenoids from seeds of Vernonia anthelmintica and their biological activities

Three new pseudo-disesquiterpenoids, vernodalidimer F – H (1 – 3), which are formed by esterification of two sesquiterpenoids, along with a new sesquiterpenoid (4) and two known sesquiterpeonids (5, 6) were isolated from the seeds of Vernonia anthelmintica. Their structures were elucidated by NMR data. The absolute configurations of 1–3 and 5 were determined by comparison of the experimental and calculated electronic circular dichroism spectra. Cytotoxicity of the isolated compounds against four human tumor cell lines were assayed. 5 exhibited strong cytotoxicity against HCT-15, PC-3, A549 and Hela cells lines with IC₅₀ values of 5.3, 5.6, 6.2, and 8.2 μM, respectively. 2 showed non-concentration dependent cytotoxicity against HCT-15, PC-3, and A549 cells lines with inhibition rate of 56.1%, 55.3%, and 50.1%, respectively. 1 and 3 showed moderate cytotoxicity against four cell lines with IC₅₀ values ranging from 12.2 ± 5.1 to 28.6 ± 2.5 μM. The influence of melanin content in B16 melanoma cells of 1, 5, and 6 were tested, and they increased melanin content by 43.6%, 28.1%, and 37.0% higher than positive control 8-methoxypsoralen.     

Calculation of substrate binding affinities for a bacterial GH78 rhamnosidase through molecular dynamics simulations

Ram2 from Pediococcus acidilactici is a rhamnosidase from the glycoside hydrolase family 78. It shows remarkable selectivity for rutinose rather than para-nitrophenyl-alpha-l-rhamnopyranoside (p-NPR). Molecular dynamics simulations were performed using a homology model of this enzyme, in complex with both substrates. Free energy calculations lead to predicted binding affinities of −34.4 and −30.6 kJ mol⁻¹ respectively, agreeing well with an experimentally estimated relative free energy of 5.4 kJ mol⁻¹. Further, the most relevant binding poses could be determined. While p-NPR preferably orients its rhamnose moiety toward the active site, rutinose interacts most strongly with its glucose moiety. A detailed hydrogen bond analysis confirms previously implicated residues in the active site (Asp217, Asp222, Trp226, Asp229 and Glu488) and quantifies the importance of individual residues for the binding. The most important amino acids are Asp229 and Phe339 which are involved in many interactions during the simulations. While Phe339 was observed in more simulations, Asp229 was involved in more persistent interactions (forming an average of at least 2 hydrogen bonds during the simulation). These analyses directly suggest mutations that could be used in a further experimental characterization of the enzyme. This study shows once more the strength of computer simulations to rationalize and guide experiments at an atomic level.     

Activation of Ni Particles into Single Ni–N Atoms for Efficient Electrochemical Reduction of CO2

Electrochemical reduction of carbon dioxide (CO₂) to fuels and value‐added industrial chemicals is a promising strategy for keeping a healthy balance between energy supply and net carbon emissions. Here, the facile transformation of residual Ni particle catalysts in carbon nanotubes into thermally stable single Ni atoms with a possible NiN₃ moiety is reported, surrounded with a porous N‐doped carbon sheath through a one‐step nanoconfined pyrolysis strategy. These structural changes are confirmed by X‐ray absorption fine structure analysis and density functional theory (DFT) calculations. The dispersed Ni single atoms facilitate highly efficient electrocatalytic CO₂ reduction at low overpotentials to yield CO, providing a CO faradaic efficiency exceeding 90%, turnover frequency approaching 12 000 h^?1, and metal mass activity reaching about 10 600 mA mg^?1, outperforming current state‐of‐the‐art single atom catalysts for CO₂ reduction to CO. DFT calculations suggest that the Ni@N₃ (pyrrolic) site favors *COOH formation with lower free energy than Ni@N₄, in addition to exothermic CO desorption, hence enhancing electrocatalytic CO₂ conversion. This finding provides a simple, scalable, and promising route for the preparation of low‐cost, abundant, and highly active single atom catalysts, benefiting future practical CO₂ electrolysis.