Reference Data Set for Circular Dichroism Spectroscopy Comprised of Validated Intrinsically Disordered Protein Models

Circular dichroism (CD) spectroscopy is an analytical technique that measures the wavelength-dependent differential absorbance of circularly polarized light and is applicable to most biologically important macromolecules, such as proteins, nucleic acids, and carbohydrates. It serves to characterize the secondary structure composition of proteins, including intrinsically disordered proteins, by analyzing their recorded spectra. Several computational tools have been developed to interpret protein CD spectra. These methods have been calibrated and tested mostly on globular proteins with well-defined structures, mainly due to the lack of reliable reference structures for disordered proteins. It is therefore still largely unclear how accurately these computational methods can determine the secondary structure composition of disordered proteins. Here, we provide such a required reference data set consisting of model structural ensembles and matching CD spectra for eight intrinsically disordered proteins. Using this set of data, we have assessed the accuracy of several published CD prediction and secondary structure estimation tools, including our own CD analysis package, SESCA. Our results show that for most of the tested methods, their accuracy for disordered proteins is generally lower than for globular proteins. In contrast, SESCA, which was developed using globular reference proteins, but was designed to be applicable to disordered proteins as well, performs similarly well for both classes of proteins. The new reference data set for disordered proteins should allow for further improvement of all published methods.

An Improved Diabatization Scheme for Computing the Electronic Circular Dichroism of Proteins

We advance the quality of first-principles calculations of protein electronic circular dichroism (CD) through an amelioration of a key deficiency of a previous procedure that involved diabatization of electronic states on the amide chromophore (to obtain interamide couplings) in a β-strand conformation of a diamide. This yields substantially improved calculated far-ultraviolet (far-UV) electronic circular dichroism (CD) spectra for β-sheet conformations. The interamide couplings from the diabatization procedure for 13 secondary structural elements (13 diamide structures) are applied to compute the CD spectra for seven example proteins: myoglobin (α helix), jacalin (β strand), concanavalin A (β type I), elastase (β type II), papain (α + β), 310-helix bundle (310-helix) and snow flea antifreeze protein (polyproline). In all cases, except concanavalin A and papain, the CD spectra computed using the interamide couplings from the diabatization procedure yield improved agreement with experiment with respect to previous first-principles calculations.

KCD: A prediction web server of knowledge-based circular dichroism

We present a web server that predicts the far-UV circular dichroism (CD) spectra of proteins by utilizing their three-dimensional (3D) structures from the Protein Data Bank (PDB). The main algorithm is based on the classical theory of optical activity together with a set of atomic complex polarizabilities, which are obtained from the analysis of a series of synchrotron radiation CD spectra and their related 3D structures from the PDB. The results of our knowledge-based CD method (KCD) are in good agreement with measured spectra that could include the effect of D-amino acids. Our method also delivers some of the most accurate predictions, in comparison with the calculated spectra from well-established models. Specifically, using a metric of closeness based on normalized absolute deviations between experimental and calculated spectra, the mean values for a series of 57 test proteins give the following figures for such models: 0.26 KCD, 0.27 PDBMD2CD, 0.30 SESCA, and 0.47 DichroCalc. From another point of view, it is worth mentioning the remarkable capabilities of the recent approaches based on artificial intelligence, which can precisely predict the native structure of proteins. The structure of proteins, however, is flexible and can be modified by a diversity of environmental factors such as interactions with other molecules, mechanical stresses, variations of temperature, pH, or ionic strength. Experimental CD spectra together with reliable predictions can be utilized to assess eventual secondary structural changes. A similar kind of evaluation can be done for the case of an incomplete protein structure that has been reconstructed by using different approaches. The KCD method can be freely accessed from: https://kcd.cinvestav.mx/.

Theoretical analysis of divalent cation effects on aptamer recognition of neurotransmitter targets

Aptamer-based sensing of small molecules such as dopamine and serotonin in the brain, requires characterization of the specific aptamer sequences in solutions mimicking the in vivo environment with physiological ionic concentrations. In particular, divalent cations (Mg2+ and Ca2+) present in brain fluid, have been shown to affect the conformational dynamics of aptamers upon target recognition. Thus, for biosensors that transduce aptamer structure switching as the signal response, it is critical to interrogate the influence of divalent cations on each unique aptamer sequence. Herein, we demonstrate the potential of molecular dynamics (MD) simulations to predict the behaviour of dopamine and serotonin aptamers on sensor surfaces. The simulations enable molecular-level visualization of aptamer conformational changes that, in some cases, are significantly influenced by divalent cations. The correlations of theoretical simulations with experimental findings validate the potential for MD simulations to predict aptamer-specific behaviors on biosensors.

Tidying up the conformational ensemble of a disordered peptide by computational prediction of spectroscopic fingerprints

