Rational peptide design for regulating liquid-liquid phase separation on the basis of residue-residue contact energy

Since liquid-liquid phase separation (LLPS) of proteins is governed by their intrinsically disordered regions (IDRs), it can be controlled by LLPS-regulators that bind to the IDRs. The artificial design of LLPS-regulators based on this mechanism can be leveraged in biological and therapeutic applications. However, the fabrication of artificial LLPS-regulators remains challenging. Peptides are promising candidates for artificial LLPS-regulators because of their ability to potentially bind to IDRs complementarily. In this study, we provide a rational peptide design methodology for targeting IDRs based on residue-residue contact energy obtained using molecular dynamics (MD) simulations. This methodology provides rational peptide sequences that function as LLPS regulators. The peptides designed with the MD-based contact energy showed dissociation constants of 35-280 nM for the N-terminal IDR of the tumor suppressor p53, which are significantly lower than the dissociation constants of peptides designed with the conventional 3D structure-based energy, demonstrating the validity of the present peptide design methodology. Importantly, all of the designed peptides enhanced p53 droplet formation. The droplet-forming peptides were converted to droplet-deforming peptides by fusing maltose-binding protein (a soluble tag) to the designed peptides. Thus, the present peptide design methodology for targeting IDRs is useful for regulating droplet formation.

Prediction and evaluation of the 3D structure of Macadamia integrifolia antimicrobial protein 2 (MiAMP2) and its interaction with palmitoleic acid or oleic acid: An integrated computational approach

This study investigated the physicochemical properties and 3D structure of Macadamia integrifolia antimicrobial protein 2 (MiAMP2) and its interaction with palmitoleic acid (POA) or oleic acid (OA) in macadamia oil. The 3D structure of MiAMP2 was constructed for the first time by ab initio modelling using the TrRosetta server. The results showed that MiAMP2 was highly hydrophilic and had seven disulfide bonds and higher α-helix and β-sheet/turn contents. Molecular simulation showed that the hydrophobic pocket of MiAMP2 created a favourable environment for the binding of POA and OA. Free energy landscape and independent gradient model (IGM) analyses revealed that hydrogen bonds and van der Waals forces were the major driving forces stabilizing complexes formed by MiAMP2 and POA or OA. The present study provides a theoretical basis and new insight for the future development and utilization of macadamia nut protein in the food industry.

Modelling of noble anaesthetic gases and high hydrostatic pressure effects in lipid bilayers

Our objective was to study molecular processes that might be responsible for inert gas narcosis and high-pressure nervous syndrome. The classical molecular dynamics trajectories (200 ns) of dioleoylphosphatidylcholine (DOPC) bilayers simulated by the Berger force field were evaluated for water and the atomic distribution of noble gases around DOPC molecules in the pressure range of 1-1000 bar and at a temperature of 310 K. Xenon and argon have been tested as model gases for general anaesthetics, and neon has been investigated for distortions that are potentially responsible for neurological tremors in hyperbaric conditions. The analysis of stacked radial pair distribution functions of DOPC headgroup atoms revealed the explicit solvation potential of the gas molecules, which correlates with their dimensions. The orientational dynamics of water molecules at the biomolecular interface should be considered as an influential factor, while excessive solvation effects appearing in the lumen of membrane-embedded ion channels could be a possible cause of inert gas narcosis. All the noble gases tested exhibit similar order parameter patterns for both DOPC acyl chains, which are opposite of the patterns found for the order parameter curve at high hydrostatic pressures in intact bilayers. This finding supports the 'critical volume' hypothesis of anaesthesia pressure reversal. The irregular lipid headgroup-water boundary observed in DOPC bilayers saturated with neon in the pressure range of 1-100 bar could be associated with the possible manifestation of neurological tremors at the atomic scale. The non-immobiliser neon also demonstrated the highest momentum impact on the normal component of the DOPC diffusion coefficient representing the monolayer undulation rate, which indicates that enhanced diffusivity rather than atomic size is the key factor.

Conformational changes of globular proteins upon adsorption on a hydrophobic surface

Coarse-grained Monte Carlo simulations are used to study thermal denaturation of small globular proteins adsorbed on a hydrophobic surface. Though helices are more stable than sheets, they are highly deformed in the adsorbed protein.

Random coil to globular thermal response of a protein (H3.1) with three knowledge-based coarse-grained potentials

The effect of temperature on the conformation of a histone (H3.1) is studied by a coarse-grained Monte Carlo simulation based on three knowledge-based contact potentials (MJ, BT, BFKV). Despite unique energy and mobility profiles of its residues, the histone H3.1 undergoes a systematic (possibly continuous) structural transition from a random coil to a globular conformation on reducing the temperature. The range over which such a systematic response in variation of the radius of gyration (R(g)) with the temperature (T) occurs, however, depends on the potential, i.e. ΔT(MJ) ≈ 0.013-0.020, ΔT(BT) ≈ 0.018-0.026, and ΔT(BFKV) ≈ 0.006-0.013 (in reduced unit). Unlike MJ and BT potentials, results from the BFKV potential show an anomaly where the magnitude of R(g) decreases on raising the temperature in a range ΔT(A) ≈ 0.015-0.018 before reaching its steady-state random coil configuration. Scaling of the structure factor, S(q) ∝ q(-1/ν), with the wave vector, q=2π/λ, and the wavelength, λ, reveals a systematic change in the effective dimension (D(e)∼1/ν) of the histone with all potentials (MJ, BT, BFKV): D(e)∼3 in the globular structure with D(e)∼2 for the random coil. Reproducibility of the general yet unique (monotonic) structural transition of the protein H3.1 with the temperature (in contrast to non-monotonic structural response of a similar but different protein H2AX) with three interaction sets shows that the knowledge-based contact potential is viable tool to investigate structural response of proteins. Caution should be exercise with the quantitative comparisons due to differences in transition regimes with these interactions.