Backbone Root Mean Square Research Articles

A Molecular Dynamics Simulation Study of the Effects of βGln114 Mutation on the Dynamic Behavior of the Catalytic Site of the Tryptophan Synthase.

L-tryptophan (l-Trp), a vital amino acid for the survival of various organisms, is synthesized by the enzyme tryptophan synthase (TS) in organisms such as eubacteria, archaebacteria, protista, fungi, and plantae. TS, a pyridoxal 5'-phosphate (PLP)-dependent enzyme, comprises α and β subunits that typically form an α2β2 tetramer. The enzyme's activity is regulated by the conformational switching of its α and β subunits between the open (T state) and closed (R state) conformations. Many microorganisms rely on TS for growth and replication, making the enzyme and the l-Trp biosynthetic pathway potential drug targets. For instance, Mycobacterium tuberculosis, Chlamydiae bacteria, Streptococcus pneumoniae, Francisella tularensis, Salmonella bacteria, and Cryptosporidium parasitic protozoa depend on l-Trp synthesis. Antibiotic-resistant salmonella strains have emerged, underscoring the need for novel drugs targeting the l-Trp biosynthetic pathway, especially for salmonella-related infections. A single amino acid mutation can significantly impact enzyme function, affecting stability, conformational dynamics, and active or allosteric sites. These changes influence interactions, catalytic activity, and protein-ligand/protein-protein interactions. This study focuses on the impact of mutating the βGln114 residue on the catalytic and allosteric sites of TS. Extensive molecular dynamics simulations were conducted on E(PLP), E(AEX1), E(A-A), and E(C3) forms of TS using the WT, βQ114A, and βQ114N versions. The results show that both the βQ114A and βQ114N mutations increase protein backbone root mean square deviation fluctuations, destabilizing all TS forms. Conformational and hydrogen bond analyses suggest the significance of βGln114 drifting away from cofactor/intermediates and forming hydrogen bonds with water molecules necessary for l-Trp biosynthesis. The βQ114A mutation creates a gap between βAla114 and cofactor/intermediates, hindering hydrogen bond formation due to short side chains and disrupting β-sites. Conversely, the βQ114N mutation positions βAsn114 closer to cofactor/intermediates, forming hydrogen bonds with O3 of cofactors/intermediates and nearby water molecules, potentially disrupting the l-Trp biosynthetic mechanism.

Fibroblast growth factor 5 (FGF5) and its missense mutant FGF5-H174 underlying trichomegaly: a molecular dynamics simulation investigation

The missense mutation Y174H of FGF5 (FGF5-H174) had been associated with trichomegaly, characterized by abnormally long and pigmented eyelashes. The amino acid tyrosine (Tyr/Y) at position 174 is conserved across many species, proposedly holding important characteristics for the functions of FGF5. Microsecond molecular dynamics simulations along with protein-protein docking and residue interacting network analysis were employed to investigate the structural dynamics and binding mode of both wild-type (FGF5-WT) and its mutated counterpart (FGF5-H174). It was found that the mutation decreased number of hydrogen bonds within the protein, sheet secondary structure, interaction of residue 174 with other residues, and number of salt-bridges. On the other hand, the mutation increased solvent accessible surface area, number of hydrogen bonds between the protein and solvent, coil secondary structure, protein C-alpha backbone root mean square deviation, protein residue root mean square fluctuations, as well as occupied conformational space. In addition, protein-protein docking integrated with molecular dynamics simulations and molecular mechanics - Poisson-Boltzmann surface area (MM/PBSA) binding energy calculation demonstrated that the mutated variant possessed stronger binding affinity towards fibroblast growth factor receptor 1 (FGFR1). However, residue interaction network analysis demonstrated that the binding mode of the FGFR1-FGF5-H174 complex was drastically different from that of the FGFR1-FGF5-WT complex. In conclusion, the missense mutation conferred more instability within itself and stronger binding affinity towards FGFR1 with distinctively altered binding mode or residue connectivity. These findings might help explain the decreased pharmacological activity of FGF5-H174 towards FGFR1, underlying trichomegaly. Communicated by Ramaswamy H. Sarma