The most advanced structure prediction methods are powerless in exploring the conformational ensemble of disordered peptides and proteins and for this reason the "protein folding problem" remains unsolved. We present a novel methodology that enables the accurate prediction of spectroscopic fingerprints (circular dichroism, infrared, Raman, and Raman optical activity), and by this allows for "tidying up" the conformational ensembles of disordered peptides and disordered regions in proteins. This concept is elaborated for and applied to a dodecapeptide, whose spectroscopic fingerprint is measured and theoretically predicted by means of enhanced-sampling molecular dynamics coupled with quantum mechanical calculations. Following this approach, we demonstrate that peptides lacking a clear propensity for ordered secondary-structure motifs are not randomly, but only conditionally disordered. This means that their conformational landscape, or phase-space, can be well represented by a basis-set of conformers including about 10 to 100 structures. The implications of this finding have profound consequences both for the interpretation of experimental electronic and vibrational spectral features of peptides in solution and for the theoretical prediction of these features using accurate and computationally expensive techniques. The here-derived methods and conclusions are expected to fundamentally impact the rationalization of so-far elusive structure-spectra relationships for disordered peptides and proteins, towards improved and versatile structure prediction methods.

DichroIDP: a method for analyses of intrinsically disordered proteins using circular dichroism spectroscopy

Intrinsically disordered proteins (IDPs) are comprised of significant numbers of residues that form neither helix, sheet, nor any other canonical type of secondary structure. They play important roles in a broad range of biological processes, such as molecular recognition and signalling, largely due to their chameleon-like ability to change structure from unordered when free in solution to ordered when bound to partner molecules. Circular dichroism (CD) spectroscopy is a widely-used method for characterising protein secondary structures, but analyses of IDPs using CD spectroscopy have suffered because the methods and reference datasets used for the empirical determination of secondary structures do not contain adequate representations of unordered structures. This work describes the creation, validation and testing of a standalone Windows-based application, DichroIDP, and a new reference dataset, IDP175, which is suitable for analyses of proteins containing significant amounts of disordered structure. DichroIDP enables secondary structure determinations of IDPs and proteins containing intrinsically disordered regions.

Therapeutic phosphorodiamidate morpholino oligonucleotides: Physical properties, solution structures, and folding thermodynamics

Elucidating the structure-function relationships for therapeutic RNA mimicking phosphorodiamidate morpholino oligonucleotides (PMOs) is challenging due to the lack of information about their structures. While PMOs have been approved by the US Food and Drug Administration for treatment of Duchenne muscular dystrophy, no structural information on these unique, charge-neutral, and stable molecules is available. We performed circular dichroism and solution viscosity measurements combined with molecular dynamics simulations and machine learning to resolve solution structures of 22-mer, 25-mer, and 30-mer length PMOs. The PMO conformational dynamics are defined by the competition between non-polar nucleobases and uncharged phosphorodiamidate groups for shielding from solvent exposure. PMO molecules form non-canonical, partially helical, stable folded structures with a small 1.4- to 1.7-nm radius of gyration, low count of three to six base pairs and six to nine base stacks, characterized by -34 to -51 kcal/mol free energy, -57 to -103 kcal/mol enthalpy, and -23 to -53 kcal/mol entropy for folding. The 4.5- to 6.2-cm³/g intrinsic viscosity and Huggins constant of 4.5-9.9 are indicative of extended and aggregating systems. The results obtained highlight the importance of the conformational ensemble view of PMO solution structures, thermodynamic stability of their non-canonical structures, and concentration-dependent viscosity properties. These principles form a paradigm to understand the structure-properties-function relationship for therapeutic PMOs to advance the design of new RNA-mimic-based drugs.

A portable elliptical dichroism spectrometer targeting secondary structural features of tumorous protein for pancreatic cancer detection

Stereochemical analysis is essential for understanding the complex function of biomolecules. Various direct and indirect approaches can be used to explore the allosteric configuration. However, the size, cost, and delicate nature of these systems limit their biomedical usage. Here, we constructed elliptical dichroism (ED) spectrometer for biomedical applications, whose performance is validated by experiment and theoretical simulation (Jones/Mueller calculus and time-dependent density-functional theory). Instead of complicated control of circular polarization, ED spectrometer adopted the absorbance of left- and right-oriented elliptically polarized light. With a simplified design, we demonstrated the potential of ED spectrometry as an alternative for secondary structural analysis of biomolecules, their conformation and chirality. It not only provides a portable, low-cost alternative to the sophisticated instruments currently used for structural analysis of biomolecules but also provides superior translational features: low sample consumption(200 μl), easy operation, and multiple working modes, for noninvasive cancer detection.

Anion-Sensing Properties of Cyclopentaphenylalanine

