A decision support system based on artificial intelligence and systems biology for the simulation of pancreatic cancer patient status.

Oncology treatments require continuous individual adjustment based on the measurement of multiple clinical parameters. Prediction tools exploiting the patterns present in the clinical data could be used to assist decision making and ease the burden associated to the interpretation of all these parameters. The goal of this study was to predict the evolution of patients with pancreatic cancer at their next visit using information routinely recorded in health records, providing a decision-support system for clinicians. We selected hematological variables as the visit's clinical outcomes, under the assumption that they can be predictive of the evolution of the patient. Multivariate models based on regression trees were generated to predict next-visit values for each of the clinical outcomes selected, based on the longitudinal clinical data as well as on molecular data sets streaming from in silico simulations of individual patient status at each visit. The models predict, with a mean prediction score (balanced accuracy) of 0.79, the evolution trends of eosinophils, leukocytes, monocytes, and platelets. Time span between visits and neutropenia were among the most common factors contributing to the predicted evolution. The inclusion of molecular variables from the systems-biology in silico simulations provided a molecular background for the observed variations in the selected outcome variables, mostly in relation to the regulation of hematopoiesis. In spite of its limitations, this study serves as a proof of concept for the application of next-visit prediction tools in real-world settings, even when available data sets are small.

On Quality Thresholds for the Clustering of Molecular Structures.

It has been recently suggested that diametral (so-called quality) similarity thresholds are superior to radial ones for the clustering of molecular three-dimensional structures (González-Alemán et al., 2020). The argument has been made for two clustering algorithms available in various software packages for the analysis of molecular structures from ensembles generated by computer simulations, attributed to Daura et al. (1999) (radial threshold) and Heyer et al. (1999) (diametral threshold). Here, we compare these two algorithms using the root-mean-squared difference (rmsd) between the Cartesian coordinates of selected atoms as pairwise similarity metric. We discuss formally the relation between these two methods and illustrate their behavior with two examples, a set of points in two dimensions and the coordinates of the tau polypeptide along a trajectory extracted from a replica-exchange molecular-dynamics simulation (Shea and Levine, 2016). We show that the two methods produce equally sized clusters as long as adequate choices are made for the respective thresholds. The real issue is not whether the threshold is radial or diametral but how to choose in either case a threshold value that is physically meaningful. We will argue that, when clustering molecular structures with the rmsd as a metric, the simplest best guess for a threshold is actually radial in nature.

Antibacterial Activity of T22, a Specific Peptidic Ligand of the Tumoral Marker CXCR4.

CXCR4 is a cytokine receptor used by HIV during cell attachment and infection. Overexpressed in the cancer stem cells of more than 20 human neoplasias, CXCR4 is a convenient antitumoral drug target. T22 is a polyphemusin-derived peptide and an effective CXCR4 ligand. Its highly selective CXCR4 binding can be exploited as an agent for the cell-targeted delivery and internalization of associated antitumor drugs. Sharing chemical and structural traits with antimicrobial peptides (AMPs), the capability of T22 as an antibacterial agent remains unexplored. Here, we have detected T22-associated antimicrobial activity and biofilm formation inhibition over Escherichia coli, Staphylococcus aureus and Pseudomonas aeruginosa, in a spectrum broader than the reference AMP GWH1. In contrast to GWH1, T22 shows neither cytotoxicity over mammalian cells nor hemolytic activity and is active when displayed on protein-only nanoparticles through genetic fusion. Under the pushing need for novel antimicrobial agents, the discovery of T22 as an AMP is particularly appealing, not only as its mere addition to the expanding catalogue of antibacterial drugs. The recognized clinical uses of T22 might allow its combined and multivalent application in complex clinical conditions, such as colorectal cancer, that might benefit from the synchronous destruction of cancer stem cells and local bacterial biofilms.

CNN-PepPred: an open-source tool to create convolutional NN models for the discovery of patterns in peptide sets-application to peptide-MHC class II binding prediction.

