Binding region of alanopine dehydrogenase predicted by unbiased molecular dynamics simulations of ligand diffusion

Aplospojaveedins A-C, unusual sulfur-containing alkaloids produced by the endophytic fungus Aplosporella javeedii using OSMAC strategy

Three sulfur-containing alkaloids aplospojaveedins A-C (1-3) with a hitherto undescribed carbon skeleton comprising octahy-dronaphthalene, α, β-unsaturated lactam and glycine-cysteine moieties were isolated from Aplosporella javeedii. Their structures were elucidated by 1D and 2D NMR spectroscopy, HR-MS, X-ray diffraction analysis, DFT-NMR and TDDFT-ECD calculations. A plausible biosynthetic pathway and putative targets are described. The blind docking suggested that 1-3 may have functional effects on several putative targets such as the GPCR cannabinoid receptor 2 or the integrin α5β1 complex.

Extracting binding energies and binding modes from biomolecular simulations of fragment binding to endothiapepsin

Fragment-based drug discovery (FBDD) aims to discover a set of small binding fragments that may be subsequently linked together. Therefore, in-depth knowledge of the individual fragments' structural and energetic binding properties is essential. In addition to experimental techniques, the direct simulation of fragment binding by molecular dynamics (MD) simulations became popular to characterize fragment binding. However, former studies showed that long simulation times and high computational demands per fragment are needed, which limits applicability in FBDD. Here, we performed short, unbiased MD simulations of direct fragment binding to endothiapepsin, a well-characterized model system of pepsin-like aspartic proteases. To evaluate the strengths and limitations of short MD simulations for the structural and energetic characterization of fragment binding, we predicted the fragments' absolute free energies and binding poses based on the direct simulations of fragment binding and compared the predictions to experimental data. The predicted absolute free energies are in fair agreement with the experiment. Combining the MD data with binding mode predictions from molecular docking approaches helped to correctly identify the most promising fragments for further chemical optimization. Importantly, all computations and predictions were done within 5 days, suggesting that MD simulations may become a viable tool in FBDD projects.

Heterologous Expression and Partial Characterization of a Putative Opine Dehydrogenase from a Metagenomic Sequence of Desulfohalobium retbaense

The aim of this research was to prove the function of the putative opine dehydrogenase from Desulfohalobium retbaense and to characterize the enzyme in terms of functional and kinetic parameters. A putative opine dehydrogenase was identified from a metagenomic library by a sequence-based technique search of the metagenomic library, and afterward was successfully heterologously produced in Escherichia coli. In order to examine its potential for applications in the synthesis of secondary amines, first the substrate specificity of the enzyme towards different amino donors and amino acceptors was determined. The highest affinity was observed towards small amino acids, preferentially L-alanine, and when it comes to α-keto acids, pyruvate proved to be a preferential amino acceptor. The highest activity was observed at pH 6.5 in the absence of salts. The enzyme showed remarkable stability in a wide range of experimental conditions, such as broad pH stability (from 6.0-11.0 after 30 min incubation in buffers at a certain pH), stability in the presence of NaCl up to 3.0 M for 24 h, it retained 80 % of the initial activity after 1 h incubation at 45 °C, and 65 % of the initial activity after 24 h incubation in 30 % dimethyl sulfoxide.

Opine dehydrogenases, an underexplored enzyme family for the enzymatic synthesis of chiral amines

Opines and opine-type chemicals are valuable natural products with diverse biochemical roles, and potential synthetic building blocks of bioactive compounds. Their synthesis involves reductive amination of ketoacids with amino acids. This transformation has high synthetic potential in producing enantiopure secondary amines. Nature has evolved opine dehydrogenases for this chemistry. To date, only one enzyme has been used as biocatalyst, however, analysis of the available sequence space suggests more enzymes to be exploited in synthetic organic chemistry. This review summarizes the current knowledge of this underexplored enzyme class, highlights key molecular, structural, and catalytic features with the aim to provide a comprehensive general description of opine dehydrogenases, thereby supporting future enzyme discovery and protein engineering studies.

Complete reconstruction of dasatinib unbinding pathway from c-Src kinase by supervised molecular dynamics simulation method; assessing efficiency and trustworthiness of the method