Cyclic pentaphenylalanine was studied as an efficient anion sensor for halides, thiocyanate and oxoanions in acetonitrile and methanol. Stability constants of the corresponding complexes were determined by means of fluorimetric, spectrophotometric, ¹H NMR, and microcalorimetric titrations. A detailed structural overview of receptor–anion complexes was obtained by classical molecular dynamics (MD) simulations. The results of ¹H NMR and MD studies indicated that the bound anions were coordinated by the amide groups of cyclopeptide, as expected. Circular dichroism (CD) titrations were also carried out in acetonitrile. To the best of our knowledge, this is the first example of the detection of anion binding by cyclopeptide using CD spectroscopy. The CD spectra were calculated from the structures obtained by MD simulations and were qualitatively in agreement with the experimental data. The stoichiometry of almost all complexes was 1:1 (receptor:anion), except for dihydrogen phosphate where the binding of dihydrogen phosphate dimer was observed in acetonitrile. The affinity of the cyclopeptide receptor was correlated with the structure of anion coordination sphere, as well as with the solvation properties of the examined solvents.

Exploring the Parallel G-Quadruplex Nucleic Acid World: A Spectroscopic and Computational Investigation on the Binding of the c-myc Oncogene NHE III1 Region by the Phytochemical Polydatin

Trans-polydatin (tPD), the 3-β-D-glucoside of the well-known nutraceutical trans-resveratrol, is a natural polyphenol with documented anti-cancer, anti-inflammatory, cardioprotective, and immunoregulatory effects. Considering the anticancer activity of tPD, in this work, we aimed to explore the binding properties of this natural compound with the G-quadruplex (G4) structure formed by the Pu22 [d(TGAGGGTGGGTAGGGTGGGTAA)] DNA sequence by exploiting CD spectroscopy and molecular docking simulations. Pu22 is a mutated and shorter analog of the G4-forming sequence known as Pu27 located in the promoter of the c-myc oncogene, whose overexpression triggers the metabolic changes responsible for cancer cells transformation. The binding of tPD with the parallel Pu22 G4 was confirmed by CD spectroscopy, which showed significant changes in the CD spectrum of the DNA and a slight thermal stabilization of the G4 structure. To gain a deeper insight into the structural features of the tPD-Pu22 complex, we performed an in silico molecular docking study, which indicated that the interaction of tPD with Pu22 G4 may involve partial end-stacking to the terminal G-quartet and H-bonding interactions between the sugar moiety of the ligand and deoxynucleotides not included in the G-tetrads. Finally, we compared the experimental CD profiles of Pu22 G4 with the corresponding theoretical output obtained using DichroCalc, a web-based server normally used for the prediction of proteins’ CD spectra starting from their “.pdb” file. The results indicated a good agreement between the predicted and the experimental CD spectra in terms of the spectral bands’ profile even if with a slight bathochromic shift in the positive band, suggesting the utility of this predictive tool for G4 DNA CD investigations.

α-Synuclein phosphorylation at serine 129 occurs after initial protein deposition and inhibits seeded fibril formation and toxicity

SignificanceConverging evidence points to the build-up of phosphorylated α-synuclein (α-syn) at residue serine 129 (pS129) in Lewy body disease, suggesting its central role in the regulation of α-syn aggregation and neuronal degeneration. However, a comprehensive understanding of the role of α-syn phosphorylation at pS129 in α-synuclenopathies pathogenesis is still lacking. Herein, we study the phosphorylation incidence and its effect on α-syn aggregation propensity and cellular toxicity. Collectively, our data suggest that pS129 occurred subsequent to initial α-syn aggregation, lessened aggregation propensity, and attenuated cytotoxicity through diverse assays. Our findings highlight major implications for a better understanding of the role of a molecular modification on protein aggregation.

Novel Design of Neuropeptide-Based Drugs with β-Sheet Breaking Potential in Amyloid-Beta Cascade: Molecular and Structural Deciphers

Our work discusses the investigation of 75 peptide-based drugs with the potential ability to break the β-sheet structures of amyloid-beta peptides from senile plaques. Hence, this study offers a unique insight into the design of neuropeptide-based drugs with β-sheet breaker potential in the amyloid-beta cascade for Alzheimer’s disease (AD). We started with five peptides (¹⁵QKLVFF²⁰, ¹⁶KLVFF²⁰, ¹⁷LVFF²⁰, ¹⁶KLVF¹⁹ and ¹⁵QKLV¹⁸), to which 14 different organic acids were attached at the N-terminal. It was necessary to evaluate the physiochemical features of these sequences due to the biological correlation with our proposal. Hence, the preliminary analysis of different pharmacological features provided the necessary data to select the peptides with the best biocompatibility for administration purposes. Our approaches demonstrated that the peptides ¹⁷LVFF²⁰, NA-¹⁷LVFF²⁰, ¹⁶KLVF¹⁹ and NA-¹⁶KLVF¹⁹ (NA-nicotinic acid) have the ability to interfere with fibril formation and hence improve the neuro and cognitive functions. Moreover, the peptide conjugate NA-¹⁶KLVF¹⁹ possesses attractive pharmacological properties, demonstrated by in silico and in vitro studies. Tandem mass spectrometry showed no fragmentation for the spectra of ¹⁶KLVF¹⁹. Such important results suggest that under the action of protease, the peptide cleavage does not occur at all. Additionally, circular dichroism confirmed docking simulations and showed that NA-¹⁶KLVF¹⁹ may improve the β-sheet breaker mechanism, and thus the entanglement process of amyloid-beta peptides can be more effective.