SUMMARY: The ability to unveil binding patterns in peptide sets has important applications in several biomedical areas, including the development of vaccines. We present an open-source tool, CNN-PepPred, that uses convolutional neural networks to discover such patterns, along with its application to peptide-HLA class II binding prediction. The tool can be used locally on different operating systems, with CPUs or GPUs, to train, evaluate, apply and visualize models. AVAILABILITY AND IMPLEMENTATION: CNN-PepPred is freely available as a Python tool with a detailed User's Guide at https://github.com/ComputBiol-IBB/CNN-PepPred. The site includes the peptide sets used in this study, extracted from IEDB (www.iedb.org). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Aggregation-prone peptides modulate activity of bovine interferon gamma released from naturally occurring protein nanoparticles.

Efficient protocols for the production of recombinant proteins are indispensable for the development of the biopharmaceutical sector. Accumulation of recombinant proteins in naturally-occurring protein aggregates is detrimental to biopharmaceutical development. In recent years, the view of protein aggregates has changed with the recognition that they are a valuable source of functional recombinant proteins. In this study, bovine interferon-gamma (rBoIFN-γ) was engineered to enhance the formation of protein aggregates, also known as protein nanoparticles (NPs), by the addition of aggregation-prone peptides (APPs) in the generally recognized as safe (GRAS) bacterial Lactococcus lactis expression system. The L6K2, HALRU and CYOB peptides were selected to assess their intrinsic aggregation capability to nucleate protein aggregation. These APPs enhanced the tendency of the resulting protein to aggregate at the expense of total protein yield. However, fine physico-chemical characterization of the resulting intracellular protein NPs, the protein released from them and the protein purified from the soluble cell fraction indicated that the compactability of protein conformations was directly related to the biological activity of variants of IFN-γ, used here as a model protein with therapeutic potential. APPs enhanced the aggregation tendency of fused rBoIFN-γ while increasing compactability of protein species. Biological activity of rBoIFN-γ was favored in more compacted conformations. Naturally-occurring protein aggregates can be produced in GRAS microorganisms as protein depots of releasable active protein. The addition of APPs to enhance the aggregation tendency has a positive impact in overall compactability and functionality of resulting protein conformers.

Interaction of camel Lactoferrin derived peptides with DNA: a molecular dynamics study.

BACKGROUND: Lactoferrampin (LFampin), Lactoferricin (LFcin), and LFchimera are three well-known antimicrobial peptides derived from Lactoferrin and proposed as alternatives for antibiotics. Although the intracellular activity of these peptides has been previously demonstrated, their mode of action is not yet fully understood. Here, we performed a molecular dynamics simulation study to understand the molecular interactions between camel Lactoferrin derived peptides, including CLFampin, CLFcin, and CLFchimera, and DNA as an important intracellular target. RESULTS: Our results indicate that all three peptides bind to DNA, albeit with different propensities, with CLFchimera showing the highest binding affinity. The secondary structures of the peptides, modeled on Lactoferrin, did not undergo significant changes during simulation, supporting their functional relevance. Main residues involved in the peptide-DNA interaction were identified based on binding free energy estimates calculated over 200 ns, which, as expected, confirmed strong electrostatic interactions between DNA phosphate groups and positively charged peptide side chains. Interaction between the different concentrations of CLFchimera and DNA revealed that after binding of four copies of CLFchimera to DNA, hydrogen bonds between the two strands of DNA start to break from one of the termini. CONCLUSIONS: Importantly, our results revealed that there is no DNA-sequence preference for peptide binding, in line with a broad antimicrobial activity. Moreover, the results showed that the strength of the interaction between DNA and CLFchimera is concentration dependent. The insight provided by these results can be used for the rational redesign of natural antimicrobial peptides targeting the bacterial DNA.

A Bowman-Birk protease inhibitor purified, cloned, sequenced and characterized from the seeds of Maclura pomifera (Raf.) Schneid.

MAIN CONCLUSION: A new BBI-type protease inhibitor with remarkable structural characteristics was purified, cloned, and sequenced from seeds of Maclura pomifera , a dicotyledonous plant belonging to the Moraceae family. In this work, we report a Bowman-Birk inhibitor (BBI) isolated, purified, cloned, and characterized from Maclura pomifera seeds (MpBBI), the first of this type from a species belonging to Moraceae family. MpBBI was purified to homogeneity by RP-HPLC, total RNA was extracted from seeds of M. pomifera, and the 3'RACE-PCR method was applied to obtain the cDNA, which was cloned and sequenced. Peptide mass fingerprinting (PMF) analysis showed correspondence between the in silico-translated protein and MpBBI, confirming that it corresponds to a new plant protease inhibitor. The obtained cDNA encoded a polypeptide of 65 residues and possesses 10 cysteine residues, with molecular mass of 7379.27, pI 6.10, and extinction molar coefficient of 9105 M-1 cm-1. MpBBI inhibits strongly trypsin with K i in the 10-10 M range and was stable in a wide array of pH and extreme temperatures. MpBBI comparative modeling was applied to gain insight into its 3D structure and highlighted some distinguishing features: (1) two non-identical loops, (2) loop 1 (CEEESRC) is completely different from any known BBI, and (3) the amount of disulphide bonds is also different from any reported BBI from dicot plants.

