Evolutionary deletions within the SARS-CoV-2 genome as signature trends for virus fitness and adaptation

Coronaviruses are large RNA viruses that can infect and spread among humans and animals. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), responsible for coronavirus disease 2019, has evolved since its first detection in December 2019. Deletions are a common occurrence in SARS-CoV-2 evolution, particularly in specific genomic sites, and may be associated with the emergence of highly competent lineages. While deletions typically have a negative impact on viral fitness, some persist and become fixed in viral populations, indicating that they may confer advantageous benefits for the virus's adaptive evolution. This work presents a literature review and data analysis on structural losses in the SARS-CoV-2 genome and the potential relevance of specific signatures for enhanced viral fitness and spread.

Seroprevalence of Dengue, Chikungunya and Zika at the epicenter of the congenital microcephaly epidemic in Northeast Brazil: A population-based survey

BACKGROUND

The four Dengue viruses (DENV) serotypes were re-introduced in Brazil's Northeast region in a couple of decades, between 1980's and 2010's, where the DENV1 was the first detected serotype and DENV4 the latest. Zika (ZIKV) and Chikungunya (CHIKV) viruses were introduced in Recife around 2014 and led to large outbreaks in 2015 and 2016, respectively. However, the true extent of the ZIKV and CHIKV outbreaks, as well as the risk factors associated with exposure to these viruses remain vague.

METHODS

We conducted a stratified multistage household serosurvey among residents aged between 5 and 65 years in the city of Recife, Northeast Brazil, from August 2018 to February 2019. The city neighborhoods were stratified and divided into high, intermediate, and low socioeconomic strata (SES). Previous ZIKV, DENV and CHIKV infections were detected by IgG-based enzyme linked immunosorbent assays (ELISA). Recent ZIKV and CHIKV infections were assessed through IgG3 and IgM ELISA, respectively. Design-adjusted seroprevalence were estimated by age group, sex, and SES. The ZIKV seroprevalence was adjusted to account for the cross-reactivity with dengue. Individual and household-related risk factors were analyzed through regression models to calculate the force of infection. Odds Ratio (OR) were estimated as measure of effect.

PRINCIPAL FINDINGS

A total of 2,070 residents' samples were collected and analyzed. The force of viral infection for high SES were lower as compared to low and intermediate SES. DENV seroprevalence was 88.7% (CI95%:87.0-90.4), and ranged from 81.2% (CI95%:76.9-85.6) in the high SES to 90.7% (CI95%:88.3-93.2) in the low SES. The overall adjusted ZIKV seroprevalence was 34.6% (CI95%:20.0-50.9), and ranged from 47.4% (CI95%:31.8-61.5) in the low SES to 23.4% (CI95%:12.2-33.8) in the high SES. The overall CHIKV seroprevalence was 35.7% (CI95%:32.6-38.9), and ranged from 38.6% (CI95%:33.6-43.6) in the low SES to 22.3% (CI95%:15.8-28.8) in the high SES. Surprisingly, ZIKV seroprevalence rapidly increased with age in the low and intermediate SES, while exhibited only a small increase with age in high SES. CHIKV seroprevalence according to age was stable in all SES. The prevalence of serological markers of ZIKV and CHIKV recent infections were 1.5% (CI95%:0.1-3.7) and 3.5% (CI95%:2.7-4.2), respectively.

CONCLUSIONS

Our results confirmed continued DENV transmission and intense ZIKV and CHIKV transmission during the 2015/2016 epidemics followed by ongoing low-level transmission. The study also highlights that a significant proportion of the population is still susceptible to be infected by ZIKV and CHIKV. The reasons underlying a ceasing of the ZIKV epidemic in 2017/18 and the impact of antibody decay in susceptibility to future DENV and ZIKV infections may be related to the interplay between disease transmission mechanism and actual exposure in the different SES.

Inhibition of 3-Hydroxykynurenine Transaminase from Aedes aegypti and Anopheles gambiae: A Mosquito-Specific Target to Combat the Transmission of Arboviruses

Arboviral infections such as Zika, chikungunya, dengue, and yellow fever pose significant health problems globally. The population at risk is expanding with the geographical distribution of the main transmission vector of these viruses, the Aedes aegypti mosquito. The global spreading of this mosquito is driven by human migration, urbanization, climate change, and the ecological plasticity of the species. Currently, there are no specific treatments for Aedes-borne infections. One strategy to combat different mosquito-borne arboviruses is to design molecules that can specifically inhibit a critical host protein. We obtained the crystal structure of 3-hydroxykynurenine transaminase (AeHKT) from A. aegypti, an essential detoxification enzyme of the tryptophan metabolism pathway. Since AeHKT is found exclusively in mosquitoes, it provides the ideal molecular target for the development of inhibitors. Therefore, we determined and compared the free binding energy of the inhibitors 4-(2-aminophenyl)-4-oxobutyric acid (4OB) and sodium 4-(3-phenyl-1,2,4-oxadiazol-5-yl)butanoate (OXA) to AeHKT and AgHKT from Anopheles gambiae, the only crystal structure of this enzyme previously known. The cocrystallized inhibitor 4OB binds to AgHKT with K _i of 300 μM. We showed that OXA binds to both AeHKT and AgHKT enzymes with binding energies 2-fold more favorable than the crystallographic inhibitor 4OB and displayed a 2-fold greater residence time τ upon binding to AeHKT than 4OB. These findings indicate that the 1,2,4-oxadiazole derivatives are inhibitors of the HKT enzyme not only from A. aegypti but also from A. gambiae.