Modeling TSH Receptor Dimerization at the Transmembrane Domain.

Biophysical studies have established that the thyrotropin (TSH) receptor (TSHR) undergoes posttranslational modifications including dimerization. Following our earlier simulation of a TSHR-transmembrane domain (TMD) monomer (called TSHR-TMD-TRIO) we have now proceeded with a molecular dynamics simulation (MD) of TSHR-TMD dimerization using this improved membrane-embedded model. The starting structure was the TMD protein with all extracellular and intracellular loops and internal waters, which was placed in the relative orientation of the model originally generated with Brownian dynamics. Furthermore, this model was embedded in a DPPC lipid bilayer further solvated with water and added salt. Data from the MD simulation studies showed that the dimeric subunits stayed in the same relative orientation and distance during the 1000 ns of study. Comparison of representative conformations of the individual monomers when dimerized with the conformations from the monomer simulation showed subtle differences as represented by the backbone root mean square deviations. Differences in the conformations of the ligand-binding sites, suggesting variable affinities for these "hot spots," were also revealed by comparing the docking scores of 46 small-molecule ligands that included known TSHR agonists and antagonists as well as their derivatives. These data add further insight into the tendency of the TSHR-TMD to form dimeric and oligomeric structures and show that the differing conformations influence small-molecule binding sites within the TMD.

Molecular modeling and dynamics studies of the synthetic small molecule agonists with GPR17 and P2Y1 receptor

The human Guanine Protein coupled membrane Receptor 17 (hGPR17), an orphan receptor that activates uracil nucleotides and cysteinyl leukotrienes is considered as a crucial target for the neurodegenerative diseases. Yet, the detailed molecular interaction of potential synthetic ligands of GPR17 needs to be characterized. Here, we have studied a comparative analysis on the interaction specificity of GPR17-ligands with hGPR17 and human purinergic G protein-coupled receptor (hP2Y1) receptors. Previously, we have simulated the interaction stability of synthetic ligands such as T0510.3657, AC1MLNKK, and MDL29951 with hGPR17 and hP2Y1 receptor in the lipid environment. In the present work, we have comparatively studied the protein-ligand interaction of hGPR17-T0510.3657 and P2Y1-MRS2500. Sequence analysis and structural superimposition of hGPR17 and hP2Y1 receptor revealed the similarities in the structural arrangement with the local backbone root mean square deviation (RMSD) value of 1.16 Å and global backbone RMSD value of 5.30 Å. The comparative receptor-ligand interaction analysis between hGPR17 and hP2Y1 receptor exposed the distinct binding sites in terms of geometrical properties. Further, the molecular docking of T0510.3657 with the hP2Y1 receptor have shown non-specific interaction. The experimental validation also revealed that Gi‐coupled activation of GPR17 by specific ligands leads to the adenylyl cyclase inhibition, while there is no inhibition upon hP2Y1 activation. Overall, the above findings suggest that T0510.3657-GPR17 binding specificity could be further explored for the treatment of numerous neuronal diseases. Communicated by Ramaswamy H. Sarma

In silico investigations of some Cyperus rotundus compounds as potential anti-inflammatory inhibitors of 5-LO and LTA4H enzymes