A Comparative Analysis of the Peptide Repertoires of HLA-DR Molecules Differentially Associated With Rheumatoid Arthritis.

OBJECTIVE: To evaluate similarity of the peptide repertoires bound to HLA-DR molecules that are differentially associated with rheumatoid arthritis (RA), and to define structural features of the shared peptides. METHODS: Peptide pools bound to HLA-DRB1*01:01, HLA-DRB1*04:01, and HLA-DRB1*10:01 (RA associated) and those bound to HLA-DRB1*15:01 (non-RA-associated) were purified and analyzed by liquid chromatography (LC) matrix-assisted laser desorption ionization-time-of-flight mass spectrometry (MS) and LC-ion-trap MS. Peptide pools from each allotype were compared in terms of size, protein origin, composition, and affinity (both theoretical and experimental with some peptides). Finally, 1 peptide sequenced from DR1, DR4, and DR10, but not from DR15, was modeled in complex with all 4 HLA-DRB1 molecules and HLA-DRB5*01:01. RESULTS: A total of 6,309 masses and 962 unique peptide sequences were compared. DR10 shared 29 peptides with DR1, 9 with DR4, and 1 with DR15; DR1 shared 6 peptides with DR4 and 9 with DR15; and DR4 and DR15 shared 4 peptides. The direct identification of peptide ligands indicated that DR1 and DR10 were the most similar molecules regarding the peptides that they could share. The peptides common to these molecules contained a high proportion of Leu at P4 and basic residues at P8 binding core positions. CONCLUSION: The degree of overlap between peptide repertoires associated with different HLA-DR molecules is low. The repertoires associated with DR1 and DR10 have the highest similarity among the molecules analyzed (∼10% overlap). Among the peptides shared between DR1 and DR10, a high proportion contained Leu(4) and basic residues at the P8 position of the binding core.

Effect of Oxidative Damage on the Stability and Dimerization of Superoxide Dismutase 1.

During their life cycle, proteins are subject to different modifications involving reactive oxygen species. Such oxidative damage to proteins may lead to the formation of insoluble aggregates and cytotoxicity and is associated with age-related disorders including neurodegenerative diseases, cancer, and diabetes. Superoxide dismutase 1 (SOD1), a key antioxidant enzyme in human cells, is particularly susceptible to such modifications. Moreover, this homodimeric metalloenzyme has been directly linked to both familial and sporadic amyotrophic lateral sclerosis (ALS), a devastating, late-onset motor neuronal disease, with more than 150 ALS-related mutations in the SOD1 gene. Importantly, oxidatively damaged SOD1 aggregates have been observed in both familial and sporadic forms of the disease. However, the molecular mechanisms as well as potential implications of oxidative stress in SOD1-induced cytotoxicity remain elusive. In this study, we examine the effects of oxidative modification on SOD1 monomer and homodimer stability, the key molecular properties related to SOD1 aggregation. We use molecular dynamics simulations in combination with thermodynamic integration to study microscopic-level site-specific effects of oxidative "mutations" at the dimer interface, including lysine, arginine, proline and threonine carbonylation, and cysteine oxidation. Our results show that oxidative damage of even single residues at the interface may drastically destabilize the SOD1 homodimer, with several modifications exhibiting a comparable effect to that of the most drastic ALS-causing mutations known. Additionally, we show that the SOD1 monomer stability decreases upon oxidative stress, which may lead to partial local unfolding and consequently to increased aggregation propensity. Importantly, these results suggest that oxidative stress may play a key role in development of ALS, with the mutations in the SOD1 gene being an additional factor.

Bottom-Up Instructive Quality Control in the Biofabrication of Smart Protein Materials.

