Dual function of Zika virus NS2B-NS3 protease.

Zika virus (ZIKV) serine protease, indispensable for viral polyprotein processing and replication, is composed of the membrane-anchored NS2B polypeptide and the N-terminal domain of the NS3 polypeptide (NS3pro). The C-terminal domain of the NS3 polypeptide (NS3hel) is necessary for helicase activity and contains an ATP-binding site. We discovered that ZIKV NS2B-NS3pro binds single-stranded RNA with a Kd of ~0.3 µM, suggesting a novel function. We tested various structural modifications of NS2B-NS3pro and observed that constructs stabilized in the recently discovered "super-open" conformation do not bind RNA. Likewise, stabilizing NS2B-NS3pro in the "closed" (proteolytically active) conformation using substrate inhibitors abolished RNA binding. We posit that RNA binding occurs when ZIKV NS2B-NS3pro adopts the "open" conformation, which we modeled using highly homologous dengue NS2B-NS3pro crystallized in the open conformation. We identified two positively charged fork-like structures present only in the open conformation of NS3pro. These forks are conserved across Flaviviridae family and could be aligned with the positively charged grove on NS3hel, providing a contiguous binding surface for the negative RNA strand exiting helicase. We propose a "reverse inchworm" model for a tightly intertwined NS2B-NS3 helicase-protease machinery, which suggests that NS2B-NS3pro cycles between open and super-open conformations to bind and release RNA enabling long-range NS3hel processivity. The transition to the closed conformation, likely induced by the substrate, enables the classical protease activity of NS2B-NS3pro.

S-nitrosylated TDP-43 triggers aggregation, cell-to-cell spread, and neurotoxicity in hiPSCs and in vivo models of ALS/FTD.

Rare genetic mutations result in aggregation and spreading of cognate proteins in neurodegenerative disorders; however, in the absence of mutation (i.e., in the vast majority of "sporadic" cases), mechanisms for protein misfolding/aggregation remain largely unknown. Here, we show environmentally induced nitrosative stress triggers protein aggregation and cell-to-cell spread. In patient brains with amyotrophic lateral sclerosis (ALS)/frontotemporal dementia (FTD), aggregation of the RNA-binding protein TDP-43 constitutes a major component of aberrant cytoplasmic inclusions. We identify a pathological signaling cascade whereby reactive nitrogen species cause S-nitrosylation of TDP-43 (forming SNO-TDP-43) to facilitate disulfide linkage and consequent TDP-43 aggregation. Similar pathological SNO-TDP-43 levels occur in postmortem human FTD/ALS brains and in cell-based models, including human-induced pluripotent stem cell (hiPSC)-derived neurons. Aggregated TDP-43 triggers additional nitrosative stress, representing positive feed forward leading to further SNO-TDP-43 formation and disulfide-linked oligomerization/aggregation. Critically, we show that these redox reactions facilitate cell spreading in vivo and interfere with the TDP-43 RNA-binding activity, affecting SNMT1 and phospho-(p)CREB levels, thus contributing to neuronal damage in ALS/FTD disorders.

S-Nitrosylation of p62 Inhibits Autophagic Flux to Promote α-Synuclein Secretion and Spread in Parkinson's Disease and Lewy Body Dementia.

Dysregulation of autophagic pathways leads to accumulation of abnormal proteins and damaged organelles in many neurodegenerative disorders, including Parkinson's disease (PD) and Lewy body dementia (LBD). Autophagy-related dysfunction may also trigger secretion and spread of misfolded proteins, such as α-synuclein (α-syn), the major misfolded protein found in PD/LBD. However, the mechanism underlying these phenomena remains largely unknown. Here, we used cell-based models, including human induced pluripotent stem cell-derived neurons, CRISPR/Cas9 technology, and male transgenic PD/LBD mice, plus vetting in human postmortem brains (both male and female). We provide mechanistic insight into this pathologic pathway. We find that aberrant S-nitrosylation of the autophagic adaptor protein p62 causes inhibition of autophagic flux and intracellular buildup of misfolded proteins, with consequent secretion resulting in cell-to-cell spread. Thus, our data show that pathologic protein S-nitrosylation of p62 represents a critical factor not only for autophagic inhibition and demise of individual neurons, but also for α-syn release and spread of disease throughout the nervous system.SIGNIFICANCE STATEMENT In Parkinson's disease and Lewy body dementia, dysfunctional autophagy contributes to accumulation and spread of aggregated α-synuclein. Here, we provide evidence that protein S-nitrosylation of p62 inhibits autophagic flux, contributing to α-synuclein aggregation and spread.