The present study aimed to experimentally identify the essential oil of Algerian Cyperus rotundus L. and to model the interaction of some known anti-inflammatory molecules with two key enzymes involved in inflammation, 5-Lypoxygenase (5-LO) and leukotriene A4 hydrolase (LTA4H). Gas chromatography/gas chromatography–mass spectrometry (GC/GC-MS) revealed that 92.7% of the essential oil contains 35 compounds, including oxygenated sesquiterpenes (44.2%), oxygenated monoterpenes (30.2%), monoterpene hydrocarbons (11.8%) and sesquiterpene hydrocarbons (6.5%). The major identified oxygenated terpenes are humulene oxide II, caryophyllene oxide, khusinol, agarospirol, spathulinol and trans-pinocarveol Myrtenol and α-terpineol are known to exhibit anti-inflammatory activities. Several complexes obtained after docking the natural terpenes with 5-LO and LTA4H have shown strong hydrogen bonding interactions. The best docking energies were found with α-terpineol, Myrtenol and khusinol. The interaction between the natural products and amino-acid residues HIS367, ILE673 and GLN363 appears to be critical for 5-LO inhibition, while the interaction with residues GLU271, HIS295, TYR383, TYR378, GLU318, GLU296 and ASP375 is critical for LTA4H inhibition. Molecular dynamics (MD) trajectories of the selected docked complexes showed stable backbone root mean square deviation (RMSD), supporting the stability of the natural product–enzyme interaction. Communicated by Ramaswamy H. Sarma

Exploring the Binding Interaction of Raf Kinase Inhibitory Protein With the N-Terminal of C-Raf Through Molecular Docking and Molecular Dynamics Simulation.

Protein-protein interactions are indispensable physiological processes regulating several biological functions. Despite the availability of structural information on protein-protein complexes, deciphering their complex topology remains an outstanding challenge. Raf kinase inhibitory protein (RKIP) has gained substantial attention as a favorable molecular target for numerous pathologies including cancer and Alzheimer’s disease. RKIP interferes with the RAF/MEK/ERK signaling cascade by endogenously binding with C-Raf (Raf-1 kinase) and preventing its activation. In the current investigation, the binding of RKIP with C-Raf was explored by knowledge-based protein-protein docking web-servers including HADDOCK and ZDOCK and a consensus binding mode of C-Raf/RKIP structural complex was obtained. Molecular dynamics (MD) simulations were further performed in an explicit solvent to sample the conformations for when RKIP binds to C-Raf. Some of the conserved interface residues were mutated to alanine, phenylalanine and leucine and the impact of mutations was estimated by additional MD simulations and MM/PBSA analysis for the wild-type (WT) and constructed mutant complexes. Substantial decrease in binding free energy was observed for the mutant complexes as compared to the binding free energy of WT C-Raf/RKIP structural complex. Furthermore, a considerable increase in average backbone root mean square deviation and fluctuation was perceived for the mutant complexes. Moreover, per-residue energy contribution analysis of the equilibrated simulation trajectory by HawkDock and ANCHOR web-servers was conducted to characterize the key residues for the complex formation. One residue each from C-Raf (Arg398) and RKIP (Lys80) were identified as the druggable “hot spots” constituting the core of the binding interface and corroborated by additional long-time scale (300 ns) MD simulation of Arg398Ala mutant complex. A notable conformational change in Arg398Ala mutant occurred near the mutation site as compared to the equilibrated C-Raf/RKIP native state conformation and an essential hydrogen bonding interaction was lost. The thirteen binding sites assimilated from the overall analysis were mapped onto the complex as surface and divided into active and allosteric binding sites, depending on their location at the interface. The acquired information on the predicted 3D structural complex and the detected sites aid as promising targets in designing novel inhibitors to block the C-Raf/RKIP interaction.

Pharmacophore Modelling-Based Drug Repurposing Approaches for SARS-CoV-2 Therapeutics.

The recent outbreak of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused a devastating effect globally with no effective treatment. The swift strategy to find effective treatment against coronavirus disease 2019 (COVID-19) is to repurpose the approved drugs. In this pursuit, an exhaustive computational method has been used on the DrugBank compounds targeting nsp16/nsp10 complex (PDB code: 6W4H). A structure-based pharmacophore model was generated, and the selected model was escalated to screen DrugBank database, resulting in three compounds. These compounds were subjected to molecular docking studies at the protein-binding pocket employing the CDOCKER module available with the Discovery Studio v18. In order to discover potential candidate compounds, the co-crystallized compound S-adenosyl methionine (SAM) was used as the reference compound. Additionally, the compounds remdesivir and hydroxycholoroquine were employed for comparative docking. The results have shown that the three compounds have demonstrated a higher dock score than the reference compounds and were upgraded to molecular dynamics simulation (MDS) studies. The MDS results demonstrated that the three compounds, framycetin, kanamycin, and tobramycin, are promising candidate compounds. They have represented a stable binding mode at the targets binding pocket with an average protein backbone root mean square deviation below 0.3 nm. Additionally, they have prompted the hydrogen bonds during the entire simulations, inferring that the compounds have occupied the active site firmly. Taken together, our findings propose framycetin, kanamycin, and tobramycin as potent putative inhibitors for COVID-19 therapeutics.