Structural Characteristics of Glycocins: Unraveling the Role of S-Linked Carbohydrates and Other Structural Elements

Glycocins are antimicrobial peptides with glycosylations, often an S-linked monosaccharide. Their recent structure elucidation has brought forth questions about their mechanisms of action as well as the impact of S-glycosylation on their structural behavior. Here, we investigated structural characteristics of glycocins using a computational approach. Depending on the peptide's class (sublancin- or glycocin F-like), the sugar changes the peptide's flexibility. Also, the presence of glycosylation is necessary for the lack of structure of Asm1. The C-terminal tail in glycocin F-like peptides influenced their structured regions, acting like a regulator. These findings corroborate the versatility of these post-translational modifications, pointing toward their potential use in molecular engineering.

Knowledge-Based Atomic Polarizabilities Used to Model Circular Dichroism Spectra of Proteins

We present a model of circular dichroism for proteins, which is mainly based on both the classical theory of optical activity and a series of effective atomic polarizabilities. Such polarizabilities are extracted from the analysis of a set of synchrotron radiation circular dichroism spectra and their corresponding three-dimensional structures from the Protein Data Bank. Each modeled spectrum is obtained from the protein atomic coordinates and the identification of its secondary structure elements. The resulting spectra are in good agreement with additional experimental data and also with the predictions of some other models. Among them, only our approach is able to describe the effect of d-amino acids. Moreover, our model is also utilized to evaluate protein reconstructions as well as structural changes.

Probing Intermolecular Interactions of Amyloidogenic Fragments of SOD1 by Site-Specific Tryptophan and Its Noncanonical Derivative

Transient amyloid intermediates are likely to be cytotoxic and play an essential role in amyloid-associated neurodegenerative diseases. Characterization of their structural and dynamic evolution is the key to elucidating the molecular mechanism of amyloid formation. Here, combining circular dichroism (CD), exciton couplet theory, and Fourier transform infrared spectroscopy with site-specific tryptophan (Trp) and its noncanonical derivative 5-cyano-tryptochan (Trp_5CN), we developed a method to monitor strand-to-strand tertiary and sheet-to-sheet quaternary interactions in the aggregation cascades of an amyloidogenic fragment from protein SOD1_28-38 (with the sequence KVKVWGSIKGL). We found that the exciton couplet generated from the B_b band of Trp can be used as a probe for side chain interactions. Its sensitivity can be further improved by four times with the incorporation of Trp_5CN. We further observed a red-shift of ∼2 cm^-1 and a broadening of ∼2 cm^-1 in the IR band generated from the CN stretch during the aggregation, which we attributed to the transition from a corkscrew-like structure to a cross-linked intermediate phase. We show here that the integration of optical methods with unique aromatic side chain-related probes is able to elucidate amyloid intermolecular interactions and even capture elusive transient intermediates on and off the amyloid assembling pathway.

Benzene, Toluene, and Monosubstituted Derivatives: Diabatic Nature of the Oscillator Strengths of S1 ← S0 Transitions

For benzene, toluene, aniline, fluorobenzene, and phenol, even sophisticated treatments of electron correlation, such as MRCI and XMS-CASPT2 calculations, show oscillator strengths typically lower than experiment. Inclusion of a simple pseudo-diabatization approach to perturb the S1 state with approximate vibronic coupling to the S2 state for each molecule results in more accurate oscillator strengths. Their absolute values agree better with experiment for all molecules except aniline. When the coupling between the S1 and S2 states is strong at the S0 geometry, the simple diabatization scheme performs less well with respect to the oscillator strengths relative to the adiabatic values. However, we expect the scheme to be useful in many cases where the coupling is weak to moderate (where the maximum component of the coupling has a magnitude less than 1.5 au). Such calculations give an insight into the effects of vibronic coupling of excited states on UV/vis spectra.

Near-Ultraviolet Circular Dichroism and Two-Dimensional Spectroscopy of Polypeptides

A fully quantitative theory of the relationship between protein conformation and optical spectroscopy would facilitate deeper insights into biophysical and simulation studies of protein dynamics and folding. In contrast to intense bands in the far-ultraviolet, near-UV bands are much weaker and have been challenging to compute theoretically. We report some advances in the accuracy of calculations in the near-UV, which were realised through the consideration of the vibrational structure of the electronic transitions of aromatic side chains.

Ancient gene duplications in RNA viruses revealed by protein tertiary structure comparisons