Quantitative profiling of protease specificity.

Proteases are an important class of enzymes, whose activity is central to many physiologic and pathologic processes. Detailed knowledge of protease specificity is key to understanding their function. Although many methods have been developed to profile specificities of proteases, few have the diversity and quantitative grasp necessary to fully define specificity of a protease, both in terms of substrate numbers and their catalytic efficiencies. We have developed a concept of "selectome"; the set of substrate amino acid sequences that uniquely represent the specificity of a protease. We applied it to two closely related members of the Matrixin family-MMP-2 and MMP-9 by using substrate phage display coupled with Next Generation Sequencing and information theory-based data analysis. We have also derived a quantitative measure of substrate specificity, which accounts for both the number of substrates and their relative catalytic efficiencies. Using these advances greatly facilitates elucidation of substrate selectivity between closely related members of a protease family. The study also provides insight into the degree to which the catalytic cleft defines substrate recognition, thus providing basis for overcoming two of the major challenges in the field of proteolysis: 1) development of highly selective activity probes for studying proteases with overlapping specificities, and 2) distinguishing targeted proteolysis from bystander proteolytic events.

Stress tensor and constant pressure simulation for polarizable Gaussian multipole model.

Our previous article has established the theory of molecular dynamics (MD) simulations for systems modeled with the polarizable Gaussian multipole (pGM) electrostatics [Wei et al., J. Chem. Phys. 153(11), 114116 (2020)]. Specifically, we proposed the covalent basis vector framework to define the permanent multipoles and derived closed-form energy and force expressions to facilitate an efficient implementation of pGM electrostatics. In this study, we move forward to derive the pGM internal stress tensor for constant pressure MD simulations with the pGM electrostatics. Three different formulations are presented for the flexible, rigid, and short-range screened systems, respectively. The analytical formulations were implemented in the SANDER program in the Amber package and were first validated with the finite-difference method for two different boxes of pGM water molecules. This is followed by a constant temperature and constant pressure MD simulation for a box of 512 pGM water molecules. Our results show that the simulation system stabilized at a physically reasonable state and maintained the balance with the externally applied pressure. In addition, several fundamental differences were observed between the pGM and classic point charge models in terms of the simulation behaviors, indicating more extensive parameterization is necessary to utilize the pGM electrostatics.

Determining the atomic charge of calcium ion requires the information of its coordination geometry in an EF-hand motif.

It is challenging to parameterize the force field for calcium ions (Ca2+) in calcium-binding proteins because of their unique coordination chemistry that involves the surrounding atoms required for stability. In this work, we observed a wide variation in Ca2+ binding loop conformations of the Ca2+-binding protein calmodulin, which adopts the most populated ternary structures determined from the molecular dynamics simulations, followed by ab initio quantum mechanical (QM) calculations on all 12 amino acids in the loop that coordinate Ca2+ in aqueous solution. Ca2+ charges were derived by fitting to the electrostatic potential in the context of a classical or polarizable force field (PFF). We discovered that the atomic radius of Ca2+ in conventional force fields is too large for the QM calculation to capture the variation in the coordination geometry of Ca2+ in its ionic form, leading to unphysical charges. Specifically, we found that the fitted atomic charges of Ca2+ in the context of PFF depend on the coordinating geometry of electronegative atoms from the amino acids in the loop. Although nearby water molecules do not influence the atomic charge of Ca2+, they are crucial for compensating for the coordination of Ca2+ due to the conformational flexibility in the EF-hand loop. Our method advances the development of force fields for metal ions and protein binding sites in dynamic environments.