The impact of cell factory quality control on material properties is a neglected but critical issue in the fabrication of protein biomaterials, which are unique in merging structure and function. The molecular chaperoning of protein conformational status is revealed here as a potent molecular instructor of the macroscopic properties of self-assembling, cell-targeted protein nanoparticles, including biodistribution upon in vivo administration.

From crystal structure to in silico epitope discovery in the Burkholderia pseudomallei flagellar hook-associated protein FlgK.

Melioidosis, caused by the Gram-negative bacterium Burkholderia pseudomallei, is a potentially fatal infection that is endemic in Southeast Asia and Northern Australia that is poorly controlled by antibiotics. Research efforts to identify antigenic components for a melioidosis vaccine have led to the identification of several proteins, including subunits forming the flagella that mediate bacterial motility, host colonization, and virulence. This study focuses on the B. pseudomallei flagellar hook-associated protein (FlgK(Bp)), and provides the first insights into the 3D structure of FlgK proteins as targets for structure-based antigen engineering. The FlgK(Bp) crystal structure (presented here at 1.8-Å resolution) reveals a multidomain fold, comprising two small ß-domains protruding from a large elongated α-helical bundle core. The evident structural similarity to flagellin, the flagellar filament subunit protein, suggests that, depending on the bacterial species, flagellar hook-associated proteins are likely to show a conserved, elongated α-helical bundle scaffold coupled to a variable number of smaller domains. Furthermore, we present immune serum recognition data confirming, in agreement with previous findings, that recovered melioidosis patients produce elevated levels of antibodies against FlgK(Bp), in comparison with seronegative and seropositive healthy subjects. Moreover, we show that FlgK(Bp) has cytotoxic effects on cultured murine macrophages, suggesting an important role in bacterial pathogenesis. Finally, computational epitope prediction methods applied to the FlgK(Bp) crystal structure, coupled with in vitro mapping, allowed us to predict three antigenic regions that locate to discrete protein domains. Taken together, our results point to FlgK(Bp) as a candidate for the design and production of epitope-containing subunits/domains as potential vaccine components.

In vivo architectonic stability of fully de novo designed protein-only nanoparticles.

The fully de novo design of protein building blocks for self-assembling as functional nanoparticles is a challenging task in emerging nanomedicines, which urgently demand novel, versatile, and biologically safe vehicles for imaging, drug delivery, and gene therapy. While the use of viruses and virus-like particles is limited by severe constraints, the generation of protein-only nanocarriers is progressively reachable by the engineering of protein-protein interactions, resulting in self-assembling functional building blocks. In particular, end-terminal cationic peptides drive the organization of structurally diverse protein species as regular nanosized oligomers, offering promise in the rational engineering of protein self-assembling. However, the in vivo stability of these constructs, being a critical issue for their medical applicability, needs to be assessed. We have explored here if the cross-molecular contacts between protein monomers, generated by end-terminal cationic peptides and oligohistidine tags, are stable enough for the resulting nanoparticles to overcome biological barriers in assembled form. The analyses of renal clearance and biodistribution of several tagged modular proteins reveal long-term architectonic stability, allowing systemic circulation and tissue targeting in form of nanoparticulate material. This observation fully supports the value of the engineered of protein building blocks addressed to the biofabrication of smart, robust, and multifunctional nanoparticles with medical applicability that mimic structure and functional capabilities of viral capsids.

Sheltering DNA in self-organizing, protein-only nano-shells as artificial viruses for gene delivery.

By recruiting functional domains supporting DNA condensation, cell binding, internalization, endosomal escape and nuclear transport, modular single-chain polypeptides can be tailored to associate with cargo DNA for cell-targeted gene therapy. Recently, an emerging architectonic principle at the nanoscale has permitted tagging protein monomers for self-organization as protein-only nanoparticles. We have studied here the accommodation of plasmid DNA into protein nanoparticles assembled with the synergistic assistance of end terminal poly-arginines (R9) and poly-histidines (H6). Data indicate a virus-like organization of the complexes, in which a DNA core is surrounded by a solvent-exposed protein layer. This finding validates end-terminal cationic peptides as pleiotropic tags in protein building blocks for the mimicry of viral architecture in artificial viruses, representing a promising alternative to the conventional use of viruses and virus-like particles for nanomedicine and gene therapy. FROM THE CLINICAL EDITOR: Finding efficient gene delivery methods still represents a challenge and is one of the bottlenecks to the more widespread application of gene therapy. The findings presented in this paper validate the application of end-terminal cationic peptides as pleiotropic tags in protein building blocks for "viral architecture mimicking" in artificial viruses, representing a promising alternative to the use of viruses and virus-like particles for gene delivery.