To date only a handful of duplicated genes have been described in RNA viruses. This shortage can be attributed to different factors, including the RNA viruses with high mutation rate that would make a large genome more prone to acquire deleterious mutations. This may explain why sequence-based approaches have only found duplications in their most recent evolutionary history. To detect earlier duplications, we performed protein tertiary structure comparisons for every RNA virus family represented in the Protein Data Bank. We present a list of thirty pairs of possible paralogs with <30 per cent sequence identity. It is argued that these pairs are the outcome of six duplication events. These include the α and β subunits of the fungal toxin KP6 present in the dsRNA Ustilago maydis virus (family Totiviridae), the SARS-CoV (Coronaviridae) nsp3 domains SUD-N, SUD-M and X-domain, the Picornavirales (families Picornaviridae, Dicistroviridae, Iflaviridae and Secoviridae) capsid proteins VP1, VP2 and VP3, and the Enterovirus (family Picornaviridae) 3C and 2A cysteine-proteases. Protein tertiary structure comparisons may reveal more duplication events as more three-dimensional protein structures are determined and suggests that, although still rare, gene duplications may be more frequent in RNA viruses than previously thought. Keywords: gene duplications; RNA viruses.

Combining Experimental Data and Computational Methods for the Non-Computer Specialist

Experimental methods are indispensable for the study of the function of biological macromolecules, not just as static structures, but as dynamic systems that change conformation, bind partners, perform reactions, and respond to different stimulus. However, providing a detailed structural interpretation of the results is often a very challenging task. While experimental and computational methods are often considered as two different and separate approaches, the power and utility of combining both is undeniable. The integration of the experimental data with computational techniques can assist and enrich the interpretation, providing new detailed molecular understanding of the systems. Here, we briefly describe the basic principles of how experimental data can be combined with computational methods to obtain insights into the molecular mechanism and expand the interpretation through the generation of detailed models.

How accurate is circular dichroism-based model validation?

Circular dichroism (CD) spectroscopy is highly sensitive to the secondary structure (SS) composition of proteins. Several methods exist to either estimate the SS composition of a protein or to validate existing structural models using its CD spectrum. The accuracy and precision of these methods depend on the quality of both the measured CD spectrum and the used reference structure. Using a large reference protein set with high-quality CD spectra and synthetic data derived from this set, we quantified deviations from both ideal spectra and reference structures due to experimental limitations. We also determined the impact of these deviations on SS estimation, CD prediction, and SS validation methods of the SESCA analysis package. With regard to the CD spectra, our results suggest intensity scaling errors and non-SS contributions as the main causes of inaccuracies. These factors also can lead to overestimated model errors during validation. The errors of the used reference structures combine non-additively with errors caused by the CD spectrum, which increases the uncertainty of model validation. We have further shown that the effects of scaling errors in the CD spectrum can be nearly eliminated by appropriate re-scaling, and that the accuracy of model validation methods can be improved by accounting for typical non-SS contributions. These improvements have now been implemented within the SESCA package and are available at: https://www.mpibpc.mpg.de/sesca .

Computed optical spectra of SARS-CoV-2 proteins

•

Calculated circular dichroism spectra in the far- and near-UV spectra.

•

Calculated infra-red (IR) spectra in the amide I region.

•

Based on experimental structures and computational models of SARS-CoV-2 proteins.

•

Near-UV CD spectra offer greatest sensitivity to conformation.

Treatment for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus that causes Covid-19, may well be predicated on knowledge of the structures of protein of this virus. However, often these cannot be determined easily or quickly. Herein, we provide calculated circular dichroism (CD) spectra in the far- and near-UV, and infra-red (IR) spectra in the amide I region for experimental structures and computational models of SARS-CoV-2 proteins. The near-UV CD spectra offer greatest sensitivity in assessing the accuracy of models.

Unfolding Dynamics of a Photoswitchable Helical Peptide

We present an atomistic force field for the azo-moiety of the photoswitchable FK-11-X peptide. We use the parameters to study the unfolding of the peptide through molecular dynamics simulations. The unfolded ensemble contains many different structures, ranging from a partially unfolded peptide to a fully unfolded structure. The averaged computed far-ultraviolet circular dichroism (CD) spectrum of the set of structures, which was simulated using the newly developed force field, agrees well with experiment. The rate of the simulated unfolding process was estimated to have a time constant of 5.80 ± 0.03 ns from the time evolution of the CD spectrum of the peptide, computed from the backbone conformations sampled over 40 simulated trajectories. Our estimated time constant is faster than, but not inconsistent with, previous experimental estimates from time-resolved infrared and optical rotatory dispersion spectroscopy.

Mitigation of peri-implantitis by rational design of bifunctional peptides with antimicrobial properties