Efficient formulation of polarizable Gaussian multipole electrostatics for biomolecular simulations.

Molecular dynamics simulations of biomolecules have been widely adopted in biomedical studies. As classical point-charge models continue to be used in routine biomolecular applications, there have been growing demands on developing polarizable force fields for handling more complicated biomolecular processes. Here, we focus on a recently proposed polarizable Gaussian Multipole (pGM) model for biomolecular simulations. A key benefit of pGM is its screening of all short-range electrostatic interactions in a physically consistent manner, which is critical for stable charge-fitting and is needed to reproduce molecular anisotropy. Another advantage of pGM is that each atom's multipoles are represented by a single Gaussian function or its derivatives, allowing for more efficient electrostatics than other Gaussian-based models. In this study, we present an efficient formulation for the pGM model defined with respect to a local frame formed with a set of covalent basis vectors. The covalent basis vectors are chosen to be along each atom's covalent bonding directions. The new local frame can better accommodate the fact that permanent dipoles are primarily aligned along covalent bonds due to the differences in electronegativity of bonded atoms. It also allows molecular flexibility during molecular simulations and facilitates an efficient formulation of analytical electrostatic forces without explicit torque computation. Subsequent numerical tests show that analytical atomic forces agree excellently with numerical finite-difference forces for the tested system. Finally, the new pGM electrostatics algorithm is interfaced with the particle mesh Ewald (PME) implementation in Amber for molecular simulations under the periodic boundary conditions. To validate the overall pGM/PME electrostatics, we conducted an NVE simulation for a small water box of 512 water molecules. Our results show that to achieve energy conservation in the polarizable model, it is important to ensure enough accuracy on both PME and induction iteration. It is hoped that the reformulated pGM model will facilitate the development of future force fields based on the pGM electrostatics for applications in biomolecular systems and processes where polarization plays crucial roles.

Interactions between motor domains in kinesin-14 Ncd - a molecular dynamics study.

Minus-end directed, non-processive kinesin-14 Ncd is a dimeric protein with C-terminally located motor domains (heads). Generation of the power-stroke by Ncd consists of a lever-like rotation of a long superhelical 'stalk' segment while one of the kinesin's heads is bound to the microtubule. The last ∼30 amino acids of Ncd head play a crucial but still poorly understood role in this process. Here, we used accelerated molecular dynamics simulations to explore the conformational dynamics of several systems built upon two crystal structures of Ncd, the asymmetrical T436S mutant in pre-stroke/post-stroke conformations of two partner subunits and the symmetrical wild-type protein in pre-stroke conformation of both subunits. The results revealed a new conformational state forming following the inward motion of the subunits and stabilized with several hydrogen bonds to residues located on the border or within the C-terminal linker, i.e. a modeled extension of the C-terminus by residues 675-683. Forming of this new, compact Ncd conformation critically depends on the length of the C-terminus extending to at least residue 681. Moreover, the associative motion leading to the compact conformation is accompanied by a partial lateral rotation of the stalk. We propose that the stable compact conformation of Ncd may represent an initial state of the working stroke.

Interaction of the cryptic fragment of myelin basic protein with mitochondrial voltage-dependent anion-selective channel-1 affects cell energy metabolism.

In demyelinating nervous system disorders, myelin basic protein (MBP), a major component of the myelin sheath, is proteolyzed and its fragments are released in the neural environment. Here, we demonstrated that, in contrast with MBP, the cellular uptake of the cryptic 84-104 epitope (MBP84-104) did not involve the low-density lipoprotein receptor-related protein-1, a scavenger receptor. Our pull-down assay, mass spectrometry and molecular modeling studies suggested that, similar with many other unfolded and aberrant proteins and peptides, the internalized MBP84-104 was capable of binding to the voltage-dependent anion-selective channel-1 (VDAC-1), a mitochondrial porin. Molecular modeling suggested that MBP84-104 directly binds to the N-terminal α-helix located midway inside the 19 ß-blade barrel of VDAC-1. These interactions may have affected the mitochondrial functions and energy metabolism in multiple cell types. Notably, MBP84-104 caused neither cell apoptosis nor affected the total cellular ATP levels, but repressed the aerobic glycolysis (lactic acid fermentation) and decreased the l-lactate/d-glucose ratio (also termed as the Warburg effect) in normal and cancer cells. Overall, our findings implied that because of its interactions with VDAC-1, the cryptic MBP84-104 peptide invoked reprogramming of the cellular energy metabolism that favored enhanced cellular activity, rather than apoptotic cell death. We concluded that the released MBP84-104 peptide, internalized by the cells, contributes to the reprogramming of the energy-generating pathways in multiple cell types.