RGD-based cell ligands for cell-targeted drug delivery act as potent trophic factors.

Integrin-binding, Arg-Gly-Asp (RGD)-containing peptides are the most widely used agents to deliver drugs, nanoparticles, and imaging agents. Although in nature, several protein-mediated signal transduction events depend on RGD motifs, the potential of RGD-empowered materials in triggering undesired cell-signaling cascades has been neglected. Using an RGD-functionalized protein nanoparticle, we show here that the RGD motif acts as a powerful trophic factor, supporting extracellular signal-regulated kinase 1/2 (ERK1/2)-linked cell proliferation and partial differentiation of PC12 cells, a neuronlike model. FROM THE CLINICAL EDITOR: This work focuses on RGD peptides, which are among the most commonly used tags for targeted drug delivery. They also promote protoneurite formation and expression of neuronal markers (MAP2) in model PC12 cells, which is an unexpected but relevant event in the functionalization of drugs and their nanocarriers.

Disulfide bond formation and activation of Escherichia coli ß-galactosidase under oxidizing conditions.

Escherichia coli ß-galactosidase is probably the most widely used reporter enzyme in molecular biology, cell biology, and biotechnology because of the easy detection of its activity. Its large size and tetrameric structure make this bacterial protein an interesting model for crystallographic studies and atomic mapping. In the present study, we investigate a version of Escherichia coli ß-galactosidase produced under oxidizing conditions, in the cytoplasm of an Origami strain. Our data prove the activation of this microbial enzyme under oxidizing conditions and clearly show the occurrence of a disulfide bond in the ß-galactosidase structure. Additionally, the formation of this disulfide bond is supported by the analysis of a homology model of the protein that indicates that two cysteines located in the vicinity of the catalytic center are sufficiently close for disulfide bond formation.

Nanoparticulate architecture of protein-based artificial viruses is supported by protein-DNA interactions.

UNLABELLED: AIM & METHODS: We have produced two chimerical peptides of 10.2 kDa, each contain four biologically active domains, which act as building blocks of protein-based nonviral vehicles for gene therapy. In solution, these peptides tend to aggregate as amorphous clusters of more than 1000 nm, while the presence of DNA promotes their architectonic reorganization as mechanically stable nanometric spherical entities of approximately 80 nm that penetrate mammalian cells through arginine-glycine-aspartic acid cell-binding domains and promote significant transgene expression levels. RESULTS & CONCLUSION: The structural analysis of the protein in these hybrid nanoparticles indicates a molecular conformation with predominance of α-helix and the absence of cross-molecular, ß-sheet-supported protein interactions. The nanoscale organizing forces generated by DNA-protein interactions can then be observed as a potentially tunable, critical factor in the design of protein-only based artificial viruses for gene therapy.

Structural and biochemical characterization of NarE, an iron-containing ADP-ribosyltransferase from Neisseria meningitidis.

NarE is a 16 kDa protein identified from Neisseria meningitidis, one of the bacterial pathogens responsible for meningitis. NarE belongs to the family of ADP-ribosyltransferases (ADPRT) and catalyzes the transfer of ADP-ribose moieties to arginine residues in target protein acceptors. Many pathogenic bacteria utilize ADP-ribosylating toxins to modify and alter essential functions of eukaryotic cells. NarE is further the first ADPRT which could be shown to bind iron through a Fe-S center, which is crucial for the catalytic activity. Here we present the NMR solution structure of NarE, which shows structural homology to other ADPRTs. Using NMR titration experiments we could identify from Chemical Shift Perturbation data both the NAD binding site, which is in perfect agreement with a consensus sequence analysis between different ADPRTs, as well as the iron coordination site, which consists of 2 cysteines and 2 histidines. This atypical iron coordination is also capable to bind zinc. These results could be fortified by site-directed mutagenesis of the catalytic region, which identified two functionally crucial residues. We could further identify a main interaction region of NarE with antibodies using two complementary methods based on antibody immobilization, proteolytic digestion, and mass spectrometry. This study combines structural and functional features of NarE providing for the first time a characterization of an iron-dependent ADPRT.

