Improved Modeling of Peptidic Foldamers Using a Quantum Chemical Parametrization Based on Torsional Minimum Energy Path Matching

The increasing interest in novel foldamer constructs demands an accurate computational treatment on an extensive timescale. However, it is still a challenge to derive a force field (FF) that can reproduce the experimentally known fold while also allowing the spontaneous exploration of other structures. Here, aiming at a realistic reproduction of backbone torsional barriers, the relevant proper dihedrals of acyclic β2‐, β3‐ and β2,3‐amino acids were added to the CHARMM FF and optimized using a novel, self‐consistent iterative procedure based on quantum chemical relaxed scans. The new FF was validated by molecular dynamics simulations on three acyclic peptides. While they resided most of the time in their preferred fold (>80 % in helices and >50 % in hairpin), they also visited other conformations. Owing to the CHARMM36m‐consistent parametrization, the proposed extension is suitable for exploring new foldamer structures and assemblies, and their interactions with diverse biomolecules

Актуальність впровадження систем газового обліку в сучасних умовах

Free energy calculation has long been an important goal for molecular dynamics simulation and force field development, but historically it has been challenged by limited performance, accuracy, and creation of topologies for arbitrary small molecules. This has made it difficult to systematically compare different sets of parameters to improve existing force fields, but in the past few years several authors have developed increasingly automated procedures to generate parameters for force fields such as Amber, CHARMM, and OPLS. Here, we present a new framework that enables fully automated generation of GROMACS topologies for any of these force fields and an automated setup for parallel adaptive optimization of high-throughput free energy calculation by adjusting lambda point placement on the fly. As a small example of this automated pipeline, we have calculated solvation free energies of 50 different small molecules using the GAFF, OPLS-AA, and CGenFF force fields and four different water models, and by including the often neglected polarization costs, we show that the common charge models are somewhat underpolarized.QC 20150505</p

Molecular Basis of Ligand Dissociation in β-Adrenergic Receptors

The important and diverse biological functions of β-adrenergic receptors (βARs) have promoted the search for compounds to stimulate or inhibit their activity. In this regard, unraveling the molecular basis of ligand binding/unbinding events is essential to understand the pharmacological properties of these G protein-coupled receptors. In this study, we use the steered molecular dynamics simulation method to describe, in atomic detail, the unbinding process of two inverse agonists, which have been recently co-crystallized with β1 and β2ARs subtypes, along four different channels. Our results indicate that this type of compounds likely accesses the orthosteric binding site of βARs from the extracellular water environment. Importantly, reconstruction of forces and energies from the simulations of the dissociation process suggests, for the first time, the presence of secondary binding sites located in the extracellular loops 2 and 3 and transmembrane helix 7, where ligands are transiently retained by electrostatic and Van der Waals interactions. Comparison of the residues that form these new transient allosteric binding sites in both βARs subtypes reveals the importance of non-conserved electrostatic interactions as well as conserved aromatic contacts in the early steps of the binding process

Ligand binding to an Allergenic Lipid Transfer Protein Enhances Conformational Flexibility resulting in an Increase in Susceptibility to Gastroduodenal Proteolysis

Non-specific lipid transfer proteins (LTPs) are a family of lipid-binding molecules that are widely distributed across flowering plant species, many of which have been identified as allergens. They are highly resistant to simulated gastroduodenal proteolysis, a property that may play a role in determining their allergenicity and it has been suggested that lipid binding may further increase stability to proteolysis. It is demonstrated that LTPs from wheat and peach bind a range of lipids in a variety of conditions, including those found in the gastroduodenal tract. Both LTPs are initially cleaved during gastroduodenal proteolysis at three major sites between residues 39–40, 56–57 and 79–80, with wheat LTP being more resistant to cleavage than its peach ortholog. The susceptibility of wheat LTP to proteolyic cleavage increases significantly upon lipid binding. This enhanced digestibility is likely to be due to the displacement of Tyr79 and surrounding residues from the internal hydrophobic cavity upon ligand binding to the solvent exposed exterior of the LTP, facilitating proteolysis. Such knowledge contributes to our understanding as to how resistance to digestion can be used in allergenicity risk assessment of novel food proteins, including GMOs

Molecular Dynamics Analysis of Apolipoprotein-D - Lipid Hydroperoxide Interactions: Mechanism for Selective Oxidation of Met-93

Background: Recent studies suggest reduction of radical-propagating fatty acid hydroperoxides to inert hydroxides by interaction with apolipoprotein-D (apoD) Met93 may represent an antioxidant function for apoD. The nature and structural consequences of this selective interaction are unknown. Methodology/Principal Findings: Herein we used molecular dynamics (MD) analysis to address these issues. Longtimescale simulations of apoD suggest lipid molecules are bound flexibly, with the molecules free to explore multiple conformations in a binding site at the entrance to the classical lipocalin ligand-binding pocket. Models of 5s- 12s- and 15s hydroperoxyeicosatetraenoic acids were created and the lipids found to wrap around Met93 thus providing a plausible mechanism by which eicosatetraenoic acids bearing hydroperoxides on different carbon atoms can interact with Met93 to yield Met93 sulfoxide (Met93SO). Simulations of glycosylated apoD indicated that a second solvent exposed Met at position 49 was shielded by a triantennerary N-glycan attached to Asn45 thereby precluding lipid interactions. MD simulations of apoD showed B-factors of the loop containing Met93SO were higher in the oxidized protein, indicating increased flexibility that is predicted to destabilize the protein and promote self-association. Conclusions/Significance: These studies provide novel insights into the mechanisms that may contribute to the antioxidant function of apoD and the structural consequences that result if Met93SO is not redox-cycled back to its native state

