Structural characterization of HDAC2-MTA1 complex

Major histocompatibility complex class I (MHC-I), a key element of the acquired immune system, plays essential roles in activating CD8+ T cells by recognizing intracellular antigens derived from pa...

Protein conformational transition from alpha-helices to beta-sheets precedes aggregation of proteins implicated in many diseases, including Alzheimer and prion diseases. Direct characterization of such transitions is often hindered by the complicated nature of the interaction network among amino acids. A recently engineered small protein-like peptide with a simple amino acid composition features a temperature-driven alpha-helix to beta-sheet conformational change. Here we studied the conformational transition of this peptide by molecular dynamics simulations. We observed a critical temperature, below which the peptide folds into an alpha-helical coiled-coil state and above which the peptide misfolds into beta-rich structures with a high propensity to aggregate. The structures adopted by this peptide during low temperature simulations have a backbone root mean square deviation less than 2 A from the crystal structure. At high temperatures, this peptide adopts an amyloid-like structure, which is mainly composed of coiled anti-parallel beta-sheets with the cross-beta-signature of amyloid fibrils. Most strikingly, we observed conformational conversions in which an alpha-helix is converted into a beta-strand by proximate stable beta-sheets with exposed hydrophobic surfaces and unsaturated hydrogen bonds. Our study suggested a possible generic molecular mechanism of the template-mediated aggregation process, originally proposed by Prusiner (Prusiner, S. B. (1998) Proc. Natl. Acad. Sci. U. S. A. 95, 13363-13383) to account for prion infectivity.

ATP binding cassette transporters are integral membrane proteins that use the energy released from ATP hydrolysis at the two nucleotide binding domains (NBDs) to translocate a wide variety of substrates through a channel at the two transmembrane domains (TMDs) across the cell membranes. MsbA from Gram-negative bacteria is a lipid and multidrug resistance ATP binding cassette exporter that can undergo large scale conformational changes between the outward-facing and the inward-facing conformations revealed by crystal structures in different states. Here, we use targeted molecular dynamics simulation methods to explore the atomic details of the conformational transition from the outward-facing to the inward-facing states of MsbA. The molecular dynamics trajectories revealed a clear spatiotemporal order of the conformational movements. The disruption of the nucleotide binding sites at the NBD dimer interface is the very first event that initiates the following conformational changes, verifying the assumption that the conformational conversion is triggered by ATP hydrolysis. The conserved x-loops of the NBDs were identified to participate in the interaction network that stabilizes the cytoplasmic tetrahelix bundle of the TMDs and play an important role in mediating the cross-talk between the NBD and TMD. The movement of the NBD dimer is transmitted through x-loops to break the tetrahelix bundle, inducing the packing rearrangements of the transmembrane helices at the cytoplasmic side and the periplasmic side sequentially. The packing rearrangement within each periplasmic wing of TMD that results in exposure of the substrate binding sites occurred at the end stage of the trajectory, preventing the wrong timing of the binding site accessibility.

Protein-ligand interactions are a necessary prerequisite for signal transduction, immunoreaction, and gene regulation. Protein-ligand interaction studies are important for understanding the mechanisms of biological regulation, and they provide a theoretical basis for the design and discovery of new drug targets. In this study, we analyzed the molecular interactions of protein-ligand which was docked by AutoDock 4.2 software. In AutoDock 4.2 software, we used a new search algorithm, hybrid algorithm of random drift particle swarm optimization and local search (LRDPSO), and the classical Lamarckian genetic algorithm (LGA) as energy optimization algorithms. The best conformations of each docking algorithm were subjected to molecular dynamic (MD) simulations to further analyze the molecular mechanisms of protein-ligand interactions. Here, we analyze the binding energy between protein receptors and ligands, the interactions of salt bridges and hydrogen bonds in the docking region, and the structural changes during complex unfolding. Our comparison of these complexes highlights differences in the protein-ligand interactions between the two docking methods. It also shows that salt bridge and hydrogen bond interactions play a crucial role in protein-ligand stability. The present work focuses on extracting the deterministic characteristics of docking interactions from their dynamic properties, which is important for understanding biological functions and determining which amino acid residues are crucial to docking interactions.

Interplay between salt bridge, hydrogen bond and anion-π interactions in thiocyanate binding