An artificial neural network model to predict structure-based protein-protein free energy of binding from Rosetta-calculated properties

The prediction of the free energy (ΔG) of binding for protein-protein complexes is of general scientific interest as it has a variety of applications in the fields of molecular and chemical biology, materials science, and biotechnology. Despite its centrality in understanding protein association phenomena and protein engineering, the ΔG of binding is a daunting quantity to obtain theoretically. In this work, we devise a novel Artificial Neural Network (ANN) model to predict the ΔG of binding for a given three-dimensional structure of a protein-protein complex with Rosetta-calculated properties. Our model was tested using two data sets, and it presented a root-mean-square error ranging from 1.67 kcal mol^-1 to 2.45 kcal mol^-1, showing a better performance compared to the available state-of-the-art tools. Validation of the model for a variety of protein-protein complexes is showcased.

GM2/GM3 controls the organizational status of CD82/Met microdomains: further studies in GM2/GM3 complexation

At cell surface gangliosides might associate with signal transducers proteins, grown factor receptors, integrins, small G-proteins and tetraspanins establishing microdomains, which play important role in cell adhesion, cell activation, motility, and growth. Previously, we reported that GM2 and GM3 form a heterodimer that interacts with the tetraspanin CD82, controlling epithelial cell mobility by inhibiting integrin-hepatocyte growth factor-induced cMet tyrosine kinase signaling. By using molecular dynamics simulations to study the molecular basis of GM2/GM3 interaction we demonstrate, here, that intracellular levels of Ca²⁺ mediate GM2/GM3 complexation via electrostatic interaction with their carboxyl groups, while hydrogen bonds between the ceramide groups likely aid stabilizing the complex. The presence of GM2/GM3 complex alters localization of CD82 on cell surface and therefore downstream signalization. These data contribute for the knowledge of how glycosylation may control signal transduction and phenotypic changes.

Identification of a Zika NS2B epitope as a biomarker for severe clinical phenotypes

The identification of specific biomarkers for Zika infection and its clinical complications is fundamental to mitigate the infection spread, which has been associated with a broad range of neurological sequelae. We present the characterization of antibody responses in serum samples from individuals infected with Zika, presenting non-severe (classical) and severe (neurological disease) phenotypes, with high-density peptide arrays comprising the Zika NS1 and NS2B proteins. The data pinpoints one strongly IgG-targeted NS2B epitope in non-severe infections, which is absent in Zika patients, where infection progressed to the severe phenotype. This differential IgG profile between the studied groups was confirmed by multivariate data analysis. Molecular dynamics simulations and circular dichroism have shown that the peptide in solution presents itself in a sub-optimal conformation for antibody recognition, which led us to computationally engineer an artificial protein able to stabilize the NS2B epitope structure. The engineered protein was used to interrogate paired samples from mothers and their babies presenting Zika-associated microcephaly and confirmed the absence of NS2B IgG response in those samples. These findings suggest that the assessment of antibody responses to the herein identified NS2B epitope is a strong candidate biomarker for the diagnosis and prognosis of Zika-associated neurological disease.

The ongoing evolution of variants of concern and interest of SARS-CoV-2 in Brazil revealed by convergent indels in the amino (N)-terminal domain of the spike protein

Mutations at both the receptor-binding domain (RBD) and the amino (N)-terminal domain (NTD) of the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Spike (S) glycoprotein can alter its antigenicity and promote immune escape. We identified that SARS-CoV-2 lineages circulating in Brazil with mutations of concern in the RBD independently acquired convergent deletions and insertions in the NTD of the S protein, which altered the NTD antigenic-supersite and other predicted epitopes at this region. Importantly, we detected the community transmission of different P.1 lineages bearing NTD indels ∆69-70 (which can impact several SARS-CoV-2 diagnostic protocols), ∆144 and ins214ANRN, and a new VOI N.10 derived from the B.1.1.33 lineage carrying three NTD deletions (∆141-144, ∆211, and ∆256-258). These findings support that the ongoing widespread transmission of SARS-CoV-2 in Brazil generates new viral lineages that might be more resistant to antibody neutralization than parental variants of concern.

Differential Frequencies of HLA-DRB1, DQA1, and DQB1 Alleles and Haplotypes Are Observed in the Arbovirus-Related Neurological Syndromes

BACKGROUND

