Density functional study of molecular interactions in secondary structures of proteins

Human Serum Albumin in the Presence of Small Platinum Nanoparticles

Noble metal materials, especially platinum nanoparticles (Pt NPs), have immense potential in nanomedicine as therapeutic agents on account of their high electron density and their high surface area. Intravenous injection is proposed as the best mode to deliver the product to patients. However, our understanding of the reaction of nanoparticles with blood components, especially proteins, is far behind the explosive development of these agents. Using synchrotron radiation circular dichroism (SRCD), we investigated the structural and stability changes of human serum albumin (HSA) upon interaction with PEG-OH coated Pt NPs at nanomolar concentrations, conditions potentially encountered for intravenous injection. There is no strong complexation found between HSA and Pt NPs. However, for the highest molar ratio of NP:HSA of 1:1, an increase of 18 °C in the thermal unfolding of HSA was observed, which is attributed to increased thermal stability of HSA generated by preferential hydration. This work proposes a new and fast method to probe the potential toxicity of nanoparticles intended for clinical use with intravenous injection.

An Innovative Methodology to Characterize, at the Molecular Scale, Interactions in Polysaccharide Aqueous Solutions

Characterizing molecular interactions at the microscopic level remains difficult and, therefore, represents a key target to better understand macromolecule and biomacromolecule behaviors in solution, alone, or in mixtures with others. Therefore, accurate characterization in liquid media, especially in aqueous solutions, without causing any perturbation of the system in which they are studied, is quite difficult. To this purpose, the present paper describes an innovative methodology based on fluorescence spectrophotometry. Two molecular fluorescent probes, namely 8-anilino-1-naphtalenesulfonic acid (ANS) and 2-benzofuryl-3-hydroxy-4(1H)-quinolone (3HQ-Bf), were selected to characterize, respectively, the dipole-dipole interactions and hydrophobic micro-domains, for the first one, and hydrogen bonding, for the second. As a support to study molecular interactions, xanthan, galactomannan, and corresponding mixtures of these substances which are well known to exhibit a synergy of interactions in well-defined mixture conditions were chosen. Once the methodology was set up, the existence of the three types of interactions in these systems was demonstrated, thus allowing the elucidation of the mechanisms of interactions at the molecular scale.

Quantum chemical studies on hydrogen bonds in helical secondary structures

We present a brief review of our recent computational studies of hydrogen bonds (H-bonds) in helical secondary structures of proteins, α-helix and 3₁₀-helix, using a Negative Fragmentation Approach with density functional theory. We found that the depolarized electronic structures of the carbonyl oxygen of the ith residue and the amide hydrogen of the (i + 4)th residue cause weaker H-bond in an α-helix than in an isolated H-bond. Our calculations showed that the H-bond energies in the 3₁₀-helix were also weaker than those of the isolated H-bonds. In the 3₁₀-helices, the adjacent N-H group at the (i + 1)th residue was closer to the C=O group of the H-bond pair than the adjacent C=O group in the 3₁₀-helices, whereas the adjacent C=O group at the (i + 1)th residue was close to the H-bond acceptor in α-helices. Therefore, the destabilization of the H-bond is attributed to the depolarization caused by the adjacent residue of the helical backbone connecting the H-bond donor and acceptor. The differences in the change in electron density revealed that such depolarizations were caused by the local electronic interactions in their neighborhood inside the helical structure and redistributed the electron density. We also present the improvements in the force field of classical molecular simulation, based on our findings.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12551-022-01034-5.

Computational biophysics and structural biology of proteins-a Special Issue in honor of Prof. Haruki Nakamura's 70th birthday

Receiving his initial training jointly in theoretical and applied physics at the University of Tokyo, Professor Haruki Nakamura has had a long and eventful scientific career, along the way helping to shape the way that biophysics is carried out in Japan. Concentrating his research efforts on the simulation of protein structure and function, he has, over his career arc, acted as director of the Institute for Protein Research (Osaka, Japan), director of the Protein Data Bank of Japan (PDBj), president of the Biophysical Society of Japan (BSJ), president of the Protein Science Society of Japan (PSSJ), and group leader and professor of Bioinformatics and Computational Structural Biology at Osaka University. In 2022, Prof. Haruki Nakamura turned 70 years old, and to mark this occasion, his scientific colleagues from around the world have combined their efforts to produce this Festschrift Issue of the IUPAB Biophysical Reviews journal around the theme of the computational biophysics and structural biology of proteins.