Over the past years, rational drug design has gained lots of attention since employing it gave the world targeted therapy and more effective treatment solutions. Structure-based drug design (SBDD) is an excellent tool in rational drug design that takes advantage of accurate methods such as unbiased molecular dynamics (UMD) simulation for designing and optimizing molecular entities by understanding the binding and unbinding pathways of the binders. Supervised molecular dynamics (SuMD) simulation is a branch of UMD in which long-duration simulations are turned into short simulations, called replica, and a specific parameter is monitored throughout the simulation. In this work, we utilized this strategy to reconstruct the unbinding pathway of the anticancer drug dasatinib from its target protein, the c-Src kinase. Several unbinding events with valuable details were achieved. Then, to assess the efficiency and trustworthiness of the SuMD method, the unbinding pathway was also reconstructed by conventional UMD simulation, which uncovered some of the limitations of this method, such as limited sampling of the active site and finding the metastable states in the unbinding pathway. Furthermore, in times like these, when the world is desperate to find treatments for the Covid-19 disease, we think these methods are of exceptional value.Communicated by Ramaswamy H. Sarma.

Biophysical and pharmacokinetic characterization of a small-molecule inhibitor of RUNX1/ETO tetramerization with anti-leukemic effects

Acute myeloid leukemia (AML) is a malignant disease of immature myeloid cells and the most prevalent acute leukemia among adults. The oncogenic homo-tetrameric fusion protein RUNX1/ETO results from the chromosomal translocation t(8;21) and is found in AML patients. The nervy homology region 2 (NHR2) domain of ETO mediates tetramerization; this oligomerization is essential for oncogenic activity. Previously, we identified the first-in-class small-molecule inhibitor of NHR2 tetramer formation, 7.44, which was shown to specifically interfere with NHR2, restore gene expression down-regulated by RUNX1/ETO, inhibit the proliferation of RUNX1/ETO-depending SKNO-1 cells, and reduce the RUNX1/ETO-related tumor growth in a mouse model. However, no biophysical and structural characterization of 7.44 binding to the NHR2 domain has been reported. Likewise, the compound has not been characterized as to physicochemical, pharmacokinetic, and toxicological properties. Here, we characterize the interaction between the NHR2 domain of RUNX1/ETO and 7.44 by biophysical assays and show that 7.44 interferes with NHR2 tetramer stability and leads to an increase in the dimer population of NHR2. The affinity of 7.44 with respect to binding to NHR2 is K_lig = 3.75 ± 1.22 µM. By NMR spectroscopy combined with molecular dynamics simulations, we show that 7.44 binds with both heteroaromatic moieties to NHR2 and interacts with or leads to conformational changes in the N-termini of the NHR2 tetramer. Finally, we demonstrate that 7.44 has favorable physicochemical, pharmacokinetic, and toxicological properties. Together with biochemical, cellular, and in vivo assessments, the results reveal 7.44 as a lead for further optimization towards targeted therapy of t(8;21) AML.

Development of a First-in-Class Small-Molecule Inhibitor of the C-Terminal Hsp90 Dimerization

Heat shock proteins 90 (Hsp90) are promising therapeutic targets due to their involvement in stabilizing several aberrantly expressed oncoproteins. In cancerous cells, Hsp90 expression is elevated, thereby exerting antiapoptotic effects, which is essential for the malignant transformation and tumor progression. Most of the Hsp90 inhibitors (Hsp90i) under investigation target the ATP binding site in the N-terminal domain of Hsp90. However, adverse effects, including induction of the prosurvival resistance mechanism (heat shock response or HSR) and associated dose-limiting toxicity, have so far precluded their clinical approval. In contrast, modulators that interfere with the C-terminal domain (CTD) of Hsp90 do not inflict HSR. Since the CTD dimerization of Hsp90 is essential for its chaperone activity, interfering with the dimerization process by small-molecule protein-protein interaction inhibitors is a promising strategy for anticancer drug research. We have developed a first-in-class small-molecule inhibitor (5b) targeting the Hsp90 CTD dimerization interface, based on a tripyrimidonamide scaffold through structure-based molecular design, chemical synthesis, binding mode model prediction, assessment of the biochemical affinity, and efficacy against therapy-resistant leukemia cells. 5b reduces xenotransplantation of leukemia cells in zebrafish models and induces apoptosis in BCR-ABL1⁺ (T315I) tyrosine kinase inhibitor-resistant leukemia cells, without inducing HSR.