We took advantage of the 2015-2016 Brazilian arbovirus outbreak (Zika [ZIKV]/dengue/chikungunya viruses) associated with neurological complications to type HLA-DRB1/DQA1/DQB1 variants in patients exhibiting neurological complications and in bone marrow donors from the same endemic geographical region.

METHODS

DRB1/DQA1/DQB1 loci were typed using sequence-specific oligonucleotides. In silico studies were performed using X-ray resolved dimer constructions.

RESULTS

The DQA1*01, DQA1*05, DQB1*02, or DQB1*06 genotypes/haplotypes and DQA1/DQB1 haplotypes that encode the putative DQA1/DQB1 dimers were overrepresented in the whole group of patients and in patients exhibiting peripheral neurological spectrum disorders (PSD) or encephalitis spectrum disorders (ESD). The DRB1*04, DRB1*13, and DQA1*03 allele groups protected against arbovirus neurological manifestation, being underrepresented in whole group of patients and ESD and PSD groups. Genetic and in silico studies revealed that DQA1/DQB1 dimers (1) were primarily associated with susceptibility to arbovirus infections; (2) can bind to a broad range of ZIKV peptides (235 of 1878 peptides, primarily prM and NS2A); and (3) exhibited hydrophilic and highly positively charged grooves when compared to the DRA1/DRB1 cleft. The protective dimer (DRA1/DRB1*04) bound a limited number of ZIKV peptides (40 of 1878 peptides, primarily prM).

CONCLUSION

Protective haplotypes may recognize arbovirus peptides more specifically than susceptible haplotypes.

Immune evasion of SARS-CoV-2 variants of concern is driven by low affinity to neutralizing antibodies

SARS-CoV-2 VOC immune evasion is mainly due to lower cross-reactivity from previously elicited class I/II neutralizing antibodies, while increased affinity to hACE2 plays a minor role. The affinity between antibodies and VOCs is impacted by remodeling of the electrostatic surface potential of the Spike RBDs. The P.3 variant is a putative VOC.

Repositioning Lopinavir, an HIV Protease Inhibitor, as a Promising Antifungal Drug: Lessons Learned from Candida albicans-In Silico, In Vitro and In Vivo Approaches

The repurposing strategy was applied herein to evaluate the effects of lopinavir, an aspartic protease inhibitor currently used in the treatment of HIV-infected individuals, on the globally widespread opportunistic human fungal pathogen Candida albicans by using in silico, in vitro and in vivo approaches in order to decipher its targets on fungal cells and its antifungal mechanisms of action. Secreted aspartic proteases (Saps) are the obviously main target of lopinavir. To confirm this hypothesis, molecular docking assays revealed that lopinavir bound to the Sap2 catalytic site of C. albicans as well as inhibited the Sap hydrolytic activity in a typically dose-dependent manner. The inhibition of Saps culminated in the inability of C. albicans yeasts to assimilate the unique nitrogen source (albumin) available in the culture medium, culminating with fungal growth inhibition (IC₅₀ = 39.8 µM). The antifungal action of lopinavir was corroborated by distinct microscopy analyses, which evidenced drastic and irreversible changes in the morphology that justified the fungal death. Furthermore, our results revealed that lopinavir was able to (i) arrest the yeasts-into-hyphae transformation, (ii) disturb the synthesis of neutral lipids, including ergosterol, (iii) modulate the surface-located molecules, such as Saps and mannose-, sialic acid- and N-acetylglucosamine-containing glycoconjugates, (iv) diminish the secretion of hydrolytic enzymes, such as Saps and esterase, (v) negatively influence the biofilm formation on polystyrene surface, (vi) block the in vitro adhesion to epithelial cells, (vii) contain the in vivo infection in both immunocompetent and immunosuppressed mice and (viii) reduce the Sap production by yeasts recovered from kidneys of infected animals. Conclusively, the exposed results highlight that lopinavir may be used as a promising repurposing drug against C. albicans infection as well as may be used as a lead compound for the development of novel antifungal drugs.

Rotational Profiler: A Fast, Automated, and Interactive Server to Derive Torsional Dihedral Potentials for Classical Molecular Simulations

Rotational Profiler provides an analytical algorithm to compute sets of classical torsional dihedral parameters by fitting an empirical energy profile to a reference one that can be obtained experimentally or by quantum-mechanical methods. The resulting profiles are compatible with the functional forms in the most widely used biomolecular force fields (e.g., GROMOS, AMBER, OPLS, and CHARMM). The linear least-squares regression method is used to generate sets of parameters that best satisfy the fitting. Rotational Profiler is free to use, analytical, and force field/package independent. The formalism is herein described, and its usage, in an interactive and automated manner, is made available as a Web server at http://rotprof.lncc.br.

NWChem: Past, present, and future