Depolarizing Effects in Hydrogen Bond Energy in 310-Helices Revealed by Quantum Chemical Analysis

Hydrogen-bond (H-bond) energies in 3₁₀-helices of short alanine peptides were systematically examined by precise DFT calculations with the negative fragmentation approach (NFA), a modified method based on the molecular tailoring approach. The contribution of each H-bond was evaluated in detail from the 3₁₀-helical conformation of total energies (whole helical model, WH_3-10 model), and the results were compared with the property of H-bond in α-helix from our previous study. The H-bond energies of the WH_3-10 model exhibited tendencies different from those exhibited by the α-helix in that they depended on the helical position of the relevant H-bond pair. H-bond pairs adjacent to the terminal H-bond pairs were observed to be strongly destabilized. The analysis of electronic structures indicated that structural characteristics cause the destabilization of the H-bond in 3₁₀-helices. We also found that the longer the helix length, the more stable the H-bond in the terminal pairs of the WH_3-10 model, suggesting the action of H-bond cooperativity.

The structure and properties of MFG-E8 and the In vitro assessment of its toxic effects on myoblast cells

Four high-molecular-weight protein fractions of milk fat globule membrane (MFGM) were isolated from bovine milk. Sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE), MALDI-TOF/TOF™ and Liquid chromatography electrospray ionization-tandem mass spectrometry (LC-ESI-MS/MS) were used to measure the molecular sizes of the MFGM. Fourier transform infrared spectroscopy (FT-IR) and circular dichroism (CD) were performed to determine the conformations of the MFGM. The results showed that the main protein (98.33%) in MFGM protein fraction 2 was Milk fat globule epidermal growth factor-VIII (MFG-E8), with a molecular weight of 47.82 kDa. The secondary structural component measurements showed that the MFG-E8 consisted of 5% helix, 70% sheet and 25% random coil, and the results matched well with the prediction by SSPro 5.1 bioinformatic analysis. The thermograms analysis revealed that Td and△H of MFG-E8 were 60.50°Cand 132.29 kJ/mol. The in vitro digestibility of MFG-E8 showed that it can be enzymatically hydrolyzed in the stomach and relatively stable in the intestinal fluid. The in vitro C₂C₁₂ and Caco2 cell activity tests indicated that MFG-E8 promoted the proliferation of C₂C₁₂ myoblast cells without cytotoxicity. The biological functional properties of MFG-E8 may be related to the fact that MFG-E8 possesses a high level of β-sheet structure. Our results suggested that MFG-E8 possesses broad prospects not only for use in functional food products but also as a source of natural anti-sarcopenia drugs.

A DFT study of structure and stability of pleated and rippled cross-β sheets with hydrophobic sidechains

The rippled cross-β sheet, a topography, in which mirror-image peptides are arranged with alternating chirality into a periodic two-dimensional network, is burgeoning as a new design principle for materials and biomedical applications. Experiments by the Schneider, Nilsson, and Raskatov labs have independently shown diverse racemic mixtures of aggregation-prone peptide of different sizes to favor the rippled over the pleated topography. Yet, systematic ab initio studies are lacking, and the field is yet to develop rules that would enable the design of new rippled cross-β frameworks from first principles. Here, DFT calculations were performed on a set of model systems, designed to begin understanding the impact that bulky, hydrophobic sidechains have upon the formation of pleated and rippled cross-β frameworks. It is hoped that this study will help stimulate the development of a predictive, general framework to enable rational design of rippled cross-β sheets in the future.

Hydrogen bond donors and acceptors are generally depolarized in α-helices as revealed by a molecular tailoring approach

Hydrogen-bond (H-bond) interaction energies in α-helices of short alanine peptides were systematically examined by precise density functional theory calculations, followed by a molecular tailoring approach. The contribution of each H-bond interaction in α-helices was estimated in detail from the entire conformation energies, and the results were compared with those in the minimal H-bond models, in which only H-bond donors and acceptors exist with the capping methyl groups. The former interaction energies were always significantly weaker than the latter energies, when the same geometries of the H-bond donors and acceptors were applied. The chemical origin of this phenomenon was investigated by analyzing the differences among the electronic structures of the local peptide backbones of the α-helices and those of the minimal H-bond models. Consequently, we found that the reduced H-bond energy originated from the depolarizations of both the H-bond donor and acceptor groups, due to the repulsive interactions with the neighboring polar peptide groups in the α-helix backbone. The classical force fields provide similar H-bond energies to those in the minimal H-bond models, which ignore the current depolarization effect, and thus they overestimate the actual H-bond energies in α-helices. © 2019 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.