Substrate binding and translocation of the serotonin transporter studied by docking and molecular dynamics simulations

The serotonin (5-HT) transporter (SERT) plays an important role in the termination of 5-HT-mediated neurotransmission by transporting 5-HT away from the synaptic cleft and into the presynaptic neuron. In addition, SERT is the main target for antidepressant drugs, including the selective serotonin reuptake inhibitors (SSRIs). The three-dimensional (3D) structure of SERT has not yet been determined, and little is known about the molecular mechanisms of substrate binding and transport, though such information is very important for the development of new antidepressant drugs. In this study, a homology model of SERT was constructed based on the 3D structure of a prokaryotic homologous leucine transporter (LeuT) (PDB id: 2A65). Eleven tryptamine derivates (including 5-HT) and the SSRI (S)-citalopram were docked into the putative substrate binding site, and two possible binding modes of the ligands were found. To study the conformational effect that ligand binding may have on SERT, two SERT–5-HT and two SERT–(S)-citalopram complexes, as well as the SERT apo structure, were embedded in POPC lipid bilayers and comparative molecular dynamics (MD) simulations were performed. Our results show that 5-HT in the SERT–5-HTB complex induced larger conformational changes in the cytoplasmic parts of the transmembrane helices of SERT than any of the other ligands. Based on these results, we suggest that the formation and breakage of ionic interactions with amino acids in transmembrane helices 6 and 8 and intracellular loop 1 may be of importance for substrate translocation

Molecular Dynamics Simulations Suggest Ligand’s Binding to Nicotinamidase/Pyrazinamidase

The research on the binding process of ligand to pyrazinamidase (PncA) is crucial for elucidating the inherent relationship between resistance of Mycobacterium tuberculosis and PncA’s activity. In the present study, molecular dynamics (MD) simulation methods were performed to investigate the unbinding process of nicotinamide (NAM) from two PncA enzymes, which is the reverse of the corresponding binding process. The calculated potential of mean force (PMF) based on the steered molecular dynamics (SMD) simulations sheds light on an optimal binding/unbinding pathway of the ligand. The comparative analyses between two PncAs clearly exhibit the consistency of the binding/unbinding pathway in the two enzymes, implying the universality of the pathway in all kinds of PncAs. Several important residues dominating the pathway were also determined by the calculation of interaction energies. The structural change of the proteins induced by NAM’s unbinding or binding shows the great extent interior motion in some homologous region adjacent to the active sites of the two PncAs. The structure comparison substantiates that this region should be very important for the ligand’s binding in all PncAs. Additionally, MD simulations also show that the coordination position of the ligand is displaced by one water molecule in the unliganded enzymes. These results could provide the more penetrating understanding of drug resistance of M. tuberculosis and be helpful for the development of new antituberculosis drugs

Towards the simulation of biomolecules: optimisation of peptide-capped glycine using FFLUX

YesThe optimisation of a peptide-capped glycine using the novel force field FFLUX is presented. FFLUX is a force field based on the machine-learning method kriging and the topological energy partitioning method called Interacting Quantum Atoms. FFLUX has a completely different architecture to that of traditional force fields, avoiding (harmonic) potentials for bonded, valence and torsion angles. In this study, FFLUX performs an optimisation on a glycine molecule and successfully recovers the target density-functional theory energy with an error of 0.89 ± 0.03 kJ mol−1. It also recovers the structure of the global minimum with a root-mean-squared deviation of 0.05 Å (excluding hydrogen atoms). We also show that the geometry of the intra-molecular hydrogen bond in glycine is recovered accurately.EPSRC Established Career Fellowship [grant number EP/K005472

Structural basis of PROTAC cooperative recognition for selective protein degradation

Inducing macromolecular interactions with small molecules to activate cellular signaling is a challenging goal. PROTACs (proteolysis-targeting chimeras) are bifunctional molecules that recruit a target protein in proximity to an E3 ubiquitin ligase to trigger protein degradation. Structural elucidation of the key ternary ligase-PROTAC-target species and its impact on target degradation selectivity remain elusive. We solved the crystal structure of Brd4 degrader MZ1 in complex with human VHL and the Brd4 bromodomain (Brd4BD2). The ligand folds into itself to allow formation of specific intermolecular interactions in the ternary complex. Isothermal titration calorimetry studies, supported by surface mutagenesis and proximity assays, are consistent with pronounced cooperative formation of ternary complexes with Brd4BD2. Structure-based-designed compound AT1 exhibits highly selective depletion of Brd4 in cells. Our results elucidate how PROTAC-induced de novo contacts dictate preferential recruitment of a target protein into a stable and cooperative complex with an E3 ligase for selective degradation

Anti-HIV-1 integrase potency of methylgallate from Alchornea cordifolia using in vitro and in silico approaches:

According to the 2018 report of the United Nations Programme on HIV/AIDS (UNAIDS), acquired immune deficiency syndrome (AIDS), a disease caused by the human immunodeficiency virus (HIV), remains a significant public health problem. The non-existence of a cure or effective vaccine for the disease and the associated emergence of resistant viral strains imply an urgent need for the discovery of novel anti-HIV drug candidates. The current study aimed to identify potential anti-retroviral compounds from Alchornea cordifolia