Effect of Acetonitrile on the Conformation of Bovine Serum Albumin.

The use of organic solvents in drug delivery systems (DDSs) either to produce albumin nanoparticles or to manipulate the binding of target molecules to albumin, a promising nanocarrier material, presents challenges due to the conformational changes induced in the protein. In this study, we investigated the alterations in the conformation of bovine serum albumin (BSA) caused by acetonitrile (ACN) in aqueous solution by using a combination of spectroscopic analysis and molecular dynamics (MD) simulations. Ultraviolet (UV) absorption, fluorescence, and infrared (IR) absorption spectroscopy were used to analyze the BSA conformation in the solutions containing 0-60 vol % ACN.

Amide-to-ester substitution as a stable alternative to N-methylation for increasing membrane permeability in cyclic peptides.

Naturally occurring peptides with high membrane permeability often have ester bonds on their backbones. However, the impact of amide-to-ester substitutions on the membrane permeability of peptides has not been directly evaluated. Here we report the effect of amide-to-ester substitutions on the membrane permeability and conformational ensemble of cyclic peptides related to membrane permeation.

CycPeptMPDB: A Comprehensive Database of Membrane Permeability of Cyclic Peptides.

Recently, cyclic peptides have been considered breakthrough drugs because they can interact with "undruggable" targets such as intracellular protein-protein interactions. Membrane permeability is an essential indicator of oral bioavailability and intracellular targeting, and the development of membrane-permeable peptides is a bottleneck in cyclic peptide drug discovery. Although many experimental data on membrane permeability of cyclic peptides have been reported, a comprehensive database is not yet available.

Realization of the structural fluctuation of biomolecules in solution: Generalized Langevin mode analysis.

A new theoretical method, referred to as Generalized Langevin Mode Analysis (GLMA), is proposed to analyze the mode of structural fluctuations of a biomolecule in solution. The method combines the two theories in the statistical mechanics, or the Generalized Langevin theory and the RISM/3D-RISM theory, to calculate the second derivative, or the Hessian matrix, of the free energy surface of a biomolecule in aqueous solution, which consists of the intramolecular interaction among atoms in the biomolecule and the solvation free energy. The method is applied to calculate the wave-number spectrum of an alanine dipeptide in water for which the optical heterodyne-detected Raman-induced spectroscopy (RIKES) spectrum is available to compare with.

Lipid Composition Is Critical for Accurate Membrane Permeability Prediction of Cyclic Peptides by Molecular Dynamics Simulations.

Cyclic peptides have attracted attention as a promising pharmaceutical modality due to their potential to selectively inhibit previously undruggable targets, such as intracellular protein-protein interactions. Poor membrane permeability is the biggest bottleneck hindering successful drug discovery based on cyclic peptides. Therefore, the development of computational methods that can predict membrane permeability and support elucidation of the membrane permeation mechanism of drug candidate peptides is much sought after.

Large-Scale Membrane Permeability Prediction of Cyclic Peptides Crossing a Lipid Bilayer Based on Enhanced Sampling Molecular Dynamics Simulations.

Membrane permeability is a significant obstacle facing the development of cyclic peptide drugs. However, membrane permeation mechanisms are poorly understood. To investigate common features of permeable (and nonpermeable) designs, it is necessary to reproduce the membrane permeation process of cyclic peptides through the lipid bilayer.

Computational Screening of a Functional Cyclodextrin Derivative for Suppressing a Side Effect of Doxorubicin.

The binding affinity of the beta-cyclodextrin (β-CyD) derivatives with Doxorubicin (Dox) is evaluated by means of the 3D-RISM/KH theory combined with the molecular dynamics simulation in order to screen the compounds for suppressing a side-effect of the cancer drug. A protocol revised for the external and conformational entropies of the host-guest system is employed to calculate the binding free energy. It is found that the direct interactions of CyD with Dox and the desolvation free-energies of the both compounds largely cancel out to leave moderate contributions to the affinity, which are comparable to those from the entropies.

Molecular Recognition and Self-Organization in Life Phenomena Studied by a Statistical Mechanics of Molecular Liquids, the RISM/3D-RISM Theory.

There are two molecular processes that are essential for living bodies to maintain their life: the molecular recognition, and the self-organization or self-assembly. Binding of a substrate by an enzyme is an example of the molecular recognition, while the protein folding is a good example of the self-organization process. The two processes are further governed by the other two physicochemical processes: solvation and the structural fluctuation.

New Protocol for Predicting the Ligand-Binding Site and Mode Based on the 3D-RISM/KH Theory.