Studies on Molecular Dynamics of Intrinsically Disordered Proteins and Their Fuzzy Complexes: A Mini-Review

The molecular dynamics (MD) method is a promising approach toward elucidating the molecular mechanisms of intrinsically disordered regions (IDRs) of proteins and their fuzzy complexes. This mini-review introduces recent studies that apply MD simulations to investigate the molecular recognition of IDRs. Firstly, methodological issues by which MD simulations treat IDRs, such as developing force fields, treating periodic boundary conditions, and enhanced sampling approaches, are discussed. Then, several examples of the applications of MD to investigate molecular interactions of IDRs in terms of the two kinds of complex formations; coupled-folding and binding and fuzzy complex. MD simulations provide insight into the molecular mechanisms of these binding processes by sampling conformational ensembles of flexible IDRs. In particular, we focused on all-atom explicit-solvent MD simulations except for studies of higher-order assembly of IDRs. Recent advances in MD methods, and computational power make it possible to dissect the molecular details of realistic molecular systems involving the dynamic behavior of IDRs.

Toward Ab Initio Protein Folding: Inherent Secondary Structure Propensity of Short Peptides from the Bioinformatics and Quantum-Chemical Perspective

By combining bioinformatics with quantum-chemical calculations, we attempt to address quantitatively some of the physical principles underlying protein folding. The former allowed us to identify tripeptide sequences in existing protein three-dimensional structures with a strong preference for either helical or extended structure. The selected representatives of pro-helical and pro-extended sequences were converted into "isolated" tripeptides-capped at N- and C-termini-and these were subjected to an extensive conformational sampling and geometry optimization (typically thousands to tens of thousands of conformers for each tripeptide). For each conformer, the QM(DFT-D3)/COSMO-RS free-energy value was then calculated, G_conf(solv). The Δ G_conf(solv) is expected to provide an objective, unbiased, and quantitatively accurate measure of the conformational preference of the particular tripeptide sequence. It has been shown that irrespective of the helical vs extended preferences of the selected tripeptide sequences in context of the protein, most of the low-energy conformers of isolated tripeptides prefer the R-helical structure. Nevertheless, pro-helical tripeptides show slightly stronger helix preference than their pro-extended counterparts. Furthermore, when the sampling is repeated in the presence of a partner tripeptide to mimic the situation in a β-sheet, pro-extended tripeptides (exemplified by the VIV) show a larger free-energy benefit than pro-helical tripeptides (exemplified by the EAM). This effect is even more pronounced in a hydrophobic solvent, which mimics the less polar parts of a protein. This is in line with our bioinformatic results showing that the majority of pro-extended tripeptides are hydrophobic. The preference for a specific secondary structure by the studied tripeptides is thus governed by the plasticity to adopt to its environment. In addition, we show that most of the "naturally occurring" conformations of tripeptide sequences, i.e., those found in existing three-dimensional protein structures, are within ∼10 kcal·mol^-1 from their global minima. In summary, our "ab initio" data suggest that complex protein structures may start to emerge already at the level of their small oligopeptidic units, which is in line with a hierarchical nature of protein folding.

Unfolding of α-helical 20-residue poly-glutamic acid analyzed by multiple runs of canonical molecular dynamics simulations

Elucidating the molecular mechanism of helix-coil transitions of short peptides is a long-standing conundrum in physical chemistry. Although the helix-coil transitions of poly-glutamic acid (PGA) have been extensively studied, the molecular details of its unfolding process still remain unclear. We performed all-atom canonical molecular dynamics simulations for a 20-residue PGA, over a total of 19 μs, in order to investigate its helix-unfolding processes in atomic resolution. Among the 28 simulations, starting with the α-helical conformation, all showed an unfolding process triggered by the unwinding of terminal residues, rather than by kinking and unwinding of the middle region of the chain. The helix-coil-helix conformation which is speculated by the previous experiments was not observed. Upon comparison between the N- and C-termini, the latter tended to be unstable and easily unfolded. While the probabilities of helix elongation were almost the same among the N-terminal, middle, and C-terminal regions of the chain, unwinding of the helix was enriched at the C-terminal region. The turn and 3₁₀-helix conformations were kinetic intermediates in the formation and deformation of α-helix, consistent with the previous computational studies for Ala-based peptides.

Computational biochemical investigation of the binding energy interactions between an estrogen receptor and its agonists

We present the energy profiles of estrogen receptor–agonist ligand interactions in atomic detail using a quantum biochemical approach.