Specialized computational chemistry packages have permanently reshaped the landscape of chemical and materials science by providing tools to support and guide experimental efforts and for the prediction of atomistic and electronic properties. In this regard, electronic structure packages have played a special role by using first-principle-driven methodologies to model complex chemical and materials processes. Over the past few decades, the rapid development of computing technologies and the tremendous increase in computational power have offered a unique chance to study complex transformations using sophisticated and predictive many-body techniques that describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, novel parallel algorithms have been able to take advantage of computational resources to address the polynomial scaling of electronic structure methods. In this paper, we briefly review the NWChem computational chemistry suite, including its history, design principles, parallel tools, current capabilities, outreach, and outlook.

Influence of directional positive Darwinian selection-driven evolution on arboviruses Dengue and Zika virulence and pathogenesis

Dengue (DENV) and Zika (ZIKV) viruses are antigenically and evolutionarily related; immunological cross-reactions between them have been associated to both cross-protection and infection-enhanced mechanisms. Here, DENV-1-4 and ZIKV were investigated through Bayesian coalescent-based approaches and selection-driven Darwinian evolution methods using robust datasets. Our findings show that both DENV and ZIKV, driven essentially by directional positive selection, have undergone evolution and diversification and that their entire polyproteins are subject to an intense directional evolution. Interestingly, positively selected codons mapped here are directly associated to DENV-1-2 virulence as well as the ZIKV burgeoning 2015-16 outbreak in the Americas, therefore, having impact on the pathogenesis of these viruses. Biochemical prediction analysis focusing on markers involved in virulence and viral transmission dynamics identified alterations in N-Glycosylation-, Phosphorylation- and Palmitoylation-sites in ZIKV sampled from different countries, hosts and isolation sources. Taking into account both DENV-ZIKV co-circulation either into and/or out of flavivirus-endemic regions, as well as recombination and quasispecies scenarios, these results indicate the action of a selection-driven evolution affecting the biology, virulence and pathogenesis of these pathogens in a non-randomized environment.

The influence of biotinylation on the ability of a computer designed protein to detect B-cells producing anti-HIV-1 2F5 antibodies

Antibodies against the HIV-1 2F5 epitope are known as one of the most powerful and broadly protective anti-HIV antibodies. Therefore, vaccine strategies that include the 2F5 epitope in their formulation require a robust method to detect specific anti-2F5 antibody production by B cells. Towards this goal, we have biotinylated a previously reported computer-designed protein carrying the HIV-1 2F5 epitope aiming the further development of a platform to detect human B-cells expressing anti-2F5 antibodies through flow cytometry. Biophysical and immunological properties of our devised protein were characterized by computer simulation and experimental methods. Biotinylation did not affect folding and improved protein stability and solubility. The biotinylated protein exhibited similar binding affinity trends compared to its unbiotinylated counterpart and was recognized by anti-HIV-1 2F5 antibodies expressed on the surface of patient-derived peripheral blood mononuclear cells. Moreover, we present a high affinity marker for the identification of epitope-specific B cells that can be used to measure the efficacy of vaccine strategies based on the HIV-1 envelope protein.

Aromatic Rings Commonly Used in Medicinal Chemistry: Force Fields Comparison and Interactions With Water Toward the Design of New Chemical Entities

The identification of lead compounds usually includes a step of chemical diversity generation. Its rationale may be supported by both qualitative (SAR) and quantitative (QSAR) approaches, offering models of the putative ligand-receptor interactions. In both scenarios, our understanding of which interactions functional groups can perform is mostly based on their chemical nature (such as electronegativity, volume, melting point, lipophilicity etc.) instead of their dynamics in aqueous, biological solutions (solvent accessibility, lifetime of hydrogen bonds, solvent structure etc.). As a consequence, it is challenging to predict from 2D structures which functional groups will be able to perform interactions with the target receptor, at which intensity and relative abundance in the biological environment, all of which will contribute to ligand potency and intrinsic activity. With this in mind, the aim of this work is to assess properties of aromatic rings, commonly used for drug design, in aqueous solution through molecular dynamics simulations in order to characterize their chemical features and infer their impact in complexation dynamics. For this, common aromatic and heteroaromatic rings were selected and received new atomic charge set based on the direction and module of the dipole moment from MP2/6-31G^* calculations, while other topological terms were taken from GROMOS53A6 force field. Afterwards, liquid physicochemical properties were simulated for a calibration set composed by nearly 40 molecules and compared to their respective experimental data, in order to validate each topology. Based on the reliance of the employed strategy, we expanded the dataset to more than 100 aromatic rings. Properties in aqueous solution such as solvent accessible surface area, H-bonds availability, H-bonds residence time, and water structure around heteroatoms were calculated for each ring, creating a database of potential interactions, shedding light on features of drugs in biological solutions, on the structural basis for bioisosterism and on the enthalpic/entropic costs for ligand-receptor complexation dynamics.

Significant Performance Enhancement of Polymer Resins by Bioinspired Dynamic Bonding

