Retrospective molecular docking study of WY-25105 ligand to β-secretase and bias of the three-dimensional structure flexibility

β-Secretase (BACE) is a very promising target in the search for a treatment for Alzheimer’s disease using a protein–ligand inhibition approach. Given the many published X-ray structures of BACE protein, structure-based drug design has been used extensively to support new inhibitor discovery programs. Due to the high flexibility and large catalytic site of this protein, sampling of the huge conformational space of the binding site is the big challenge to overcome and is the main limitation of the most widely used docking programs. Incorrect treatment of these pitfalls can introduce bias into ligand docking and could affect the results. This is especially the case with the WY-25105 compound reported by the Wyeth Corporation as a BACE ligand that did not fit into any of the known crystal structures. In the present retrospective study, a set of available X-ray enzyme structures was selected and molecular dynamics simulations were conducted to generate more diverse representative BACE protein conformations. These conformations were then used for a docking study of the WY-25105 compound. The results confirmed the need to use an ensemble of structures in protein–ligand docking for identification of new binding modes in structure-based drug design of BACE inhibitors.

A flexible-protein molecular docking study of the binding of ruthenium complex compounds to PIM1, GSK-3β, and CDK2/Cyclin A protein kinases

We employ ensemble docking simulations to characterize the interactions of two enantiomeric forms of a Ru-complex compound (1-R and 1-S) with three protein kinases, namely PIM1, GSK-3β, and CDK2/cyclin A. We show that our ensemble docking computational protocol adequately models the structural features of these interactions and discriminates between competing conformational clusters of ligand-bound protein structures. Using the determined X-ray crystal structure of PIM1 complexed to the compound 1-R as a control, we discuss the importance of including the protein flexibility inherent in the ensemble docking protocol, for the accuracy of the structure prediction of the bound state. A comparison of our ensemble docking results suggests that PIM1 and GSK-3β bind the two enantiomers in similar fashion, through two primary binding modes: conformation I, which is very similar to the conformation presented in the existing PIM1/compound 1-R crystal structure; conformation II, which represents a 180° flip about an axis through the NH group of the pyridocarbazole moiety, relative to conformation I. In contrast, the binding of the enantiomers to CDK2 is found to have a different structural profile including a suggested bound conformation, which lacks the conserved hydrogen bond between the kinase and the ligand (i.e., ATP, staurosporine, Ru-complex compound). The top scoring conformation of the inhibitor bound to CDK2 is not present among the top-scoring conformations of the inhibitor bound to either PIM1 or GSK-3β and vice-versa. Collectively, our results help provide atomic-level insights into inhibitor selectivity among the three kinases.

Stereoselectivity of chalcone isomerase with chalcone derivatives: a computational study

Chalcone isomerase (CHI) catalyzes the intramolecular cyclization of chalcones into flavonoids. The activity of CHI is essential for the biosynthesis of flavonoids precursors of floral pigments and phenylpropanoid plant defense compounds. In the present study, we explored the detailed binding structures and binding free energies for two different active site conformations of CHI with s-cis/s-trans conformers of three chalcone compounds by performing molecular dynamics (MD) simulations and binding free energy calculations. The computational results indicate that s-cis/s-trans conformers of chalcone compounds are orientated in the similar binding position in the active site of CHI and stabilized by the different first hydrogen bond network and the same second hydrogen bond network. The first hydrogen bond network results in much lower binding affinity of s-trans conformer of chalcone compound with CHI than that of s-cis conformer. The conformational change of the active site residue T48 from indirectly interacting with the substrate via the second hydrogen bond network to directly forming the hydrogen bond with the substrates cannot affect the binding mode of both conformers of chalcone compounds, but remarkably improves the binding affinity. These results show that CHI has a strong stereoselectivity. The calculated binding free energies for three chalcone compounds with CHI are consistent with the experimental activity data. In addition, several valuable insights are suggested for future rational design and discovery of high-efficiency mutants of CHI.

Entropy versus aromaticity in the conformational dynamics of aromatic rings