Reconstruction of the binding pathway of an anti-HIV drug, Indinavir, in complex with the HTLV-1 protease using unaggregated unbiased molecular dynamics simulation

Retroviruses are a growing concern for the health of human beings, and one of the dangerous members of this family is the Human T-cell Leukemia Virus 1 (HTLV-1) virus. It has affected more than 20 million people so far, and since there are no registered treatments against it yet, urgent treatment solutions are needed. One of the most promising drug targets to fight this virus is the protease enzyme of the virus's protein machinery. In this study, by utilizing a computational method called Unaggregated Unbiased Molecular Dynamics (UUMD), we reconstructed the binding pathway of a HTLV-1 protease inhibitor, Indinavir, to find the details of the binding pathway, the influential residues, and also the stable states of the binding pathway. We achieved the native conformation of the inhibitor in 6 rounds, 360 replicas by performing over 4 micro-seconds of UMD simulations. We found 3 Intermediate states between the solvated state and the native conformation state in the binding pathway. We also discovered that aromatic residues such as Trp98 and Trp98', catalytic residues Asp32 and Asp32', and the flap region's residues have the most influential roles in the binding pathway and also have the most contribution to the total interaction energies. We believe that the details found in this study would be a great guide for developing new treatment solutions against the HTLV-1 virus by inhibiting the HTLV-1 protease.

Comparative analysis of the unbinding pathways of antiviral drug Indinavir from HIV and HTLV1 proteases by supervised molecular dynamics simulation

Determining the unbinding pathways of potential small molecule compounds from their target proteins is of great significance for designing efficacious treatment solutions. One of these potential compounds is the approved HIV-1 protease inhibitor, Indinavir, which has a weak effect on the HTLV-1 protease. In this work, by employing the SuMD method, we reconstructed the unbinding pathways of Indinavir from HIV and HTLV-1 proteases to compare and understand the mechanism of the unbinding and to discover the reasons for the lack of inhibitory activity of Indinavir against the HTLV-1 protease. We achieved multiple unbinding events from both HIV and HTLV-1 proteases in which the RMSD values of Indinavir reached over 40 Å. Also, we found that the mobility and fluctuations of the flap region are higher in the HTLV-1 protease, making the drug less stable. We realized that critically positioned aromatic residues such as Trp98/Trp98' and Phe67/Phe67' in the HTLV-1 protease could make strong π-Stacking interactions with Indinavir in the unbinding pathway, which are unfavorable for the stability of Indinavir in the active site. The details found in this study can make a reasonable explanation for the lack of inhibitory activity of this drug against HTLV-1 protease. We believe the details discovered in this work can help design more effective and selective inhibitors for the HTLV-1 protease.

Unraveling the unbinding pathways of SARS-CoV-2 Papain-like proteinase known inhibitors by Supervised Molecular Dynamics simulation

The COVID-19 disease has infected and killed countless people all over the world since its emergence at the end of 2019. No specific therapy for COVID-19 is not currently available, and urgent treatment solutions are needed. Recent studies have found several potential molecular targets, and one of the most critical proteins of the SARS-CoV-2 virus work machine is the Papain-like protease (Plpro). Potential inhibitors are available, and their X-ray crystallographic structures in complex with this enzyme have been determined recently. However, their activities against this enzyme are insufficient and need to be characterized and improved to be of clinical values. Therefore, in this work, by utilizing the Supervised Molecular Dynamics (SuMD) simulation method, we achieved multiple unbinding events of Plpro inhibitors, GRL0617, and its derivates, and captured and understood the details of the unbinding pathway. We found that residues of the BL2 loop, such as Tyr268 and Gln269, play major roles in the unbinding pathways, but the most important contributing factor is the natural movements and behavior of the BL2 loop, which can control the entire process. We believe that the details found in this study can be used to refine and optimize potential inhibitors like GRL0617 and design more efficacious inhibitors as a treatment for the SARS-CoV-2 virus.

Promiscuous Esterases Counterintuitively Are Less Flexible than Specific Ones