Marine mussels use catechol-rich interfacial mussel foot proteins (mfps) as primers that attach to mineral surfaces via hydrogen, metal coordination, electrostatic, ionic, or hydrophobic bonds, creating a secondary surface that promotes bonding to the bulk mfps. Inspired by this biological adhesive primer, it is shown that a ≈1 nm thick catecholic single-molecule priming layer increases the adhesion strength of crosslinked polymethacrylate resin on mineral surfaces by up to an order of magnitude when compared with conventional primers such as noncatecholic silane- and phosphate-based grafts. Molecular dynamics simulations confirm that catechol groups anchor to a variety of mineral surfaces and shed light on the binding mode of each molecule. Here, a ≈50% toughness enhancement is achieved in a stiff load-bearing polymer network, demonstrating the utility of mussel-inspired bonding for processing a wide range of polymeric interfaces, including structural, load-bearing materials.

Mapping Putative B-Cell Zika Virus NS1 Epitopes Provides Molecular Basis for Anti-NS1 Antibody Discrimination between Zika and Dengue Viruses

B-cell epitope sequences from Zika virus (ZIKV) NS1 protein have been identified using epitope prediction tools. Mapping these sequences onto the NS1 surface reveals two major conformational epitopes and a single linear one. Despite an overall average sequence identity of ca. 55% between the NS1 from ZIKV and the four dengue virus (DENV) serotypes, epitope sequences were found to be highly conserved. Nevertheless, nonconserved epitope-flanking residues are responsible for a dramatically divergent electrostatic surface potential on the epitope regions of ZIKV and DENV2 serotypes. These findings suggest that strategies for differential diagnostics on the basis of short linear NS1 sequences are likely to fail due to immunological cross-reactions. Overall, results provide the molecular basis of differential discrimination between Zika and DENVs by NS1 monoclonal antibodies.

Amino acid polymorphisms in the fibronectin-binding repeats of fibronectin-binding protein A affect bond strength and fibronectin conformation

The Staphylococcus aureus cell surface contains cell wall-anchored proteins such as fibronectin-binding protein A (FnBPA) that bind to host ligands (e.g. fibronectin; Fn) present in the extracellular matrix of tissue or coatings on cardiac implants. Recent clinical studies have found a correlation between cardiovascular infections caused by S. aureus and nonsynonymous SNPs in FnBPA. Atomic force microscopy (AFM), surface plasmon resonance (SPR), and molecular simulations were used to investigate interactions between Fn and each of eight 20-mer peptide variants containing amino acids Ala, Asn, Gln, His, Ile, and Lys at positions equivalent to 782 and/or 786 in Fn-binding repeat-9 of FnBPA. Experimentally measured bond lifetimes (1/k_off) and dissociation constants (K_d = k_off/k_on), determined by mechanically dissociating the Fn·peptide complex at loading rates relevant to the cardiovascular system, varied from the lowest-affinity H782A/K786A peptide (0.011 s, 747 μm) to the highest-affinity H782Q/K786N peptide (0.192 s, 15.7 μm). These atomic force microscopy results tracked remarkably well to metadynamics simulations in which peptide detachment was defined solely by the free-energy landscape. Simulations and SPR experiments suggested that an Fn conformational change may enhance the stability of the binding complex for peptides with K786I or H782Q/K786I (K_d^app = 0.2-0.5 μm, as determined by SPR) compared with the lowest-affinity double-alanine peptide (K_d^app = 3.8 μm). Together, these findings demonstrate that amino acid substitutions in Fn-binding repeat-9 can significantly affect bond strength and influence the conformation of Fn upon binding. They provide a mechanistic explanation for the observation of nonsynonymous SNPs in fnbA among clinical isolates of S. aureus that cause endovascular infections.

Conformational stability of the epidermal growth factor (EGF) receptor as influenced by glycosylation, dimerization and EGF hormone binding

The epidermal growth factor receptor (EGFR) is an important transmembrane glycoprotein kinase involved the initiation or perpetuation of signal transduction cascades within cells. These processes occur after EGFR binds to a ligand [epidermal growth factor (EGF)], thus inducing its dimerization and tyrosine autophosphorylation. Previous publications have highlighted the importance of glycosylation and dimerization for promoting proper function of the receptor and conformation in membranes; however, the effects of these associations on the protein conformational stability have not yet been described. Molecular dynamics simulations were performed to characterize the conformational preferences of the monomeric and dimeric forms of the EGFR extracellular domain upon binding to EGF in the presence and absence of N-glycan moieties. Structural stability analyses revealed that EGF provides the most conformational stability to EGFR, followed by glycosylation and dimerization, respectively. The findings also support that EGF-EGFR binding takes place through a large-scale induced-fitting mechanism. Proteins 2017; 85:561-570. © 2016 Wiley Periodicals, Inc.

Molecularly Smooth Self-Assembled Monolayer for High-Mobility Organic Field-Effect Transistors