Comparison of the results of Car-Parrinello molecular dynamics simulations of isolated benzene, pyrimidine and 1,2,4-triazine molecules reveals that the unusually low population of planar geometry of the benzene ring is caused by entropy effects despite its high aromaticity. The decrease in symmetry of the molecule results in smaller changes in entropy and Gibbs free energy due to out-of-plane deformations of the ring, leading to an increase in the population of planar geometry of the ring. This leads to differences in the topology of potential energy and Gibbs free energy surfaces.

Coarse-grained simulations for organic molecular liquids based on Gay-Berne and electric multipole potentials

Coarse-grained studies of CH₃SH, CH₃CHO and CHCl₃ liquids, based on anisotropic Gay-Berne (GB) and electric multipole potentials (EMP), demonstrate that the coarse-grained model is able to qualitatively reproduce the results obtained from the atomistic model (AMOEBA polarizable force field) and allows for significant saving in computation time. It should be pointed out that the accuracy of the coarse-grained model is very sensitive to how well the anisotropic GB particle is defined and how satisfactorily the EMP sites are chosen.

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.

Molecular dynamics analysis of a series of 22 potential farnesyltransferase substrates containing a CaaX-motif

Protein farnesyltransferase (FTase) is an important target in many research fields, more markedly so in cancer investigation since several proteins known to be involved in human cancer development are thought to serve as substrates for FTase and to require farnesylation for proper biological activity. Several FTase inhibitors (FTIs) have advanced into clinical testing. Nevertheless, despite the progress in the field several functional and mechanistic doubts on the FTase catalytic activity have persisted. This work provides some crucial information on this important enzyme by describing the application of molecular dynamics simulations using specifically designed molecular mechanical parameters for a variety of 22 CaaX peptides known to work as natural substrates or inhibitors for this enzyme. The study involves a comparative analysis of several important molecular aspects, at the mechanistic level, of the behavior of substrates and inhibitors at the dynamic level, including the behavior of the enzyme and peptides, as well as their interaction, together with the effect of the solvent. Properties evaluated include the radial distribution function of the water molecules around the catalytically important zinc metal atom and cysteine sulfur of CaaX, the conformations of the substrate and inhibitor and the corresponding RMSF values, critical hydrogen bonds, and several catalytically relevant distances. These results are discussed in light of recent experimental and computational evidence that provides new insights into the activity of this enzyme.

Insights into the structural determinants for selective inhibition of nitric oxide synthase isoforms

Selective inhibition of the nitric oxide synthase isoforms (NOS) is a promising approach for the treatment of various disorders. However, given the high active site conservation among all NOS isoforms, the design of selective inhibitors is a challenging task. Analysis of the X-ray crystal structures of the NOS isoforms complexed with known inhibitors most often gives no clues about the structural determinants behind the selective inhibition since the inhibitors share the same binding conformation. Aimed at a better understanding of the structural factors responsible for selective inhibition of NOS isoforms we have performed MD simulations for iNOS, nNOS and eNOS complexed with N^ω-NO₂-L-Arg (1), and with the aminopyridine derivatives 2 and 3. The slightly better selectivity of 1 for nNOS may be assigned to the presence of extra charge–charge interactions due to its “extended” conformation. While the high affinity of 2 for iNOS can be explained by the formation of an iNOS-specific subpocket upon binding, the lack of affinity for eNOS is associated to a conformational change in Glu363. The strong van der Waals and electrostatic interactions between 3 and the active site of nNOS are most likely responsible for its higher affinity for this isoform. Owing to the elongated and narrow binding pocket of iNOS, the correct positioning of 3 over the heme group is difficult, which may account for its lower affinity toward this isoform. Brought together, our results might help to rationalize the design of selective NOS inhibitors.

Partition decoupling for multi-gene analysis of gene expression profiling data

Background

Molecular dynamics (MD) simulations are powerful tools to investigate the conformational dynamics of proteins that is often a critical element of their function. Identification of functionally relevant conformations is generally done clustering the large ensemble of structures that are generated. Recently, Self-Organising Maps (SOMs) were reported performing more accurately and providing more consistent results than traditional clustering algorithms in various data mining problems. We present a novel strategy to analyse and compare conformational ensembles of protein domains using a two-level approach that combines SOMs and hierarchical clustering.