The integration of molecular and cell biology with materials science has led to strategies to improve the interface between dental implants with the surrounding soft and hard tissues in order to replace missing teeth and restore mastication. More than 3 million implants have been placed in the US alone and this number is rising by 500,000/year. Peri-implantitis, an inflammatory response to oral pathogens growing on the implant surface threatens to reduce service life leading to eventual implant failure, and such an outcome will have adverse impact on public health and create significant health care costs. Here we report a predictive approach to peptide design, which enabled us to engineer a bifunctional peptide to combat bacterial colonization and biofilm formation, reducing the adverse host inflammatory immune response that destroys the tissue surrounding implants and shortens their lifespans. This bifunctional peptide contains a titanium-binding domain that recognizes and binds with high affinity to titanium implant surfaces, fused through a rigid spacer domain with an antimicrobial domain. By varying the antimicrobial peptide domain, we were able to predict the properties of the resulting bifunctional peptides in their entirety by analyzing the sequence-structure-function relationship. These bifunctional peptides achieve: 1) nearly 100% surface coverage within minutes, a timeframe suitable for their clinical application to existing implants; 2) nearly 100% binding to a titanium surface even in the presence of contaminating serum protein; 3) durability to brushing with a commercially available electric toothbrush; and 4) retention of antimicrobial activity on the implant surface following bacterial challenge. A bifunctional peptide film can be applied to both new implants and/or repeatedly applied to previously placed implants to control bacterial colonization mitigating peri-implant disease that threatens dental implant longevity.

SESCA: Predicting Circular Dichroism Spectra from Protein Molecular Structures

Circular dichroism (CD) spectroscopy is a highly sensitive but low-resolution technique to study the structure of proteins. Combined with molecular modeling or other complementary techniques, CD spectroscopy can provide essential information at higher resolution. To this end, we introduce a new computational method to calculate the electronic circular dichroism spectra of proteins from a structural model or ensemble using the average secondary structure composition and a precalculated set of basis spectra. The method is designed for model validation to estimate the error of a given protein structural model based on the measured CD spectrum. We compared the predictive power of our method to that of existing algorithms, namely, DichroCalc and PDB2CD, and found that it predicts CD spectra more accurately. Our results indicate that the derived basis sets are robust to both experimental errors in the reference spectra and the choice of the secondary structure classification algorithm. For over 80% of the globular reference proteins, our basis sets accurately predict the experimental spectrum solely from their secondary structure composition. For the remaining 20%, correcting for intensity normalization considerably improves the prediction power. Additionally, we show that the predictions for short peptides and an example complex of intrinsically disordered proteins strongly benefit from accounting for side-chain contributions and structural flexibility.

Enzyme-responsive progelator cyclic peptides for minimally invasive delivery to the heart post-myocardial infarction

Injectable biopolymer hydrogels have gained attention for use as scaffolds to promote cardiac function and prevent negative left ventricular (LV) remodeling post-myocardial infarction (MI). However, most hydrogels tested in preclinical studies are not candidates for minimally invasive catheter delivery due to excess material viscosity, rapid gelation times, and/or concerns regarding hemocompatibility and potential for embolism. We describe a platform technology for progelator materials formulated as sterically constrained cyclic peptides which flow freely for low resistance injection, and rapidly assemble into hydrogels when linearized by disease-associated enzymes. Their utility in vivo is demonstrated by their ability to flow through a syringe and gel at the site of MI in rat models. Additionally, synthetic functionalization enables these materials to flow through a cardiac injection catheter without clogging, without compromising hemocompatibility or cytotoxicity. These studies set the stage for the development of structurally dynamic biomaterials for therapeutic hydrogel delivery to the MI.

Similar but Not the Same: First Kinetic and Structural Analyses of a Methanol Dehydrogenase Containing a Europium Ion in the Active Site

Since the discovery of the biological relevance of rare earth elements (REEs) for numerous different bacteria, questions concerning the advantages of REEs in the active sites of methanol dehydrogenases (MDHs) over calcium(II) and of why bacteria prefer light REEs have been a subject of debate. Here we report the cultivation and purification of the strictly REE-dependent methanotrophic bacterium Methylacidiphilum fumariolicum SolV with europium(III), as well as structural and kinetic analyses of the first methanol dehydrogenase incorporating Eu in the active site. Crystal structure determination of the Eu-MDH demonstrated that overall no major structural changes were induced by conversion to this REE. Circular dichroism (CD) measurements were used to determine optimal conditions for kinetic assays, whereas inductively coupled plasma mass spectrometry (ICP-MS) showed 70 % incorporation of Eu in the enzyme. Our studies explain why bacterial growth of SolV in the presence of Eu3+ is significantly slower than in the presence of La3+ /Ce3+ /Pr3+ : Eu-MDH possesses a decreased catalytic efficiency. Although REEs have similar properties, the differences in ionic radii and coordination numbers across the series significantly impact MDH efficiency.

Genetic identification and molecular modeling characterization reveal a novel PROM1 mutation in Stargardt4-like macular dystrophy