Matrix metalloproteinases - From the cleavage data to the prediction tools and beyond.

Understanding the physiological role of any protease requires identification of both its cleavage substrates and their relative cleavage efficacy as compared with other substrates and other proteinases. Our review manuscript is focused on the cleavage preferences of the individual matrix metalloproteinases (MMPs) and the cleavage similarity and distinction that exist in the human MMP family. The recent in-depth analysis of MMPs by us and many others greatly increased knowledge of the MMP biology and structural-functional relationships among this protease family members. A better knowledge of cleavage preferences of MMPs has led us to the development of the prediction tools that are now capable of the high throughput reliable prediction and ranking the MMP cleavage sites in the peptide sequences in silico. Our software unifies and consolidates volumes of the pre-existing data. Now this prediction-ranking in silico tool is ready to be used by others. The software we developed may facilitate both the identification of the novel proteolytic regulatory pathways and the discovery of the previously uncharacterized substrates of the individual MMPs. Because now the MMP research may be based on the mathematical probability parameters rather than on either random luck or common sense alone, the researchers armed with this novel in silico tool will be better equipped to fine-tune or, at least, to sharply focus their wet chemistry experiments. This article is part of a Special Issue entitled: Matrix Metalloproteinases edited by Rafael Fridman.

Effect of phosphorylation and single nucleotide polymorphisms on caspase substrates processing.

Posttranslational modifications that involve either reversible covalent modification of proteins or irreversible proteolysis are central to the regulation of key cellular mechanisms, including apoptosis, cell-cycle regulation and signal transduction. There is mounting evidence suggesting cross-talk between proteases and kinases. For instance: caspases, a class of proteases involved in programmed cell death-apoptosis, cleave a large set of various types of proteins. Simultaneously, kinases restrict caspase activity by phosphorylating their protein substrates in the vicinity of cleavage site. In addition, the caspase cleavage pattern in target proteins may be modified as a result of single nucleotide polymorphisms (SNPs) in the coding gene. This may either create a novel cleavage site, or increase/decrease the cleavage efficiency of a substrate. Such point mutations are often associated with the onset of disease. In this study, we predicted how phosphorylation and SNPs affect known human caspase proteolytic events collected in the CASBAH and Degrabase databases by applying Random Forest caspases' substrates prediction method, as implemented in the CaspDB, and the molecular dynamics free energy simulations approach. Our analysis confirms several experimental observations. Phosphorylation could have both positive or negative regulatory effects depending on its position with respect to the caspase cleavage site. For instance, we demonstrate that phosphorylation at P1' is the most detrimental for proteolytic efficiency of caspases. Phosphorylation at the P2 and P2' positions also negatively affect the cleavage events. In addition, we uncovered SNPs in 11 caspase substrates capable of completely abolishing the cleavage site due to polymorphism at the P1 position. The findings presented here may be useful for determining the link between aberrant proteolysis and disease.

Induced polarization restricts the conformational distribution of a light-harvesting molecular triad in the ground state.