Understanding mechanisms of promiscuity is increasingly important from a fundamental and application point of view. As to enzyme structural dynamics, more promiscuous enzymes generally have been recognized to also be more flexible. However, examples for the opposite received much less attention. Here, we exploit comprehensive experimental information on the substrate promiscuity of 147 esterases tested against 96 esters together with computationally efficient rigidity analyses to understand the molecular origin of the observed promiscuity range. Unexpectedly, our data reveal that promiscuous esterases are significantly less flexible than specific ones, are significantly more thermostable, and have a significantly increased specific activity. These results may be reconciled with a model according to which structural flexibility in the case of specific esterases serves for conformational proofreading. Our results signify that an esterase sequence space can be screened by rigidity analyses for promiscuous esterases as starting points for further exploration in biotechnology and synthetic chemistry.

A boosted unbiased molecular dynamics method for predicting ligands binding mechanisms: probing the binding pathway of dasatinib to Src-kinase

SUMMARY

Small molecules such as metabolites and drugs play essential roles in biological processes and pharmaceutical industry. Knowing their interactions with biomacromolecular targets demands a deep understanding of binding mechanisms. Dozens of papers have suggested that discovering of the binding event by means of conventional unbiased molecular dynamics (MD) simulation urges considerable amount of computational resources, therefore, only one who holds a cluster or a supercomputer can afford such extensive simulations. Thus, many researchers who do not own such resources are reluctant to take the benefits of running unbiased MD simulation, in full atomistic details, when studying a ligand binding pathway. Many researchers are impelled to be content with biased MD simulations which seek its validation due to its intrinsic preconceived framework. In this work, we have presented a workable stratagem to encourage everyone to perform unbiased (unguided) MD simulations, in this case a protein-ligand binding process, by typical desktop computers and so achieve valuable results in nanosecond time scale. Here, we have described a dynamical binding's process of an anticancer drug, the dasatinib, to the c-Src kinase in full atomistic details for the first time, without applying any biasing force or potential which may lead the drug to artificial interactions with the protein. We have attained multiple independent binding events which occurred in the nanosecond time scales, surprisingly as little as ∼30 ns. Both the protonated and deprotonated forms of the dasatinib reached the crystallographic binding mode without having any major intermediate state during induction.

AVAILABILITY AND IMPLEMENTATION

The links of the tutorial and technical documents are accessible in the article.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

TopSuite Web Server: A Meta-Suite for Deep-Learning-Based Protein Structure and Quality Prediction

Proteins carry out the most fundamental processes of life such as cellular metabolism, regulation, and communication. Understanding these processes at a molecular level requires knowledge of their three-dimensional structures. Experimental techniques such as X-ray crystallography, NMR spectroscopy, and cryogenic electron microscopy can resolve protein structures but are costly and time-consuming and do not work for all proteins. Computational protein structure prediction tries to overcome these problems by predicting the structure of a new protein using existing protein structures as a resource. Here we present TopSuite, a web server for protein model quality assessment (TopScore) and template-based protein structure prediction (TopModel). TopScore provides meta-predictions for global and residue-wise model quality estimation using deep neural networks. TopModel predicts protein structures using a top-down consensus approach to aid the template selection and subsequently uses TopScore to refine and assess the predicted structures. The TopSuite Web server is freely available at https://cpclab.uni-duesseldorf.de/topsuite/.

Rational drug repurposing for cancer by inclusion of the unbiased molecular dynamics simulation in the structure-based virtual screening approach: Challenges and breakthroughs

Managing cancer is now one of the biggest concerns of health organizations. Many strategies have been developed in drug discovery pipelines to help rectify this problem and two of the best ones are drug repurposing and computational methods. The combination of these approaches can have immense impact on the course of drug discovery. In silico drug repurposing can significantly reduce the time, the cost and the effort of drug development. Computational methods such as structure-based drug design (SBDD) and virtual screening can predict the potentials of small molecule binders, such as drugs, for having favorable effect on a particular molecular target. However, the demand for accuracy and efficiency of SBDD requires more sophisticated and complicated approaches such as unbiased molecular dynamics (UMD) simulation that has been recently introduced. As a complementary strategy, the knowledge acquired from UMD simulations can increase the chance of finding the right candidates and the pipeline of its administration is introduced and discussed in this review. An elaboration of this pipeline is also made by detailing an example, the binding and unbinding pathways of dasatinib-c-Src kinase complex, which shows that how influential this method can be in rational drug repurposing in cancer treatment.

Quantitative assessment of the determinant structural differences between redox-active and inactive glutaredoxins