Despite the need for molecularly smooth self-assembled monolayers (SAMs) on silicon dioxide surfaces (the most common dielectric surface), current techniques are limited to nonideal silane grafting. Here, we show unique bioinspired zwitterionic molecules forming a molecularly smooth and uniformly thin SAM in "water" in <1 min on various dielectric surfaces, which enables a dip-coating process that is essential for organic electronics to become reality. This monomolecular layer leads to high mobility of organic field-effect transistors (OFETs) based on various organic semiconductors and source/drain electrodes. A combination of experimental and computational techniques confirms strong adsorption (W_ad > 20 mJ m^-2), uniform thickness (∼0.5 or ∼1 nm) and orientation (all catechol head groups facing the oxide surface) of the "monomolecular" layers. This robust (strong adsorption), rapid, and green SAM represents a promising advancement toward the next generation of nanofabrication compared to the current nonuniform and inconsistent polysiloxane-based SAM involving toxic chemicals, long processing time (>10 h), or heat (>80 °C).

Assessing protein conformational sampling and structural stability via de novo design and molecular dynamics simulations

Molecular dynamics and de novo techniques, associated to quality parameter sets, have excelled at determining the structure of small proteins with high accuracy. To achieve a detailed description of protein conformations, these methods must critically assess the thermodynamic features of the molecular ensembles. Here, a comparison of the conformational ensemble generated by molecular dynamics and de novo techniques were carried out for six Top7-based proteins carrying gp41 HIV-1 epitopes. The native Top7, a highly stable computationally designed protein, was used as benchmark. Structural stability, flexibility, and secondary structure content were assessed. The consistency of the latter was compared to experimental circular dichroism spectra for all proteins. While both methods are capable to identify the stable from unstable chimeric proteins, the sampled conformational space and flexibility differ significantly in both methods. Molecular dynamics simulations seem to better describe secondary structure content and identify regions responsible for conformational instability. The de novo method, as implemented in Rosetta-a prime tool for protein design, overestimates secondary structure content. On the other hand, its empirical energy function is capable to predict the threshold for protein stability.

Extension and validation of the GROMOS 53A6(GLYC) parameter set for glycoproteins

An extension of the GROMOS 53A6GLYC force field for carbohydrates to encompass glycoprotein linkages is presented. The set includes new atomic charges and incorporates adequate torsional potential parameters for N-, S-, C-, P-, and O-glycosydic linkages, offering compatibility with the GROMOS force field family for proteins. Validation included the description of glycosydic linkage geometries between amino acid and monosaccharide residues, comparison of NMR-derived protein-carbohydrate and carbohydrate-carbohydrate nuclear overhauser effect (NOE) signals for glycoproteins and the effects of glycosylation on protein flexibility and dynamics.

Desorption/ionization efficiencies of triacylglycerols and phospholipids via EASI-MS

Knowledge of the major effects governing desorption/ionization efficiency is required for the development and application of ambient mass spectrometry. Although all triacylglycerols (TAG) have the same favorable protonation and cationization sites, their desorption/ionization efficiencies can vary dramatically during easy ambient sonic-spray ionization because of structural differences in the carbon chain. To quantify this somewhat surprising and drastic effect, we have performed a systematic investigation of desorption/ionization efficiencies as a function of unsaturation and length for TAG as well as for diacylglycerols, monoacylglycerols and several phospholipids (PL). Affinities for Na(+) as a function of unsaturation level have also been assayed via comprehensive metadynamics calculations to understand the influence of this phenomenon on the ionization efficiency. The results suggest that dipole-dipole interactions within a carbon chain tuned by unsaturation sites govern ionization efficiency of TAG and PL.

MDWiZ: a platform for the automated translation of molecular dynamics simulations

A variety of popular molecular dynamics (MD) simulation packages were independently developed in the last decades to reach diverse scientific goals. However, such non-coordinated development of software, force fields, and analysis tools for molecular simulations gave rise to an array of software formats and arbitrary conventions for routine preparation and analysis of simulation input and output data. Different formats and/or parameter definitions are used at each stage of the modeling process despite largely contain redundant information between alternative software tools. Such Babel of languages that cannot be easily and univocally translated one into another poses one of the major technical obstacles to the preparation, translation, and comparison of molecular simulation data that users face on a daily basis. Here, we present the MDWiZ platform, a freely accessed online portal designed to aid the fast and reliable preparation and conversion of file formats that allows researchers to reproduce or generate data from MD simulations using different setups, including force fields and models with different underlying potential forms. The general structure of MDWiZ is presented, the features of version 1.0 are detailed, and an extensive validation based on GROMACS to LAMMPS conversion is presented. We believe that MDWiZ will be largely useful to the molecular dynamics community. Such fast format and force field exchange for a given system allows tailoring the chosen system to a given computer platform and/or taking advantage of a specific capabilities offered by different software engines.

IMSPeptider: A computational peptide collision cross-section area calculator based on a novel molecular dynamics simulation protocol