Stargardt disease-4 (STGD4) is an autosomal dominant complex, genetically heterogeneous macular degeneration/dystrophy (MD) disorder. In this paper, we used targeted next generation sequencing and multiple molecular dynamics analyses to identify and characterize a disease-causing genetic variant in four generations of a Chinese family with STGD4-like MD. We found a novel heterozygous missense mutation, c.734T>C (p.L245P) in the PROM1 gene. Structurally, this mutation most likely impairs PROM1 protein stability, flexibility, and amino acid interaction network after changing the amino acid residue Leucine into Proline in the basic helix-loop-helix leucine zipper domain. Molecular dynamic simulation and principal component analysis provide compelling evidence that this PROM1 mutation contributes to disease causativeness or susceptibility variants in patients with STGD4-like MD. Thus, this finding defines new approaches in genetic characterization, accurate diagnosis, and prevention of STGD4-like MD.

Copper, differently from zinc, affects the conformation, oligomerization state and activity of bradykinin

The sole role of bradykinin (BK) as an inflammatory mediator is controversial, as recent data also support an anti-inflammatory role for BK in Alzheimer's disease (AD). The involvement of two different receptors (B1R and B2R) could be a key to understand this issue. However, although copper and zinc dyshomeostasis has been demonstrated to be largely involved in the development of AD, a detailed study of the interaction of BK with these two metal ions has never been addressed. In this work, we have applied mass spectrometry, circular dichroism as well as computational methods in order to assess if copper and zinc have the ability to modulate the conformation and oligomerization of BK. In addition, we have correlated the chemical data with the effect of metals on the activity of BK analyzed in cell cultures by biochemical procedures. The biochemical analyses on monocyte/macrophage cell culture (THP-1 Cell Line human) in line with the effect of metals on the conformation of BK showed that the presence of copper can affect the signaling cascade mediated by the BK receptors. The results obtained show a further role of metal ions, particularly copper, in the development and outcome of neuroinflammatory diseases. The possible implications in AD are discussed.

T-dependent B cell responses to Plasmodium induce antibodies that form a high-avidity multivalent complex with the circumsporozoite protein

The repeat region of the Plasmodium falciparum circumsporozoite protein (CSP) is a major vaccine antigen because it can be targeted by parasite neutralizing antibodies; however, little is known about this interaction. We used isothermal titration calorimetry, X-ray crystallography and mutagenesis-validated modeling to analyze the binding of a murine neutralizing antibody to Plasmodium falciparum CSP. Strikingly, we found that the repeat region of CSP is bound by multiple antibodies. This repeating pattern allows multiple weak interactions of single F_AB domains to accumulate and yield a complex with a dissociation constant in the low nM range. Because the CSP protein can potentially cross-link multiple B cell receptors (BCRs) we hypothesized that the B cell response might be T cell independent. However, while there was a modest response in mice deficient in T cell help, the bulk of the response was T cell dependent. By sequencing the BCRs of CSP-repeat specific B cells in inbred mice we found that these cells underwent somatic hypermutation and affinity maturation indicative of a T-dependent response. Last, we found that the BCR repertoire of responding B cells was limited suggesting that the structural simplicity of the repeat may limit the breadth of the immune response.

Vaccines aim to protect by inducing the immune system to make molecules called antibodies that can recognize molecules on the surface of invading pathogens. In the case of malaria, our most advanced vaccine candidates aim to promote the production of antibodies that recognize the circumsporozoite protein (CSP) molecule on the surface of the invasive parasite stage called the sporozoite. In this report we use X-ray crystallography to determine the structure of CSP-binding antibodies at the atomic level. We use other techniques such as isothermal titration calorimetry and structural modeling to examine how this antibody interacts with the CSP molecule. Strikingly, we found that each CSP molecule could bind 6 antibodies. This finding has implications for the immune response and may explain why high titers of antibody are needed for protection. Moreover, because the structure of the CSP repeat is quite simple we determined that the number of different kinds of antibodies that could bind this molecule are quite small. However a high avidity interaction between those antibodies and CSP can result from a process called affinity maturation that allows the body to learn how to make improved antibodies specific for pathogen molecules. These data show that while it is challenging for the immune system to recognize and neutralize CSP, it should be possible to generate viable vaccines targeting this molecule.

Modulation of polypeptide conformation through donor-acceptor transformation of side-chain hydrogen bonding ligands

Synthetic polypeptides have received increasing attention due to their ability to form higher ordered structures similar to proteins. The control over their secondary structures, which enables dynamic conformational changes, is primarily accomplished by tuning the side-chain hydrophobic or ionic interactions. Herein we report a strategy to modulate the conformation of polypeptides utilizing donor-acceptor interactions emanating from side-chain H-bonding ligands. Specifically, 1,2,3-triazole groups, when incorporated onto polypeptide side-chains, serve as both H-bond donors and acceptors at neutral pH and disrupt the α-helical conformation. When protonated, the resulting 1,2,3-triazolium ions lose the ability to act as H-bond acceptors, and the polypeptides regain their α-helical structure. The conformational change of triazole polypeptides in response to the donor-acceptor pattern was conclusively demonstrated using both experimental-based and simulation-based methods. We further showed the utility of this transition by designing smart, cell-penetrating polymers that undergo acid-activated endosomal escape in living cells.Hydrogen bonding plays a major role in determining the tridimensional structure of biopolymers. Here, the authors show that control over a polypeptide conformation can be achieved by altering the donor-acceptor properties of side-chain triazole units via protonation-deprotonation.

