Molecular dynamics simulations of lysozyme–lipid systems: probing the early steps of protein aggregation

Using the molecular dynamics simulation, the role of lipids in the lysozyme transition into the aggregation-competent conformation has been clarified. Analysis of the changes of lysozyme secondary structure upon its interactions with the model bilayer membranes composed of phosphatidylcholine and its mixtures with phosphatidylglycerol (10, 40, and 80 mol%) within the time interval of 100 ns showed that lipid-bound protein is characterized by the increased content of β-structures. Along with this, the formation of protein–lipid complexes was accompanied by the increase in the gyration radius and the decrease in RMSD of polypeptide chain. The results obtained were interpreted in terms of the partial unfolding of lysozyme molecule on the lipid matrix, with the magnitude of this effect being increased with increasing the fraction of anionic lipids. Based on the results of molecular dynamics simulation, a hypothetical model of the nucleation of lysozyme amyloid fibrils in a membrane environment was suggested.     

Molecular docking and dynamics simulations on the interaction of cationic porphyrin–anthraquinone hybrids with DNA G-quadruplexes

A series of cationic porphyrin–anthraquinone hybrids bearing either pyridine, imidazole, or pyrazole rings at the meso-positions have been investigated for their interaction with DNA G-quadruplexes by employing molecular docking and molecular dynamics simulations. Three types of DNA G-quadruplexes were utilized, which comprise parallel, antiparallel, and mixed hybrid topologies. The porphyrin hybrids have a preference to bind with parallel and mixed hybrid structures compared to the antiparallel structure. This preference arises from the end stacking of porphyrin moiety following G-stem and loop binding of anthraquinone tail, which is not found in the antiparallel due to the presence of diagonal and lateral loops that crowd the G-quartet. The binding to the antiparallel, instead, occurred with poorer affinity through both the loop and wide groove. All sites of porphyrin binding were confirmed by 6 ns molecular dynamics simulation, as well as by the negative value of the total binding free energies that were calculated using the MMPBSA method. Free energy analysis shows that the favorable contribution came from the electrostatic term, which supposedly originated from the interaction of either cationic pyridinium, pyrazole, or imidazole groups and the anionic phosphate backbone, and also from the van der Waals energy, which primarily contributed through end stacking interaction.     

Molecular glues to stabilise protein–protein interactions

Targeting protein–protein interactions (PPIs) has become a common approach to tackle various diseases whose pathobiology is driven by their mis-regulation in important signalling pathways. Modulating PPIs has tremendous untapped therapeutic potential and different approaches can be used to modulate PPIs. Initially, therapeutic effects were mostly sought by inhibiting PPIs. However, by gaining insight in the mode of action of certain therapeutic compounds, it became clear that stabilising (i.e. enhancing) PPIs can also be useful. The latter strategy is recently gaining a lot of attention, as stabilising physiologic, or even inducing novel interactions of a target protein with E3 ubiquitin ligases forms the basis of the targeted protein degradation (TPD) approach. An emerging additional example for drug discovery based on PPI stabilisation are the 14-3-3 proteins, a family of regulatory proteins, which engages in many protein–protein interactions, some of which might become therapeutical targets.     

High-throughput single-molecule studies of protein–DNA interactions

Fluorescence and force-based single-molecule studies of protein–nucleic acid interactions continue to shed critical insights into many aspects of DNA and RNA processing. As single-molecule assays are inherently low-throughput, obtaining statistically relevant datasets remains a major challenge. Additionally, most fluorescence-based single-molecule particle-tracking assays are limited to observing fluorescent proteins that are in the low-nanomolar range, as spurious background signals predominate at higher fluorophore concentrations. These technical limitations have traditionally limited the types of questions that could be addressed via single-molecule methods. In this review, we describe new approaches for high-throughput and high-concentration single-molecule biochemical studies. We conclude with a discussion of outstanding challenges for the single-molecule biologist and how these challenges can be tackled to further approach the biochemical complexity of the cell.     

Molecular dynamics simulations and contact angle of surfactant at the coal–water interface

