Search

Structural and charge interaction analyses of the binding site of the aspartate 95 mutants of the Desulfovibrio vulgaris flavodoxin reveals the role played by the hydrogen bonding in the modulation of the redox potentials of the bound flavodoxin, particularly the carboxylate group of the mutant.

Researchers describe the structural changes that occur as the lipase from Thermomyces lanuginosa is activated in the presence of lipids. During this process the enzyme exists in a low activity form, an activated form, and a fully active form.

The dynamics of the enzyme ribonuclease-Sa as a function of temperature is explored through a series of molecular dynamics simulations. The enzyme is preorganized dynamically, regions with low plasticity confer specificity while the more flexible regions have the fluidity to facilitate energetically inexpensive conformational rearrangements such as those required to achieve the transition state.

The rotation of a water molecule buried inside a protein bovine pancreatic trypsin inhibitor (BPTI) was examined. Reaction path calculations using conjugate peak refinement were employed to compute the transition state and activation energy for the rotation in both BPTI and ice. The behavior of the structural water, which results to the interchange of two water hydrogens, was found to be similar in both ice and BPTI. The complex process involves two successive rotations around the orthogonal axes.

NMR structure determination and molecular dynamics simulations are performed to understand the critical interactions between subunits B and F in the [A.sub.1][A.sub.O] ATP synthase of Methanosarcina mazei [G.sup.0]1. Movements of both helices of [B.sub.Mm] are observed to accommodate the incoming C-terminus of [F.sub.Mm] and connect the events of ion pumping and nucleotide binding in the [A.sub.1][A.sub.O] ATP synthase.

The molecular dynamics simulations along with experimental information were used to reproduce experimental trends in binding affinities of the transcription factor p53 peptides with the E3-ubiquitin ligase MDM2. The results offer insight into the origin of ensemble protein-peptide complex and importance of the adaptive surface of MDM2 for drug design.

The studies of papain, the member of the cysteine proteinase superfamily are carried out to demonstrate the role of noncovalent interactions in enzyme active center chemistry. The findings provide insight into the generation of substantial nucleophilic reactivity in the (Cys25)-[S.sup.-]/(His159)-[Im.sup.+]H ion pair by decreasing the level of mutual solvation of the anionic and cationic partners.

The solution studies and crystal structure analysis of N-lithocholyl insulin, in which residue Lys(super B29) has been acylated with a bile acid derivative, are presented. The unit cell in the crystal consists of an insulin hexamer containing two zinc ions, which have two m-cresol molecules, bound at each dimer-dimer interface stabilizing an R6 conformation.

A comparative analysis of properties extracted from two independent molecular dynamics (MD) simulations and from structures determined by Nuclear magnetic resonance (NMR) spectroscopy and X-ray diffraction crystallography (XRD) is presented. The simulations present a view of the dynamics of Rnase Sa that is different but complementary to NMR and X-ray.

The vibrational entropy change on the burial of a crystallographically well-ordered water molecule in bovine pancreatic trypsin inhibitor (BPTI) using normal-mode analysis is dissected. A detailed analysis is performed of how the translational and rotational degrees of freedom of the isolated water molecule are transformed into vibrational modes in the complex.

Reactivation of tumor-suppressing activity of p53 protein by targeting its negative regulator MDM2/MDMX has been pursued as a potential anticancer strategy. A promising dual inhibitor of MDM2/MDMX that has been developed and is currently in clinical trials is the stapled peptide ATSP-7041. The activity of this molecule is reported to be modulated in the presence of serum. Albumin is the most abundant protein in serum and is known to bind reversibly to several molecules. To study this interaction, we develop a protocol combining molecular modeling, docking, and simulations. Exhaustive docking of the peptide with representative simulated structures of human serum albumin led to the identification of probable binding sites on the surface of the protein, including both known canonical and novel binding sites. Sequence differences at putative peptide-binding sites in human and mouse albumin result in differing interaction energies with the peptide and enable us to rationalize the observed differences in vivo. In general, the findings should help in guiding the design of features in such peptides that may affect their distribution and cell permeability, opening a new window in structure-guided design strategies. [ABSTRACT FROM AUTHOR]

The addition of a phenol to a hexameric insulin solution leads to the formation of a stable hexamer whose stability arises from a rearrangement in which residues B1-B8 form an (alpha)-helix in the R-state. This R-state, in turn, is partly stabilized by the nonpolar interactions between the phenolic molecule and residue B5 His at the dimer-dimer surface. To test if a tyrosine side chain could mimic the effect of phenol binding and induce the R-state, a B5 His -> Tyr mutant human insulin was constructed. The results showed that it could not since the mutation destabilized the T-state.

The article discusses a study which demonstrated that the C-terminal end of the transactivation domain of p53 modulates binding to MDM2. This led researchers to isolate several p53 peptide variants of which the W23L and P27S mutants displayed the lowest and highest affinities for MDM2, with binding free energies 3.0 and 2.3 kcal/mol weaker and stronger than wild type p53.