Crystal structure and mutational analysis of the human TRIM7 B30.2 domain provide insights into the molecular basis of its binding to glycogenin-1

Tripartite motif (TRIM)7 is an E3 ubiquitin ligase that was first identified through its interaction with glycogenin-1 (GN1), the autoglucosyltransferase that initiates glycogen biosynthesis. A growing body of evidence indicates that TRIM7 plays an important role in cancer development, viral pathogenesis, and atherosclerosis and, thus, represents a potential therapeutic target. TRIM family proteins share a multidomain architecture with a conserved N-terminal TRIM and a variable C-terminal domain. Human TRIM7 contains the canonical TRIM motif and a B30.2 domain at the C terminus. To contribute to the understanding of the mechanism of action of TRIM7, we solved the X-ray crystal structure of its B30.2 domain (TRIM7B30.2) in two crystal forms at resolutions of 1.6 Å and 1.8 Å. TRIM7B30.2 exhibits the typical B30.2 domain fold, consisting of two antiparallel β-sheets of seven and six strands, arranged as a distorted β-sandwich. Furthermore, two long loops partially cover the concave face of the β-sandwich defined by the β-sheet of six strands, thus forming a positively charged cavity. We used sequence conservation and mutational analyses to provide evidence of a putative binding interface for GN1. These studies showed that Leu423, Ser499, and Cys501 of TRIM7B30.2 and the C-terminal 33 amino acids of GN1 are critical for this binding interaction. Molecular dynamics simulations also revealed that hydrogen bond and hydrophobic interactions play a major role in the stability of a modeled TRIM7B30.2-GN1 C-terminal peptide complex. These data provide useful information that could be used to target this interaction for the development of potential therapeutic agents.

Computational identification of fertility functions of bovine Reprimo gene

Improvement in fertility is possible through gene assisted selection. Cattle fertility and genes underlying it should be thoroughly studied and exploited to find solution to declining cattle fertility. Reprimo (RPRM) gene is a pleiotropic gene involved in suppression of cancer, regulation of mitotic cell cycle, cell cycle arrest and regulation of survival. Comparison of protein tertiary structures is important in inferring functional characteristics of new proteins. This study used computational approach to identify some fertility functions of bovine RPRM gene using motif prediction and protein structure comparison. Amino acid sequences of bovine RPRM gene and some other cattle fertility genes were retrieved from GenBank. Motifs in the amino acid sequence of bovine RPRM gene were predicted using PROSITE software. The domain structure of bovine RPRM protein was predicted using simple modular architecture research tool (SMART). Protein tertiary structures (3D structures) of bovine RPRM gene and other cattle fertility genes were predicted with Phyre2 software. To have structural and functional similarity, it has been found that protein structure after superimposition should have Root Mean Square Deviation (RMSD) value less than or equal to 2Å. The protein 3D structures of other cattle fertility genes were superposed against the protein 3D structure of bovine RPRM gene using SuperPose web server and the proteins with RMSD value of 2Å or less were predicted as proteins with similar functions and structures as bovine RPRM gene. The predicted motifs (N-glycosylation site, N-myristoylation site, and cAMP and cGMP-dependent protein kinase phosphorylation site) and protein structure comparison revealed that, bovine RPRM gene and bovine growth hormone gene have the same fertility function with alpha carbon and backbone root mean square deviations of 1.94 Å and 1.81 Å, respectively. It follows therefore that other fertility functions of bovine RPRM gene included sexual maturation, steroidogenesis, gametogenesis, gonadal differentiation and gonadotrophin secretion which are the functions of growth hormone gene.