The light-harvesting molecular triad consisting of carotenoid polyene (C), diaryl-porphyrin (P) and pyrrole-fullerene (C60) is a donor-acceptor molecule capable of absorbing incident light in the visible range. Its ability to convert solar energy to electrical excitation and charge separation energy suggests a great potential in real-world applications. The ensemble of its conformations under ambient conditions varies widely according to its electronic state. In previous work, we applied a non-polarizable model to study the conformational distribution of the molecular triad in the ground and charge separated states. However, due to the lack of polarization, which imparts subtle changes in the charge distribution on atoms, molecular simulations fail to produce accurate average dipole moments. We developed the first polarizable model for a molecular triad to investigate the structural and dynamic properties of a molecular triad in the ground state in an explicit organic solvent, tetrahydrofuran (THF). We performed first-principles electronic structure calculations of the individual components in the triad as well as THF and then fit the partial atomic charges to the electrostatic potential using the i-RESP methodology. We validated these force field parameters by comparing the thermodynamic and dynamic properties obtained from molecular dynamics simulations with those from experiments. We enhanced the sampling of the triad conformations with replica exchange molecular dynamics simulations. We characterized the effects of induced polarization on the structural stability of the triad by analyzing the free energy landscapes constructed with polarizable force fields. Furthermore, by using principal component analysis, we found that the molecular triad conformations adopted a small range of torsional angles with induced polarization. The triad conformation solvated in polar solvent with a polarizable force field qualitatively agrees with that obtained from nuclear magnetic resonance spectroscopy.

Basis for substrate recognition and distinction by matrix metalloproteinases.

Genomic sequencing and structural genomics produced a vast amount of sequence and structural data, creating an opportunity for structure-function analysis in silico [Radivojac P, et al. (2013) Nat Methods 10(3):221-227]. Unfortunately, only a few large experimental datasets exist to serve as benchmarks for function-related predictions. Furthermore, currently there are no reliable means to predict the extent of functional similarity among proteins. Here, we quantify structure-function relationships among three phylogenetic branches of the matrix metalloproteinase (MMP) family by comparing their cleavage efficiencies toward an extended set of phage peptide substrates that were selected from ∼ 64 million peptide sequences (i.e., a large unbiased representation of substrate space). The observed second-order rate constants [k(obs)] across the substrate space provide a distance measure of functional similarity among the MMPs. These functional distances directly correlate with MMP phylogenetic distance. There is also a remarkable and near-perfect correlation between the MMP substrate preference and sequence identity of 50-57 discontinuous residues surrounding the catalytic groove. We conclude that these residues represent the specificity-determining positions (SDPs) that allowed for the expansion of MMP proteolytic function during evolution. A transmutation of only a few selected SDPs proximal to the bound substrate peptide, and contributing the most to selectivity among the MMPs, is sufficient to enact a global change in the substrate preference of one MMP to that of another, indicating the potential for the rational and focused redesign of cleavage specificity in MMPs.

Matrix Metalloproteinase (MMP) Proteolysis of the Extracellular Loop of Voltage-gated Sodium Channels and Potential Alterations in Pain Signaling.

Congenital insensitivity to pain (CIP) or congenital analgesia is a rare monogenic hereditary condition. This disorder is characterized by the inability to perceive any form of pain. Nonsense mutations in Nav.1.7, the main pain signaling voltage-gated sodium channel, lead to its truncations and, consequently, to the inactivation of the channel functionality. However, a non-truncating homozygously inherited missense mutation in a Bedouin family with CIP (Nav1.7-R907Q) has also been reported. Based on our currently acquired in-depth knowledge of matrix metalloproteinase (MMP) cleavage preferences, we developed the specialized software that predicts the presence of the MMP cleavage sites in the peptide sequences. According to our in silico predictions, the peptide sequence of the exposed extracellular unstructured region linking the S5-S6 transmembrane segments in the DII domain of the human Nav1.7 sodium channel is highly sensitive to MMP-9 proteolysis. Intriguingly, the CIP R907Q mutation overlaps with the predicted MMP-9 cleavage site sequence. Using MMP-9 proteolysis of the wild-type, CIP, and control peptides followed by mass spectrometry of the digests, we demonstrated that the mutant sequence is severalfold more sensitive to MMP-9 proteolysis relative to the wild type. Because of the substantial level of sequence homology among sodium channels, our data also implicate MMP proteolysis in regulating the cell surface levels of the Nav1.7, Nav1.6, and Nav1.8 channels, but not Nav1.9. It is likely that the aberrantly accelerated MMP-9 proteolysis during neurogenesis is a biochemical rational for the functional inactivation in Nav1.7 and that the enhanced cleavage of the Nav1.7-R907Q mutant is a cause of CIP in the Bedouin family.