Wettability of nonylphenol ethoxylate with four ethylene oxide groups (NP-4) on a subbituminous coal was carried out. As the concentration of NP-4 gradually increases, the contact angle firstly increases and then decreases with maximum contact angle at about critical micelle concentration (CMC) of NP-4. The monolayer adsorption behaviour of NP-4 on the model surface of Hatcher subbituminous coal was investigated by means of molecular dynamics simulations. The surfactant molecules could be detected at the water–coal interface. The water molecules are repelled and stronger hydrophobicity of the coal is obtained in the presence of NP-4, which are consistent with contact angle results at low concentration. The aggregated structure of the surfactant molecules on the coal surface in terms of head group and tail group density profiles along the perpendicular direction shows that the ethoxylate groups of the surfactant are attached at the solid surfaces. The negative interaction energy between NP-4 and the subbituminous coal surface calculated suggests that adsorption process is spontaneous. The self-diffusion coefficients results indicate that the presence of NP-4 causes higher water mobility meaning improving the hydrophobicity of low-rank coal, which is consistent with the experimental results of contact angle.     

Molecular dynamics simulations of GPCR–cholesterol interaction: An emerging paradigm

G protein-coupled receptors (GPCRs) are the largest class of molecules involved in signal transduction across cell membranes and represent major targets in the development of novel drug candidates. Membrane cholesterol plays an important role in GPCR structure and function. Molecular dynamics simulations have been successful in exploring the effect of cholesterol on the receptor and a general consensus molecular view is emerging. We review here recent molecular dynamics studies at multiple resolutions highlighting the main features of cholesterol-GPCR interaction. Several cholesterol interaction sites have been identified on the receptor that are reminiscent of nonannular sites. These cholesterol hot-spots are highly dynamic and have a microsecond time scale of exchange with the bulk lipids. A few consensus sites (such as the CRAC site) have been identified that correspond to higher cholesterol interaction. Interestingly, high plasticity is observed in the modes of cholesterol interaction and several sites have been suggested to have high cholesterol occupancy. We therefore believe that these cholesterol hot-spots are indicative of ‘high occupancy sites’ rather than ‘binding sites’. The results suggest that the energy landscape of cholesterol association with GPCRs corresponds to a series of shallow minima interconnected by low barriers. These specific interactions, along with general membrane effects, have been observed to modulate GPCR organization. Membrane cholesterol effects on receptor structure and organization, that in turn influences receptor cross-talk and drug efficacy, represent a new frontier in GPCR research. This article is part of a Special Issue entitled: Lipid-protein interactions. Guest Editors: Amitabha Chattopadhyay and Jean-Marie Ruysschaert.     

Landscape of π–π and sugar–π contacts in DNA–protein interactions

There were 1765 contacts identified between DNA nucleobases or deoxyribose and cyclic (W, H, F, Y) or acyclic (R, E, D) amino acids in 672 X-ray structures of DNA–protein complexes. In this first study to compare π-interactions between the cyclic and acyclic amino acids, visual inspection was used to categorize amino acid interactions as nucleobase π–π (according to biological edge) or deoxyribose sugar–π (according to sugar edge). Overall, 54% of contacts are nucleobase π–π interactions, which involve all amino acids, but are more common for Y, F, and R, and involve all DNA nucleobases with similar frequencies. Among binding arrangements, cyclic amino acids prefer more planar (stacked) π-systems than the acyclic counterparts. Although sugar–π interactions were only previously identified with the cyclic amino acids and were found to be less common (38%) than nucleobase–cyclic amino acid contacts, sugar–π interactions are more common than nucleobase π–π contacts for the acyclic series (61% of contacts). Similar to DNA–protein π–π interactions, sugar–π contacts most frequently involve Y and R, although all amino acids adopt many binding orientations relative to deoxyribose. These DNA–protein π-interactions stabilize biological systems, by up to approximately ?40 kJ mol^?1 for neutral nucleobase or sugar–amino acid interactions, but up to approximately ?95 kJ mol^?1 for positively or negatively charged contacts. The high frequency and strength, despite variation in structure and composition, of these π-interactions point to an important function in biological systems.     

Evaluation of protein–ligand affinity prediction using steered molecular dynamics simulations