Mitogen-activated protein kinase-interacting kinases 1 and 2 (MNK1/2) represent a central class of enzymes that are activated by extracellular signal-regulated kinase (ERK) or p38 mitogen-activated protein (MAP) kinases. MNK1 and MNK2 coordinate cellular signaling, control production of inflammatory chemokines, and regulate cell proliferation and survival. MNK1/2 are referred to as serine/threonine kinases as they phosphorylate serine or threonine residues on their substrates. Upon activation, MNK1/2 phosphorylate eukaryotic translation initiation factor 4E (eIF4E) at Ser209, which in turn initiates ribosome assembly and protein translation. Deleterious overexpression of MNK1/2 and/or eIF4E have been reported in several diseases including cancers, neurological disorders, autism, and inflammation. Recently, there have been intense efforts toward the development of potent and selective inhibitors of MNK1/2 in both academia and industry. Herein, we review the current understanding of the structural and biological aspects of MNK1/2 and provide an update of pharmacological inhibitors of MNK1/2 including candidates in clinical trials.

The evolution of antibiotic-resistant bacteria is an ongoing and troubling development that has increased the number of diseases and infections that risk going untreated. There is an urgent need to develop alternative strategies and treatments to address this issue. One class of molecules that is attracting significant interest is that of antimicrobial peptides (AMPs). Their design and development has been aided considerably by the applications of molecular models, and we review these here. These methods include the use of tools to explore the relationships between their structures, dynamics, and functions and the increasing application of machine learning and molecular dynamics simulations. This review compiles resources such as AMP databases, AMP-related web servers, and commonly used techniques, together aimed at aiding researchers in the area toward complementing experimental studies with computational approaches.

The increasing realization of the prevalence of liquid-liquid phase separation (LLPS) across multiple length scales of biological constructs, from intracellular membraneless organelles to extracellular load-bearing tissues, has raised intriguing questions about intermolecular interactions regulating LLPS at the atomic level. Squid-beak derived histidine (His)- and tyrosine (Tyr)-rich peptides (HB peps ) have recently emerged as suitable short model peptides to precisely assess the roles of peptide motifs and single residues on the phase behavior and material properties of microdroplets obtained by LLPS. In this study, by systematically introducing single mutations in an HB pep , we have identified specific sticker residues that attract peptide chains together. We find that His and Tyr residues located near the sequence termini drive phase separation, forming interaction nodes that stabilize microdroplets. Combining quantum chemistry simulations with NMR studies, we predict atomic-level bond geometries and uncover inter-residue supramolecular interactions governing LLPS. These results are subsequently used to propose possible topological arrangements of the peptide chains, which upon expansion can help explain the three-dimensional network of microdroplets. The stability of the proposed topologies carried out through all-atom molecular dynamics simulations predicts chain topologies that are more likely to stabilize the microdroplets. Overall, this study provides useful guidelines for the de novo design of peptide coacervates with tunable phase behavior and material properties. In addition, the analysis of nanoscale topologies may pave the way to understand how client molecules can be trapped within microdroplets, with direct implications for the encapsulation and controlled release of therapeutics for drug delivery applications.

During host cell infection by flaviviruses such as dengue and Zika, acidic pH within the endosome triggers a conformational change in the envelope protein on the outer surface of the virion. This results in exposure of the ∼15 residue fusion peptide (FP) region, freeing it to induce fusion between the viral and endosomal membranes. A better understanding of the conformational dynamics of the FP in the presence of membranes, and the basis for its selectivity for anionic lipid species present within the endosome, would facilitate its therapeutic targeting with antiviral drugs and antibodies. In this work, multiscale modeling, simulations, and free energy calculations (including a total of ∼75 μs of atomic-resolution sampling), combined with imaging total internal reflection fluorescence correlation spectroscopy experiments, were employed to investigate the mechanisms of interaction of FP variants with lipid bilayers. Wild-type FPs (in the presence or absence of a fluorescein isothiocyanate tag) were shown to possess a funneled conformational landscape governing their exit from solvent and penetration into the lipid phase and to exhibit an electrostatically favored >2-fold affinity for membranes containing anionic species over purely zwitterionic ones. Conversely, the landscape was abolished in a nonfunctional point mutant, leading to a 2-fold drop in host membrane affinity. Collectively, our data reveal how the highly conserved flavivirus FP has evolved to funnel its conformational space toward a maximally fusogenic state anchored within the endosomal membrane. Therapeutically targeting the accessible ensemble of FP conformations may represent a new, rational strategy for blocking viral infection.