Sequence-derived structural features driving proteolytic processing.

Proteolytic signaling, or regulated proteolysis, is an essential part of many important pathways such as Notch, Wnt, and Hedgehog. How the structure of the cleaved substrate regions influences the efficacy of proteolytic processing remains underexplored. Here, we analyzed the relative importance in proteolysis of various structural features derived from substrate sequences using a dataset of more than 5000 experimentally verified proteolytic events captured in CutDB. Accessibility to the solvent was recognized as an essential property of a proteolytically processed polypeptide chain. Proteolytic events were found nearly uniformly distributed among three types of secondary structure, although with some enrichment in loops. Cleavages in α-helices were found to be relatively abundant in regions apparently prone to unfolding, while cleavages in ß-structures tended to be located at the periphery of ß-sheets. Application of the same statistical procedures to proteolytic events divided into separate sets according to the catalytic classes of proteases proved consistency of the results and confirmed that the structural mechanisms of proteolysis are universal. The estimated prediction power of sequence-derived structural features, which turned out to be sufficiently high, presents a rationale for their use in bioinformatic prediction of proteolytic events.

Substrate cleavage profiling suggests a distinct function of Bacteroides fragilis metalloproteinases (fragilysin and metalloproteinase II) at the microbiome-inflammation-cancer interface.

Enterotoxigenic anaerobic Bacteroides fragilis is a significant source of inflammatory diarrheal disease and a risk factor for colorectal cancer. Two distinct metalloproteinase types (the homologous 1, 2, and 3 isoforms of fragilysin (FRA1, FRA2, and FRA3, respectively) and metalloproteinase II (MPII)) are encoded by the B. fragilis pathogenicity island. FRA was demonstrated to be important to pathogenesis, whereas MPII, also a potential virulence protein, remained completely uncharacterized. Here, we, for the first time, extensively characterized MPII in comparison with FRA3, a representative of the FRA isoforms. We employed a series of multiplexed peptide cleavage assays to determine substrate specificity and proteolytic characteristics of MPII and FRA. These results enabled implementation of an efficient assay of MPII activity using a fluorescence-quenched peptide and contributed to structural evidence for the distinct substrate cleavage preferences of MPII and FRA. Our data imply that MPII specificity mimics the dibasic Arg↓Arg cleavage motif of furin-like proprotein convertases, whereas the cleavage motif of FRA (Pro-X-X-Leu-(Arg/Ala/Leu)↓) resembles that of human matrix metalloproteinases. To the best of our knowledge, MPII is the first zinc metalloproteinase with the dibasic cleavage preferences, suggesting a high level of versatility of metalloproteinase proteolysis. Based on these data, we now suggest that the combined (rather than individual) activity of MPII and FRA is required for the overall B. fragilis virulence in vivo.

Assessment of Amino Acid Electrostatic Parametrizations of the Polarizable Gaussian Multipole Model.

Accurate parametrization of amino acids is pivotal for the development of reliable force fields for molecular modeling of biomolecules such as proteins. This study aims to assess amino acid electrostatic parametrizations with the polarizable Gaussian Multipole (pGM) model by evaluating the performance of the pGM-perm (with atomic permanent dipoles) and pGM-ind (without atomic permanent dipoles) variants compared to the traditional RESP model. The 100-conf-combterm fitting strategy on tetrapeptides was adopted, in which (1) all peptide bond atoms (-CO-NH-) share identical set of parameters and (2) the total charges of the two terminal N-acetyl (ACE) and N-methylamide (NME) groups were set to neutral. The accuracy and transferability of electrostatic parameters across peptides with varying lengths and real-world examples were examined. The results demonstrate the enhanced performance of the pGM-perm model in accurately representing the electrostatic properties of amino acids. This insight underscores the potential of the pGM-perm model and the 100-conf-combterm strategy for the future development of the pGM force field.