Interactions of transcriptional corepressors such as histone deacetylase 7 (HDAC7), a class IIa HDAC, with myocyte enhancer factor-2 (MEF2) regulate MEF2 activity. Despite previous investigations exploring interactions between HDAC7 and MEF2, a detailed characterization of the HDAC7-MEF2 functional complex is still lacking. Herein, we first modeled the structure of the HDAC7-MEF2A complex and investigated the inter-protein interactions using all-atom molecular dynamics (MD) simulations. We identified specific amino acids within HDAC7 and MEF2A that participate in interactions such as salt bridges, hydrogen bonds, and hydrophobic interactions. Our results reveal a salt bridge formed between LYS96(HDAC7) and ASP63(MEF2A). Our analysis also predicted formations of reliable hydrogen bonds between SER82(HDAC7) and ASP63(MEF2A) as well as LYS96(HDAC7) and ASP63(MEF2A). In addition, clustering of hydrophobic residues at the interface contributes in stabilizing the HDAC7-MEF2A complex. Results from multiple sequence alignment show that most of the HDAC7 residues that are predicted to associate with MEF2A are conserved in at least three class IIa HDACs and all predicted residues in MEF2A are conserved in MEF2s. We also found that the association of DNA to MEF2A has no significant effect on HDAC7-MEF2A interactions. Our results may also provide useful insights into the interactions between other class IIa HDACs and MEF2s.

Oculocutaneous albinism type 3 (OCA3) is an autosomal recessive disorder caused by mutations in the TYRP1 gene. Tyrosinase-related protein 1 (Tyrp1) is involved in eumelanin synthesis, catalyzing the oxidation of 5,6-dihydroxyindole-2-carboxylic acid oxidase (DHICA) to 5,6-indolequinone-2-carboxylic acid (IQCA). Here, for the first time, four OCA3-causing mutations of Tyrp1, C30R, H215Y, D308N, and R326H, were investigated computationally to understand Tyrp1 protein stability and catalytic activity. Using the Tyrp1 crystal structure (PDB:5M8L), global mutagenesis was conducted to evaluate mutant protein stability. Consistent with the foldability parameter, C30R and H215Y should exhibit greater instability, and two other mutants, D308N and R326H, are expected to keep a native conformation. SDS-PAGE and Western blot analysis of the purified recombinant proteins confirmed that the foldability parameter correctly predicted the effect of mutations critical for protein stability. Further, the mutant variant structures were built and simulated for 100 ns to generate free energy landscapes and perform docking experiments. Free energy landscapes formed by Y362, N378, and T391 indicate that the binding clefts of C30R and H215Y mutants are larger than the wild-type Tyrp1. In docking simulations, the hydrogen bond and salt bridge interactions that stabilize DHICA in the active site remain similar among Tyrp1, D308N, and R326H. However, the strengths of these interactions and stability of the docked ligand may decrease proportionally to mutation severity due to the larger and less well-defined natures of the binding clefts in mutants. Mutational perturbations in mutants that are not unfolded may result in allosteric alterations to the active site, reducing the stability of protein-ligand interactions.

Discovery of benzo[d]isothiazole derivatives as novel scaffold inhibitors targeting the programmed cell death-1/programmed cell death-ligand 1 (PD-1/PD-L1) interaction through “ring fusion” strategy

Studying the recognition mechanism of TcaR and ssDNA using molecular dynamic simulations

ObjectivesThree strong interactions between amino acid side chains (salt bridge, cation-π, and amide bridge) are studied that are stronger than (or comparable to) the common hydrogen bond interactions, and play important roles in protein-protein interactions.MethodsQuantum chemical methods MP2 and CCSD(T) are used in calculations of interaction energies and structural optimizations.ResultsThe energies of three types of amino acid side chain interactions in gaseous phase and in aqueous solutions are calculated using high level quantum chemical methods and basis sets. Typical examples of amino acid salt bridge, cation-π, and amide bridge interactions are analyzed, including the inhibitor design targeting neuraminidase (NA) enzyme of influenza A virus, and the ligand binding interactions in the HCV p7 ion channel. The inhibition mechanism of the M2 proton channel in the influenza A virus is analyzed based on strong amino acid interactions.Conclusion(1) The salt bridge interactions between acidic amino acids (Glu- and Asp-) and alkaline amino acids (Arg+, Lys+ and His+) are the strongest residue-residue interactions. However, this type of interaction may be weakened by solvation effects and broken by lower pH conditions. (2) The cation- interactions between protonated amino acids (Arg+, Lys+ and His+) and aromatic amino acids (Phe, Tyr, Trp and His) are 2.5 to 5-fold stronger than common hydrogen bond interactions and are less affected by the solvation environment. (3) The amide bridge interactions between the two amide-containing amino acids (Asn and Gln) are three times stronger than hydrogen bond interactions, which are less influenced by the pH of the solution. (4) Ten of the twenty natural amino acids are involved in salt bridge, or cation-, or amide bridge interactions that often play important roles in protein-protein, protein-peptide, protein-ligand, and protein-DNA interactions.