Introduction of ion mobility mass spectrometry (IMS/MS) into the proteomic workflow provides an orthogonal separation to the widely used LC-MS platforms. IMS also provides structural information that could facilitate peptide identification. However, the lack of tools capable of predictive power in a high-throughput fashion makes peptide global profiling quite challenging. To target this issue, a computational workflow was developed based on biophysical principles to predict the collision cross-section area (CCS) of peptides as measured from IMS/MS experiments. Hosted on a web server, it allows the user to input a primary sequence (query) and retrieve information on peptide structure, sequence, and corresponding CCS. The current version is designed to identify peptide sequences up to 23 residues in length, in its higher charge state, based on a match of the molecule m/z and CCS. The protocol was validated against a 128-sequences-dataset and CCS predicted within 2.8% average error.

GROMOS 53A6GLYC, an Improved GROMOS Force Field for Hexopyranose-Based Carbohydrates

An improved parameter set for explicit-solvent simulations of carbohydrates (referred to as GROMOS 53A6GLYC) is presented, allowing proper description of the most stable conformation of all 16 possible aldohexopyranose-based monosaccharides. This set includes refinement of torsional potential parameters associated with the determination of hexopyranose rings conformation by fitting to their corresponding quantum-mechanical profiles. Other parameters, as the rules for third and excluded neighbors, are taken directly from the GROMOS 53A6 force field. Comparisons of the herein presented parameter set to our previous version (GROMOS 45A4), the GLYCAM06 force field, and available NMR data are presented in terms of ring puckering free energies, conformational distribution of the hydroxymethyl group, and glycosidic linkage geometries for 16 selected monosaccharides and eight disaccharides. The proposed parameter modifications have shown a significant improvement for the above-mentioned quantities over the two tested force fields, while retaining full compatibility with the GROMOS 53A6 and 54A7 parameter sets for other classes of biomolecules.

Conformational variability of organophosphorus hydrolase upon soman and paraoxon binding

The bacterial enzyme organophosphorus hydrolase (OPH) exhibits both catalytic and substrate promiscuity. It hydrolyzes bonds in a variety of phosphotriester (P-O), phosphonothioate (P-S), phosphofluoridate (P-F), and phosphonocyanate (F-CN) compounds. However, its catalytic efficiency varies markedly for different substrates, limiting the broad-range application of OPH as catalyst in the bioremediation of pesticides and chemical war agents. In the present study, pK(a) calculations and multiple explicit-solvent molecular dynamics (MD) simulations were performed to characterize and contrast the structural dynamics of OPH bound to two substrates hydrolyzed with very distinct catalytic efficiencies: the nerve agent soman (O-pinacolylmethylphosphonofluoridate) and the pesticide paraoxon (diethyl p-nitrophenyl phosphate). pK(a) calculations for the substrate-bound and unbound enzyme showed a significant pK(a) shift from standard values (ΔpK(a) = ±3 units) for residues His254 and Arg275. MD simulations of protonated His254 revealed a dynamic hydrogen bond network connecting the catalytic residue Asp301 via His254 to Asp232, Asp233, Arg275, and Asp235, and is consistent with a previously postulated proton relay mechanism to ferry protons away from the active site with substrates that do not require activation of the leaving group. Hydrogen bonds between Asp301 and His254 were persistent in the OPH-paraoxon complex but not in the OPH-soman one, suggesting a potential role for such interaction in the more efficient hydrolysis of paraoxon over soman by OPH. These results are in line with previous mutational studies of residue His254, which led to an increase of the catalytic efficiency of OPH over soman yet decreased its efficiency for paraoxon. In addition, comparative analysis of the molecular trajectories for OPH bound to soman and paraoxon suggests that binding of the latter facilitates the conformational transition of OPH from the open to the closed substate promoting a tighter binding of paraoxon.

Chitosan molecular structure as a function of N-acetylation

Molecular dynamics simulations have been carried out to characterize the structure and solubility of chitosan nanoparticle-like structures as a function of the deacetylation level (0, 40, 60, and 100%) and the spatial distribution of the N-acetyl groups in the particles. The polysaccharide chains of highly N-deacetylated particles where the N-acetyl groups are uniformly distributed present a high flexibility and preference for the relaxed two-fold helix and five-fold helix motifs. When these groups are confined to a given region of the particle, the chains adopt preferentially a two-fold helix with ϕ and ψ values close to crystalline chitin. Nanoparticles with up to 40% acetylation are moderately soluble, forming stable aggregates when the N-acetyl groups are unevenly distributed. Systems with 60% or higher N-acetylation levels are insoluble and present similar degrees of swelling regardless the distribution of their N-acetyl groups. Overall particle solvation is highly affected by electrostatic forces resulting from the degree of acetylation. The water mobility and orientation around the polysaccharide chains affects the stability of the intramolecular O3-HO3((n)) ···O5((n +) (1)) hydrogen bond, which in turn controls particle aggregation.

The role of nonbonded interactions in the conformational dynamics of organophosphorous hydrolase adsorbed onto functionalized mesoporous silica surfaces