Class I glutaredoxins are enzymatically active, glutathione-dependent oxidoreductases, whilst class II glutaredoxins are typically enzymatically inactive, Fe-S cluster-binding proteins. Enzymatically active glutaredoxins harbor both a glutathione-scaffold site for reacting with glutathionylated disulfide substrates and a glutathione-activator site for reacting with reduced glutathione. Here, using yeast ScGrx7 as a model protein, we comprehensively identified and characterized key residues from four distinct protein regions, as well as the covalently bound glutathione moiety, and quantified their contribution to both interaction sites. Additionally, we developed a redox-sensitive GFP2-based assay, which allowed the real-time assessment of glutaredoxin structure-function relationships inside living cells. Finally, we employed this assay to rapidly screen multiple glutaredoxin mutants, ultimately enabling us to convert enzymatically active and inactive glutaredoxins into each other. In summary, we have gained a comprehensive understanding of the mechanistic underpinnings of glutaredoxin catalysis and have elucidated the determinant structural differences between the two main classes of glutaredoxins.

Glutaredoxins play a central role in numerous biological processes including cellular redox homeostasis and Fe-S cluster biogenesis. Here the authors establish the molecular basis for glutaredoxin redox catalysis through comprehensive biochemical and structural analyses.

Resolving dynamics and function of transient states in single enzyme molecules

We use a hybrid fluorescence spectroscopic toolkit to monitor T4 Lysozyme (T4L) in action by unraveling the kinetic and dynamic interplay of the conformational states. In particular, by combining single-molecule and ensemble multiparameter fluorescence detection, EPR spectroscopy, mutagenesis, and FRET-positioning and screening, and other biochemical and biophysical tools, we characterize three short-lived conformational states over the ns-ms timescale. The use of 33 FRET-derived distance sets, to screen available T4L structures, reveal that T4L in solution mainly adopts the known open and closed states in exchange at 4 µs. A newly found minor state, undisclosed by, at present, more than 500 crystal structures of T4L and sampled at 230 µs, may be actively involved in the product release step in catalysis. The presented fluorescence spectroscopic toolkit will likely accelerate the development of dynamic structural biology by identifying transient conformational states that are highly abundant in biology and critical in enzymatic reactions.

T4 Lysozyme (T4L) is a model protein whose structure is extensively studied. Here the authors combine single-molecule and ensemble FRET measurements, FRET-positioning and screening and EPR spectroscopy to study the structural dynamics of T4L and describe its conformational landscape during the catalytic cycle by an extended Michaelis–Menten mechanism and identify an excited conformational state of the enzyme.

Structural Model of the ETR1 Ethylene Receptor Transmembrane Sensor Domain

The structure, mechanism of action and copper stoichiometry of the transmembrane sensor domain of the plant ethylene receptor ETR1 and homologs have remained elusive, hampering the understanding on how the perception of the plant hormone ethylene is transformed into a downstream signal. We generated the first structural model of the transmembrane sensor domain of ETR1 by integrating ab initio structure prediction and coevolutionary information. To refine and independently validate the model, we determined protein-related copper stoichiometries on purified receptor preparations and explored the helix arrangement by tryptophan scanning mutagenesis. All-atom molecular dynamics simulations of the dimeric model reveal how ethylene can bind proximal to the copper ions in the receptor, illustrating the initial stages of the ethylene perception process.

Targeting HSP90 dimerization via the C terminus is effective in imatinib-resistant CML and lacks the heat shock response

Heat shock protein 90 (HSP90) stabilizes many client proteins, including the BCR-ABL1 oncoprotein. BCR-ABL1 is the hallmark of chronic myeloid leukemia (CML) in which treatment-free remission (TFR) is limited, with clinical and economic consequences. Thus, there is an urgent need for novel therapeutics that synergize with current treatment approaches. Several inhibitors targeting the N-terminal domain of HSP90 are under investigation, but side effects such as induction of the heat shock response (HSR) and toxicity have so far precluded their US Food and Drug Administration approval. We have developed a novel inhibitor (aminoxyrone [AX]) of HSP90 function by targeting HSP90 dimerization via the C-terminal domain. This was achieved by structure-based molecular design, chemical synthesis, and functional preclinical in vitro and in vivo validation using CML cell lines and patient-derived CML cells. AX is a promising potential candidate that induces apoptosis in the leukemic stem cell fraction (CD34+CD38-) as well as the leukemic bulk (CD34+CD38+) of primary CML and in tyrosine kinase inhibitor (TKI)-resistant cells. Furthermore, BCR-ABL1 oncoprotein and related pro-oncogenic cellular responses are downregulated, and targeting the HSP90 C terminus by AX does not induce the HSR in vitro and in vivo. We also probed the potential of AX in other therapy-refractory leukemias. Therefore, AX is the first peptidomimetic C-terminal HSP90 inhibitor with the potential to increase TFR in TKI-sensitive and refractory CML patients and also offers a novel therapeutic option for patients with other types of therapy-refractory leukemia because of its low toxicity profile and lack of HSR.