Multistep mutational transformation of a protein fold through structural intermediates.

New protein folds may evolve from existing folds through metamorphic evolution involving a dramatic switch in structure. To mimic pathways by which amino acid sequence changes could induce a change in fold, we designed two folded hybrids of Xfaso 1 and Pfl 6, a pair of homologous Cro protein sequences with ~40% identity but different folds (all-α vs. α + β, respectively). Each hybrid, XPH1 or XPH2, is 85% identical in sequence to its parent, Xfaso 1 or Pfl 6, respectively; 55% identical to its noncognate parent; and ~70% identical to the other hybrid. XPH1 and XPH2 also feature a designed hybrid chameleon sequence corresponding to the C-terminal region, which switched from α-helical to β-sheet structure during Cro evolution. We report solution nuclear magnetic resonance (NMR) structures of XPH1 and XPH2 at 0.3 Å and 0.5 Å backbone root mean square deviation (RMSD), respectively. XPH1 retains a global fold generally similar to Xfaso 1, and XPH2 retains a fold similar to Pfl 6, as measured by TM-align scores (~0.7), DALI Z-scores (7-9), and backbone RMSD (2-3 Å RMSD for the most ordered regions). However, these scores also indicate significant deviations in structure. Most notably, XPH1 and XPH2 have different, and intermediate, secondary structure content relative to Xfaso 1 and Pfl 6. The multistep progression in sequence, from Xfaso 1 to XPH1 to XPH2 to Pfl 6, thus involves both abrupt and gradual changes in folding pattern. The plasticity of some protein folds may allow for "polymetamorphic" evolution through intermediate structures.

Hydrophobicity-driven unfolding of Trp-cage encapsulated between graphene sheets.

Understanding the interaction between proteins and graphene not only helps elucidate the behaviors of proteins in confined geometries, but is also imperative to the development of a plethora of graphene-based biotechnologies, such as the graphene liquid cell transmission electron microscopy. To discuss the overall geometrical-thermal effects on proteins, we performed molecular dynamics simulations of hydrated Trp-cage miniprotein sandwiched between two graphene sheets and in the bulk environment at the temperatures below and above its unfolding temperature. The structural fluctuations of Trp-cage were characterized using the backbone root mean square displacement and the radius of gyration, from which the free energy landscape of Trp-cage was further constructed. We observed that at both temperatures the confined protein became adsorbed to the graphene surfaces and exhibited unfolded structures. Residue-specific analyses clearly showed the preference for the graphene to interact with the hydrophobic regions of Trp-cage. These results suggested that the conformation space accessible to the protein results from the competition between the thermodynamic driving forces and the geometrical restraints. While confinement usually tends to restrict the conformation of proteins by volume exclusion, it may also induce the unfolding of proteins by hydrophobic interactions.

Temperature effect on the structure and conformational fluctuations in two zinc knuckles from the mouse mammary tumor virus

Zinc fingers are small protein domains in which zinc plays a structural role, contributing to the stability of the zinc–peptide complex. Zinc fingers are structurally diverse and are present in proteins that perform a broad range of functions in various cellular processes, such as replication and repair, transcription and translation, metabolism and signaling, cell proliferation, and apoptosis. Zinc fingers typically function as interaction modules and bind to a wide variety of compounds, such as nucleic acids, proteins, and small molecules. In this study, we investigated the structural properties, in solution, of the proximal and distal zinc knuckles of the nucleocapsid (NC) protein from the mouse mammary tumor virus (MMTV) (MMTV NC). For this purpose, we performed a series of molecular dynamics simulations in aqueous solution at 300 K, 333 K, and 348 K. The temperature effect was evaluated in terms of root mean square deviation of the backbone atoms and root mean square fluctuation of the coordinating residue atoms. The stability of the zinc coordination sphere was analyzed based upon the time profile of the interatomic distances between the zinc ions and the chelator atoms. The results indicate that the hydrophobic character of the proximal zinc finger is dominant at 333 K. The low mobility of the coordinating residues suggests that the strong electrostatic effect exerted by the zinc ion on its coordinating residues is not influenced by the increase in temperature. The evolution of the structural parameters of the coordination sphere of the distal zinc finger at 300 K gives us a reasonable picture of the unfolding pathway, as proposed by Bombarda and coworkers (Bombarda et al., 2005), which can predict the binding order of the four conserved ligand-binding residues. Our results support the conclusion that the structural features can vary significantly between the two zinc knuckles of MMTV NC.