The enzyme organophosphorous hydrolase (OPH) catalyzes the hydrolysis of a wide variety of organophosphorous compounds with high catalytic efficiency and broad substrate specificity. The immobilization of OPH in functionalized mesoporous silica (FMS) surfaces increases significantly its catalytic specific activity, as compared to the enzyme in solution, with important applications for the detection and decontamination of insecticides and chemical warfare agents. Experimental measurements of immobilization efficiency as a function of the charge and coverage percentage of different functional groups have been interpreted as electrostatic forces being the predominant interactions underlying the adsorption of OPH onto FMS surfaces. Explicit solvent molecular dynamics simulations have been performed for OPH in bulk solution and adsorbed onto two distinct interaction potential models of the FMS functional groups to investigate the relative contributions of nonbonded interactions to the conformational dynamics and adsorption of the protein. Our results support the conclusion that electrostatic interactions are responsible for the binding of OPH to the FMS surface. However, these results also show that van der Waals forces are detrimental for interfacial adhesion. In addition, it is found that OPH adsorption onto the FMS models favors a protein conformation whose active site is fully accessible to the substrate, in contrast to the unconfined protein.

In vitro evolution of a peptide with a hematite binding motif that may constitute a natural metal-oxide binding archetype

Phage-display technology was used to evolve peptides that selectively bind to the metal-oxide hematite (Fe2O3) from a library of approximately 3 billion different polypeptides. The sequences of these peptides contained the highly conserved amino acid motif, Ser/Thr-hydrophobic/aromatic-Ser/Thr-Pro-Ser/Thr. To better understand the nature of the peptide-metal oxide binding demonstrated by these experiments, molecular dynamics simulations were carried out for Ser-Pro-Ser at a hematite surface. These simulations show that hydrogen bonding occurs between the two serine amino acids and the hydroxylated hematite surface and that the presence of proline between the hydroxide residues restricts the peptide flexibility, thereby inducing a structural-binding motif. A search of published sequence data revealed that the binding motif (Ser/Thr-Pro-Ser/Thr) is adjacent to the terminal heme-binding domain of both OmcA and MtrC, which are outer membrane cytochromes from the metal-reducing bacterium Shewanella oneidensis MR-1. The entire five amino acid consensus sequence (Ser/Thr-hydrophobic/ aromatic-Ser/Thr-Pro-Ser/Thr) was also found as multiple copies in the primary sequences of metal-oxide binding proteins Sil1 and Sil2 from Thalassiosira pseudonana. We suggest that this motif constitutes a natural metal-oxide binding archetype that could be exploited in enzyme-based biofuel cell design and approaches to synthesize tailored metal-oxide nanostructures.

Protein under pressure: Molecular dynamics simulation of the arc repressor

Experimental nuclear magnetic resonance results for the Arc Repressor have shown that this dimeric protein dissociates into a molten globule at high pressure. This structural change is accompanied by a modification of the hydrogen-bonding pattern of the intermolecular beta-sheet: it changes its character from intermolecular to intramolecular with respect to the two monomers. Molecular dynamics simulations of the Arc Repressor, as a monomer and a dimer, at elevated pressure have been performed with the aim to study this hypothesis and to identify the major structural and dynamical changes of the protein under such conditions. The monomer appears less stable than the dimer. However, the complete dissociation has not been seen because of the long timescale needed to observe this phenomenon. In fact, the protein structure altered very little when increasing the pressure. It became slightly compressed and the dynamics of the side-chains and the unfolding process slowed down. Increasing both, temperature and pressure, a tendency of conversion of intermolecular into intramolecular hydrogen bonds in the beta-sheet region has been detected, supporting the mentioned hypothesis. Also, the onset of denaturation of the separated chains was observed.

Variational particle number approach for rational compound design

Within density functional theory, a variational particle number approach for rational compound design (RCD) is presented. An expression for RCD is obtained in terms of minimization of a suitably defined energy penalty functional whose gradients are the nuclear and the electronic chemical potential. Using combined quantum and molecular mechanics, a nonpeptidic anticancer drug candidate is designed.

An improved nucleic acid parameter set for the GROMOS force field

Over the past decades, the GROMOS force field for biomolecular simulation has primarily been developed for performing molecular dynamics (MD) simulations of polypeptides and, to a lesser extent, sugars. When applied to DNA, the 43A1 and 45A3 parameter sets of the years 1996 and 2001 produced rather flexible double-helical structures, in which the Watson-Crick hydrogen-bonding content was more limited than expected. To improve on the currently available parameter sets, the nucleotide backbone torsional-angle parameters and the charge distribution of the nucleotide bases are reconsidered based on quantum-chemical data. The new 45A4 parameter set resulting from this refinement appears to perform well in terms of reproducing solution NMR data and canonical hydrogen bonding. The deviation between simulated and experimental observables is now of the same order of magnitude as the uncertainty in the experimental values themselves.