Computing infrared spectra of proteins using the exciton model

The ability to compute from first principles the infrared spectrum of a protein in solution phase representing a biological system would provide a useful connection to atomistic models of protein structure and dynamics. Indeed, such calculations are a vital complement to 2DIR experimental measurements, allowing the observed signals to be interpreted in terms of detailed structural and dynamical information. In this article, we have studied nine structurally and spectroscopically well-characterized proteins, representing a range of structural types. We have simulated the equilibrium conformational dynamics in an explicit point charge water model. Using the resulting trajectories based on MD simulations, we have computed the one and two dimensional infrared spectra in the Amide I region, using an exciton approach, in which a local mode basis of carbonyl stretches is considered. The role of solvent in shifting the Amide I band (by 30 to 50 cm-1 ) is clearly evident. Similarly, the conformational dynamics contribute to the broadening of peaks in the spectrum. The inhomogeneous broadening in both the 1D and 2D spectra reflects the significant conformational diversity observed in the simulations. Through the computed 2D cross-peak spectra, we show how different pulse schemes can provide additional information on the coupled vibrations. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

Quantitative first principles calculations of protein circular dichroism in the near-ultraviolet

Vibrational structure in the near-UV circular dichroism (CD) spectra of proteins is an important source of information on protein conformation and can be exploited to study structure and folding. A fully quantitative theory of the relationship between protein conformation and optical spectroscopy would facilitate deeper interpretation of and insight into biophysical and simulation studies of protein dynamics and folding. We have developed new models of the aromatic side chain chromophores toluene, p-cresol and 3-methylindole, which incorporate ab initio calculations of the Franck-Condon effect into first principles calculations of CD using an exciton approach. The near-UV CD spectra of 40 proteins are calculated with the new parameter set and the correlation between the computed and the experimental intensity from 270 to 290 nm is much improved. The contribution of individual chromophores to the CD spectra has been calculated for several mutants and in many cases helps rationalize changes in their experimental spectra. Considering conformational flexibility by using families of NMR structures leads to further improvements for some proteins and illustrates an informative level of sensitivity to side chain conformation. In several cases, the near-UV CD calculations can distinguish the native protein structure from a set of computer-generated misfolded decoy structures.

Understanding nucleic acid structural changes by comparing wide-angle x-ray scattering (WAXS) experiments to molecular dynamics simulations

Wide-angle x-ray scattering (WAXS) is emerging as a powerful tool for increasing the resolution of solution structure measurements of biomolecules. Compared to its better known complement, small angle x-ray scattering (SAXS), WAXS targets higher scattering angles and can enhance structural studies of molecules by accessing finer details of solution structures. Although the extension from SAXS to WAXS is easy to implement experimentally, the computational tools required to fully harness the power of WAXS are still under development. Currently, WAXS is employed to study structural changes and ligand binding in proteins; however, the methods are not as fully developed for nucleic acids. Here, we show how WAXS can qualitatively characterize nucleic acid structures as well as the small but significant structural changes driven by the addition of multivalent ions. We show the potential of WAXS to test all-atom molecular dynamics (MD) simulations and to provide insight into understanding how the trivalent ion cobalt(III) hexammine (CoHex) affects the structure of RNA and DNA helices. We find that MD simulations capture the RNA structural change that occurs due to addition of CoHex.

Engineering Recombinant Protein Sensors for Quantifying Histone Acetylation

H3K14ac (acetylation of lysine 14 of histone H3) is one of the most important epigentic modifications. Aberrant changes in H3K14ac have been associated with various diseases, including cancers and neurological disorders. Tools that enable detection and quantification of H3K14ac levels in cell extracts and in situ are thus of critical importance to reveal its role in various biological processes. Current detection techniques of specific histone modifications, however, are constrained by tedious sample pretreatments, lack of quantitative accuracy, and reliance on high quality antibodies. To address this issue, we engineered recombinant sensors that are suitable for probing histone acetylation levels using various biological samples. The protein sensor contains recongition domain(s) with sequences derived from the bromodomain of human polybromo-1 (PB1), a natural H3K14ac reader domain. Various sensor designs were tested using nuclear extracts and live cells. The sensor containing dimeric repeats of bromodomain was found most effective in quantifying H3K14ac level in both in vitro and in situ assays. The sensor has a linear detection range of 0.5-50 nM when mixed with nuclear extracts. The sensor colocalizes with H3K14ac antibodies in situ when transfected into human embryonic kidney 293T (HEK293T) cells and is thus capable of providing spatial details of histone modification within the nucleus. Corrected nuclear fluorescence intensity was used to quantify the modification level in situ and found to correlate well with our in vitro assays. Our sensor offers a novel tool to characterize the histone modification level using nuclear extracts and probe histone modification change in live cells.