Recent advances in anion–π interactions

Over the past 10 years, anion–π interaction has been recognized as an important weak force that may occur between anionic systems and electron- deficient aromatics. Lately, this supramolecular contact has experienced a rapidly growing interest, as reflected by numerous recent literature reports. The present paper highlights the tremendous progress achieved in the field by emphasizing three important studies involving anion–π interactions published in 2010. In addition, a pioneering search of the Protein Data Bank (PDB) reveals short anion–π contacts in some protein structures

Prediction using step-wise L1, L2 regularization and feature selection for small data sets with large number of features

<p>Abstract</p> <p>Background</p> <p>Machine learning methods are nowadays used for many biological prediction problems involving drugs, ligands or polypeptide segments of a protein. In order to build a prediction model a so called training data set of molecules with measured target properties is needed. For many such problems the size of the training data set is limited as measurements have to be performed in a wet lab. Furthermore, the considered problems are often complex, such that it is not clear which molecular descriptors (features) may be suitable to establish a strong correlation with the target property. In many applications all available descriptors are used. This can lead to difficult machine learning problems, when thousands of descriptors are considered and only few (e.g. below hundred) molecules are available for training.</p> <p>Results</p> <p>The CoEPrA contest provides four data sets, which are typical for biological regression problems (few molecules in the training data set and thousands of descriptors). We applied the same two-step training procedure for all four regression tasks. In the first stage, we used optimized L1 regularization to select the most relevant features. Thus, the initial set of more than 6,000 features was reduced to about 50. In the second stage, we used only the selected features from the preceding stage applying a milder L2 regularization, which generally yielded further improvement of prediction performance. Our linear model employed a soft loss function which minimizes the influence of outliers.</p> <p>Conclusions</p> <p>The proposed two-step method showed good results on all four CoEPrA regression tasks. Thus, it may be useful for many other biological prediction problems where for training only a small number of molecules are available, which are described by thousands of descriptors.</p

Connectivity independent protein-structure alignment: a hierarchical approach

BACKGROUND: Protein-structure alignment is a fundamental tool to study protein function, evolution and model building. In the last decade several methods for structure alignment were introduced, but most of them ignore that structurally similar proteins can share the same spatial arrangement of secondary structure elements (SSE) but differ in the underlying polypeptide chain connectivity (non-sequential SSE connectivity). RESULTS: We perform protein-structure alignment using a two-level hierarchical approach implemented in the program GANGSTA. On the first level, pair contacts and relative orientations between SSEs (i.e. α-helices and β-strands) are maximized with a genetic algorithm (GA). On the second level residue pair contacts from the best SSE alignments are optimized. We have tested the method on visually optimized structure alignments of protein pairs (pairwise mode) and for database scans. For a given protein structure, our method is able to detect significant structural similarity of functionally important folds with non-sequential SSE connectivity. The performance for structure alignments with strictly sequential SSE connectivity is comparable to that of other structure alignment methods. CONCLUSION: As demonstrated for several applications, GANGSTA finds meaningful protein-structure alignments independent of the SSE connectivity. GANGSTA is able to detect structural similarity of protein folds that are assigned to different superfamilies but nevertheless possess similar structures and perform related functions, even if these proteins differ in SSE connectivity

Protonation State-Dependent Communication in Cytochrome c Oxidase

Proton transfer in cytochrome c oxidase from the cellular inside to the binuclear redox center (BNC) can occur through two distinct pathways, the D- and K-channels. For the protein to function as both a redox enzyme and a proton pump, proton transfer into the protein toward the BNC or toward a proton loading site (and ultimately through the membrane) must be highly regulated. The PR → F transition is the first step in a catalytic cycle that requires proton transfer from the bulk at the N-side to the BNC. Molecular dynamics simulations of the PR → F intermediate of this transition, with 16 different combinations of protonation states of key residues in the D- and K-channel, show the impact of the K-channel on the D-channel to be protonation-state dependent. Strength as well as means of communication, correlations in positions, or communication along the hydrogen-bonded network depends on the protonation state of the K-channel residue K362. The conformational and hydrogen-bond dynamics of the D-channel residue N139 is regulated by an interplay of protonation in the D-channel and K362. N139 thus assumes a gating function by which proton passage through the D-channel toward E286 is likely facilitated for states with protonated K362 and unprotonated E286. In contrast, proton passage through the D-channel is hindered by N139’s preference for a closed conformation in situations with protonated E286

Insights from surface enhanced resonance Raman spectroscopy and QM/MM calculations