In computational drug design, ranking a series of compound analogs in a manner that is consistent with experimental affinities remains a challenge. In this study, we evaluated the prediction of protein–ligand binding affinities using steered molecular dynamics simulations. First, we investigated the appropriate conditions for accurate predictions in these simulations. A conic harmonic restraint was applied to the system for efficient sampling of work values on the ligand unbinding pathway. We found that pulling velocity significantly influenced affinity predictions, but that the number of collectable trajectories was less influential. We identified the appropriate pulling velocity and collectable trajectories for binding affinity predictions as 1.25 Å/ns and 100, respectively, and these parameters were used to evaluate three target proteins (FK506 binding protein, trypsin, and cyclin-dependent kinase 2). For these proteins using our parameters, the accuracy of affinity prediction was higher and more stable when Jarzynski’s equality was employed compared with the second-order cumulant expansion equation of Jarzynski’s equality. Our results showed that steered molecular dynamics simulations are effective for predicting the rank order of ligands; thus, they are a potential tool for compound selection in hit-to-lead and lead optimization processes.     

Assessment and significance of protein–protein interactions during development of protein biopharmaceuticals

Early development of protein biotherapeutics using recombinant DNA technology involved progress in the areas of cloning, screening, expression and recovery/purification. As the biotechnology industry matured, resulting in marketed products, a greater emphasis was placed on development of formulations and delivery systems requiring a better understanding of the chemical and physical properties of newly developed protein drugs. Biophysical techniques such as analytical ultracentrifugation, dynamic and static light scattering, and circular dichroism were used to study protein–protein interactions during various stages of development of protein therapeutics. These studies included investigation of protein self-association in many of the early development projects including analysis of highly glycosylated proteins expressed in mammalian CHO cell cultures. Assessment of protein–protein interactions during development of an IgG1 monoclonal antibody that binds to IgE were important in understanding the pharmacokinetics and dosing for this important biotherapeutic used to treat severe allergic IgE-mediated asthma. These studies were extended to the investigation of monoclonal antibody–antigen interactions in human serum using the fluorescent detection system of the analytical ultracentrifuge. Analysis by sedimentation velocity analytical ultracentrifugation was also used to investigate competitive binding to monoclonal antibody targets. Recent development of high concentration protein formulations for subcutaneous administration of therapeutics posed challenges, which resulted in the use of dynamic and static light scattering, and preparative analytical ultracentrifugation to understand the self-association and rheological properties of concentrated monoclonal antibody solutions.     

Peptides and peptidomimetics as regulators of protein–protein interactions

1. Download : Download high-res image (114KB)
2. Download : Download full-size image

     

Rational design of the survivin/CDK4 complex by combining protein–protein docking and molecular dynamics simulations

Survivin, the smallest inhibitor of apoptosis protein (IAP), is a valid target for cancer research. It mediates both the apoptosis pathway and the cell cycle and has been proposed to form a complex with the cyclin-dependent kinase protein CDK4. The resulting complex transports CDK4 from the cytosol to the nucleus, where CDK4 participates in cell division. Survivin has been recognized as a node protein that interacts with several partners; disruption of the formed complexes can lead to new anticancer compounds. We propose a rational model of the survivin/CDK4 complex that fulfills the experimental evidence and that can be used for structure-based design of inhibitors modifying its interface recognition. In particular, the suggested complex involves the alpha helical domain of survivin and resembles the mode of binding of survivin in the survivin/borealin X-ray structure. The proposed model has been obtained by combining protein–protein docking, fractal-based shape complementarity, electrostatics studies and extensive molecular dynamics simulations.

Systems for the detection and analysis of protein–protein interactions

The analysis of protein–protein interactions is important for developing a better understanding of the functional annotations of proteins that are involved in various biochemical reactions in vivo. The discovery that a protein with an unknown function binds to a protein with a known function could provide a significant clue to the cellular pathway concerning the unknown protein. Therefore, information on protein–protein interactions obtained by the comprehensive analysis of all gene products is available for the construction of interactive networks consisting of individual protein–protein interactions, which, in turn, permit elaborate biological phenomena to be understood. Systems for detecting protein–protein interactions in vitro and in vivo have been developed, and have been modified to compensate for limitations. Using these novel approaches, comprehensive and reliable information on protein–protein interactions can be determined. Systems that permit this to be achieved are described in this review.K. Kuroda, M. Kato and J. Mima contributed equally to this work.