Docking and molecular dynamics simulations of the Fyn‐SH3 domain with free and phospholipid bilayer‐associated 18.5‐kDa myelin basic protein (MBP)—Insights into a noncanonical and fuzzy interaction

The molecular details of the association between the human Fyn‐SH3 domain, and the fragment of 18.5‐kDa myelin basic protein (MBP) spanning residues S38–S107 (denoted as xα2‐peptide, murine sequence numbering), were studied in silico via docking and molecular dynamics over 50‐ns trajectories. The results show that interaction between the two proteins is energetically favorable and heavily dependent on the MBP proline‐rich region (P93‐P98) in both aqueous and membrane environments. In aqueous conditions, the xα2‐peptide/Fyn‐SH3 complex adopts a “sandwich”"‐like structure. In the membrane context, the xα2‐peptide interacts with the Fyn‐SH3 domain via the proline‐rich region and the β‐sheets of Fyn‐SH3, with the latter wrapping around the proline‐rich region in a form of a clip. Moreover, the simulations corroborate prior experimental evidence of the importance of upstream segments beyond the canonical SH3‐ligand. This study thus provides a more‐detailed glimpse into the context‐dependent interaction dynamics and importance of the β‐sheets in Fyn‐SH3 and proline‐rich region of MBP. Proteins 2017; 85:1336–1350. © 2017 Wiley Periodicals, Inc.     

Reproducible polypeptide folding and structure prediction using molecular dynamics simulations

The folding of a polypeptide from an extended state to a well-defined conformation is studied using microsecond classical molecular dynamics (MD) simulations and replica exchange molecular dynamics (REMD) simulations in explicit solvent and in vacuo. It is shown that the solvated peptide folds many times in the REMD simulations but only a few times in the conventional simulations. From the folding events in the classical simulations we estimate an approximate folding time of 1-2 micros. The REMD simulations allow enough sampling to deduce a detailed Gibbs free energy landscape in three dimensions. The global minimum of the energy landscape corresponds to the native state of the peptide as determined previously by nuclear magnetic resonance (NMR) experiments. Starting from an extended state it takes about 50 ns before the native structure appears in the REMD simulations, about an order of magnitude faster than conventional MD. The calculated melting curve is in good qualitative agreement with experiment. In vacuo, the peptide collapses rapidly to a conformation that is substantially different from the native state in solvent.     

Molecular dynamics simulations of bovine rhodopsin: influence of protonation states and different membrane-mimicking environments

G-protein coupled receptors (GPCRs) are a protein family of outstanding pharmaceutical interest. GPCR homology models, based on the crystal structure of bovine rhodopsin, have been shown to be valuable tools in the drug-design process. The initial model is often refined by molecular dynamics (MD) simulations, a procedure that has been recently discussed controversially. We therefore analyzed MD simulations of bovine rhodopsin in order to identify contacts that could serve as constraints in the simulation of homology models. Additionally, the effect of an N-terminal truncation, the nature of the membrane mimic, the influence of varying protonation states of buried residues and the importance of internal water molecules was analyzed. All simulations were carried out using the program-package GROMACS. While N-terminal truncation negatively influenced the overall protein stability, a stable simulation was possible in both solvent environments. As regards the protonation state of titratable sites, the experimental data could be reproduced by the program UHBD (University of Houston Brownian Dynamics), suggesting its application for studying homology models of GPCRs. A high flexibility was observed for internal water molecules at some sites. Finally, interhelical hydrogen-bonding interactions could be derived, which can now serve as constraints in the simulations of GPCR homology models.     

Computational exploration of polymorphisms in 5-hydoxytryptamine 5-HT₁A and 5-HT₂A receptors associated with psychiatric disease

The huge polymorphic data have been prioritized towards a specific disease based on sequence and structure homology tools to a large extent. In this study, we have explored the potential non-synonymous Single Nucleotide Polymorphism (nsSNP) in serotonin (5-HT) receptors involved in psychotic syndromes and their response pathway. The most damaging point mutations were screened from 12 classes of serotonin receptors comprising 7743 variants. In 5HT(1A) receptor, two alleles were found to be highly deleterious located at ligand binding extracellular-2 and one at intracellular loop-3 domains. Similarly, we found two alleles predicted to be highly damaging in 5HT(2A) residing at N and C-Terminal domains. The above alleles were further confirmed based on their flexibility and stability difference using the molecular dynamic simulation analysis. Integrating these results appeared promising for being able to filter out potential non-synonymous Single Nucleotide Polymorphisms for neuropsychiatric disorders.     