Constrained Cyclic Coordinate Descent for Cryo-EM Images at Medium Resolutions: Beyond the Protein Loop Closure Problem.

The cyclic coordinate descent (CCD) method is a popular loop closure method in protein structure modeling. It is a robotics algorithm originally developed for inverse kinematic applications. We demonstrate an effective method of building the backbone of protein structure models using the principle of CCD and a guiding trace. For medium-resolution 3-dimensional (3D) images derived using cryo-electron microscopy (cryo-EM), it is possible to obtain guiding traces of secondary structures and their skeleton connections. Our new method, constrained cyclic coordinate descent (CCCD), builds α-helices, β-strands, and loops quickly and fairly accurately along predefined traces. We show that it is possible to build the entire backbone of a protein fairly accurately when the guiding traces are accurate. In a test of 10 proteins, the models constructed using CCCD show an average of 3.91Å of backbone root mean square deviation (RMSD). When the CCCD method is incorporated in a simulated annealing framework to sample possible shift, translation, and rotation freedom, the models built with the true topology were ranked high on the list, with an average backbone RMSD100 of 3.76Å. CCCD is an effective method for modeling atomic structures after secondary structure traces and skeletons are extracted from 3D cryo-EM images.

PEPstrMOD: structure prediction of peptides containing natural, non-natural and modified residues.

BackgroundIn the past, many methods have been developed for peptide tertiary structure prediction but they are limited to peptides having natural amino acids. This study describes a method PEPstrMOD, which is an updated version of PEPstr, developed specifically for predicting the structure of peptides containing natural and non-natural/modified residues.ResultsPEPstrMOD integrates Forcefield_NCAA and Forcefield_PTM force field libraries to handle 147 non-natural residues and 32 types of post-translational modifications respectively by performing molecular dynamics using AMBER. AMBER was also used to handle other modifications like peptide cyclization, use of D-amino acids and capping of terminal residues. In addition, GROMACS was used to implement 210 non-natural side-chains in peptides using SwissSideChain force field library. We evaluated the performance of PEPstrMOD on three datasets generated from Protein Data Bank; i) ModPep dataset contains 501 non-natural peptides, ii) ModPep16, a subset of ModPep, and iii) CyclicPep contains 34 cyclic peptides. We achieved backbone Root Mean Square Deviation between the actual and predicted structure of peptides in the range of 3.81–4.05 Å.ConclusionsIn summary, the method PEPstrMOD has been developed that predicts the structure of modified peptide from the sequence/structure given as input. We validated the PEPstrMOD application using a dataset of peptides having non-natural/modified residues. PEPstrMOD offers unique advantages that allow the users to predict the structures of peptides having i) natural residues, ii) non-naturally modified residues, iii) terminal modifications, iv) post-translational modifications, v) D-amino acids, and also allows extended simulation of predicted peptides. This will help the researchers to have prior structural information of modified peptides to further design the peptides for desired therapeutic property. PEPstrMOD is freely available at http://osddlinux.osdd.net/raghava/pepstrmod/.ReviewersThis article was reviewed by Prof Michael Gromiha, Dr. Bojan Zagrovic and Dr. Zoltan Gaspari.Electronic supplementary materialThe online version of this article (doi:10.1186/s13062-015-0103-4) contains supplementary material, which is available to authorized users.

Insights on the structural perturbations in human MTHFR Ala222Val mutant by protein modeling and molecular dynamics