Understanding the coupling between heme reduction and proton translocation in cytochrome c oxidase (CcO) is still an open problem. The propionic acids of heme a3 have been proposed to act as a proton loading site (PLS) in the proton pumping pathway, yet this proposal could not be verified by experimental data so far. We have set up an experiment where the redox states of the two hemes in CcO can be controlled via external electrical potential. Surface enhanced resonance Raman (SERR) spectroscopy was applied to simultaneously monitor the redox state of the hemes and the protonation state of the heme propionates. Simulated spectra based on QM/MM calculations were used to assign the resonant enhanced CH2 twisting modes of the propionates to the protonation state of the individual heme a and heme a3 propionates respectively. The comparison between calculated and measured H2OD2O difference spectra allowed a sound band assignment. In the fully reduced enzyme at least three of the four heme propionates were found to be protonated whereas in the presence of a reduced heme a and an oxidized heme a3 only protonation of one heme a3 propionates was observed. Our data supports the postulated scenario where the heme a3 propionates are involved in the proton pathway

Quantification of Local Electric Field Changes at the Active Site of Cytochrome c Oxidase by Fourier Transform Infrared Spectroelectrochemical Titrations

Cytochrome c oxidase (CcO) is a transmembrane protein complex that reduces molecular oxygen to water while translocating protons across the mitochondrial membrane. Changes in the redox states of its cofactors trigger both O2 reduction and vectorial proton transfer, which includes a proton-loading site, yet unidentified. In this work, we exploited carbon monoxide (CO) as a vibrational Stark effect (VSE) probe at the binuclear center of CcO from Rhodobacter sphaeroides. The CO stretching frequency was monitored as a function of the electrical potential, using Fourier transform infrared (FTIR) absorption spectroelectrochemistry. We observed three different redox states (R4CO, R2CO, and O), determined their midpoint potential, and compared the resulting electric field to electrostatic calculations. A change in the local electric field strength of +2.9 MV/cm was derived, which was induced by the redox transition from R4CO to R2CO. We performed potential jump experiments to accumulate the R2CO and R4CO species and studied the FTIR difference spectra in the protein fingerprint region. The comparison of the experimental and computational results reveals that the key glutamic acid residue E286 is protonated in the observed states, and that its hydrogen-bonding environment is disturbed upon the redox transition of heme a3. Our experiments also suggest propionate A of heme a3 changing its protonation state in concert with the redox state of a second cofactor, heme a. This supports the role of propionic acid side chains as part of the proton-loading site

The redox-coupled proton-channel opening in cytochrome c oxidase

Cytochrome c oxidase (CcO), a redox-coupled proton pump, catalyzes the reduction of molecular oxygen to water, thereby establishing the transmembrane proton gradient that fuels ATP synthesis. CcO employs two channels for proton uptake, the D- and the K-channel. In contrast to the D-channel, the K-channel does not constitute a continuous pathway of H-bonds for proton conduction and is only active in the reductive phase rendering its proton transport mechanism enigmatic. Theoretical studies have suggested selective hydration changes within the K-channel to become activated and being essential for vectorial proton transport. Here, we unravel a previously unidentified mechanism for transient proton channel activation by combining computational studies with site-directed nano-environmental probing of protonation, structural changes, and water dynamics. We show that electrostatic changes at the binuclear center lead to long-range conformational changes propagating to the K-channel entrance as evidenced by time-resolved fluorescence depolarization experiments and molecular dynamics (MD) simulations. These redox-induced long-range structural rearrangements affect the H-bond network at the K-channel's protein surface as shown by pKa-shift analysis of a local probe in experiment and simulation. Concomitantly, selective channel hydration at the K-channel entrance was revealed by dipolar relaxation studies to be associated with channel opening. We propose that instead of a singular change, it is the intricate interplay of these individual redox-triggered changes in the cause–effect relationship that defines the mechanism for transient proton conduction of the K-channel

Supplementary data for the article: Dragelj, J. Lj.; Stanković, I. M.; Božinovski, D. M.; Meyer, T.; Veljković, D. Z.; Medaković, V. B.; Knapp, E.-W.; Zarić, S. D. C-H/O Interactions of Aromatic CH Donors within Proteins: A Crystallographic Study. Crystal Growth and Design 2016, 16 (4), 1948–1957. https://doi.org/10.1021/acs.cgd.5b01543

Supplementary material for: [https://doi.org/10.1021/acs.cgd.5b01543]Related to published version: [http://cherry.chem.bg.ac.rs/handle/123456789/1910

Protein dynamics with off-lattice Monte Carlo moves

A Monte Carlo method for dynamics simulation of all-atom protein models is introduced, to reach long times not accessible to conventional molecular dynamics. The considered degrees of freedom are the dihedrals at C$_\alpha$-atoms. Two Monte Carlo moves are used: single rotations about torsion axes, and cooperative rotations in windows of amide planes, changing the conformation globally and locally, respectively. For local moves Jacobians are used to obtain an unbiased distribution of dihedrals. A molecular dynamics energy function adapted to the protein model is employed. A polypeptide is folded into native-like structures by local but not by global moves.Comment: 10 pages, 4 Postscript figures, uses epsf.sty and a4.sty; scheduled tentatively for Phys.Rev.E issue of 1 March 199