Beta-lactamase inhibitor protein (BLIP) binds a variety of class A beta-lactamases with affinities ranging from micromolar to picomolar. Whereas the TEM-1 and SHV-1 beta-lactamases are almost structurally identical, BLIP binds TEM-1 approximately 1000-fold tighter than SHV-1. Determining the underlying source of this affinity difference is important for understanding the molecular basis of beta-lactamase inhibition and mechanisms of protein-protein interface specificity and affinity. Here we present the 1.6A resolution crystal structure of SHV-1.BLIP. In addition, a point mutation was identified, SHV D104E, that increases SHV.BLIP binding affinity from micromolar to nanomolar. Comparison of the SHV-1.BLIP structure with the published TEM-1.BLIP structure suggests that the increased volume of Glu-104 stabilizes a key binding loop in the interface. Solution of the 1.8A SHV D104K.BLIP crystal structure identifies a novel conformation in which this binding loop is removed from the interface. Using these structural data, we evaluated the ability of EGAD, a program developed for computational protein design, to calculate changes in the stability of mutant beta-lactamase.BLIP complexes. Changes in binding affinity were calculated within an error of 1.6 kcal/mol of the experimental values for 112 mutations at the TEM-1.BLIP interface and within an error of 2.2 kcal/mol for 24 mutations at the SHV-1.BLIP interface. The reasonable success of EGAD in predicting changes in interface stability is a promising step toward understanding the stability of the beta-lactamase.BLIP complexes and computationally assisted design of tight binding BLIP variants.

Allosteric Effects of Sodium Ion Binding on Activation of the M3 Muscarinic G-Protein-Coupled Receptor

All-atom molecular dynamics simulations are used with the highly mobile membrane mimetic method to study the α6-α7 peptide of the critical yeast Osh4 peripheral membrane protein. This research focuses on the impact of 1-palmitoyl-2-oleoyl-sn-glycero-phosphatidylinoside 4,5-bisphosphate (PIP2) lipids and 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-l-serine on the protein's ability to bind to the membrane. Details of the binding mechanism are described qualitatively and quantitatively by measuring the position of the deepest residues, angle of the peptide during binding, root mean square deviation of the atomic positions within the peptide, and interaction energy, while changing variables, such as the force field used and the presence of the PIP2 lipids. The negatively charged PIP2 has a large head group that is a few Ångstroms above the main membrane phosphates enabling the PIP2 lipids to interact with the peptide before it binds deeper into the membrane. The PIP2 lipids can alter the position of the peptide during binding by recruiting charged residues on the α7 helix, such as R344 and R347. Residues R347 and R344 are unusual because they are slightly out of the reach of the main membrane phosphates but optimally positioned to interact with the PIP2 lipids. The salt-bridge interactions can also typically occur between cationic peptide residues such as R314, K325, and K336. The force field interaction effect on peptide binding was also investigated by changing the standard CHARMM36m to an improved description between some amino acids and lipid moieties (Phys. Chem. Chem. Phys. 20, 8432-8449). This resulted in the total number of salt bridges and hydrogen bonds being drastically reduced, the interaction energy was also reduced, and there was more balance between electrostatic and nonpolar interactions, but the general bound structure is maintained. This work is an important initial step to understand the effect of the Osh4 protein on the membrane binding and to quantify the effect of PIP2 lipids on this domain.

Structural conservation in the glutathione binding in Sphingomonas sp. glutaredoxin Grx3 and variations for cold adaptation

Six helix surface positions of protein G (Gbeta1) were redesigned using a computational protein design algorithm, resulting in the five fold mutant Gbeta1m2. Gbeta1m2 is well folded with a circular dichroism spectrum nearly identical to that of Gbeta1, and a melting temperature of 91 degrees C, approximately 6 degrees C higher than that of Gbeta1. The crystal structure of Gbeta1m2 was solved to 2.0 A resolution by molecular replacement. The absence of hydrogen bond or salt bridge interactions between the designed residues in Gbeta1m2 suggests that the increased stability of Gbeta1m2 is due to increased helix propensity and more favorable helix dipole interactions.