Results

The conformational dynamics of the α-spectrin SH3 protein domain and six single mutants were analysed by MD simulations. The Cα's Cartesian coordinates of conformations sampled in the essential space were used as input data vectors for SOM training, then complete linkage clustering was performed on the SOM prototype vectors. A specific protocol to optimize a SOM for structural ensembles was proposed: the optimal SOM was selected by means of a Taguchi experimental design plan applied to different data sets, and the optimal sampling rate of the MD trajectory was selected. The proposed two-level approach was applied to single trajectories of the SH3 domain independently as well as to groups of them at the same time. The results demonstrated the potential of this approach in the analysis of large ensembles of molecular structures: the possibility of producing a topological mapping of the conformational space in a simple 2D visualisation, as well as of effectively highlighting differences in the conformational dynamics directly related to biological functions.

Conclusions

The use of a two-level approach combining SOMs and hierarchical clustering for conformational analysis of structural ensembles of proteins was proposed. It can easily be extended to other study cases and to conformational ensembles from other sources.     

Exploration of conformational changes in lactose permease upon sugar binding and proton transfer through coarse‐grained simulations

Escherichia coli lactose permease (LacY) actively transports lactose and other galactosides across cell membranes through lactose/H⁺ symport process. Lactose/H⁺ symport is a highly complex process that involves sugar translocation, H⁺ transfer, and large‐scale protein conformational changes. The complete picture of lactose/H⁺ symport is largely unclear due to the complexity and multiscale nature of the process. In this work, we develop the force field for sugar molecules compatible with PACE, a hybrid and coarse‐grained force field that couples the united‐atom protein models with the coarse‐grained MARTINI water/lipid. After validation, we implement the new force field to investigate the binding of a ‐d ‐galactopyranosyl‐1‐thio‐ ‐d ‐galactopyranoside (TDG) molecule to a wild‐type LacY. Results show that the local interactions between TDG and LacY at the binding pocket are consistent with the X‐ray experiment. Transitions from inward‐facing to outward‐facing conformations upon TDG binding and protonation of Glu269 have been achieved from ～5.5 µs simulations. Both the opening of the periplasmic side and closure of the cytoplasmic side of LacY are consistent with double electron–electron resonance and thiol cross‐linking experiments. Our analysis suggests that the conformational changes of LacY are a cumulative consequence of interdomain H‐bonds breaking at the periplasmic side, interdomain salt‐bridge formation at the cytoplasmic side, and the TDG orientational changes during the transition.     

Monte carlo study of the percolation in two-dimensional polymer systems

The structure of a two-dimensional film formed by adsorbed polymer chains was studied by means of Monte Carlo simulations. The polymer chains were represented by linear sequences of lattice beads and positions of these beads were restricted to vertices of a two-dimensional square lattice. Two different Monte Carlo methods were employed to determine the properties of the model system. The first was the random sequential adsorption (RSA) and the second one was based on Monte Carlo simulations with a Verdier-Stockmayer sampling algorithm. The methodology concerning the determination of the percolation thresholds for an infinite chain system was discussed. The influence of the chain length on both thresholds was presented and discussed. It was shown that the RSA method gave considerably lower thresholds for longer chains. This behavior can be explained by a different pool of chain conformations used in the calculations in both methods under consideration.

Effects of organic solvents and substrate binding on trypsin in acetonitrile and hexane media

In this work, we used molecular dynamic (MD) simulation to study trypsin with and without a six-amino-acid peptide bound in three different solvents (water, acetonitrile and hexane) in order to provide molecular information for well understanding the structure and function of enzymes in non-aqueous media. The results show that the enzyme is more compact and less native-like in hexane than in the other two polar solvents. The substrate could stabilize the native protein structure in the two polar media, but not in the non-polar hexane. There are no significant differences in the conformation of the S1 pocket upon the substrate binding in water and acetonitrile media while a reverse behavior is observed in hexane media, implying a possible induced fit binding mechanism in the non-polar media. The substrate binding enhances the stability of catalytic H-bond network since it could expel the solvent molecules from the active site. The enzyme and the substrate appear to be more appropriate to the reactive conformation in the organic solvents compared with aqueous solution. There is much greater substrate binding strength in hexane media than the water and acetonitrile ones since the polar solvent significantly weakens electrostatic interactions, which are observed to be the main driving force to the binding. In addition, some residues of the S1 pocket could remain favorable contribution to the binding despite the solvent change, but with differences in the contribution extent, the number and the type of residues between the three media.