Methylenetetrahydrofolate reductase (MTHFR) protein catalyzes the only biochemical reaction which produces methyltetrahydrofolate, the active form of folic acid essential for several molecular functions. The Ala222Val polymorphism of human MTHFR encodes a thermolabile protein associated with increased risk of neural tube defects and cardiovascular disease. Experimental studies have shown that the mutation does not affect the kinetic properties of MTHFR, but inactivates the protein by increasing flavin adenine dinucleotide (FAD) loss. The lack of completely solved crystal structure of MTHFR is an impediment in understanding the structural perturbations caused by the Ala222Val mutation; computational modeling provides a suitable alternative. The three-dimensional structure of human MTHFR protein was obtained through homology modeling, by taking the MTHFR structures from Escherichia coli and Thermus thermophilus as templates. Subsequently, the modeled structure was docked with FAD using Glide, which revealed a very good binding affinity, authenticated by a Glide XP score of −10.3983 (kcal mol−1). The MTHFR was mutated by changing Alanine 222 to Valine. The wild-type MTHFR-FAD complex and the Ala222Val mutant MTHFR-FAD complex were subjected to molecular dynamics simulation over 50 ns period. The average difference in backbone root mean square deviation (RMSD) between wild and mutant variant was found to be ~.11 Å. The greater degree of fluctuations in the mutant protein translates to increased conformational stability as a result of mutation. The FAD-binding ability of the mutant MTHFR was also found to be significantly lowered as a result of decreased protein grip caused by increased conformational flexibility. The study provides insights into the Ala222Val mutation of human MTHFR that induces major conformational changes in the tertiary structure, causing a significant reduction in the FAD-binding affinity.

A Disulfide Stabilized β-Sandwich Defines the Structure of a New Cysteine Framework M-Superfamily Conotoxin.

The structure of a new cysteine framework (-C-CC-C-C-C-) "M"-superfamily conotoxin, Mo3964, shows it to have a β-sandwich structure that is stabilized by inter-sheet cross disulfide bonds. Mo3964 decreases outward K(+) currents in rat dorsal root ganglion neurons and increases the reversal potential of the NaV1.2 channels. The structure of Mo3964 (PDB ID: 2MW7 ) is constructed from the disulfide connectivity pattern, i.e., 1-3, 2-5, and 4-6, that is hitherto undescribed for the "M"-superfamily conotoxins. The tertiary structural fold has not been described for any of the known conus peptides. NOE (549), dihedral angle (84), and hydrogen bond (28) restraints, obtained by measurement of (h3)JNC' scalar couplings, were used as input for structure calculation. The ensemble of structures showed a backbone root mean square deviation of 0.68 ± 0.18 Å, with 87% and 13% of the backbone dihedral (ϕ, ψ) angles lying in the most favored and additional allowed regions of the Ramachandran map. The conotoxin Mo3964 represents a new bioactive peptide fold that is stabilized by disulfide bonds and adds to the existing repertoire of scaffolds that can be used to design stable bioactive peptide molecules.

Comparative Modeling and Molecular Dynamics Simulation of Substrate Binding in Human Fatty Acid Synthase: Enoyl Reductase and β-Ketoacyl Reductase Catalytic Domains.

Fatty acid synthase (FASN, EC 2.3.1.85), is a multi-enzyme dimer complex that plays a critical role in lipogenesis. This lipogenic enzyme has gained importance beyond its physiological role due to its implications in several clinical conditions-cancers, obesity, and diabetes. This has made FASN an attractive pharmacological target. Here, we have attempted to predict the theoretical models for the human enoyl reductase (ER) and β-ketoacyl reductase (KR) domains based on the porcine FASN crystal structure, which was the structurally closest template available at the time of this study. Comparative modeling methods were used for studying the structure-function relationships. Different validation studies revealed the predicted structures to be highly plausible. The respective substrates of ER and KR domains-namely, trans-butenoyl and β-ketobutyryl-were computationally docked into active sites using Glide in order to understand the probable binding mode. The molecular dynamics simulations of the apo and holo states of ER and KR showed stable backbone root mean square deviation trajectories with minimal deviation. Ramachandran plot analysis showed 96.0% of residues in the most favorable region for ER and 90.3% for the KR domain, respectively. Thus, the predicted models yielded significant insights into the substrate binding modes of the ER and KR catalytic domains and will aid in identifying novel chemical inhibitors of human FASN that target these domains.