An efficient algorithm to find the binding position and mode of small ligands bound at an active site of protein is proposed based on the spatial distribution function (SDF) obtained from the three-dimensional reference interaction site model (3D-RISM) theory with the Kovalenko-Hirata (KH) closure relation. The ligand examined includes hydrophobic, acidic, and basic molecules and zwitterions. Eighteen different types of proteins, which serve as targets for those ligands, are selected to examine the robustness of the algorithm.

Predicting the Binding Mode of 2-Hydroxypropyl-β-cyclodextrin to Cholesterol by Means of the MD Simulation and the 3D-RISM-KH Theory.

It has been found that a cyclodextrin derivative, 2-hydroxypropyl-β-cyclodextrin (HPβCD), has reasonable therapeutic effect on Niemann-Pick disease type C, which is caused by abnormal accumulation of unesterified cholesterol and glycolipids in the lysosomes and shortage of esterified cholesterol in other cellular compartments. We study the binding affinity and mode of HPβCD with cholesterol to elucidate the possible mechanism of HPβCD for removing cholesterol from the lysosomes. The dominant binding mode of HPβCD with cholesterol is found based on the molecular dynamics simulation and a statistical mechanics theory of liquids, or the three-dimensional reference interaction site model theory with Kovalenko-Hirata closure relation.

Perspective: Structural fluctuation of protein and Anfinsen's thermodynamic hypothesis.

The thermodynamics hypothesis, casually referred to as "Anfinsen's dogma," is described theoretically in terms of a concept of the structural fluctuation of protein or the first moment (average structure) and the second moment (variance and covariance) of the structural distribution. The new theoretical concept views the unfolding and refolding processes of protein as a shift of the structural distribution induced by a thermodynamic perturbation, with the variance-covariance matrix varying. Based on the theoretical concept, a method to characterize the mechanism of folding (or unfolding) is proposed.

A Systematic Analysis of the Binding Affinity between the Pim-1 Kinase and Its Inhibitors Based on the MM/3D-RISM/KH Method.

A systematic study of the binding affinities of 16 lead compounds targeting the Pim-1 kinase based on the 3D-RISM/KH theory and MD simulations is reported. The results show a correlation coefficient R = 0.69 between the theoretical and experimental values of the binding free energy.

Predicting the binding free energy of the inclusion process of 2-hydroxypropyl-β-cyclodextrin and small molecules by means of the MM/3D-RISM method.

A protocol to calculate the binding free energy of a host-guest system is proposed based on the MM/3D-RISM method, taking cyclodextrin derivatives and their ligands as model systems. The protocol involves the procedure to identify the most probable binding mode (MPBM) of receptors and ligands by means of the umbrella sampling method. The binding free energies calculated by the MM/3D-RISM method for the complexes of the seven ligands with the MPBM of the cyclodextrin, and with the fluctuated structures around it, are in agreement with the corresponding experimental data in a semi-quantitative manner.

Topological and sequence information predict that foldons organize a partially overlapped and hierarchical structure.

It has been suggested that proteins have substructures, called foldons, which can cooperatively fold into the native structure. However, several prior investigations define foldons in various ways, citing different foldon characteristics, thereby making the concept of a foldon ambiguous. In this study, we perform a Gō model simulation and analyze the characteristics of substructures that cooperatively fold into the native-like structure.

Implication of the cause of differences in 3D structures of proteins with high sequence identity based on analyses of amino acid sequences and 3D structures.

Background: Proteins that share a high sequence homology while exhibiting drastically different 3D structures are investigated in this study. Recently, artificial proteins related to the sequences of the GA and IgG binding GB domains of human serum albumin have been designed. These artificial proteins, referred to as GA and GB, share 98% amino acid sequence identity but exhibit different 3D structures, namely, a 3α bundle versus a 4β + α structure.

Analyses of the folding properties of ferredoxin-like fold proteins by means of a coarse-grained Gō model: relationship between the free energy profiles and folding cores.

The folding mechanisms of proteins with multi-state transitions, the role of the intermediate states, and the precise mechanism how each transition occurs are significant on-going research issues. In this study, we investigate ferredoxin-like fold proteins which have a simple topology and multi-state transitions. We analyze the folding processes by means of a coarse-grained Gō model.

Incorporating into a Cα Go model the effects of geometrical restriction on Cα atoms caused by side chain orientations.

Coarse-grained Go models have been widely used for studying protein-folding mechanisms. Despite the simplicity of the model, these can reproduce the essential features of the folding process of a protein. However, it is also known that side chains significantly contribute to the folding mechanism.