Comparison of the structural characteristics of Cu2+-bound and unbound α-syn12 peptide obtained in simulations using different force fields

The effects of Cu²⁺ binding and the utilization of different force fields when modeling the structural characteristics of α-syn12 peptide were investigated. To this end, we performed extensive temperature replica exchange molecular dynamics (T-REMD) simulations on Cu²⁺-bound and unbound α-syn12 peptide using the GROMOS 43A1, OPLS-AA, and AMBER03 force fields. Each replica was run for 300 ns. The structural characteristics of α-syn12 peptide were studied based on backbone dihedral angle distributions, free-energy surfaces obtained with different reaction coordinates, favored conformations, the formation of different Turn structures, and the solvent exposure of the hydrophobic residues. The findings show that AMBER03 prefers to sample helical structures for the unbound α-syn12 peptide and does not sample any β-hairpin structure for the Cu²⁺-bound α-syn12 peptide. In contrast, the central structure of the major conformational clusters for the Cu²⁺-bound and unbound α-syn12 peptide according to simulations performed using the GROMOS 43A1 and OPLS-AA force fields is a β-hairpin with Turn_9-6. Cu²⁺ can also promote the formation of the β-hairpin and increase the solvent exposure of hydrophobic residues, which promotes the aggregation of α-syn12 peptide. This study can help us to understand the mechanisms through which Cu²⁺ participates in the fibrillation of α-syn12 peptide at the atomic level, which in turn represents a step towards elucidating the nosogenesis of Parkinson’s disease.

Structural adaptation of extreme halophilic proteins through decrease of conserved hydrophobic contact surface

Background

Disrupting protein-protein interactions by small organic molecules is nowadays a promising strategy employed to block protein targets involved in different pathologies. However, structural changes occurring at the binding interfaces make difficult drug discovery processes using structure-based drug design/virtual screening approaches. Here we focused on two homologous calcium binding proteins, calmodulin and human centrin 2, involved in different cellular functions via protein-protein interactions, and known to undergo important conformational changes upon ligand binding.

Results

In order to find suitable protein conformations of calmodulin and centrin for further structure-based drug design/virtual screening, we performed in silico structural/energetic analysis and molecular docking of terphenyl (a mimicking alpha-helical molecule known to inhibit protein-protein interactions of calmodulin) into X-ray and NMR ensembles of calmodulin and centrin. We employed several scoring methods in order to find the best protein conformations. Our results show that docking on NMR structures of calmodulin and centrin can be very helpful to take into account conformational changes occurring at protein-protein interfaces.

Conclusions

NMR structures of protein-protein complexes nowadays available could efficiently be exploited for further structure-based drug design/virtual screening processes employed to design small molecule inhibitors of protein-protein interactions.     

Molecular dynamics simulation of temperature induced unfolding of animal prion protein

To elucidate the structural stability and the unfolding dynamics of the animal prion protein, the temperature induced structural evolution of turtle prion protein (tPrPc) and bank vole prion protein (bvPrPc) have been performed with molecular dynamics (MD) simulation. The unfolding behaviors of secondary structures showed that the α-helix was more stable than β-sheet. Extension and disruption of β-sheet commonly appeared in the temperature induced unfolding process. The conversion of α-helix to π-helix occurred more readily at the elevating temperature. Furthermore, it was suggested in this work that the unfolding of prion protein could be regulated by the temperature.

Ligand-induced conformational changes in a thermophilic ribose-binding protein

Background