High‐resolution modeling of antibody structures by a combination of bioinformatics, expert knowledge, and molecular simulations

In the second antibody modeling assessment, we used a semiautomated template-based structure modeling approach for 11 blinded antibody variable region (Fv) targets. The structural modeling method involved several steps, including template selection for framework and canonical structures of complementary determining regions (CDRs), homology modeling, energy minimization, and expert inspection. The submitted models for Fv modeling in Stage 1 had the lowest average backbone root mean square deviation (RMSD) (1.06 Å). Comparison to crystal structures showed the most accurate Fv models were generated for 4 out of 11 targets. We found that the successful modeling in Stage 1 mainly was due to expert-guided template selection for CDRs, especially for CDR-H3, based on our previously proposed empirical method (H3-rules) and the use of position specific scoring matrix-based scoring. Loop refinement using fragment assembly and multicanonical molecular dynamics (McMD) was applied to CDR-H3 loop modeling in Stage 2. Fragment assembly and McMD produced putative structural ensembles with low free energy values that were scored based on the OSCAR all-atom force field and conformation density in principal component analysis space, respectively, as well as the degree of consensus between the two sampling methods. The quality of 8 out of 10 targets improved as compared with Stage 1. For 4 out of 10 Stage-2 targets, our method generated top-scoring models with RMSD values of less than 1 Å. In this article, we discuss the strengths and weaknesses of our approach as well as possible directions for improvement to generate better predictions in the future.

Inhibitory Mechanism of an Allosteric Antibody Targeting the Glucagon Receptor

Elevated glucagon levels and increased hepatic glucagon receptor (GCGR) signaling contribute to hyperglycemia in type 2 diabetes. We have identified a monoclonal antibody that inhibits GCGR, a class B G-protein coupled receptor (GPCR), through a unique allosteric mechanism. Receptor inhibition is mediated by the binding of this antibody to two distinct sites that lie outside of the glucagon binding cleft. One site consists of a patch of residues that are surface-exposed on the face of the extracellular domain (ECD) opposite the ligand-binding cleft, whereas the second binding site consists of residues in the αA helix of the ECD. A docking model suggests that the antibody does not occlude the ligand-binding cleft. We solved the crystal structure of GCGR ECD containing a naturally occurring G40S mutation and found a shift in the register of the αA helix that prevents antibody binding. We also found that alterations in the αA helix impact the normal function of GCGR. We present a model for the allosteric inhibition of GCGR by a monoclonal antibody that may form the basis for the development of allosteric modulators for the treatment of diabetes and other class B GPCR-related diseases.

Ab initio structure prediction of the antibody hypervariable H3 loop

Antibodies have the capability of binding a wide range of antigens due to the diversity of the six loops constituting the complementarity determining region (CDR). Among the six loops, the H3 loop is the most diverse in structure, length, and sequence identity. Prediction of the three-dimensional structures of antibodies, especially the CDR loops, is an important step in the computational design and engineering of novel antibodies for improved affinity and specificity. Although it has been demonstrated that the conformation of the five non-H3 loops can be accurately predicted by comparing their sequences against databases of canonical loop conformations, no such connection has been established for H3 loops. In this work, we present the results for ab initio structure prediction of the H3 loop using conformational sampling and energy calculations with the program Prime on a dataset of 53 loops ranging in length from 4 to 22 residues. When the prediction is performed in the crystal environment and including symmetry mates, the median backbone root mean square deviation (RMSD) is 0.5 Å to the crystal structure, with 91% of cases having an RMSD of less than 2.0 Å. When the prediction is performed in a noncrystallographic environment, where the scaffold is constructed by swapping the H3 loops between homologous antibodies, 70% of cases have an RMSD below 2.0 Å. These results show promise for ab initio loop predictions applied to modeling of antibodies.