Recognition motif and mechanism of ripening inhibitory peptides in plant hormone receptor ETR1

Synthetic peptides derived from ethylene-insensitive protein 2 (EIN2), a central regulator of ethylene signalling, were recently shown to delay fruit ripening by interrupting protein-protein interactions in the ethylene signalling pathway. Here, we show that the inhibitory peptide NOP-1 binds to the GAF domain of ETR1 - the prototype of the plant ethylene receptor family. Site-directed mutagenesis and computational studies reveal the peptide interaction site and a plausible molecular mechanism for the ripening inhibition.

Chlorflavonin Targets Acetohydroxyacid Synthase Catalytic Subunit IlvB1 for Synergistic Killing of Mycobacterium tuberculosis

The flavonoid natural compound chlorflavonin was isolated from the endophytic fungus Mucor irregularis, which was obtained from the Cameroonian medicinal plant Moringa stenopetala. Chlorflavonin exhibited strong growth inhibitory activity in vitro against Mycobacterium tuberculosis (MIC₉₀ 1.56 μM) while exhibiting no cytotoxicity toward the human cell lines MRC-5 and THP-1 up to concentrations of 100 μM. Mapping of resistance-mediating mutations employing whole-genome sequencing, chemical supplementation assays, and molecular docking studies as well as enzymatic characterization revealed that chlorflavonin specifically inhibits the acetohydroxyacid synthase catalytic subunit IlvB1, causing combined auxotrophies to branched-chain amino acids and to pantothenic acid. While exhibiting a bacteriostatic effect in monotreatment, chlorflavonin displayed synergistic effects with the first-line antibiotic isoniazid and particularly with delamanid, leading to a complete sterilization in liquid culture in combination treatment. Using a fluorescent reporter strain, intracellular activity of chlorflavonin against Mycobacterium tuberculosis inside infected macrophages was demonstrated and was superior to streptomycin treatment.

Determinants of FIV and HIV Vif sensitivity of feline APOBEC3 restriction factors

Background

Feline immunodeficiency virus (FIV) is a global pathogen of Felidae species and a model system for Human immunodeficiency virus (HIV)-induced AIDS. In felids such as the domestic cat (Felis catus), APOBEC3 (A3) genes encode for single-domain A3Z2s, A3Z3 and double-domain A3Z2Z3 anti-viral cytidine deaminases. The feline A3Z2Z3 is expressed following read-through transcription and alternative splicing, introducing a previously untranslated exon in frame, encoding a domain insertion called linker. Only A3Z3 and A3Z2Z3 inhibit Vif-deficient FIV. Feline A3s also are restriction factors for HIV and Simian immunodeficiency viruses (SIV). Surprisingly, HIV-2/SIV Vifs can counteract feline A3Z2Z3.

Results

To identify residues in feline A3s that Vifs need for interaction and degradation, chimeric human–feline A3s were tested. Here we describe the molecular direct interaction of feline A3s with Vif proteins from cat FIV and present the first structural A3 model locating these interaction regions. In the Z3 domain we have identified residues involved in binding of FIV Vif, and their mutation blocked Vif-induced A3Z3 degradation. We further identified additional essential residues for FIV Vif interaction in the A3Z2 domain, allowing the generation of FIV Vif resistant A3Z2Z3. Mutated feline A3s also showed resistance to the Vif of a lion-specific FIV, indicating an evolutionary conserved Vif–A3 binding. Comparative modelling of feline A3Z2Z3 suggests that the residues interacting with FIV Vif have, unlike Vif-interacting residues in human A3s, a unique location at the domain interface of Z2 and Z3 and that the linker forms a homeobox-like domain protruding of the Z2Z3 core. HIV-2/SIV Vifs efficiently degrade feline A3Z2Z3 by possible targeting the linker stretch connecting both Z-domains.