Chronic myeloid leukemia (CML) is a myeloproliferative disease caused by bcr-abl1, a constitutively active tyrosine kinase fusion gene responsible for an abnormal proliferation of leukemic stem cells (LSCs). Inhibition of BCR-ABL1 kinase activity offers long-term relief to CML patients. However, for a proportion of them, BCR-ABL1 inhibition will become ineffective at treating the disease, and CML will progress to blast crisis (BC) CML with poor prognosis. BC-CML is often associated with excessive phosphorylated eukaryotic translation initiation factor 4E (eIF4E), which renders LSCs capable of proliferating via self-renewal, oblivious to BCR-ABL1 inhibition. In vivo, eIF4E is exclusively phosphorylated on Ser209 by MNK1/2. Consequently, a selective inhibitor of MNK1/2 should reduce the level of phosphorylated eIF4E and re-sensitize LSCs to BCR-ABL1 inhibition, thus hindering the proliferation of BC LSCs. We report herein the structure-activity relationships and pharmacokinetic properties of a selective MNK1/2 inhibitor clinical candidate, ETC-206, which in combination with dasatinib prevents BC-CML LSC self-renewal in vitro and enhances dasatinib antitumor activity in vivo.

When using non-natural amino acids in computational simulations of proteins, it is necessary to ensure appropriate parameterization of the new amino acids toward the creation of appropriate input files. In particular, the charges on the atoms may have to be derived de novo and ad hoc for the new species. As there are many variables in the charge derivation process, an investigation was devised to compare different approaches and determine their effect on simulations. This was done with the purpose to identify the methods which produced results compatible with the existing parameters. It was found in this study that all analyzed charge derivation methods reproduce with sufficient accuracy the literature values and can be used with confidence when parameterizing novel species., Competing Interests: The authors declare no competing financial interest.

Membrane-active antibiotics are of great interest in fighting bacterial resistance. α-Mangostin is a membrane-active molecule, but there are no details of its mechanism of action at the atomistic level. We have employed free-energy simulations and microsecond-long conventional molecular dynamics simulations to study the mode of interaction of α-mangostin with a model bacterial membrane and compare it with the mechanisms of three hydrophobic molecules (ciprofloxacin, xanthone, and tetracycline). We find that α-mangostin is thermodynamically more favored to insert into the membrane compared to the other three molecules. Apart from tetracycline, which is largely hydrophilic, the other three molecules aggregate in water; however, only α-mangostin can penetrate into the lipid tail region of the membrane. When it reaches a high concentration in the lipid tail region, α-mangostin can form tubular clusters that span the two head group regions of the membrane, resulting in a large number of water translocations along the transmembrane aggregates. Structure-activity relationship analysis revealed two structural properties that characterize α-mangostin, namely, the two isoprenyl groups and the polar groups present in the aromatic rings, which result in "disruptive amphiphilicity" and hence its excellent membrane activity., Competing Interests: The authors declare no competing financial interest.

Overexpression of the eukaryotic initiation factor 4E (eIF4E) is linked to a variety of cancers. Both mitogen-activated protein kinases-interacting kinases 1 and 2 (Mnk1/2) activate the oncogene eIF4E through posttranslational modification (phosphorylating it at the conserved Ser209). Inhibition of Mnk prevents eIF4E phosphorylation, making the Mnk-eIF4E axis a potential therapeutic target for oncology. Recently, the design and synthesis of a series of novel potent compounds inhibiting the Mnk1/2 kinases were carried out in-house. Here, we describe computational models of the interactions between Mnk1/2 kinases and these inhibitors. Molecular modeling combined with free energy calculations show that these compounds bind to the inactive forms of the kinases. All compounds adopt similar conformations in the catalytic sites of both kinases, stabilized by hydrogen bonds with the hinge regions and with the catalytic Lys78 (Mnk1) and Lys113 (Mnk2). These hydrogen bond interactions clearly play a critical role in determining the conformational stability and potency of the compounds. We also find that van der Waals interactions with an allosteric pocket are key to their binding and potency. Two distinct hydration sites that appear to further stabilize the ligand binding/interactions were observed. Critically, the inclusion of explicit water molecules in the calculations results in improving the agreement between calculated and experimental binding free energies., Competing Interests: The authors declare no competing financial interest.

A computational ligand-mapping approach to detect protein surface pockets that interact with hydrophobic moieties is presented. In this method, we incorporated benzene molecules into explicit solvent molecular dynamics simulations of various protein targets. The benzene molecules successfully identified the binding locations of hydrophobic hot-spot residues and all-hydrocarbon cross-links from known peptidic ligands. They also unveiled cryptic binding sites that are occluded by side chains and the protein backbone. Our results demonstrate that ligand-mapping molecular dynamics simulations hold immense promise to guide the rational design of peptidic modulators of protein-protein interactions, including that of stapled peptides, which show promise as an exciting new class of cell-penetrating therapeutic molecules.

By using a phage display derived peptide as an initial template, compounds have been developed that are highly specific against Mdm2/Mdm4. These compounds exhibit greater potency in p53 activation and protein-protein interaction assays than a compound derived from the p53 wild-type sequence. Unlike Nutlin, a small molecule inhibitor of Mdm2/Mdm4, the phage derived compounds can arrest cells resistant to p53 induced apoptosis over a wide concentration range without cellular toxicity, suggesting they are highly suitable for cyclotherapy.