Members of the periplasmic binding protein (PBP) superfamily are involved in transport and signaling processes in both prokaryotes and eukaryotes. Biological responses are typically mediated by ligand-induced conformational changes in which the binding event is coupled to a hinge-bending motion that brings together two domains in a closed form. In all PBP-mediated biological processes, downstream partners recognize the closed form of the protein. This motion has also been exploited in protein engineering experiments to construct biosensors that transduce ligand binding to a variety of physical signals. Understanding the mechanistic details of PBP conformational changes, both global (hinge bending, twisting, shear movements) and local (rotamer changes, backbone motion), therefore is not only important for understanding their biological function but also for protein engineering experiments.

Results

Here we present biochemical characterization and crystal structure determination of the periplasmic ribose-binding protein (RBP) from the hyperthermophile Thermotoga maritima in its ribose-bound and unliganded state. The T. maritima RBP (tmRBP) has 39% sequence identity and is considerably more resistant to thermal denaturation ( ^app T _m value is 108°C) than the mesophilic Escherichia coli homolog (ecRBP) ( ^app T _m value is 56°C). Polar ligand interactions and ligand-induced global conformational changes are conserved among ecRBP and tmRBP; however local structural rearrangements involving side-chain motions in the ligand-binding site are not conserved.

Conclusion

Although the large-scale ligand-induced changes are mediated through similar regions, and are produced by similar backbone movements in tmRBP and ecRBP, the small-scale ligand-induced structural rearrangements differentiate the mesophile and thermophile. This suggests there are mechanistic differences in the manner by which these two proteins bind their ligands and are an example of how two structurally similar proteins utilize different mechanisms to form a ligand-bound state.     

Effects of protein binding on a lipid bilayer containing local anesthetic articaine, and the potential of mean force calculation: a molecular dynamics simulation approach

Articaine, as a local anesthetic drug has been simulated in neutral and charged forms, and its interaction with the dimyristoylphosphatidylcholine (DMPC) lipid bilayer membrane is investigated by molecular dynamics simulation using GROMACS software. In order to obtain the optimum location of the drug molecules, as they penetrate into the membrane, umbrella sampling is applied and the free energy is calculated. The effect of protein binding to DMPC membrane on the process of drug diffusion through the membrane is considered. Five simulation systems are designed and by applying the potential of mean force, the molecular dynamics simulation on the system is performed. In light of the obtained results, the electrostatic potential, variation of lipid bilayer’s order parameter and the diffusion coefficient of drug are discussed.

Interaction Between Umami Peptide and Taste Receptor T1R1/T1R3

The umami taste receptor is a heterodimer composed of two members of the T1R taste receptor family: T1R1 and T1R3. The homology models of the ligand binding domains of the human umami receptor have been constructed based on crystallographic structures of the taste receptor of the central nervous system. Furthermore, the molecular simulations of the ligand binding domain show that the likely conformation was that T1R1 protein exists in the closed conformation, and T1R3 in the open conformation in the heterodimer. The molecular docking study of T1R1 and T1R3 in complex with four peptides, including Lys–Gly–Asp–Glu–Ser–Leu–Leu–Ala, Ser–Glu–Glu, G1u–Ser, and Asp–Glu–Ser, displayed that the amino acid residue of SER146 and Glu277 in T1R3 may play great roles in the synergism of umami taste. This docking result further validated the robustness of the model. In the paper, binding of umami peptide and the T1R1/T1R3 receptor was first described and the interaction is the base of umami activity theory.     

Why tazobactam and sulbactam have different intermediates population with SHV-1 β-lactamase: a molecular dynamics study

The imine intermediates of tazobactam and sulbactam bound to SHV-1 β-lactamase were investigated by molecular dynamics (MD) simulation respectively. Hydrogen bond networks around active site were found different between tazobactam and sulbactam acyl-enzymes. In tazobactam imine intermediate, it was observed that the triazolyl ring formed stable hydrogen bonds with Asn170 and Thr167. The results suggest that conformation of imine determined the population of intermediates. In imine intermediate of tazobactam, the triazolyl ring is trapped in Thr_Asn pocket, and it restricts the rotation of C5-C6 bond so that tazobactam can only generate trans enamine intermediate. Further, conformational cluster analyses are performed to substantiate the results. These findings provide an explanation for the corresponding experimental results, and will be potentially useful in the development of new inhibitors.