Protein nanodisk assembling and intracellular trafficking powered by an arginine-rich (R9) peptide.

AIMS: Arginine(R)-rich cationic peptides are powerful tools in drug delivery since, alone or when associated with polyplexes, proteins or chemicals, they confer DNA condensation, membrane translocation and blood-brain barrier crossing abilities. The unusual stability and high in vivo performance of their associated drugs suggest a particulate organization or R(n) complexes, which this study aimed to explore. MATERIALS & METHODS: We have analyzed the particulate organization and biological performance in DNA delivery of a model, R9-containing green fluorescent protein by dynamic light scattering, transmission electron microscopy, atomic force microscopy, single cell confocal microscopy and flow cytometry. RESULTS: A deep nanoscale examination of R9-powered constructs reveals a novel and promising feature of R9, that when fused to a scaffold green fluorescent protein, promote its efficient self-assembling as highly stable, regular disk-shaped nanoparticles of 20 x 3 nm. These constructs are efficiently internalized in mammalian cells and rapidly migrate through the cytoplasm towards the nucleus in a fully bioactive form. Besides, such particulate platforms accommodate, condense and deliver plasmid DNA to the nucleus and promote plasmid-driven transgene expression. CONCLUSION: The architectonic properties of arginine-rich peptides at the nanoscale reveal a new category of protein nanoparticles, namely nanodisks, and provide novel strategic concepts and architectonic tools for the tailored construction of new-generation artificial viruses for gene therapy and drug delivery.

Polypeptide folding on a conformational-space network: dependence of network topology on the structural discretization procedure.

Mapping the conformational space of a polypeptide onto a network of conformational states involves a number of subjective choices, mostly in relation to the definition of conformation and its discrete nature in a network framework. Here, we evaluate the robustness of the topology of conformational-space networks derived from Molecular Dynamics (MD) simulations with respect to the use of different discretization (clustering) methods, variation of their parameters, simulation length and analysis time-step, and removing high-frequency motions from the coordinate trajectories. In addition, we investigate the extent to which polypeptide dynamics can be reproduced on the resulting networks when assuming Markovian behavior. The analysis is based on eight 500 ns and eight 400 ns MD simulations in explicit water of two 10-residue peptides. Three clustering algorithms were used, two of them based on the pair-wise root-mean-square difference between structures and one on dihedral-angle patterns. A short characteristic path length and a power-law behavior of the probability distribution of the node degree are obtained irrespective of the clustering method or the value of any of the tested parameters. The average cliquishness is consistently one or two orders of magnitude larger than that of a random realization of a network of corresponding size and connectivity. The cliquishness as function of node degree and the kinetic properties of the networks are found to be most dependent on clustering method and/or parameters. Although Markovian simulations on the networks reproduce cluster populations accurately, their kinetic properties most often differ from those observed in the MD simulations.

Molecular dynamics simulations of a reversibly folding beta-heptapeptide in methanol: influence of the treatment of long-range electrostatic interactions.

Eight 100-ns molecular dynamics simulations of a beta-heptapeptide in methanol at 340 K (within cubic periodic computational boxes of about 6-nm edge) are reported and compared. These simulations were performed with three different charge-state combinations at the peptide termini, one of them with or without a neutralizing chloride counterion, and using either the lattice-sum (LS) or reaction-field (RF) scheme to handle electrostatic interactions. The choice of the electrostatic scheme has essentially no influence on the folding-unfolding equilibrium when the peptide termini are uncharged and only a small influence when the peptide is positively charged at its N-terminus (with or without inclusion of a neutralizing chloride counterion). However, when the peptide is zwitterionic, the LS scheme leads to preferential sampling of the high-dipole folded helical state, whereas the RF scheme leads to preferential sampling of a low-dipole unfolded salt-bridged state. A continuum electrostatics analysis based on the sampled configurations (zwitterionic case) suggests that the LS scheme stabilizes the helical state through artificial periodicity, but that the magnitude of this perturbation is essentially negligible (compared to the thermal energy) for the large box size and relatively polar solvent considered. The results thus provide clear evidence (continuum electrostatics analysis) for the absence of LS artifacts and some indications (still not definitive because of the limited sampling of the folding-unfolding transition) for the presence of RF artifacts in this specific system.