PCMRESP: A Method for Polarizable Force Field Parameter Development and Transferability of the Polarizable Gaussian Multipole Models Across Multiple Solvents.

The transferability of force field parameters is a crucial aspect of high-quality force fields. Previous investigations have affirmed the transferability of electrostatic parameters derived from polarizable Gaussian multipole models (pGMs) when applied to water oligomer clusters, polypeptides across various conformations, and different sequences. In this study, we introduce PCMRESP, a novel method for electrostatic parametrization in solution, intended for the development of polarizable force fields. We utilized this method to assess the transferability of three models: a fixed charge model and two variants of pGM models. Our analysis involved testing these models on 377 small molecules and 100 tetra-peptides in five representative dielectric environments: gas, diethyl ether, dichloroethane, acetone, and water. Our findings reveal that the inclusion of atomic polarization significantly enhances transferability and the incorporation of permanent atomic dipoles, in the form of covalent bond dipoles, leads to further improvements. Moreover, our tests on dual-solvent strategies demonstrate consistent transferability for all three models, underscoring the robustness of the dual-solvent approach. In contrast, an evaluation of the traditional HF/6-31G* method indicates poor transferability for the pGM-ind and pGM-perm models, suggesting the limitations of this conventional approach.

R.E.D. Server: a web service for deriving RESP and ESP charges and building force field libraries for new molecules and molecular fragments.

R.E.D. Server is a unique, open web service, designed to derive non-polarizable RESP and ESP charges and to build force field libraries for new molecules/molecular fragments. It provides to computational biologists the means to derive rigorously molecular electrostatic potential-based charges embedded in force field libraries that are ready to be used in force field development, charge validation and molecular dynamics simulations. R.E.D. Server interfaces quantum mechanics programs, the RESP program and the latest version of the R.E.D. tools. A two step approach has been developed. The first one consists of preparing P2N file(s) to rigorously define key elements such as atom names, topology and chemical equivalencing needed when building a force field library. Then, P2N files are used to derive RESP or ESP charges embedded in force field libraries in the Tripos mol2 format. In complex cases an entire set of force field libraries or force field topology database is generated. Other features developed in R.E.D. Server include help services, a demonstration, tutorials, frequently asked questions, Jmol-based tools useful to construct PDB input files and parse R.E.D. Server outputs as well as a graphical queuing system allowing any user to check the status of R.E.D. Server jobs.

Transferability of the Electrostatic Parameters of the Polarizable Gaussian Multipole Model.

Accuracy and transferability are the two highly desirable properties of molecular mechanical force fields. Compared with the extensively used point-charge additive force fields that apply fixed atom-centered point partial charges to model electrostatic interactions, polarizable force fields are thought to have the advantage of modeling the atomic polarization effects. Previous works have demonstrated the accuracy of the recently developed polarizable Gaussian multipole (pGM) models. In this work, we assessed the transferability of the electrostatic parameters of the pGM models with (pGM-perm) and without (pGM-ind) atomic permanent dipoles in terms of reproducing the electrostatic potentials surrounding molecules/oligomers absent from electrostatic parameterizations. Encouragingly, both the pGM-perm and pGM-ind models show significantly improved transferability than the additive model in the tests (1) from water monomer to water oligomer clusters; (2) across different conformations of amino acid dipeptides and tetrapeptides; (3) from amino acid tetrapeptides to longer polypeptides; and (4) from nucleobase monomers to Watson-Crick base pair dimers and tetramers. Furthermore, we demonstrated that the double-conformation fittings using amino acid tetrapeptides in the αR and ß conformations can result in good transferability not only across different tetrapeptide conformations but also from tetrapeptides to polypeptides with lengths ranging from 1 to 20 repetitive residues for both the pGM-ind and pGM-perm models. In addition, the observation that the pGM-ind model has significantly better accuracy and transferability than the point-charge additive model, even though they have an identical number of parameters, strongly suggest the importance of intramolecular polarization effects. In summary, this and previous works together show that the pGM models possess both accuracy and transferability, which are expected to serve as foundations for the development of next-generation polarizable force fields for modeling various polarization-sensitive biological systems and processes.