The recognition of membrane‐bound PtdIns3P by PX domains

Phox‐homology (PX) domains target proteins to the organelles of the secretary and endocytic systems by binding to phosphatidylinositol phospholipids (PIPs). Among all the structures of PX domains that have been solved, only three have been solved in a complex with the main physiological ligand: PtdIns3P. In this work, molecular dynamic simulations have been used to explore the structure and dynamics of the p40^phox–PX domain and the SNX17–PX domain and their interaction with membrane‐bound PtdIns3P. In the simulations, both PX domains associated spontaneously with the membrane‐bound PtdIns3P and formed stable complexes. The interaction between the p40^phox–PX domain and PtdIns3P in the membrane was found to be similar to the crystal structure of the p40^phox–PX–PtdIns3P complex that is available. The interaction between the SNX17–PX domain and PtdIns3P was similar to that observed in the p40^phox–PX–PtdIns3P complex; however, some residues adopted different orientations. The simulations also showed that nonspecific interactions between the β1–β2 loop and the membrane play an important role in the interaction of membrane bound PtdIns3P and different PX domains. The behaviour of unbound PtdIns3P within a 2‐oleoyl‐1‐palmitoyl‐sn‐glycero‐3‐phosphocholine (POPC) membrane environment was also examined and compared to the available experimental data and simulation studies of related molecules. Proteins 2014; 82:2332–2342. © 2014 Wiley Periodicals, Inc.     

Homology modelling and molecular dynamics simulations of a protein serine/threonine phosphatase stp1 in Staphylococcus aureus N315: a potential drug target

The prevalence of methicillin-resistant Staphylococcus aureus and vanomycin intermediate S. aureus infections is on the rise, globally. This poses a huge challenge due to limited therapeutic options and the limited number of bacterial-specific drug targets available for due conservation with the human host. A serine/threonine phosphatase/kinase stp1/stk1 phospho-signalling system in S. aureus, which is just beginning to be understood, has been shown to be of importance in virulence and susceptibility to glycopeptide antibiotics. In this study, 3D structure of stp1 (clinical strain of S. aureus N315) was predicted using a homology modelling tool MODELLER. The validation of the predicted model was done using various tools such as PROCHECK, ERRAT, VERIFY-3D and ProSA. Molecular dynamics (MD) study was carried out using GROMACS to refine the least energy model generated from MODELLER9v11 and it was compared with the template. The template used was the crystal structure of serine/threonine phosphatase stp1 in Streptoccocus agalactiae (Protein Data Bank ID: 2PK0) with 38% identity with the query. Various validation tools showed the quality of the model generated using MODELLER. PROCHECK predicted 100% residues in the allowed region, ERRAT with overall quality factor of 76.47, VERIFY-3D with average score of >0.2 in 81.78% of residues, WHATIF with packaging quality score of > ? 5 for all residues and ProSA with Z-score of ? 7.02. MD simulation of the protein showed some fluctuations in the aqueous environment and changes in the ligand binding residues after simulation. The availability of the 3D-structural information of a viable drug target in S. aureus stp1 is expected to facilitate structure–activity relationship and interactions with proteins.     

Structure of RecX protein complex with the presynaptic RecA filament: Molecular dynamics simulations and small angle neutron scattering

Using molecular modeling techniques we have built the full atomic structure and performed molecular dynamics simulations for the complexes formed by Escherichia coli RecX protein with a single-stranded oligonucleotide and with RecA presynaptic filament. Based on the modeling and SANS experimental data a sandwich-like filament structure formed two chains of RecX monomers bound to the opposite sides of the single stranded DNA is proposed for RecX::ssDNA complex. The model for RecX::RecA::ssDNA include RecX binding into the grove of RecA::ssDNA filament that occurs mainly via Coulomb interactions between RecX and ssDNA. Formation of RecX::RecA::ssDNA filaments in solution was confirmed by SANS measurements which were in agreement with the spectra computed from the molecular dynamics simulations.