Conclusions

Here we identified in feline A3s residues important for binding of FIV Vif and a unique protein domain insertion (linker). To understand Vif evolution, a structural model of the feline A3 was developed. Our results show that HIV Vif binds human A3s differently than FIV Vif feline A3s. The linker insertion is suggested to form a homeo-box domain, which is unique to A3s of cats and related species, and not found in human and mouse A3s. Together, these findings indicate a specific and different A3 evolution in cats and human.

Electronic supplementary material

The online version of this article (doi:10.1186/s12977-016-0274-9) contains supplementary material, which is available to authorized users.

Assessing the potential of atomistic molecular dynamics simulations to probe reversible protein-protein recognition and binding

Protein-protein recognition and binding are governed by diffusion, noncovalent forces and conformational flexibility, entangled in a way that only molecular dynamics simulations can dissect at high resolution. Here we exploited ubiquitin's noncovalent dimerization equilibrium to assess the potential of atomistic simulations to reproduce reversible protein-protein binding, by running submicrosecond simulations of systems with multiple copies of the protein at millimolar concentrations. The simulations essentially fail because they lead to aggregates, yet they reproduce some specificity in the binding interfaces as observed in known covalent and noncovalent ubiquitin dimers. Following similar observations in literature we hint at electrostatics and water descriptions as the main liable force field elements, and propose that their optimization should consider observables relevant to multi-protein systems and unfolded proteins. Within limitations, analysis of binding events suggests salient features of protein-protein recognition and binding, to be retested with improved force fields. Among them, that specific configurations of relative direction and orientation seem to trigger fast binding of two molecules, even over 50 Å distances; that conformational selection can take place within surface-to-surface distances of 10 to 40 Å i.e. well before actual intermolecular contact; and that establishment of contacts between molecules further locks their conformations and relative orientations.

Peptide Binding to a PDZ Domain by Electrostatic Steering via Nonnative Salt Bridges

We have captured the binding of a peptide to a PDZ domain by unbiased molecular dynamics simulations. Analysis of the trajectories reveals on-pathway encounter complex formation, which is driven by electrostatic interactions between negatively charged carboxylate groups in the peptide and positively charged side chains surrounding the binding site. In contrast, the final stereospecific complex, which matches the crystal structure, features completely different interactions, namely the burial of the hydrophobic side chain of the peptide C-terminal residue and backbone hydrogen bonds. The simulations show that nonnative salt bridges stabilize kinetically the encounter complex during binding. Unbinding follows the inverse sequence of events with the same nonnative salt bridges in the encounter complex. Thus, in contrast to protein folding, which is driven by native interactions, the binding of charged peptides can be steered by nonnative interactions, which might be a general mechanism, e.g., in the recognition of histone tails by bromodomains.

Molecular dynamics simulations and structure-guided mutagenesis provide insight into the architecture of the catalytic core of the ectoine hydroxylase

Many bacteria amass compatible solutes to fend-off the detrimental effects of high osmolarity on cellular physiology and water content. These solutes also function as stabilizers of macromolecules, a property for which they are referred to as chemical chaperones. The tetrahydropyrimidine ectoine is such a compatible solute and is widely synthesized by members of the Bacteria. Many ectoine producers also synthesize the stress protectant 5-hydroxyectoine from the precursor ectoine, a process that is catalyzed by the ectoine hydroxylase (EctD). The EctD enzyme is a member of the non-heme-containing iron(II) and 2-oxoglutarate-dependent dioxygenase superfamily. A crystal structure of the EctD protein from the moderate halophile Virgibacillus salexigens has previously been reported and revealed the coordination of the iron catalyst, but it lacked the substrate ectoine and the co-substrate 2-oxoglutarate. Here we used this crystal structure as a template to assess the likely positioning of the ectoine and 2-oxoglutarate ligands within the active site by structural comparison, molecular dynamics simulations, and site-directed mutagenesis. Collectively, these approaches suggest the positioning of the iron, ectoine, and 2-oxoglutarate ligands in close proximity to each other and with a spatial orientation that will allow the region-selective and stereo-specific hydroxylation of (4S)-ectoine to (4S,5S)-5-hydroxyectoine. Our study thus provides a view into the catalytic core of the ectoine hydroxylase and suggests an intricate network of interactions between the three ligands and evolutionarily highly conserved residues in members of the EctD protein family.